update

.gitignore

@@ -1,2 +1,3 @@
 *.g2
 *.log
+*.*Store

high/agora/config_agora_sigma00.dice

@@ -3,7 +3,7 @@
 # to re-use the previous computation,
 # and thus reduces the global execution time
 #                                   x[kpc]  y[kpc]  z[kpc]  vx[km/s]  vy[km/s]  vz[km/s]  spin[deg]  incl[deg]
-Galaxy         params/agora_sigma30.params  0.      0.      0.      0.        0.        0.        0.        0.
+Galaxy         params/agora_sigma00.params  0.      0.      0.      0.        0.        0.        0.        0.
 # Number of OpenMP threads
 Nthreads                16
 # Compute the mean inter-particle distance
@@ -15,7 +15,7 @@ MaxCompNumber           64
 OutputRz                0
 # Name of the final IC file
 # The file extension wil be added by DICE
-Filename                out/agora_sigma30
+Filename                out/agora_sigma00
 # File format of the IC file [Gadget1,Gadget2]
 ICformat      Gadget2
 # Cosmological parameters for Virial quantities
@@ -23,11 +23,11 @@ ICformat      Gadget2
 # Redshift
 Redshift      1.0
 # Hubble constant
-H0            71.0            
+H0            73.0            
 # Dark energy density parameter
-OmegaL            0.70
+OmegaL            0.6911
 # Matter density parameter
-OmegaM            0.30
+OmegaM            0.3089
 # Spatial curvature density paramter
 OmegaK            0.00
 # Normalize the sum of mass fractions to one

high/agora/config_agora_sigma20.dice

@@ -3,7 +3,7 @@
 # to re-use the previous computation,
 # and thus reduces the global execution time
 #                                   x[kpc]  y[kpc]  z[kpc]  vx[km/s]  vy[km/s]  vz[km/s]  spin[deg]  incl[deg]
-Galaxy         params/agora_sigma10.params  0.      0.      0.      0.        0.        0.        0.        0.
+Galaxy         params/agora_sigma20.params  0.      0.      0.      0.        0.        0.        0.        0.
 # Number of OpenMP threads
 Nthreads                16
 # Compute the mean inter-particle distance
@@ -15,7 +15,7 @@ MaxCompNumber           64
 OutputRz                0
 # Name of the final IC file
 # The file extension wil be added by DICE
-Filename                out/agora_sigma10
+Filename                out/agora_sigma20
 # File format of the IC file [Gadget1,Gadget2]
 ICformat      Gadget2
 # Cosmological parameters for Virial quantities

high/agora/nml/base.nml

@@ -11,7 +11,7 @@ rt=.false.
 /
 
 &OUTPUT_PARAMS
-delta_tout=0.25
+delta_tout=1
 tend=50.0
 /
 
@@ -22,7 +22,7 @@ tendmov=50.
 movie=.true.
 nw_frame=1080
 nh_frame=1080
-levelmax_frame=14
+levelmax_frame=16
 movie_vars_txt='stars','dens','temp'
 proj_axis='zx'
 xcentre_frame=160,0.0,0.0,0.0 160,0.0,0.0,0.0
@@ -36,7 +36,7 @@ deltaz_frame=30,0.0 30,0.0
 
 &INIT_PARAMS
 filetype='dice'
-initfile(1)='/bulk1/teyssier/hlasco/dice_ics/med/agora/out'
+initfile(1)='/users/hlascomb/dice_ics/high/agora/out'
 /
 
 &UNITS_PARAMS
@@ -46,9 +46,9 @@ units_length  = 0.308567758128200E+22 ! kpc
 /
 
 &AMR_PARAMS
-levelmin=7
+levelmin=8
-levelmax=14
+levelmax=16
-ngridmax=3000000
+ngridmax=1000000
 npartmax=5000000
 boxlen=320.0
 /
@@ -83,8 +83,8 @@ n_star=0.01                         ! SF density threshold. [0.01, 10]
 eps_star=1.0                        ! SF efficency. Should be 1 with sf_virial=.true.       
 T2_star=10.                         ! ISM polytropic temperature
 g_star=1.0                          ! ISM polytropic index
-m_star=10.0                         ! Star particle mass 
+m_star=1.0                          ! Star particle mass 
-sf_log_properties=.false.           ! do not set this to .true.
+sf_log_properties=.true.
 /
 
 &FEEDBACK_PARAMS
@@ -107,14 +107,14 @@ riemann='hllc'
 /
 
 &REFINE_PARAMS
-mass_sph=1e-6
+mass_sph=1e-7
 m_refine=10*8.
 interpol_type=0
 interpol_var=1
 /
 
 &DICE_PARAMS
-ic_file='agora_sigma20.g2'          ! Name of the initial conditions file
+ic_file='agora_sigma20.g2'
 ic_nfile=1                    ! If greater than one, look for files with name matching ic_file//'.n'
 ic_ifout=1                    ! Change ramses output index for restarts
 ic_format='Gadget2'           ! Format of the initial conditions. 'Gadget1' or 'Gadget2'

high/agora/params/agora_sigma00.params

@@ -49,11 +49,11 @@ hydro_eq_niter          3
 # Fraction of the virial mass in the component 1
 mass_frac1              0.9578424465115449
 # Number of particles for the component 1
-npart1                  100000
+npart1                  1000000
 # Target mass of individual particles [Msol]
 #part_mass1              2e6
 # Number of particles for the potential computation
-npart_pot1              200000
+npart_pot1              2000000
 # Target mass of individual particles for potential computation [Msol]
 #part_mass_pot1          1e6
 # Component 1 density model
@@ -84,7 +84,7 @@ rcore1                  0.
 # the concentration parameter of an NFW halo with a mass M200*mass_frac
 concentration1          10
 # Gravitational softening for the poisson solver [kpc]
-softening1              0.6
+softening1              0.3
 # Component 1 radial density cut [kpc] 
 cut1                    0.0
 # Component 1 thickness parameters
@@ -132,8 +132,8 @@ accept_max1             0.99
 # Component 2: Thin stellar disk
 ################################
 mass_frac2              0.03066166043860979
-npart2                  100000
+npart2                  1000000
-npart_pot2              200000
+npart_pot2              2000000
 model2                  3
 # If the value is zero and the particle type is not 1 the size is determined using the spin conservation
 # (Fitting formula from Mo, Mao & White 1998)
