{'fraction_non_alphanumeric': 0.02740266140956136, 'fraction_numerical': 0.022556267455232464, 'mean_word_length': 4.291749978266539, 'pattern_counts': {'":': 0, '<': 3, '<?xml version=': 0, '>': 5, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Gas Stripping in Simulated Galaxies with a Multiphase ISM

Cluster galaxies moving through the intracluster medium (ICM) are expected to lose some of their interstellar medium (ISM) through ISM-ICM interactions. We perform high resolution (40 pc) three-dimensional hydrodynamical simulations of a galaxy undergoing ram pressure stripping including radiative cooling in order to investigate stripping of a multiphase medium. The clumpy, multiphase ISM is self-consistently produced by the inclusion of radiative cooling, and spans six orders of magnitude in gas density. We find no large variations in the amount of gas lost whether or not cooling is involved, although the gas in the multiphase galaxy is stripped more quickly and to a smaller radius. We also see significant differences in the morphology of the stripped disks. This occurs because the multiphase medium naturally includes high density clouds set inside regions of lower density. We find that the lower density gas is stripped quickly from any radius of the galaxy, and the higher density gas can then be ablated. If high density clouds survive, through interaction with the ICM they lose enough angular momentum to drift towards the center of the galaxy where they are no longer stripped. Finally, we find that low ram pressure values compress gas into high density clouds that could lead to enhanced star formation, while high ram pressure leads to a smaller amount of high-density gas.


INTRODUCTION
It has long been known that a galaxy's environment has an impact on its morphology: spirals dominate in the field, and ellipticals and S0s dominate in dense cluster environments (Hubble & Humason 1931). This has been quantified in the density-morphology relation (Oemler 1974;Dressler 1980), but this relation alone does not determine whether the cluster environment affects the formation or the evolution of galaxies. According to the Butcher-Oemler effect, clusters at z ≥ 0.2 are bluer and contain more spirals than nearby clusters, indicating that cluster environments may be correlated with the evolution of galaxies from spirals to earlier types (Butcher & Oemler 1978). More recently, by examining a specific merging cluster, Tran et al. (2005) presents a strong relationship between the Butcher-Oemler effect and infalling galaxies.
The environment-driven evolution of spiral galaxies into S0s is caused by mechanisms that physically remove gas and/or induce star formation until almost no gas remains. Galaxy-intracluster medium (ICM) interactions are ones in which the gas in a galaxy interacts with the ambient intracluster hot gas. One such process is ram pressure stripping (RPS), which removes the interstellar medium (ISM) of a galaxy as it moves through the ICM (Gunn & Gott 1972). Ram pressure could also compress the gas within a galaxy to cause a burst of star formation that would consume all of the remaining gas (Fujita & Nagashima 1999). A galaxy-ICM interaction does not affect the stellar component of a galaxy. The cluster potential can also strip both gas and stars or compress gas to cause an increased star formation rate (Byrd & Valtonen 1990). These galaxy-cluster interactions affect both the stellar and gas components of galaxies. There are also galaxy-galaxy interactions that can take place between galaxies. These include mergers between galaxies with low relative velocities, and galaxy harassment: high-speed interactions between cluster galaxies (Hashimoto et al. 1998;Bekki 1999; Barnes & Hernquist 1991;Bekki 1998;Moore et al. 1996). These interactions can cause an increased star formation rate and will also affect both the stellar and gas component of galaxies.
In our earlier work (Tonnesen, Bryan & van Gorkom 2007), we studied cluster galaxies that had evolved in a cosmological simulation. This simulation included dark matter and stellar particles, and followed gas through an adaptive mesh refinement (AMR) code. It included radiative cooling and star formation. This work was unique because it allowed the cluster to form in a cosmological simulation, but had the resolution to follow individual galaxies (3 kpc). We found that most galaxies that lost all their gas did so without losing stellar mass, indicating that ram pressure stripping was the most effective evolutionary mechanism. This convinced us to look more closely at ram pressure stripping of galaxies.
While a substantial body of literature exists on ram pressure stripping (e.g. Roediger & Hensler 2005;Roediger & Brüggen 2006;Vollmer, Hutchmeier & van Driel 2005; see references in van Gorkom 2004), the effect of a multiphase ISM has only begun to be explored in simulations. Schulz and Struck (2001) included cooling and found that it caused a disk instability that led to an inflow of dense gas, and that including cooling allows the disk to remain cool and thin despite ICM heating. Stripping of a multiphase ISM was also studied at higher resolution by Quilis, Moore & Bower (2000). They did not include radiative cooling, but instead randomly placed holes and overdensities in their ISM disk in an attempt to model an inhomogeneous ISM. They find that stripping is much more complete as ICM streams through the holes in their galaxy and prevents stripped material from falling back onto the galaxy as earlier simulations had shown. More recently, Kronberger et al (2008) have modeled cooling in the context of star formation. They find increased star formation in the center of their galaxy and star formation in the stripped material (although it does not escape the galactic halo).
Observations have been inconclusive about the effect of stripping on dense molecular clouds. In a survey of Virgo cluster galaxies, Kenney & Young (1986, 1989 found that CO remains in HI deficient galaxies, indicating that molecular clouds are not being stripped entirely. Whether this is because molecular clouds reside in the center of galaxies where even atomic gas is rarely stripped, or because molecular clouds are too dense to be stripped is unclear. The likelihood of molecular gas being stripped varies with each individual case. Clemens et al. (2000Clemens et al. ( , 2001 have examined in detail the interacting galaxy pair NGC 4490 and NGC 4485. They find evidence that in moving through the extended HI envelope surrounding the two galaxies, NGC 4485 has had atomic gas, molecular gas, and dust removed via ram pressure stripping. Crowl et al. (2005) have studied the Virgo cluster galaxy NGC 4402, and observe ram pressure stripping and dense cloud ablation. They find that it is possible to strip the lower density ISM without stripping dense clouds, but that the surviving clouds will eventually be ablated by the ICM wind. They also observe an HII region which has likely formed within the stripped gas, so molecular gas was either stripped or formed from stripped material. However, in NGC 4848, Vollmer et al (2001) find that the displaced molecular gas they observe is formed from re-accreted lower density gas, and that the molecular gas was not stripped at all. NGC 4438 has been observed to have both stripped and surviving molecular gas, determined by the clouds' positions in the galaxy . Another galaxy in the Virgo cluster, NGC 4522, has no HII regions in the outer disk, but extraplanar HII regions (Kenney & Koopmann 1999;Kenney, van Gorkom & Vollmer 2004).
As discussed in Kenney et al (2004), there are four possibilities for stripping of molecular clouds. (1) Molecular clouds can be directly stripped from the disk, although there is probably an upper limit to the surface density for which this is possible. (2) Molecular clouds, which are thought to have a lifetime of about 10 7 years (Blitz & Shu 1980;Larson 2003, Hartmann 2003, can also evolve via star formation into lower density gas, which can then be stripped directly. (3) Another possibility is that the molecular clouds are not directly stripped out of the galaxy, but are instead ablated by the ICM wind on a longer timescale (Nulsen 1982). (4) Finally, dense clouds may be coupled to the rest of the ISM by magnetic fields or some other mechanism, so that the entire ISM is stripped out together.
In this paper, we run a set of high resolution simulations (40 pc resolution, which is small enough to marginally resolve Giant Molecular Clouds) to understand how a multiphase ISM could affect ram pressure stripping of a galaxy's gas. It is important to recognize that we do not attempt to include all of the physics involved in the ISM. We use these simulations to understand how the density fluctuations that are observed in the multiphase ISM of galaxies will affect gas loss. We vary ram pressure strength to investigate the result of different wind strengths. Using our simulations, we are able to look in detail at possibilities (1) and (3). We also consider whether including molecular clouds affects the large scale stripping properties such as the total amount of gas stripped and the radius to which gas is stripped.
In an introduction to our code, we provide the general characteristics of our simulation ( §2.1-3). We then discuss our suite of runs and the measurements we use to compare them ( §2.4-5). In §3 we discuss how cooling affects a galaxy evolving in a static ICM. We then ( §4) compare our simulations. In §5 we briefly discuss the limitations and implications of our results. Finally, we conclude in §6 with a summary of our results and the overall picture of ram pressure stripping of a galaxy with a multiphase ISM.

Simulation
We have simulated a galaxy using the adaptive mesh refinement (AMR) code Enzo. Although this code can follow stellar and dark matter particles, we chose to use a fixed stellar and dark matter potential for ease of computation, although the self-gravity of the gas is selfconsistently computed. To follow the gas, we employ an adaptive mesh for solving the fluid equations including gravity (Bryan 1999;Norman & Bryan 1999;O'Shea et al. 2004). The code begins with a fixed, static grid and automatically adds refined grids as required in order to resolve important features in the flow as defined by enhancements in the gas density.
The simulation includes radiative cooling using the Sarazin & White (1987) cooling curve. We allow cooling to 8,000 K, which results in gas with neutral hydrogen temperatures without overcooling a large fraction of our gas. With this cutoff, we find that we still can form gas clouds with densities typical of molecular clouds, although it is clear that the internal structure of such clouds is not reproduced in detail (see Figure 7 for the density mass distribution in the gas disks of our galaxies). We discuss this choice in greater detail in §5.5.
Our box is 311 kpc on a side. Our coarsest resolution is 2.5 kpc. We allow 6 levels of refinement for our runs that include radiative cooling, for a best resolution of 40 pc. The runs without radiative cooling have 5 levels of refinement. We refine our simulation using baryon mass, with our minimum overdensity for refinement set as 51.2, which we found immediately refined the entire galactic disk to our best resolution. We discuss resolution effects in more detail in §5.4. Our galaxy is placed at (0.5,0.5,0.22) or (155.5,155.5,68.42) kpc, so that we can follow the stripped gas as far as possible. The ICM wind flows along the z-axis in the positive direction, with the lower z boundary set for inflow and higher z boundary set as outflow. The x and y boundaries are set to reflecting.

The Galaxy
We model a massive spiral galaxy with a flat rotation curve of 200 km s −1 . It consists of a gas disk that is followed using adaptive mesh refinement (including selfgravity of the gas -a crucial ingredient required to form self-gravitating molecular clouds), as well as the static potentials of a stellar disk, a stellar bulge, and a dark matter halo. We directly follow Roediger & Brüggen (2006) in our modeling of the stellar and dark matter potential and gas disk. Briefly, we model the stellar disk using a Plummer-Kuzmin disk (see Miyamoto & Nagai (Hernquist 1993), and the dark matter halo using the spherical model of Burkert (1995). This dark matter halo model is compatible with observed rotation curves (Burkert 1995;Trachternach et al. 2008). The equation for the analytic potential is in Mori & Burkert (2000). The gas is described as a softened exponential disk: Given this gas density distribution in the disk, the gas temperature and pressure are calculated to maintain the disk in hydrostatic equilibrium with the surrounding ICM in the z direction. The gas disk's rotational velocity is set so that the combination of the centrifugal force and the pressure gradient of the disk balances the radial gravitational force. We cut the gas disk smoothly by multiplying the gas density distribution by 0.5(1+cos(π(R−20 kpc)/26(21) kpc)) for 20 kpc < R ≤ 26(21) kpc. See our galaxy parameters in Tables 1 and 2. We use different initial gas disk radii in order to more easily compare disks when the ICM wind hits the galaxy. We discuss this in greater detail in §3.

ICM Conditions
Our galaxies evolve in a pressurized ICM. Because a galaxy moving through the ICM can be more easily simulated by modeling a fixed galaxy within a moving ICM, our galaxy remains in the same place in our simulated box (at least the stars and dark matter do). The galaxy initially evolves in a static ICM, to examine the stability in the static ICM and to allow cool, dense gas to form. Later, we trigger a constant ICM inflow along the z-axis, which is always face-on to the galaxy. See the Appendix for exactly when the wind hits each galaxy. Briefly, the wind is triggered so that it hits the galaxy while gas is still distributed out to a large radius, but enough time has passed so that high density gas clouds have formed. The exact time was chosen to be when the galaxy in the corresponding RCNW case had gas collapsed to ρ ≥ 10 −20 g cm −3 , typical of densities found in molecular clouds. The measurements of the ICM density and difference in velocity between the galaxies and the ICM for: the cluster in the large simulation we examined in an earlier paper (Tonnesen, Bryan & van Gorkom 2007, as blue triangles), the mean values we used for our three detailed simulations (black diamonds), and the values used by Roediger and Brüggen (2006) (black asterisks). To guide the eye, we also plot lines of constant ram pressure at 1 × 10 −11 dynes cm −2 , 1 × 10 −12 dynes cm −2 , and 1 × 10 −13 dynes cm −2 . See §2.3 for discussion.
We choose to study three ram pressure strengths using realistic values for the ICM density and velocity. The three ram pressure strengths we choose are P ram = ρv 2 ICM = 6.4 × 10 −12 dynes cm −2 , 1 × 10 −12 dynes cm −2 , and 6.4 × 10 −13 dynes cm −2 (see Table 3). We choose the maximum and minimum ram pressure strengths to match two of the values used by Roediger and Brüggen (2006), and we pick 1 × 10 −12 dynes cm −2 because it is the highest ram pressure value at the virial radius of our cluster simulation (Tonnesen et al. 2007).
To choose the corresponding ICM parameters, ρ and v ICM , we use the results from our earlier cluster simulation to find the mean density, velocity, and temperature of the ICM at the three ram pressures (Tonnesen, Bryan & van Gorkom 2007). This is shown in Figure 1, where for comparison we plot the values used by Roediger & Brüggen (2006). As discussed earlier, we set the inflow boundary condition in order to model the ICM wind. However, because we need to allow our galaxies to cool and evolve before the ICM wind begins to strip them, the galaxies are initialized in a static ICM and then, after a multiphase ISM has developed, the boundary values are set to generate a wind with the desired characteristics. To get the parameters for this initial ICM we backtrack from the density, temperature and velocity of the ICM wind using the Rankine-Hugoniot jump conditions, assuming a Mach number of 3.5. We choose this Mach number because it gives us a supersonic shock (which leads to cleaner behavior at the inflow boundary), while allowing our initial static ICM to have a similar density and temperature as the ICM wind. This procedure is important because we compare the evolution of our stripped galaxy to one that only cools in a static medium that has the same density and temperature as the flowing ICM.
We use the same density and temperature as in the moving ICM for the cases in which the ICM always remains static, because the stripped galaxies spend more time in the simulation in the moving ICM.

Suite of Simulations
We perform nine full runs because while attempting to understand how cooling affects galactic gas stripping, we also need to understand how cooling alone affects a galaxy. We examine how three different ram pressure strengths affect galaxies that include and do not include radiative cooling. We choose a naming convention for the simulations which indicates the parameters used for that simulation (see also Table 3). The three ram pressure strengths also correspond to the three ICM thermal pressures, and are denoted PH, PM, PL for high, medium, and low. If cooling is turned on, the next identifier is RC, if not, NC. The six runs with an ICM wind end with W. In addition to these six runs, we also study three runs of galaxies evolving with radiative cooling in a static ICM with the same temperature and density as the moving ICM. All these runs end with NW, for no wind. For details on these nine simulations, see Table  3, as well as the appendix. Because we are interested primarily in the impact of radiative cooling on gas stripping, we only consider cases with a face-on wind. These cases are likely to be the most affected by an ICM wind because the gas density fluctuations in the disk, caused by radiative cooling, are most apparent in a face-on view.

The Measurements
We look closely at three types of measurements: 1) the total amount of gas in the galaxy, 2) the radius of the dense gas in the galaxy, and 3) the amount of gas mass at different densities.
We first measure the total amount of gas that remains in the galaxy. We measure this in two ways: the amount of gas mass in a cylinder centered on the galaxy with a height of 10 kpc and a radius of 27 kpc, and the total amount of bound gas mass in the entire box. Mass is bound if the thermal and kinetic energy are smaller than the gravitational energy from the static potential of the stars and dark matter. We can compare these values to an analytic estimate of how much gas a galaxy should lose using two analytic mass loss estimates, corresponding to two stages of gas stripping. We will follow Roediger & Brüggen (2007) by naming the initial stage ram pressure pushing and the later stage continuous stripping. We can predict the amount of gas loss from ram pressure pushing by measuring the amount of gas mass outside the analytic Gunn & Gott (1972) stripping radius, the radius at which the restoring force per unit area from the galaxy is equal to the ram pressure acting on the galaxy, or where ρv 2 = 2πGΣ star Σ gas . We calculate the radius using the initial gas distribution and measure the mass outside this radius immediately before the ICM wind hits the galaxy. Once initial gas stripping has occurred, there is a slower continuous stripping phase which is caused by the Kelvin-Helmholtz instability. Nulsen (1982) derived an equation for the mass loss rate of a spherical cloud, which Roediger & Brüggen (2007) adapted for a disk: M KH = 0.5πR 2 ρ ICM v gal . Because we use a constant density and velocity ICM wind, only the changing radius of the galaxy will change this rate of gas loss. -The radii of two galaxies with no cooling and in a static ICM at two different pressures. The radii fluctuate slightly at the epicyclic frequency. We believe this is because of the precipitous pressure decrease at the edge of the galaxy.
We measure the stripping radius in the simulation by finding the radius of gas with a density ρ > 10 −26 g cm −3 (or approximately a number density of n ∼ 10 −2 cm −3 ). To mitigate the effect of asymmetry on this measurement, we find the radius in 12 azimuthal slices around the disk and take the average. Because cooling causes both overdensities and underdensities, we measure this radius by finding the distance to the furthest cell with gas above our threshold density. As may be expected, this increases the radius we measure for the galaxies that include cooling, and has no effect on the radius of the galaxies without cooling. We compare this to the analytic stripping radius.
We find that -as in the galaxies without cooling -there does seem to be a radial epicyclic oscillation. However, in the higher ICM pressure cases (PHRCNW and PMR-CNW), low density and low pressure regions in the disk allow non-rotating ICM gas to enter and interact with the rotating denser clouds. These clouds therefore slow and lose angular momentum, leading the disk to settle towards the center. This drop in disk radius can be seen in the bottom panels of Figure 6. It is notable that in the lowest pressure ICM there is very little entrainment of the ICM gas, and the radius of the gas is dominated by epicyclic oscillations. We also find evidence of this process, and how it varies with ICM pressure, by following the total amount of gas in these galaxies. For example, PHRCNW gains 3.8 × 10 8 M ⊙ over the 950 Myr run, while PLRCNW gains 2.1 × 10 8 M ⊙ , because of the different amounts of gas entrained through cooling of the disk.
Because the radius of the gas disk decreases as the gas cools (see the bottom panels of Figure 6), we lower the gas radius of the galaxies without cooling to the radius that the radiatively cooled galaxies reach due to their entrainment of the ICM (as described in the previous paragraph). Therefore, although the radius of the galaxies with cooling begins at 25 kpc, the radius of the galaxies without cooling start at 21 kpc. Despite this different initial radius, we set the radial gas scale length, a gas , so that the galaxies with and without cooling have the same amount of mass. This is all done in order to ensure a fair comparison between the radiatively cooling and non-cooling runs. Myr, 250 Myr, 500 Myr, and 750 Myr after the wind has hit the galaxy. The projections are shaded to reflect surface density of the gas with a logarithmic stretch, and are 103.7 kpc across. There are a number of important differences in the morphology of the gas disks that are reflected in our quantitative comparisons. First, the gas in the galaxy with radiative cooling (PHRCW) forms small clouds with radii of about 100 pc near the center of the disk (within 7 kpc) that are of higher density than anywhere in the gas disk without cooling (PHNCW). Only these central clouds attain molecular densities, even in the galaxy that is never affected by an ICM wind. Further, at any radius of the disk in PHRCW, there is both higher and lower density gas than in the disk of PHNCW.

THE EFFECTS OF RAM PRESSURE ON A MULTIPHASE
This density structure has an important impact on the way that gas is stripped from the disk. First, it allows ICM gas to create holes in the disk, which leads to ablation of material from a wide range of radii, rather than just at the outer edge of the disk, as in the no-cooling run (PHNCW). As we will show in more detail below, this allows the gas to be stripped more quickly (see also Quilis et al 2000). In addition, the gas that is stripped tends to stream directly behind the disk, rather than flaring to the side, and therefore creates a narrower wake. We will examine the wake in more detail in a future paper.
The density fluctuations in the disk in the PHRCW galaxy create a ring of higher density gas from about 10 kpc to 16 kpc. This gas never collapses to the highest densities. In fact, in the PHRCNW galaxy, gas with ρ > 10 22 g cm −3 only persists for ∼100 Myr in the middle of the simulation. Although the vast majority of gas in this ring is always of a low enough density for direct stripping to occur, allowing the simulation to run until this ring had fully fragmented might result in a different stripping scenario. It would also allow gas more time to drift towards the center of the galactic potential, making it more difficult to strip. However, because little molecular gas is found outside the center of galaxies in both these simulations and in observations (e.g. Leroy et al. 2008), we don't expect the ring structure to make much difference in how much gas can be stripped.
Although in this paper we will not discuss the stripped gas in any detail, it is notable that the stripped gas in PHRCW has more density substructure than that of PH-NCW. Higher density gas leaves the disk in PHRCW than in PHNCW (although not gas with molecular densities), and stripped clumps survive for at least 100 kpc. We begin our quantitative analysis of the impact of cooling by comparing the amount of gas mass that is lost in the cases with and without cooling ( Figure 5). In general, the amount of gas mass lost by the galaxies that include cooling is about the same as the gas mass lost by the galaxies that do not include cooling. However, the timescale for stripping is much shorter for the radiatively cooled galaxies. These differences are illustrated in Table   Fig. 4.-A face-on and edge-on view of gas surface density simulation PHNCW, at 40, 250, 500, and 750 Myr after the wind has hit the galaxy. Note the continued smooth distribution of gas both in the disk and in the stripped gas throughout the simulation. Both the box size and color scale are the same as in Figure 3 4. Note that the last column in this table is the time from when the wind first hits the galaxy to when the galaxy has lost 85% of the total amount of gas lost.

Gas Mass Loss
In the highest ram pressure runs (PHNCW and PHRCW), after the wind has been hitting the galaxy for 800 Myr, both the galaxies with and without cooling have lost similar amounts of mass: 5.845 × 10 9 M ⊙ and 5.86 × 10 9 M ⊙ in the cooling and no cooling cases, respectively. We can compare these values to an analytic prediction for how much mass the galaxies should lose based on the amount of gas mass outside the analytic Gunn & Gott (1972) stripping radius measured immediately before the wind hits the galaxy. We find that the predictions are similar, although smaller, than the measured values: the galaxy with cooling should lose 5.46 × 10 9 M ⊙ and the galaxy without cooling should lose 5.22 × 10 9 M ⊙ . In the case with cooling, almost the entirety (86%) of this mass loss occurs within 150 Myr, while the same amount of mass loss takes about 650 Myr in the PHNCW case. Also, both of these galaxies have a fallback event (gas that is initially stripped from the disk region, although still gravitationally bound, falling back onto the disk) about 250 Myr after the wind hits the galaxy, although in PHNCW, the fallback is more drastic and seems to be followed by a few minor fallback events. Although the PHRCW case seems to stop losing gas despite the continued ICM wind, losing less than 5×10 8 M ⊙ in the final 300 Myr we follow the galaxy, the PHNCW case does not slow its gas loss by nearly as much, losing twice as much gas in the same amount of time. This difference is fit well by the difference in radius between the two cases and the dependence of Kelvin-Helmholtz stripping on radius squared.
In PMRCW, our run with the middle ram pressure value, after 900 Myr of the wind hitting the galaxy, it has lost about 1.3 × 10 9 M ⊙ of gas mass. In fact, it loses almost this entire amount of gas (95%) within 340 Myr of the wind first hitting the galaxy. On the other hand, in PMNCW, the galaxy has only lost 9.5×10 8 M ⊙ in the 900 Myr since the ICM wind has hit the galaxy. Unlike the highest ram pressure case, the analytic prediction for the amount of gas loss is larger than the measured amount: 1.7×10 9 M ⊙ for PMRCW and 1.6×10 9 M ⊙ for PMNCW. Again, the flatness of the PMNCW gas mass curve and the steeper grade of the PMRCW gas mass curve may be attributed to Kelvin-Helmholtz stripping effecting only the galaxy without cooling. Also, there does not seem to be any steep drop in the PMNCW galaxy's gas that is the hallmark of ram pressure pushing. Even more than in the highest ram pressure case, the galaxy with cooling seems to be only instantaneously stripped, without much longer-term Kelvin-Helmholtz stripping.
The lowest ram pressure runs (PLRCW and PLNCW) are very similar to the middle ram pressure runs, except that only 8.5 × 10 8 and 5.8 × 10 8 M ⊙ of gas are lost in the galaxies with and without cooling, respectively. Like the middle ram pressure case, the analytic prediction for the amount of gas loss is more than the amount measured, 1.0 × 10 9 (PLRCW) and 9.3 × 10 8 (PLNCW) M ⊙ . Even when the ram pressure is only 6.4 × 10 −13 dynes cm −2 , the galaxy with radiative cooling appears to lose gas through ram pressure pushing.

Runs
M lost M G&G r strip r G&G t 85% PHRCW 5.845 × 10 9 5.46 × 10 9 11 11.46 148 PHNCW 5.86 × 10 9 5.22 × 10 9 7 11.46 625 PMRCW 1.3 × 10 9 1.7 × 10 9 11 16.8 288 PMNCW 9.5 × 10 8 1.6 × 10 9 19.5 16.8 810 PLRCW 8.5 × 10 8 1 × 10 9 15 18.22 316 PLNCW 5.8 × 10 8 9.3 × 10 8 20? 18.22 754 Although similar amounts of gas mass are lost whether or not the galaxy is allowed to radiatively cool, there are a few trends. First, in the middle and low ram pressure cases, the amount of gas lost by the galaxy with radiative cooling is more than that lost by the galaxy without radiative cooling. Although this qualitatively agrees with the analytic predictions for the amount of gas mass loss indicating that the radiatively cooled galaxy will lose more gas, the difference in measured mass loss between the RC and NC cases is greater than the difference in the analytic prediction. Thus, it may be that for lower ram pressures, a multiphase ISM is easier to strip. Also, galaxies that radiatively cool lose their gas more quickly. There is slightly more mass in total outside the stripping radius in the radiatively cooling galaxies, and there is more mass at higher density. A plausible reason for the quick stripping of gas mass, then, is that a lumpy gas distribution is more quickly and easily stripped than gas with a smoother profile. Later continuous stripping seems to mainly affect the galaxies without cooling, perhaps because their larger cohesive gas disk is more affected by the Kelvin-Helmholtz instability.

Stripping Radius
In general, the galaxies that include radiative cooling are stripped to a smaller radius than predicted using the Gunn & Gott (1972) analytic equation, and end at a smaller radius than stripped galaxies without cooling and at a smaller radius than galaxies with radiative cooling that are not affected by an ICM wind ( Figure 6). Again, see Table 4 for the stripping radius values.
In the highest ram pressure cases, the analytic stripping radius is 11.5 kpc, very similar to the stripping radius of the galaxy in PHNCW. The oscillations in the radius correspond to fallback episodes in Figure 5. The final radius of the PHRCW case is about 7 kpc. The final radius of the PHRCNW galaxy is just above 11 kpc. This means that including cooling without a wind causes galactic gas to spiral towards the center of the disk, most likely because the entrainment of nonrotating ICM into the disk creates drag and angular momentum loss. The entrainment of gas into the disk means that the average gas density of the resulting disk would be higher than before, making it more difficult to strip. In fact, comparing the amount of gas at different densities in PHRCW and PHNCW about 15 Myr before the ICM wind hits the galaxy, we find that PHRCW has more high density gas outside 10 kpc and a higher fraction of the gas in the disk region is high density than in PHNCW. Therefore, the holes in the disk must allow for easier stripping down to smaller radii. The holes may also create more drag against surviving gas, allowing it to spiral towards the center of the disk more quickly.
In the middle ram pressure cases, the analytic stripping radius is 16.8 kpc. The final radius of the galaxy in PMNCW is between 19 and 20 kpc. However, the initial radius drop is to about 17 kpc, while a static galaxy without cooling, only undergoing epicyclic variations, has a smallest radius of about 19 kpc. The final average radius of the galaxy in PMRCNW is steady at 15 kpc, although it is very asymmetric. The galaxy in PMRCW is stripped down to a radius of 11 kpc. Again, a mixture of stripping and cooling results in the smallest radius.
In the lowest ram pressure cases, the analytic stripping radius is 18.2 kpc. It is difficult to say whether the galaxy in PLNCW is stripped to a smaller radius because of the large oscillation in radius over time. PLRCW is stripped to about 15 kpc. In this case more than in the previous two, it is clear that it must be the holes punched by the ICM wind in the gas disk that allow the gas to be stripped to a small galactic radius. Neither PLNCW nor PLRCNW have much change in the radius of the gas disk, so it must be a combination of the ICM wind and the clumpiness of the gas in the disk that result in the lowest gas radius.
By examining the stripping radius of gas in these galaxies, we find that allowing radiative cooling results in a smaller radius when hit by an ICM wind. Again, the most likely explanation is that the clumpiness of the disk gas allows for stripping to a smaller radius.

Densest Gas
As gas in the galaxy cools, it becomes more dense, eventually attaining densities greater than 1000 cm −3 , certainly in the range of molecular gas. Although with a resolution of 40 pc, we only marginally resolve giant molecular clouds, we can begin to answer whether ram pressure has any effect on the densest gas in a galaxy (and hence on star formation).
When comparing the three cases with radiative cooling and no ICM wind (bottom panels in Figure 7), it is clear that a higher pressure ICM results in a faster formation of high density gas and the formation of more high density gas. However, as we will see, including an ICM wind can either strip the surrounding lower density gas that would accrete onto clouds, or create more high density gas when the wind pressure is low.
We find that with the highest ram pressure of 6.4 × 10 −12 dynes cm −2 , the amount of higher density gas (above 10 −22 g cm −3 ) is less than if the galaxy were to radiatively cool without an ICM wind (see Figure  7). This is true even when (as in Figure 8), we shift the PHRCW line showing how much gas has cooled as a function of time by 200 Myr in order to account for the fact that the galaxy in PHRCW has been in the higher pressure ICM for almost 200 Myr less than the galaxy in PHRCNW. When the wind hits the galaxy, less dense gas throughout the disk is accelerated out of the galaxy, leaving less nearby lower density gas that would later either be accreted onto nearby clouds or cool and collapse into new clouds. However, gas with densities above 10 −22 g cm −3 is never seen outside of the disk (always within 1 kpc of the disk plane), so molecular clouds are not being stripped in their entirety. Examining the snapshots indicates that higher density clouds do leave the disk, leading us to hypothesize that dense clouds can be ablated until they are of just low enough column density to be stripped from the disk.
The scenario seems less clear in the medium ram pressure case (1 × 10 −12 dynes cm −2 ). Although the amount of higher density gas increases more slowly in PMRCW than in PMRCNW, both end their runs with similar amounts of high density gas (see Figure 8).
In the lowest ram pressure case (6.4 × 10 −13 dynes cm −2 ), the galaxy that is hit by the ICM wind has more of the highest density gas (above 10 −20 g cm −3 ) than the galaxy that evolves in the post-shock ICM. In this case, the added ram pressure from the wind appears to compress the gas, leading to more higher density gas.
In summary, molecular clouds do not get stripped, even if the ram pressure is 6.4 × 10 −12 dynes cm −2 . They do not grow within a high ram pressure ICM wind, however, because lower density gas is accelerated out of even the inner disk (inside the stripping radius of 7 kpc). However, with low ram pressure, 6.4×10 −13 dynes cm −2 , more high density gas is formed, likely because the ram pressure, while not enough to strip much gas from the disk, is enough to compress gas to a higher density.

Comparison to Observations
A galaxy is considered HI deficient if its deficiency value is 0.3 or greater. Deficiency is defined as DEF = < log X > T,D -log X obs (Giovanelli & Haynes 1985), where X obs is the observed HI mass of the cluster galaxy and < log X > T,D is averaged over a sample of field galaxies of morphological type T and optical diameter D. Given this definition, the highest ram pressure case is almost but not quite HI deficient. However, gas is truncated to less than half the original gas radius in the largest ram pressure case, and is down to less than 75% in the the lowest ram pressure case. If we were to roughly classify these galaxies using the scheme of Koopmann & Kenney (2004), assuming that HII regions will follow our gas distribution, all three of our runs result in truncated disks.
We can also compare the amount of high density gas we find in our galaxies to the amount of molecular gas observed using CO. Kenney & Young (1986, 1989 found that HI deficient galaxies in Virgo have undepleted CO, and Nakanishi et al. (2006) found that the molecular fraction in three galaxies near the center of Virgo is high. We find that in our high ram pressure run, the fraction of high density (i.e. molecular) gas actually remains similar to the case with no wind. On the other hand, in the low ram pressure case, the fraction of high density gas does increase as the ICM wind compresses the ISM. This case better matches the observations of high molecular gas fractions (Nakanishi et al. 2006). It is possible that galaxies will move through low ram pressure environments on their way to higher ram pressure, so the initial compression followed by strong stripping could result in HI depleted disks with unaffected, or increased, amounts of molecular gas. However, we do note that these results are somewhat sensitive to the resolution and details of the cooling (discussed in more detail below).
We do not include star formation in our simulations, but we do have an enhanced amount of high density gas in our lowest ram pressure simulation, while our highest ram pressure simulation has less high density gas in comparison to their respective no-wind comparison runs. If we assume that higher density gas corresponds to a higher rate of star formation, the lowest ram pressure case would result in an enhanced star formation rate, while the highest ram pressure case would result in a lower star formation rate. Placing this into a cluster environment, we would expect to find enhanced star formation at or outside the virial radius in the lower density ICM, and post-star-forming galaxies closer to the cluster center. In Coma, Poggianti et al (2004) found that post-starburst or post-star-forming galaxies were found at the edge of ICM substructure, which seems to follow our claim that starbursts will occur outside of the densest ICM, while high ram pressure will quench star formation.
We can also use morphological cues to determine how closely these simulations match observations. A perfect example of a galaxy showing both the strength and weakness of our simulations that include a multiphase medium is the Virgo spiral NGC 4402 (Crowl et al 2005). The galaxy has upturned edges in it's dust and gas distribution, similar to the earliest projection of our galaxy without cooling. However, similar to our galaxy with radiative cooling, there is evidence that two dense clouds have survived longer than the lower density ISM and are being ablated by the wind.

The Evolving ISM
A molecular cloud may survive between a few and a few tens of Myr (Blitz & Shu 1980;Larson 2003;Hartmann 2003). Because we do not include star formation, our dense clouds of gas do not evolve into less dense pockets of gas that could be stripped. In the highest ram pressure case, holes are punched around the disk (as close as 6 kpc to the disk center) between 6 -11 Myr after the ICM wind first hits the galaxy, which indicates that as long as gas stays low density for 10 Myr it can be stripped by the wind. If our dense regions of gas could turn into very low density regions through star formation (ρ ≤ 10 −23 g cm −3 ), our galaxies that include radiative cooling could be entirely stripped of gas in the highest ram pressure case. In the lowest ram pressure case, no gas appears to be stripped from the central region of the gas disk.  (190 Myr,274 Myr,and 271 Myr,respectively). High ram pressure results in less high density gas, while low ram pressure increases the amount of high density gas.

Magnetic Fields and Cosmic Rays
Our simulation does not include magnetic fields or cosmic rays. In general in the Milky Way, the magnetic field is measured to be about 2-3 microGauss (Men et al 2008). A field of this strength would not have enough magnetic pressure (2 − 5 × 10 −13 dynes cm −2 ) to significantly lessen the direct impact of ram pressure on the galactic gas. However, there are measurements of the magnetic field in molecular clouds (Crutcher 1991 and references therein) of up to 10 3 microGauss, which can certainly withstand ram pressures of order 10 −12 dynes cm −2 . These magnetic fields could help protect molecular clouds from being ablated by the ICM wind, so could help them survive. Also, a magnetic field could attach a molecular cloud to less dense gas, either anchoring the surrounding gas or helping to drag out the molecular cloud. Using this simulation we cannot conjecture on the likely impact of magnetic fields on the total amount of gas stripped.
We also do not include cosmic rays in our simulation, which would add energy to the galactic gas and increase the vertical distribution (e.g. Joung & Mac Low 2006). It is also possible that the increase in gas pressure would lower the amount of entrainment into the disk from the lower pressure ICM, which would in turn lessen the number of very low density regions and make the disk more cohesive. This would make stripping more difficult because fewer holes could be punched through the disk by the ICM wind. However, the energy input from cosmic rays would slow high density gas formation, possibly allowing for more extended stripping if the gas is not locked in very dense clouds. As with magnetic fields, the addition of cosmic rays into our simulation does not have an easily predicted effect.

Resolution
As has been discussed in detail in previous works (e.g. Tonnesen et al. 2007 andBrüggen 2006), resolution does not have a large effect on our runs without radiative cooling. The runs without radiative cooling that we discuss in this paper have a cell size of 80 pc, and we ran a comparison run with a cell size of 40 pc to check this. We found very similar results in the amount of gas and the radius of the disk. In addition, our results are similar to those in Roediger & Brggen (2006), where even lower resolution was used.
However, including radiative cooling leads to results that are very dependent on resolution. The runs discussed in this paper have a cell size of 40 pc, which leads to gas collapsing to molecular densities. For example, the surface density in the center of the disk immediately Fig. 9.-The gas mass and gas radius over time for PHRCW, PHRCW low resolution, and PHNCW. This figure shows how the low resolution run with radiative cooling is a cross between PHRCW and PHNCW, because the gas density is smoother than PHRCW, but less smooth than PHNCW. See §5.4 for discussion.
before the wind hits the galaxy is 200 M ⊙ kpc −2 , with values above 10 M ⊙ kpc −2 for much of the inner disk (this surface density corresponds to the transition between HI and molecular gas observationally found in Bigiel et al 2008).
To examine the impact of resolution, we have run a comparison simulation with the same parameters as PHRCW, but with a cell size of 80 pc. We find that gas density does increase, but the highest density gas never forms, and there is very little gas in the middle density range. Because the lack of cloud formation results in a smoother gas distribution in the galaxy disk, the ICM wind does not have the same effect on this lower resolution disk as it does in the higher resolution cases. Figure  9 shows that although the low resolution disk with cooling still gets stripped slightly more quickly than the no cooling case, it does not get stripped as quickly as the higher resolution case. This occurs because the corresponding low density regions in the disk do not form, so the ICM wind does not punch holes throughout the disk. The low resolution run also has much less fallback than PHNCW, but because very dense clouds are not formed (see Figure 10), continuous stripping occurs and gas loss does not cease after a few hundred Myr.
We also ran a pair of comparison runs between a galaxy with a cell size of 40 pc and one with a cell size of 20 pc. Because of the computing time needed to run the high resolution simulation, we only had output for 160 Myr of evolution. However, within this time, the maximum surface density of the higher resolution run is more than double that of the regular resolution run. Further, after 160 Myr the amount of gas with densities above 10 −20 g cm −3 is already more than 25% of the total galaxy gas mass and growing at a rate of 40 M ⊙ yr −1 , while gas with those densities has not yet formed in the regular resolution run (at 160 Myr, ρ max ≥ 10 −22 g cm −3 ). In fact, in the highest resolution run, 5 × 10 8 M ⊙ of gas with densities above 10 −19 g cm −3 (n > 10 5 cm −3 ) has already formed by 160 Myr into the simulation, with no sign of slowing down much higher gas density formation. A visual inspection shows that the cloud radius is set by the minimum cell size in the simulation. We conclude that there are physical processes that we do not include in the simulation, such as star formation, magnetic fields, cosmic rays, and feedback from supernovae, that would all result in lower density gas than is found in our simulations, either by adding energy to the ISM or locking the ISM into place and hindering collapse. In this case, higher resolution is not necessarily better, because our omissions may result in unphysical gas densities. In the regular resolution run, the larger cells balance, to a certain extent, the smaller-scale physics we omit.

Radiative Cooling Floor
As mentioned in §2.1, we include radiative cooling using the Sarazin & White (1987) cooling curve. We allow cooling to 8,000 K, as this results in gas with neutral hydrogen temperatures without overcooling a large fraction of our gas. With this cutoff, we find that we form clouds with densities and sizes typical of molecular clouds, al- though we make no attempt to reproduce the internal structure of the clouds. For a more detailed discussion of the density distribution in the ISM for similar simulations, but without a wind, see Tasker & Bryan (2006.
To examine the impact of our adopted cooling floor, we perform a run with no wind and cooling to 300 K which shows that after 800 Myr the density distribution of disk gas was nearly identical to the case with cooling to 8000 K, although the case with cooling to 300 K reached the distribution more quickly and leveled out. The highest density gas is less than a factor of three more dense in the 300 K case than in the 8000 K case.
To better compare how an ICM wind would affect a galaxy with a lower cooling floor, we ran a simulation which is identical to PHRCW, but with radiative cooling to 300 K. We find a few interesting results, shown in Figure 11. First, more gas is lost than in either PHRCW or PHNCW. Gas is initially lost more slowly than in the other two cases, but does not slow down with fallback as in PHNCW. Also, the final radius of the gas disk is similar to PHRCW, but it reaches this radius more slowly. Finally, the amount of gas in each of the three density ranges is similar to that in PHRCW, but the two lower density regimes have a little less gas. The highest density gas is the same in both runs throughout the simulations. We also show ( Figure 12) an image of each of the two cooling cases 150 Myr after the wind has hit the galaxy.
All of the differences in these two cases are explained by the fact that the galaxy with the lower radiative cooling floor fragments more than our standard cooled galaxy. Because more of the gas is in dense knots, the initial onslaught of the wind does not strip gas as quickly. Further, because gas that began in the outer regions of the disk survives, we see epicyclic variation in the radius, similar to our other simulations, resulting in an increase in the radius. This is unlike PHRCW, in which all the gas that began in the outskirts is quickly stripped. Because the clouds in the 300 K cooling case move to large radii, there is more time for them to be ablated as they drift back towards the center of the galaxy. This explains why gas takes longer to be stripped, why more gas is stripped, and why the amount of mid-density gas decreases in this galaxy. In the inner region of the disk, the larger amount of fragmentation means that more gas is stripped at early times from the central regions of the disk, because the gas disk is not smooth and cannot deflect the ICM wind. In fact, we do not see gas fallback, as in the other high ram pressure cases, presumably because the wind pushes gas at all radii. The longer period of gas stripping results both from ablation of the dense clouds in the center of the disk that are more exposed to the wind because of the larger amount of fragmentation, and because the dense clouds at large radii can be ablated more easily than their counterparts in the 8000 K case, which are closer to the galactic center.
It is unclear what level of fragmentation is more realistic in observed galaxies. Because of the lower Jeans length in the lower temperature floor runs, we resolve the resulting fragmentation less well and so prefer the runs with a higher minimum temperature. In addition, it is not clear if the lower minimum temperature runs are physically more realistic as we do not include some effects such as small-scale turbulence, cosmic rays and magnetic fields, all of which may provide a source of effective pressure in low temperature regions. While our simulations represent a sizable step forward in realism, it is clear that more work is required.

CONCLUSIONS
We have run a set of detailed galaxy simulations including radiative cooling to understand how a multiphase ISM interacts with an oncoming ICM wind. We find: 1. The total amount of gas mass loss is similar when comparing galaxies with and without radiative cooling. However, the time over which the gas is stripped is shorter in the galaxies with radiative cooling. The continued stripping of galaxies without radiative cooling could be Kelvin-Helmholtz stripping acting on the larger disk.
2. The stripping radius is significantly smaller for the galaxy simulations with radiative cooling. This is even true when accounting for the angular momentum loss in the ISM gas due to the interaction between the ISM and ICM gas (see section 3).
3. The morphology of the stripped disks are considerably different, with the cooling case producing a multiphase medium which results in many holes in the disk, while the case without cooling produces a coherent disk. These holes allow ICM gas to stream through the galaxy disk, affecting the evolution of the dense clouds which remain.
4. The amount of high density gas is affected differently depending on the amount of ram pressure. When the ram pressure is strong, there is less gas at all densities, including the highest density. This is because the low density gas is stripped, even from the inner regions of the galaxy, so there is less gas to feed molecular clouds. With a weaker ram pressure, there is more of the highest density gas, although there is also less of the lower density gas. A weaker wind strips less of the lower density gas, and increases the pressure on disk gas to create more high density gas.
This leads us to a picture in which radiative cooling re- Fig. 11.-Comparison between PHRCW and the case with cooling to 300 K. In both cases, the ICM wind hits the galaxies at 190 Myr. Note how extended the stripping is in the 300 K cooling case, and note that more gas is stripped. These facts are linked with the larger gas radius for much of the simulation. See §5.5 for discussion. sults in both overdensities and underdensities at all radii in the galactic gas disk. The underdensities are clearly stripped quickly throughout the disk, and the overdense regions that then have the wind streaming around them through the holes in the disk are stripped more rapidly than without a multiphase medium. However, if a dense cloud is able to survive this initial onslaught, it loses angular momentum through its interaction with the nonrotating ICM wind and therefore settles toward a lower radius. This is how the runs with cooling and a wind (the RCW cases) are stripped quickly and to a small radius while not losing much more gas than the galaxies without radiative cooling. Also, the galaxies with radiative cooling seem to stop losing gas while those without cooling continue to lose gas through the entire simulation. This may be because the Kelvin-Helmholtz instability is inefficient on small dense clouds but can affect the large cohesive disk in the no-cooling runs. Allowing cooling to a lower temperature supports our picture that dense clouds cannot be directly stripped but can be ablated over time. Although these simulations are arguably the most realistic stripping simulations performed to date, we have discussed in detail the limitations involved in our work in §5.
We acknowledge support from NSF grants AST-05-07161, AST-05-47823, and AST-06-06959, as well as computational resources from the National Center for Supercomputing Applications. We thank Jacqueline van Gorkom, Jeff Kenney, Hugh Crowl, Airee Chung, Tomer Tal, and Anne Abramson for useful discussions, as well as Elizabeth Tasker for invaluable help setting up the initial conditions. We would also like to thank our referee for editing and advice that greatly improved our paper. Fig. 12.-A tilted view of the high ram pressure cases with radiative cooling to 8000 K and 300 K, respectively, 150 Myr after the wind has hit the galaxy. Note that the galaxy with cooling to 300 K has formed more high density clouds in the disk outskirts that are not stripped. See §5.5 for discussion.

PHRCNW
This galaxy begins with a gas radius of 25 kpc (as do all of the galaxies in simulations that include radiative cooling). Throughout this run the surrounding ICM has the density and temperature of the post-shock ICM when there is a wind in the high ram pressure cases. We chose the hotter ICM because of it's longer cooling time. Initially the radius seems to fluctuate with the epicyclic period (∼ 500 Myr), however, cooling creates both over densities and underdensities in the disk. This allows some of the non-rotating ICM to fall into holes in the disk and to slow the disk's rotation, leading to a drop in the disk radius through the loss of angular momentum. This gas infall is also seen by considering the total amount of gas in the disk region: the galaxy gains 3.8 × 10 8 M ⊙ over the 950 Myr run. High density gas forms about 100 Myr into the galaxy's evolution, with the highest density gas, above 10 −20 g cm −3 , forming closer to 150 Myr. Throughout the run the galaxy gains gas with densities above 10 −22 g cm −3 , loses less dense gas as it turns into more dense gas. Although we would need to extend this run to be sure the radius would stay constant, the radius settles between 11 and 12 kpc in the final 200 Myr of the run.

PHRCW
As mentioned above, because this galaxy includes radiative cooling, it begins with a radius of 25 kpc. The surrounding density and temperature are the same as run PHNCW, shown in Table 3. After about 190 Myr, the same ICM wind hits this galaxy as hits the galaxy in PHNCW. Very little gas is lost for a few tens of Myrs, at which point 5 × 10 9 M ⊙ is lost from the galaxy in under 150 Myr. 3 × 10 9 M ⊙ is immediately gravitationally unbound from the galaxy, and almost all of the remaining 2 × 10 9 M ⊙ leaves the halo within 400 Myr. A small amount of this gas falls back onto the galactic plane. After the gas immediately removed from the galactic disk by ram pressure pushing is no longer gravitationally bound, about 550 Myr after the wind first hits the galaxy, gas loss precedes at about 1 M ⊙ yr −1 for the final 250 Myr we continue the run. The total gas mass loss is about 5.8 × 10 9 M ⊙ . As in the PHRCNW run, the radius of the disk falls before the ICM wind hits the galaxy. Once the wind hits the galaxy the gas is stripped to a radius below 9 kpc within 250 Myr. The radius after 800 Myr of ram pressure stripping is at or below 7 kpc. As in PHRCNW, the galaxy loses lower density gas (below 10 −22 g cm −3 ) throughout the run. Although the galaxy does increase all it's higher density gas, gas with densities between 10 −20 and 10 −22 g cm −3 remains constant after the ICM wind has been hitting the galaxy for 350 Myr. Also, the highest density gas does not appear until 350 Myr into the simulation, about 150 Myr after the ICM wind hits the galaxy.

PMNCW
This galaxy begins with a radius of 21 kpc. The surrounding ICM density and temperature is listed in Table 3. After the galaxy has evolved for about 80 Myr, an ICM wind hits the galaxy with a ram pressure of 1 × 10 −12 dynes cm −2 , and density and temperature listed in Table 3. The total amount of gas loss in the more than 900 Myr in which the ICM wind hits the galaxy is less than 9.5 × 10 8 M ⊙ . The mass loss rate for the first 100 Myr is less than 0.5 M ⊙ yr −1 , and increases until in the final 100 Myr it is a little over 2 M ⊙ yr −1 . The radius fluctuates, but seems to settle at about 19 kpc. In this case, all the gas that is stripped is less dense than 10 −24 g cm −3 .

PMRCNW
This galaxy begins with a radius of 21 kpc. This surrounding ICM density and temperature are the middle postshock values listed in Table 3. Again the radius seems to vary due to a combination of epicyclic variation and mixing from non-rotating ICM gas. Less ICM gas is accreted in this case than in PHRCNW, 2.1 × 10 8 M ⊙ , which can be expected since the ICM pressure is lower. As in the higher ICM pressure case, PHRCNW, throughout the run the galaxy gains gas with densities above 10 −22 g cm −3 , and loses gas less dense than 10 −22 g cm −3 . However, there is less gas with densities greater than 10 −20 g cm −3 in this run than in PHRCNW, and it does not form until about 250 Myr. The radius remains at about 15.5 kpc for the final 400 Myr of this run.

PMRCW
This galaxy begins with a radius of 25 kpc. The surrounding ICM density and temperature values are the same as in PMNCW. After the galaxy has evolved for about 275 Myr, the ICM wind hits the galaxy with a ram pressure of 1 × 10 −12 dynes cm −2 . For 100 Myr little gas is lost, and then in less than 350 Myr after the wind first hits the galaxy the majority of the gas is stripped, 1.23 × 10 9 M ⊙ of a total loss of 1.3 × 10 9 M ⊙ . The radius falls to 21 kpc before the ICM wind hits the galaxy, and then begins to fall again about 230 Myr after the wind first hits the galaxy. The final radius, just below 11 kpc, is reached in slightly less than 900 Myr after the wind hits the galaxy and is constant for the final 200 Myr of the simulation. The galaxy loses lower density gas (below 10 −22 g cm −3 ) throughout the run. Although the galaxy does increase it's higher density gas, gas with densities between 10 −20 and 10 −22 g cm −3 remains constant after the ICM wind has been hitting the galaxy for 700 Myr. The highest density gas, above 10 −20 g cm −3 , is first formed 450 Myr into the simulation, or 175 Myr after the higher pressure ICM hits the galaxy.

PLNCW
This galaxy begins with a radius of 21 kpc. The surrounding ICM density and temperature are listed in the bottom pre-shock rows of Table 3. After the galaxy has evolved for about 77 Myr, the ICM wind hits the galaxy with the post-shock values in Table 3 and a ram pressure of 6.4 × 10 −13 dynes cm −2 . The total gas loss is 5.8 × 10 8 M ⊙ . The