Changeset 4052 for trunk/LMDZ.MARS

Timestamp:

Feb 6, 2026, 12:39:37 PM (4 days ago)

Author:

emillour

Message:

Mars PCM:
Add posibilty for user to input lifted dust effective radius in callphys.def
as "reff_lift_dust" (default is the previously hard-coded value of 3.0e6 m)
EM

Location:

trunk/LMDZ.MARS

Files:

• changelog.txt (modified) (1 diff)
• libf/phymars/initracer.F90 (modified) (5 diffs)

trunk/LMDZ.MARS/changelog.txt

@@ -5042 +5042 @@
 improvedclouds: from now on, the tendencies from previous physics are added
 all-at-once rather than distributed over the microphysics substeps.
+
+== 06/02/2026 == EM
+Add posibilty for user to input lifted dust effective radius in callphys.def
+as "reff_lift_dust" (default is the previously hard-coded value of 3.0e6 m)

trunk/LMDZ.MARS/libf/phymars/initracer.F90

@@ -14 +14 @@
                           callnlte, nltemodel, dustinjection,    &
                           co2clouds, co2useh2o, meteo_flux
+use ioipsl_getin_p_mod, only : getin_p
 
 implicit none
@@ -618 +619 @@
     ! '''''''''''
     reff_lift = 3.0e-6 !3.e-6 !Effective radius of lifted dust (m)
+    call getin_p("reff_lift_dust",reff_lift)
+    
     alpha_devil(igcm_dust_mass)=9.e-9   !  dust devil lift mass coeff
     ! alpha_lift(igcm_dust_mass)=3.0e-15  !  Lifted mass coeff
@@ -629 +632 @@
     alpha_lift(igcm_dust_mass)=1.e-6 !1.e-6 !Lifted mass coeff
     if (dustinjection >= 1) then
-        reff_lift = 3.0e-6 ! Effective radius of lifted dust (m)
+        !reff_lift = 3.0e-6 ! Effective radius of lifted dust (m)! already set
         alpha_lift(igcm_dust_mass)=(4/3.)*reff_lift*rho_dust/2.4
     endif
@@ -649 +652 @@
     ! normal dust
     if (rdstorm) then
-        reff_storm=3.e-6 ! reff_lift !3.e-6
+        reff_storm=reff_lift !3.e-6
         r0_storm=reff_storm/ref_r0
         rho_q(igcm_stormdust_mass)=rho_dust
@@ -672 +675 @@
         rho_q(igcm_topdust_mass)=rho_dust
         rho_q(igcm_topdust_number)=rho_dust
-        radius(igcm_topdust_mass) = 3.e-6
-        radius(igcm_topdust_number) = 3.e-6
+        radius(igcm_topdust_mass) = reff_lift !3.e-6
+        radius(igcm_topdust_number) = reff_lift !3.e-6
     end if !(topflows)
     ! ----------------------------------------------------------------------