@@ -183,9 +183,9 @@ metal_gradient3         0
 # Temperature of the gas particles [K]
 t_init3                 1e4
 # Turbulent velocity dispersion [km/s]
-turb_sigma3             10.0
+turb_sigma3             0.0
 # Turbulence injection scale [kpc]
-turb_scale_inj3         0.34
+turb_scale_inj3         1.0
 # Turbulence dissipation scale [kpc]
 turb_scale_diss3        0.01
 # Seed for the turbulent gaussian field
@@ -203,8 +203,8 @@ gamma_poly3             1.0
 # Component 4: Stellar bulge
 ############################
 mass_frac4              0.0038322616318995422
-npart4                  12500
+npart4                  125000
-npart_pot4              25000
+npart_pot4              250000
 model4                  4
 scale_length4           0.343218
 flatz4                  0.80

high/agora/params/agora_sigma20.params

@@ -49,11 +49,11 @@ hydro_eq_niter          3
 # Fraction of the virial mass in the component 1
 mass_frac1              0.9578424465115449
 # Number of particles for the component 1
-npart1                  100000
+npart1                  1000000
 # Target mass of individual particles [Msol]
 #part_mass1              2e6
 # Number of particles for the potential computation
-npart_pot1              200000
+npart_pot1              2000000
 # Target mass of individual particles for potential computation [Msol]
 #part_mass_pot1          1e6
 # Component 1 density model
@@ -84,7 +84,7 @@ rcore1                  0.
 # the concentration parameter of an NFW halo with a mass M200*mass_frac
 concentration1          10
 # Gravitational softening for the poisson solver [kpc]
-softening1              0.6
+softening1              0.3
 # Component 1 radial density cut [kpc] 
 cut1                    0.0
 # Component 1 thickness parameters
@@ -132,8 +132,8 @@ accept_max1             0.99
 # Component 2: Thin stellar disk
 ################################
 mass_frac2              0.03066166043860979
-npart2                  100000
+npart2                  1000000
-npart_pot2              200000
+npart_pot2              2000000
 model2                  3
 # If the value is zero and the particle type is not 1 the size is determined using the spin conservation
 # (Fitting formula from Mo, Mao & White 1998)
@@ -183,11 +183,11 @@ metal_gradient3         0
 # Temperature of the gas particles [K]
 t_init3                 1e4
 # Turbulent velocity dispersion [km/s]
-turb_sigma3             30.0
+turb_sigma3             20.0
 # Turbulence injection scale [kpc]
 turb_scale_inj3         0.34
 # Turbulence dissipation scale [kpc]
-turb_scale_diss3        0.01
+turb_scale_diss3        0.005
 # Seed for the turbulent gaussian field
 turb_seed3              1234
 # Compute hydrostatic equilibrium
@@ -203,8 +203,8 @@ gamma_poly3             1.0
 # Component 4: Stellar bulge
 ############################
 mass_frac4              0.0038322616318995422
-npart4                  12500
+npart4                  125000
-npart_pot4              25000
+npart_pot4              250000
 model4                  4
 scale_length4           0.343218
 flatz4                  0.80

low/agora/nml/agora_sigma00.nml

@@ -1,145 +0,0 @@
-&RUN_PARAMS 
-cosmo=.false.
-pic=.true.
-poisson=.true.
-hydro=.true.
-nrestart=0
-nremap=25
-nsubcycle=1,1,2
-ncontrol=1
-rt=.false.
-/
-
-&OUTPUT_PARAMS
-delta_tout=0.5
-tend=20.0
-/
-
-&MOVIE_PARAMS
-imovout=2000
-imov=0
-tendmov=50.
-movie=.false.
-nw_frame=1080
-nh_frame=1080
-levelmax_frame=12
-movie_vars_txt='stars','dens','temp'
-proj_axis='zx'
-xcentre_frame=120,0.0,0.0,0.0 120,0.0,0.0,0.0
-ycentre_frame=120,0.0,0.0,0.0 120,0.0,0.0,0.0
-zcentre_frame=120,0.0,0.0,0.0 120,0.0,0.0,0.0
-deltax_frame=30,0.0 30,0.0
-deltay_frame=30,0.0 30,0.0
-deltaz_frame=30,0.0 30,0.0
-! zoom_only_frame=.true.
-/
-
-&INIT_PARAMS
-filetype='dice'
-initfile(1)='/bulk1/teyssier/hlasco/dice_ics/low/agora/out'
-/
-
-&UNITS_PARAMS
-units_density = 0.677025430198932E-22 ! 1e9 Msol/kpc^3
-units_time    = 0.470430312423675E+15 ! G=1
-units_length  = 0.308567758128200E+22 ! kpc
-/
-
-&AMR_PARAMS
-levelmin=7
-levelmax=12
-ngridmax=3000000
-npartmax=5000000
-boxlen=320.0
-/
-
-&POISSON_PARAMS 
-epsilon=1e-4
-/
-
-&BOUNDARY_PARAMS
-nboundary=6
-bound_type= 2, 2,  2,  2,  2,  2
-ibound_min=-1, 1, -1, -1, -1, -1
-ibound_max=-1, 1,  1,  1,  1,  1
-jbound_min= 0, 0, -1,  1, -1, -1
-jbound_max= 0, 0, -1,  1,  1,  1
-kbound_min= 0, 0,  0,  0, -1,  1
-kbound_max= 0, 0,  0,  0, -1,  1
-no_inflow=.true.
-/
-
-&COOLING_PARAMS
-cooling=.true.
-metal=.true.
-z_ave=1.0                           ! To explore explore [0.001, 1]
-self_shielding=.true.
-/
-
-&SF_PARAMS
-sf_virial=.true.                    ! Enable turbulent SF
-sf_model=5                          ! See mom2/star_formation.f90  
-n_star=0.01                         ! SF density threshold. [0.01, 10]
-eps_star=1.0                        ! SF efficency. Should be 1 with sf_virial=.true.       
-T2_star=10.                         ! ISM polytropic temperature
-g_star=1.0                          ! ISM polytropic index
-m_star=10.0                         ! Star particle mass 
-sf_log_properties=.false.           ! do not set this to .true.
-/
-
-&FEEDBACK_PARAMS
-momentum_feedback=1                 ! SN feedback model (from patch/mom2)
-delayed_cooling=.false.             ! Cooling starts from t=0
-yield=0.177703                      ! SN metal yield
-eta_sn=0.162637                     ! SN efficiency
-t_sne=3.0                           ! time for a star to turn into SN
-f_w=0.
-/
-
-&HYDRO_PARAMS
-gamma=1.6666
-courant_factor=0.5
-slope_type=1
-pressure_fix=.true.                 ! Scheme to prevent negative T
-beta_fix=0.5                        ! Controls threshold at which energy is 
-                                    ! truncated
-riemann='hllc'
-/
-
-&REFINE_PARAMS
-mass_sph=1e-5
-m_refine=10*8.
-interpol_type=0
-interpol_var=1
-/
-
-&DICE_PARAMS
-ic_file='agora_sigma00.g2'          ! Name of the initial conditions file
-ic_nfile=1                    ! If greater than one, look for files with name matching ic_file//'.n'
-ic_ifout=1                    ! Change ramses output index for restarts
-ic_format='Gadget2'           ! Format of the initial conditions. 'Gadget1' or 'Gadget2'
-ic_center=0.0,0.0,0.0         ! Shift center parameter. ICs are automatically shifted with boxlen/2
-ic_scale_pos=1.0              ! Scaling factor for the position vector 
-ic_scale_vel=1.0              ! Scaling factor for the velocity vector
-ic_scale_mass=1.0             ! Scaling factor for the mass
-ic_scale_u=1.0                ! Scaling factor for the internal energy
-ic_scale_age=1.0              ! Scaling factor for the particles age
-ic_scale_metal=1.0            ! Scaling factor for the metallicity
-ic_head_name='HEAD'           ! Name of the Header datablock (Gadget2 format only)
-ic_pos_name='POS '            ! Name of the position vector datablock (Gadget2 format only)
-ic_vel_name='VEL '            ! Name of the velocity vector datablock (Gadget2 format only)
-ic_mass_name='MASS'           ! Name of the mass datablock (Gadget2 format only)
-ic_id_name='ID  '             ! Name of the particle identifier datablock (Gadget2 format only)
-ic_u_name='U   '              ! Name of the internal energy datablock (Gadget2 format only)
-ic_metal_name='Z   '          ! Name of the metallicity datablock (Gadget2 format only)
-ic_age_name='AGE '            ! Name of the particle age datablock (Gadget2 format only)
-IG_rho=1.0D-6                 ! Intergalactic gas density
-IG_T2=1.0D7                   ! Intergalactic gas temperature
-IG_metal=0.0                  ! Intergalactic gas metallicity
-amr_struct=.false.            ! Reproduce the AMR structure of the Gadget2 file resulting from a ramses to gadget conversion
-gadget_scale_l=3.085677581282D21 ! Gadget file length unit
-gadget_scale_v=1.0D5             ! Gadget file velocity unit
-gadget_scale_m=1.9891D43         ! Gadget file mass unit
-gadget_scale_t=3.15360e+13       ! Gadget file time unit
-ic_skip_type        = -1      ! Skip specific particle type
-/

