Changeset 3848 for trunk/LMDZ.MARS/libf/phymars/dyn1d/init_testphys1d_mod.F90

Timestamp:

Jul 15, 2025, 12:13:09 PM (3 weeks ago)

Author:

jbclement

Message:

Mars PCM 1D:
In the 1D model, if no input profile found, the 'dust_mass' and 'dust_number' tracers are initialized according to "Conrath dust" in line with what it's done for CO2 and HDO, instead of being set to 0 which caused bugs.
JBC

File:

• trunk/LMDZ.MARS/libf/phymars/dyn1d/init_testphys1d_mod.F90 (modified) (4 diffs)

trunk/LMDZ.MARS/libf/phymars/dyn1d/init_testphys1d_mod.F90

@@ -319 +319 @@
 write(*,*) " psurf = ",psurf
 
+! Compute pressures and altitudes of atmospheric levels
+! -----------------------------------------------------
+! Vertical Coordinates
+! """"""""""""""""""""
+hybrid = .true.
+write(*,*)'Hybrid coordinates?'
+call getin("hybrid",hybrid)
+write(*,*) " hybrid = ", hybrid
+
+! Reference pressures
+pa = 20.     ! Transition pressure (for hybrid coord.)
+preff = 610. ! Reference surface pressure
+
+call disvert_noterre
+! Now that disvert has been called, initialize module vertical_layers_mod
+call init_vertical_layers(nlayer,preff,scaleheight,ap,bp,aps,bps,presnivs,pseudoalt)
+
+plev = ap + psurf*bp
+play = aps + psurf*bps
+zlay = -200.*r*log(play/plev(1))/g
+
 ! Coefficient to control the surface pressure change
 ! To be defined in "callphys.def" so that the PEM can read it asa well
@@ -331 +352 @@
     error stop 'Please, specify a correct value!'
 endif
-
-! Reference pressures
-pa = 20.     ! Transition pressure (for hybrid coord.)
-preff = 610. ! Reference surface pressure
 
 ! Aerosol properties
@@ -512 +529 @@
 write(*,*) "init_testphys1d: initializing tracers"
 if (.not. therestart1D) then
-    call read_profile(nq,nlayer,qsurf(1,:,1),q)
+    call read_profile(nq,nlayer,qsurf(1,:,1),q,play)
     do iq = 1,nq
         qsurf(1,iq,:) = qsurf(1,iq,1)
@@ -631 +648 @@
     endif
 endif !(.not. therestartfi)
-
-! Compute pressures and altitudes of atmospheric levels
-! -----------------------------------------------------
-! Vertical Coordinates
-! """"""""""""""""""""
-hybrid = .true.
-write(*,*)'Hybrid coordinates?'
-call getin("hybrid",hybrid)
-write(*,*) " hybrid = ", hybrid
-
-call disvert_noterre
-! Now that disvert has been called, initialize module vertical_layers_mod
-call init_vertical_layers(nlayer,preff,scaleheight,ap,bp,aps,bps,presnivs,pseudoalt)
-
-plev = ap + psurf*bp
-play = aps + psurf*bps
-zlay = -200.*r*log(play/plev(1))/g
 
 ! Initialize temperature profile