Adding -expand_z option

src/init_mcfost.f90

@@ -200,6 +200,7 @@ subroutine set_default_variables()
   lrandomize_gap = .false.
   lrandomize_outside_gap = .false.
   lcentre_on_sink = .false.
+  lexpand_z = .false.
   lwrite_column_density = .false.
   lwrite_mol_column_density = .false.
   lwrite_velocity = .false.
@@ -1385,6 +1386,12 @@ subroutine initialisation_mcfost()
         call get_command_argument(i_arg,s)
         read(s,*) gap_factor
         i_arg = i_arg + 1
+     case("-expand_z","-expand-z","-ez")
+        i_arg = i_arg + 1
+        lexpand_z = .true.
+        call get_command_argument(i_arg,s)
+        read(s,*) expand_z_factor
+        i_arg = i_arg + 1
      case("-cd","-column_density")
         i_arg = i_arg + 1
         lwrite_column_density = .true.
@@ -2047,6 +2054,7 @@ subroutine display_help()
   write(*,*) "        : -no_vz : force the vertical velocities to be 0"
   write(*,*) "        : -vphi_Kep : force the azimuthal velocities to be Keplerian"
   write(*,*) "        : -centre_on_sink <number> : centre the model on the sink particle"
+  write(*,*) "        : -expand_z <factor> : multipky all z values by factor"
   write(*,*) "        : -SPH_amin <size> [mum] : force the grain size that follow the gas"
   write(*,*) "        : -SPH_amax <size> [mum] : force the grain size that follow the dust"
   write(*,*) "                                   (only works with 1 grain size dump)"

src/mcfost_env.f90

@@ -5,7 +5,7 @@ module mcfost_env
   implicit none
 
   real, parameter :: mcfost_version = 4.1
-  character(8), parameter :: mcfost_release = "4.1.05"
+  character(8), parameter :: mcfost_release = "4.1.06"
   real, parameter :: required_utils_version = 4.0
 
   !character(len=128) :: web_server    = "https://ipag.osug.fr/public/pintec/"

src/mess_up_SPH.f90

@@ -5,7 +5,7 @@ module mess_up_SPH
   implicit none
 
   ! mask is set to 1 for hide particle == opacity set to 0, and set to 2 to delete particle before tesselation
-  public :: mask_Hill_sphere, randomize_azimuth, randomize_gap, mask_inside_rsph, mask_outside_rsph, mask_above_theta
+  public :: mask_Hill_sphere, randomize_azimuth, randomize_gap, mask_inside_rsph, mask_outside_rsph, mask_above_theta, expand_z
 
   private
 
@@ -242,6 +242,25 @@ end subroutine randomize_azimuth
 
   !*********************************************************
 
+  subroutine expand_z(np, xyzh, vxyzu, factor)
+
+    integer, intent(in) :: np
+    real(kind=dp), dimension(:,:), intent(inout) :: xyzh, vxyzu
+    real(kind=dp) :: factor
+
+    integer :: i
+
+    do i=1, np
+       xyzh(3,i) = xyzh(3,i) * factor
+       vxyzu(3,i) = vxyzu(3,i) * factor
+    enddo
+
+    return
+
+  end subroutine expand_z
+
+  !*********************************************************
+
   subroutine randomize_gap(np, nptmass, xyzh, vxyzu, xyzmh_ptmass ,udist, factor, inside)
     ! Randomly rotate all the particles inside or outside cylinder
     ! a width +/- factor * r_Hill of the planet

src/parameters.f90

@@ -159,9 +159,9 @@ module parametres
   ! Phantom
   logical :: ldudt_implicit, lscale_length_units, lscale_mass_units, lignore_dust
   logical :: ldelete_Hill_sphere, lrandomize_Voronoi, lrandomize_azimuth, lrandomize_gap, lrandomize_outside_gap, lcentre_on_sink
-  logical :: ldelete_inside_rsph, ldelete_outside_rsph, ldelete_above_theta
+  logical :: ldelete_inside_rsph, ldelete_outside_rsph, ldelete_above_theta, lexpand_z
   real(kind=dp) :: ufac_implicit,scale_length_units_factor,scale_mass_units_factor,correct_density_factor_elongated_cells
-  real(kind=dp) :: SPH_amin, SPH_amax, fluffyness, gap_factor, rsph_min, rsph_max
+  real(kind=dp) :: SPH_amin, SPH_amax, fluffyness, gap_factor, rsph_min, rsph_max, expand_z_factor
   logical :: lupdate_velocities, lno_vr, lno_vz, lvphi_Kep, lfluffy, lnot_random_Voronoi, lignore_sink
   integer :: isink_centre
 

src/read_phantom.f90

@@ -753,6 +753,8 @@ subroutine modify_dump(np, nptmass, xyzh, vxyzu, xyzmh_ptmass, ulength, mask)
        enddo
     endif
 
+    if (lexpand_z) call expand_z(np, xyzh, vxyzu, expand_z_factor)
+
     return
 
   end subroutine modify_dump