low/agora/nml/agora_sigma10.nml

@@ -1,145 +0,0 @@
-&RUN_PARAMS 
-cosmo=.false.
-pic=.true.
-poisson=.true.
-hydro=.true.
-nrestart=0
-nremap=25
-nsubcycle=1,1,2
-ncontrol=1
-rt=.false.
-/
-
-&OUTPUT_PARAMS
-delta_tout=0.5
-tend=50.0
-/
-
-&MOVIE_PARAMS
-imovout=2000
-imov=0
-tendmov=50.
-movie=.false.
-nw_frame=1080
-nh_frame=1080
-levelmax_frame=12
-movie_vars_txt='stars','dens'
-proj_axis='zx'
-xcentre_frame=160,0.0,0.0,0.0 160,0.0,0.0,0.0
-ycentre_frame=160,0.0,0.0,0.0 160,0.0,0.0,0.0
-zcentre_frame=160,0.0,0.0,0.0 160,0.0,0.0,0.0
-deltax_frame=30,0.0 30,0.0
-deltay_frame=30,0.0 30,0.0
-deltaz_frame=30,0.0 30,0.0
-! zoom_only_frame=.true.
-/
-
-&INIT_PARAMS
-filetype='dice'
-initfile(1)='/bulk1/teyssier/hlasco/dice_ics/low/agora/out'
-/
-
-&UNITS_PARAMS
-units_density = 0.677025430198932E-22 ! 1e9 Msol/kpc^3
-units_time    = 0.470430312423675E+15 ! G=1
-units_length  = 0.308567758128200E+22 ! kpc
-/
-
-&AMR_PARAMS
-levelmin=7
-levelmax=12
-ngridmax=3000000
-npartmax=5000000
-boxlen=320.0
-/
-
-&POISSON_PARAMS 
-epsilon=1e-4
-/
-
-&BOUNDARY_PARAMS
-nboundary=6
-bound_type= 2, 2,  2,  2,  2,  2
-ibound_min=-1, 1, -1, -1, -1, -1
-ibound_max=-1, 1,  1,  1,  1,  1
-jbound_min= 0, 0, -1,  1, -1, -1
-jbound_max= 0, 0, -1,  1,  1,  1
-kbound_min= 0, 0,  0,  0, -1,  1
-kbound_max= 0, 0,  0,  0, -1,  1
-no_inflow=.true.
-/
-
-&COOLING_PARAMS
-cooling=.true.
-metal=.true.
-z_ave=1.0                           ! To explore explore [0.001, 1]
-self_shielding=.true.
-/
-
-&SF_PARAMS
-sf_virial=.true.                    ! Enable turbulent SF
-sf_model=5                          ! See mom2/star_formation.f90  
-n_star=0.01                         ! SF density threshold. [0.01, 10]
-eps_star=1.0                        ! SF efficency. Should be 1 with sf_virial=.true.       
-T2_star=10.                         ! ISM polytropic temperature
-g_star=1.0                          ! ISM polytropic index
-m_star=10.0                         ! Star particle mass 
-sf_log_properties=.false.           ! do not set this to .true.
-/
-
-&FEEDBACK_PARAMS
-momentum_feedback=1                 ! SN feedback model (from patch/mom2)
-delayed_cooling=.false.             ! Cooling starts from t=0
-yield=0.177703                      ! SN metal yield
-eta_sn=0.162637                     ! SN efficiency
-t_sne=3.0                           ! time for a star to turn into SN
-f_w=0.
-/
-
-&HYDRO_PARAMS
-gamma=1.6666
-courant_factor=0.5
-slope_type=1
-pressure_fix=.true.                 ! Scheme to prevent negative T
-beta_fix=0.5                        ! Controls threshold at which energy is 
-                                    ! truncated
-riemann='hllc'
-/
-
-&REFINE_PARAMS
-mass_sph=1e-5
-m_refine=10*8.
-interpol_type=0
-interpol_var=1
-/
-
-&DICE_PARAMS
-ic_file='agora_sigma10.g2'          ! Name of the initial conditions file
-ic_nfile=1                    ! If greater than one, look for files with name matching ic_file//'.n'
-ic_ifout=1                    ! Change ramses output index for restarts
-ic_format='Gadget2'           ! Format of the initial conditions. 'Gadget1' or 'Gadget2'
-ic_center=0.0,0.0,0.0         ! Shift center parameter. ICs are automatically shifted with boxlen/2
-ic_scale_pos=1.0              ! Scaling factor for the position vector 
-ic_scale_vel=1.0              ! Scaling factor for the velocity vector
-ic_scale_mass=1.0             ! Scaling factor for the mass
-ic_scale_u=1.0                ! Scaling factor for the internal energy
-ic_scale_age=1.0              ! Scaling factor for the particles age
-ic_scale_metal=1.0            ! Scaling factor for the metallicity
-ic_head_name='HEAD'           ! Name of the Header datablock (Gadget2 format only)
-ic_pos_name='POS '            ! Name of the position vector datablock (Gadget2 format only)
-ic_vel_name='VEL '            ! Name of the velocity vector datablock (Gadget2 format only)
-ic_mass_name='MASS'           ! Name of the mass datablock (Gadget2 format only)
-ic_id_name='ID  '             ! Name of the particle identifier datablock (Gadget2 format only)
-ic_u_name='U   '              ! Name of the internal energy datablock (Gadget2 format only)
-ic_metal_name='Z   '          ! Name of the metallicity datablock (Gadget2 format only)
-ic_age_name='AGE '            ! Name of the particle age datablock (Gadget2 format only)
-IG_rho=1.0D-6                 ! Intergalactic gas density
-IG_T2=1.0D7                   ! Intergalactic gas temperature
-IG_metal=0.0                  ! Intergalactic gas metallicity
-amr_struct=.false.            ! Reproduce the AMR structure of the Gadget2 file resulting from a ramses to gadget conversion
-gadget_scale_l=3.085677581282D21 ! Gadget file length unit
-gadget_scale_v=1.0D5             ! Gadget file velocity unit
-gadget_scale_m=1.9891D43         ! Gadget file mass unit
-gadget_scale_t=3.15360e+13       ! Gadget file time unit
-ic_skip_type        = -1      ! Skip specific particle type
-/

low/agora/nml/agora_sigma30.nml

@@ -1,145 +0,0 @@
-&RUN_PARAMS 
-cosmo=.false.
-pic=.true.
-poisson=.true.
-hydro=.true.
-nrestart=0
-nremap=25
-nsubcycle=1,1,2
-ncontrol=1
-rt=.false.
-/
-
-&OUTPUT_PARAMS
-delta_tout=0.5
-tend=50.0
-/
-
-&MOVIE_PARAMS
-imovout=2000
-imov=0
-tendmov=50.
-movie=.false.
-nw_frame=1080
-nh_frame=1080
-levelmax_frame=12
-movie_vars_txt='stars','dens'
-proj_axis='zx'
-xcentre_frame=160,0.0,0.0,0.0 160,0.0,0.0,0.0
-ycentre_frame=160,0.0,0.0,0.0 160,0.0,0.0,0.0
-zcentre_frame=160,0.0,0.0,0.0 160,0.0,0.0,0.0
-deltax_frame=30,0.0 30,0.0
-deltay_frame=30,0.0 30,0.0
-deltaz_frame=30,0.0 30,0.0
-! zoom_only_frame=.true.
-/
-
-&INIT_PARAMS
-filetype='dice'
-initfile(1)='/bulk1/teyssier/hlasco/dice_ics/low/agora/out'
-/
-
-&UNITS_PARAMS
-units_density = 0.677025430198932E-22 ! 1e9 Msol/kpc^3
-units_time    = 0.470430312423675E+15 ! G=1
-units_length  = 0.308567758128200E+22 ! kpc
-/
-
-&AMR_PARAMS
-levelmin=7
-levelmax=12
-ngridmax=3000000
-npartmax=5000000
-boxlen=320.0
-/
-
-&POISSON_PARAMS 
-epsilon=1e-4
-/
-
-&BOUNDARY_PARAMS
-nboundary=6
-bound_type= 2, 2,  2,  2,  2,  2
-ibound_min=-1, 1, -1, -1, -1, -1
-ibound_max=-1, 1,  1,  1,  1,  1
-jbound_min= 0, 0, -1,  1, -1, -1
-jbound_max= 0, 0, -1,  1,  1,  1
-kbound_min= 0, 0,  0,  0, -1,  1
-kbound_max= 0, 0,  0,  0, -1,  1
-no_inflow=.true.
-/
-
-&COOLING_PARAMS
-cooling=.true.
-metal=.true.
-z_ave=1.0                           ! To explore explore [0.001, 1]
-self_shielding=.true.
-/
-
-&SF_PARAMS
-sf_virial=.true.                    ! Enable turbulent SF
-sf_model=5                          ! See mom2/star_formation.f90  
-n_star=0.01                         ! SF density threshold. [0.01, 10]
-eps_star=1.0                        ! SF efficency. Should be 1 with sf_virial=.true.       
-T2_star=10.                         ! ISM polytropic temperature
-g_star=1.0                          ! ISM polytropic index
-m_star=10.0                         ! Star particle mass 
-sf_log_properties=.false.           ! do not set this to .true.
-/
-
-&FEEDBACK_PARAMS
-momentum_feedback=1                 ! SN feedback model (from patch/mom2)
-delayed_cooling=.false.             ! Cooling starts from t=0
-yield=0.177703                      ! SN metal yield
-eta_sn=0.162637                     ! SN efficiency
-t_sne=3.0                           ! time for a star to turn into SN
-f_w=0.
-/
-
-&HYDRO_PARAMS
-gamma=1.6666
-courant_factor=0.5
-slope_type=1
-pressure_fix=.true.                 ! Scheme to prevent negative T
-beta_fix=0.5                        ! Controls threshold at which energy is 
-                                    ! truncated
-riemann='hllc'
-/
-
-&REFINE_PARAMS
-mass_sph=1e-5
-m_refine=10*8.
-interpol_type=0
-interpol_var=1
-/
-
-&DICE_PARAMS
-ic_file='agora_sigma30.g2'          ! Name of the initial conditions file
-ic_nfile=1                    ! If greater than one, look for files with name matching ic_file//'.n'
-ic_ifout=1                    ! Change ramses output index for restarts
-ic_format='Gadget2'           ! Format of the initial conditions. 'Gadget1' or 'Gadget2'
-ic_center=0.0,0.0,0.0         ! Shift center parameter. ICs are automatically shifted with boxlen/2
-ic_scale_pos=1.0              ! Scaling factor for the position vector 
-ic_scale_vel=1.0              ! Scaling factor for the velocity vector
-ic_scale_mass=1.0             ! Scaling factor for the mass
-ic_scale_u=1.0                ! Scaling factor for the internal energy
-ic_scale_age=1.0              ! Scaling factor for the particles age
-ic_scale_metal=1.0            ! Scaling factor for the metallicity
-ic_head_name='HEAD'           ! Name of the Header datablock (Gadget2 format only)
-ic_pos_name='POS '            ! Name of the position vector datablock (Gadget2 format only)
-ic_vel_name='VEL '            ! Name of the velocity vector datablock (Gadget2 format only)
-ic_mass_name='MASS'           ! Name of the mass datablock (Gadget2 format only)
-ic_id_name='ID  '             ! Name of the particle identifier datablock (Gadget2 format only)
-ic_u_name='U   '              ! Name of the internal energy datablock (Gadget2 format only)
-ic_metal_name='Z   '          ! Name of the metallicity datablock (Gadget2 format only)
-ic_age_name='AGE '            ! Name of the particle age datablock (Gadget2 format only)
-IG_rho=1.0D-6                 ! Intergalactic gas density
-IG_T2=1.0D7                   ! Intergalactic gas temperature
-IG_metal=0.0                  ! Intergalactic gas metallicity
-amr_struct=.false.            ! Reproduce the AMR structure of the Gadget2 file resulting from a ramses to gadget conversion
-gadget_scale_l=3.085677581282D21 ! Gadget file length unit
-gadget_scale_v=1.0D5             ! Gadget file velocity unit
-gadget_scale_m=1.9891D43         ! Gadget file mass unit
-gadget_scale_t=3.15360e+13       ! Gadget file time unit
-ic_skip_type        = -1      ! Skip specific particle type
-/

low/agora/nml/base.nml

@@ -1,4 +1,4 @@
-&RUN_PARAMS 
+&RUN_PARAMS
 cosmo=.false.
 pic=.true.
 poisson=.true.
@@ -11,7 +11,7 @@ rt=.false.
 /
 
 &OUTPUT_PARAMS
-delta_tout=0.5
+delta_tout=1
 tend=50.0
 /
 
@@ -19,7 +19,7 @@ tend=50.0
 imovout=2000
 imov=0
 tendmov=50.
-movie=.false.
+movie=.true.
 nw_frame=1080
 nh_frame=1080
 levelmax_frame=12
@@ -36,7 +36,7 @@ deltaz_frame=30,0.0 30,0.0
 
 &INIT_PARAMS
 filetype='dice'
-initfile(1)='/bulk1/teyssier/hlasco/dice_ics/low/agora/out'
+initfile(1)='/users/hlascomb/dice_ics/low/agora/out'
 /
 
 &UNITS_PARAMS
@@ -48,12 +48,12 @@ units_length  = 0.308567758128200E+22 ! kpc
 &AMR_PARAMS
 levelmin=7
 levelmax=12
-ngridmax=3000000
+ngridmax=300000
 npartmax=5000000
 boxlen=320.0
 /
 
-&POISSON_PARAMS 
+&POISSON_PARAMS
 epsilon=1e-4
 /
 
@@ -78,13 +78,13 @@ self_shielding=.true.
 
 &SF_PARAMS
 sf_virial=.true.                    ! Enable turbulent SF
-sf_model=5                          ! See mom2/star_formation.f90  
+sf_model=5                          ! See mom2/star_formation.f90
 n_star=0.01                         ! SF density threshold. [0.01, 10]
-eps_star=1.0                        ! SF efficency. Should be 1 with sf_virial=.true.       
+eps_star=1.0                        ! SF efficency. Should be 1 with sf_virial=.true.
 T2_star=10.                         ! ISM polytropic temperature
 g_star=1.0                          ! ISM polytropic index
-m_star=10.0                         ! Star particle mass 
+m_star=1.0                          ! Star particle mass
-sf_log_properties=.false.           ! do not set this to .true.
+sf_log_properties=.true.            ! do not set this to .true.
 /
 
 &FEEDBACK_PARAMS
@@ -101,7 +101,7 @@ gamma=1.6666
 courant_factor=0.5
 slope_type=1
 pressure_fix=.true.                 ! Scheme to prevent negative T
-beta_fix=0.5                        ! Controls threshold at which energy is 
+beta_fix=0.5                        ! Controls threshold at which energy is
                                     ! truncated
 riemann='hllc'
 /
@@ -114,12 +114,12 @@ interpol_var=1
 /
 
 &DICE_PARAMS
-ic_file='agora_sigma20.g2'          ! Name of the initial conditions file
+ic_file='agora_sigma20.g2'    ! Name of the initial conditions file
 ic_nfile=1                    ! If greater than one, look for files with name matching ic_file//'.n'
 ic_ifout=1                    ! Change ramses output index for restarts
 ic_format='Gadget2'           ! Format of the initial conditions. 'Gadget1' or 'Gadget2'
 ic_center=0.0,0.0,0.0         ! Shift center parameter. ICs are automatically shifted with boxlen/2
-ic_scale_pos=1.0              ! Scaling factor for the position vector 
+ic_scale_pos=1.0              ! Scaling factor for the position vector
 ic_scale_vel=1.0              ! Scaling factor for the velocity vector
 ic_scale_mass=1.0             ! Scaling factor for the mass
 ic_scale_u=1.0                ! Scaling factor for the internal energy

med/agora/config_agora_sigma10.dice

@@ -1,39 +0,0 @@
-# List of galaxy parameters file
-# Putting identical galaxies one beside the other allows DICE
-# to re-use the previous computation,
-# and thus reduces the global execution time
-#                                   x[kpc]  y[kpc]  z[kpc]  vx[km/s]  vy[km/s]  vz[km/s]  spin[deg]  incl[deg]
-Galaxy         params/agora_sigma10.params  0.      0.      0.      0.        0.        0.        0.        0.
-# Number of OpenMP threads
-Nthreads                16
-# Compute the mean inter-particle distance
-# Useful to define the softening length in SPH codes
-MeanPartDist            0
-# Maximum number of component for a galaxy model
-MaxCompNumber           64
-# Output quantities in the rz plane
-OutputRz                0
-# Name of the final IC file
-# The file extension wil be added by DICE
-Filename                out/agora_sigma10
-# File format of the IC file [Gadget1,Gadget2]
-ICformat      Gadget2
-# Cosmological parameters for Virial quantities
-# Default values are taken from Planck cosmology
-# Redshift
-Redshift      1.0
-# Hubble constant
-H0            71.0            
-# Dark energy density parameter
-OmegaL            0.70
-# Matter density parameter
-OmegaM            0.30
-# Spatial curvature density paramter
-OmegaK            0.00
-# Normalize the sum of mass fractions to one
-NormMassFact            0
-# Size of the GSL integration space
-GslWorkspaceSize        10000
-# Integration scheme: 1=GSL QAG (slow) / 2=GSL QAGIU (very slow) / 3= GSL QNG (fast)
-GslIntegrationScheme  3
-

med/agora/config_agora_sigma30.dice

@@ -1,39 +0,0 @@
-# List of galaxy parameters file
-# Putting identical galaxies one beside the other allows DICE
-# to re-use the previous computation,
-# and thus reduces the global execution time
-#                                   x[kpc]  y[kpc]  z[kpc]  vx[km/s]  vy[km/s]  vz[km/s]  spin[deg]  incl[deg]
-Galaxy         params/agora_sigma30.params  0.      0.      0.      0.        0.        0.        0.        0.
-# Number of OpenMP threads
-Nthreads                16
-# Compute the mean inter-particle distance
-# Useful to define the softening length in SPH codes
-MeanPartDist            0
-# Maximum number of component for a galaxy model
-MaxCompNumber           64
-# Output quantities in the rz plane
-OutputRz                0
-# Name of the final IC file
-# The file extension wil be added by DICE
-Filename                out/agora_sigma30
-# File format of the IC file [Gadget1,Gadget2]
-ICformat      Gadget2
-# Cosmological parameters for Virial quantities
-# Default values are taken from Planck cosmology
-# Redshift
-Redshift      1.0
-# Hubble constant
-H0            71.0            
-# Dark energy density parameter
-OmegaL            0.70
-# Matter density parameter
-OmegaM            0.30
-# Spatial curvature density paramter
-OmegaK            0.00
-# Normalize the sum of mass fractions to one
-NormMassFact            0
-# Size of the GSL integration space
-GslWorkspaceSize        10000
-# Integration scheme: 1=GSL QAG (slow) / 2=GSL QAGIU (very slow) / 3= GSL QNG (fast)
-GslIntegrationScheme  3
-

med/agora/nml/agora_sigma00.nml

@@ -1,145 +0,0 @@
-&RUN_PARAMS 
-cosmo=.false.
-pic=.true.
-poisson=.true.
-hydro=.true.
-nrestart=0
-nremap=25
-nsubcycle=1,1,2
-ncontrol=1
-rt=.false.
-/
-
-&OUTPUT_PARAMS
-delta_tout=0.25
-tend=50.0
-/
-
-&MOVIE_PARAMS
-imovout=2000
-imov=0
-tendmov=50.
-movie=.true.
-nw_frame=1080
-nh_frame=1080
-levelmax_frame=14
-movie_vars_txt='stars','dens','temp'
-proj_axis='zx'
-xcentre_frame=120,0.0,0.0,0.0 120,0.0,0.0,0.0
-ycentre_frame=120,0.0,0.0,0.0 120,0.0,0.0,0.0
-zcentre_frame=120,0.0,0.0,0.0 120,0.0,0.0,0.0
-deltax_frame=30,0.0 30,0.0
-deltay_frame=30,0.0 30,0.0
-deltaz_frame=30,0.0 30,0.0
-! zoom_only_frame=.true.
-/
-
-&INIT_PARAMS
-filetype='dice'
-initfile(1)='/bulk1/teyssier/hlasco/dice_ics/med/agora/out'
-/
-
-&UNITS_PARAMS
-units_density = 0.677025430198932E-22 ! 1e9 Msol/kpc^3
-units_time    = 0.470430312423675E+15 ! G=1
-units_length  = 0.308567758128200E+22 ! kpc
-/
-
-&AMR_PARAMS
-levelmin=7
-levelmax=14
-ngridmax=3000000
-npartmax=5000000
-boxlen=320.0
-/
-
-&POISSON_PARAMS 
-epsilon=1e-4
-/
-
-&BOUNDARY_PARAMS
-nboundary=6
-bound_type= 2, 2,  2,  2,  2,  2
-ibound_min=-1, 1, -1, -1, -1, -1
-ibound_max=-1, 1,  1,  1,  1,  1
-jbound_min= 0, 0, -1,  1, -1, -1
-jbound_max= 0, 0, -1,  1,  1,  1
-kbound_min= 0, 0,  0,  0, -1,  1
-kbound_max= 0, 0,  0,  0, -1,  1
-no_inflow=.true.
-/
-
-&COOLING_PARAMS
-cooling=.true.
-metal=.true.
-z_ave=1.0                           ! To explore explore [0.001, 1]
-self_shielding=.true.
-/
-
-&SF_PARAMS
-sf_virial=.true.                    ! Enable turbulent SF
-sf_model=5                          ! See mom2/star_formation.f90  
-n_star=0.01                         ! SF density threshold. [0.01, 10]
-eps_star=1.0                        ! SF efficency. Should be 1 with sf_virial=.true.       
-T2_star=10.                         ! ISM polytropic temperature
-g_star=1.0                          ! ISM polytropic index
-m_star=10.0                         ! Star particle mass 
-sf_log_properties=.false.           ! do not set this to .true.
-/
-
-&FEEDBACK_PARAMS
-momentum_feedback=1                 ! SN feedback model (from patch/mom2)
-delayed_cooling=.false.             ! Cooling starts from t=0
-yield=0.177703                      ! SN metal yield
-eta_sn=0.162637                     ! SN efficiency
-t_sne=3.0                           ! time for a star to turn into SN
-f_w=0.
-/
-
-&HYDRO_PARAMS
-gamma=1.6666
-courant_factor=0.5
-slope_type=1
-pressure_fix=.true.                 ! Scheme to prevent negative T
-beta_fix=0.5                        ! Controls threshold at which energy is 
-                                    ! truncated
-riemann='hllc'
-/
-
-&REFINE_PARAMS
-mass_sph=1e-6
-m_refine=10*8.
-interpol_type=0
-interpol_var=1
-/
-
-&DICE_PARAMS
-ic_file='agora_sigma00.g2'          ! Name of the initial conditions file
-ic_nfile=1                    ! If greater than one, look for files with name matching ic_file//'.n'
-ic_ifout=1                    ! Change ramses output index for restarts
-ic_format='Gadget2'           ! Format of the initial conditions. 'Gadget1' or 'Gadget2'
-ic_center=0.0,0.0,0.0         ! Shift center parameter. ICs are automatically shifted with boxlen/2
-ic_scale_pos=1.0              ! Scaling factor for the position vector 
-ic_scale_vel=1.0              ! Scaling factor for the velocity vector
-ic_scale_mass=1.0             ! Scaling factor for the mass
-ic_scale_u=1.0                ! Scaling factor for the internal energy
-ic_scale_age=1.0              ! Scaling factor for the particles age
-ic_scale_metal=1.0            ! Scaling factor for the metallicity
-ic_head_name='HEAD'           ! Name of the Header datablock (Gadget2 format only)
-ic_pos_name='POS '            ! Name of the position vector datablock (Gadget2 format only)
-ic_vel_name='VEL '            ! Name of the velocity vector datablock (Gadget2 format only)
-ic_mass_name='MASS'           ! Name of the mass datablock (Gadget2 format only)
-ic_id_name='ID  '             ! Name of the particle identifier datablock (Gadget2 format only)
-ic_u_name='U   '              ! Name of the internal energy datablock (Gadget2 format only)
-ic_metal_name='Z   '          ! Name of the metallicity datablock (Gadget2 format only)
-ic_age_name='AGE '            ! Name of the particle age datablock (Gadget2 format only)
-IG_rho=1.0D-6                 ! Intergalactic gas density
-IG_T2=1.0D7                   ! Intergalactic gas temperature
-IG_metal=0.0                  ! Intergalactic gas metallicity
-amr_struct=.false.            ! Reproduce the AMR structure of the Gadget2 file resulting from a ramses to gadget conversion
-gadget_scale_l=3.085677581282D21 ! Gadget file length unit
-gadget_scale_v=1.0D5             ! Gadget file velocity unit
-gadget_scale_m=1.9891D43         ! Gadget file mass unit
-gadget_scale_t=3.15360e+13       ! Gadget file time unit
-ic_skip_type        = -1      ! Skip specific particle type
-/

med/agora/nml/agora_sigma10.nml

@@ -1,145 +0,0 @@
-&RUN_PARAMS 
-cosmo=.false.
-pic=.true.
-poisson=.true.
-hydro=.true.
-nrestart=0
-nremap=25
-nsubcycle=1,1,2
-ncontrol=1
-rt=.false.
-/
-
-&OUTPUT_PARAMS
-delta_tout=0.25
-tend=50.0
-/
-
-&MOVIE_PARAMS
-imovout=2000
-imov=0
-tendmov=50.
-movie=.true.
-nw_frame=1080
-nh_frame=1080
-levelmax_frame=14
-movie_vars_txt='stars','dens','temp'
-proj_axis='zx'
-xcentre_frame=160,0.0,0.0,0.0 160,0.0,0.0,0.0
-ycentre_frame=160,0.0,0.0,0.0 160,0.0,0.0,0.0
-zcentre_frame=160,0.0,0.0,0.0 160,0.0,0.0,0.0
-deltax_frame=30,0.0 30,0.0
-deltay_frame=30,0.0 30,0.0
-deltaz_frame=30,0.0 30,0.0
-! zoom_only_frame=.true.
-/
-
-&INIT_PARAMS
-filetype='dice'
-initfile(1)='/bulk1/teyssier/hlasco/dice_ics/med/agora/out'
-/
-
-&UNITS_PARAMS
-units_density = 0.677025430198932E-22 ! 1e9 Msol/kpc^3
-units_time    = 0.470430312423675E+15 ! G=1
-units_length  = 0.308567758128200E+22 ! kpc
-/
-
-&AMR_PARAMS
-levelmin=7
-levelmax=14
-ngridmax=3000000
-npartmax=5000000
-boxlen=320.0
-/
-
-&POISSON_PARAMS 
-epsilon=1e-4
-/
-
-&BOUNDARY_PARAMS
-nboundary=6
-bound_type= 2, 2,  2,  2,  2,  2
-ibound_min=-1, 1, -1, -1, -1, -1
-ibound_max=-1, 1,  1,  1,  1,  1
-jbound_min= 0, 0, -1,  1, -1, -1
-jbound_max= 0, 0, -1,  1,  1,  1
-kbound_min= 0, 0,  0,  0, -1,  1
-kbound_max= 0, 0,  0,  0, -1,  1
-no_inflow=.true.
-/
-
-&COOLING_PARAMS
-cooling=.true.
-metal=.true.
-z_ave=1.0                           ! To explore explore [0.001, 1]
-self_shielding=.true.
-/
-
-&SF_PARAMS
-sf_virial=.true.                    ! Enable turbulent SF
-sf_model=5                          ! See mom2/star_formation.f90  
-n_star=0.01                         ! SF density threshold. [0.01, 10]
-eps_star=1.0                        ! SF efficency. Should be 1 with sf_virial=.true.       
-T2_star=10.                         ! ISM polytropic temperature
-g_star=1.0                          ! ISM polytropic index
-m_star=10.0                         ! Star particle mass 
-sf_log_properties=.false.           ! do not set this to .true.
-/
-
-&FEEDBACK_PARAMS
-momentum_feedback=1                 ! SN feedback model (from patch/mom2)
-delayed_cooling=.false.             ! Cooling starts from t=0
-yield=0.177703                      ! SN metal yield
-eta_sn=0.162637                     ! SN efficiency
-t_sne=3.0                           ! time for a star to turn into SN
-f_w=0.
-/
-
-&HYDRO_PARAMS
-gamma=1.6666
-courant_factor=0.5
-slope_type=1
-pressure_fix=.true.                 ! Scheme to prevent negative T
-beta_fix=0.5                        ! Controls threshold at which energy is 
-                                    ! truncated
-riemann='hllc'
-/
-
-&REFINE_PARAMS
-mass_sph=1e-6
-m_refine=10*8.
-interpol_type=0
-interpol_var=1
-/
-
-&DICE_PARAMS
-ic_file='agora_sigma10.g2'          ! Name of the initial conditions file
-ic_nfile=1                    ! If greater than one, look for files with name matching ic_file//'.n'
-ic_ifout=1                    ! Change ramses output index for restarts
-ic_format='Gadget2'           ! Format of the initial conditions. 'Gadget1' or 'Gadget2'
-ic_center=0.0,0.0,0.0         ! Shift center parameter. ICs are automatically shifted with boxlen/2
-ic_scale_pos=1.0              ! Scaling factor for the position vector 
-ic_scale_vel=1.0              ! Scaling factor for the velocity vector
-ic_scale_mass=1.0             ! Scaling factor for the mass
-ic_scale_u=1.0                ! Scaling factor for the internal energy
-ic_scale_age=1.0              ! Scaling factor for the particles age
-ic_scale_metal=1.0            ! Scaling factor for the metallicity
-ic_head_name='HEAD'           ! Name of the Header datablock (Gadget2 format only)
-ic_pos_name='POS '            ! Name of the position vector datablock (Gadget2 format only)
-ic_vel_name='VEL '            ! Name of the velocity vector datablock (Gadget2 format only)
-ic_mass_name='MASS'           ! Name of the mass datablock (Gadget2 format only)
-ic_id_name='ID  '             ! Name of the particle identifier datablock (Gadget2 format only)
-ic_u_name='U   '              ! Name of the internal energy datablock (Gadget2 format only)
-ic_metal_name='Z   '          ! Name of the metallicity datablock (Gadget2 format only)
-ic_age_name='AGE '            ! Name of the particle age datablock (Gadget2 format only)
-IG_rho=1.0D-6                 ! Intergalactic gas density
-IG_T2=1.0D7                   ! Intergalactic gas temperature
-IG_metal=0.0                  ! Intergalactic gas metallicity
-amr_struct=.false.            ! Reproduce the AMR structure of the Gadget2 file resulting from a ramses to gadget conversion
-gadget_scale_l=3.085677581282D21 ! Gadget file length unit
-gadget_scale_v=1.0D5             ! Gadget file velocity unit
-gadget_scale_m=1.9891D43         ! Gadget file mass unit
-gadget_scale_t=3.15360e+13       ! Gadget file time unit
-ic_skip_type        = -1      ! Skip specific particle type
-/

med/agora/nml/base.nml

@@ -1,4 +1,4 @@
-&RUN_PARAMS 
+&RUN_PARAMS
 cosmo=.false.
 pic=.true.
 poisson=.true.
@@ -11,7 +11,7 @@ rt=.false.
 /
 
 &OUTPUT_PARAMS
-delta_tout=0.25
+delta_tout=1
 tend=50.0
 /
 
@@ -36,7 +36,7 @@ deltaz_frame=30,0.0 30,0.0
 
 &INIT_PARAMS
 filetype='dice'
-initfile(1)='/bulk1/teyssier/hlasco/dice_ics/med/agora/out'
+initfile(1)='/users/hlascomb/dice_ics/med/agora/out'
 /
 
 &UNITS_PARAMS
@@ -48,12 +48,12 @@ units_length  = 0.308567758128200E+22 ! kpc
 &AMR_PARAMS
 levelmin=7
 levelmax=14
-ngridmax=3000000
+ngridmax=500000
-npartmax=5000000
+npartmax=3000000
 boxlen=320.0
 /
 
-&POISSON_PARAMS 
+&POISSON_PARAMS
 epsilon=1e-4
 /
 
@@ -78,13 +78,13 @@ self_shielding=.true.
 
 &SF_PARAMS
 sf_virial=.true.                    ! Enable turbulent SF
-sf_model=5                          ! See mom2/star_formation.f90  
+sf_model=5                          ! See mom2/star_formation.f90
 n_star=0.01                         ! SF density threshold. [0.01, 10]
-eps_star=1.0                        ! SF efficency. Should be 1 with sf_virial=.true.       
+eps_star=1.0                        ! SF efficency. Should be 1 with sf_virial=.true.
 T2_star=10.                         ! ISM polytropic temperature
 g_star=1.0                          ! ISM polytropic index
-m_star=10.0                         ! Star particle mass 
+m_star=1.0                         ! Star particle mass
-sf_log_properties=.false.           ! do not set this to .true.
+sf_log_properties=.true.           ! do not set this to .true.
 /
 
 &FEEDBACK_PARAMS
@@ -101,7 +101,7 @@ gamma=1.6666
 courant_factor=0.5
 slope_type=1
 pressure_fix=.true.                 ! Scheme to prevent negative T
-beta_fix=0.5                        ! Controls threshold at which energy is 
+beta_fix=0.5                        ! Controls threshold at which energy is
                                     ! truncated
 riemann='hllc'
 /
@@ -114,12 +114,12 @@ interpol_var=1
 /
 
 &DICE_PARAMS
-ic_file='agora_sigma30.g2'          ! Name of the initial conditions file
+ic_file='agora_sigma20.g2'
 ic_nfile=1                    ! If greater than one, look for files with name matching ic_file//'.n'
 ic_ifout=1                    ! Change ramses output index for restarts
 ic_format='Gadget2'           ! Format of the initial conditions. 'Gadget1' or 'Gadget2'
 ic_center=0.0,0.0,0.0         ! Shift center parameter. ICs are automatically shifted with boxlen/2
-ic_scale_pos=1.0              ! Scaling factor for the position vector 
+ic_scale_pos=1.0              ! Scaling factor for the position vector
 ic_scale_vel=1.0              ! Scaling factor for the velocity vector
 ic_scale_mass=1.0             ! Scaling factor for the mass
 ic_scale_u=1.0                ! Scaling factor for the internal energy