Changeset 3146

Timestamp:

Dec 2, 2023, 3:01:42 PM (13 months ago)

Author:

llange

Message:

Mars PCM
Bug fix in init_testphys and testphys. Index for water vapor and ice were not initialized
(but used in testphys1d when atm_wat_profile > 0). It is fixed by calling
initracer in init_testphys1d_mod.
LL

Location:

trunk/LMDZ.MARS

Files:

• changelog.txt (modified) (1 diff)
• libf/phymars/dyn1d/init_testphys1d_mod.F90 (modified) (2 diffs)

trunk/LMDZ.MARS/changelog.txt

@@ -4369 +4369 @@
 Remember to use full-layers model configuration (64x48x73) when you try to call this scheme because we need to evaluate the saturation altitude of each monochromatic wave.
 More than ~75% of the waves saturated at altitude above 120 km. 
+
+== 02/12/2023 == LL
+Bug fix in init_testphys and testphys. Index for water vapor and ice were not initialized
+(but used in testphys1d when atm_wat_profile > 0). It is fixed by calling
+initracer in init_testphys1d_mod.  

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

@@ -388 +388 @@
 call ini_fillgeom(1,latitude,longitude,(/1.0/))
 call conf_phys(1,llm,nq)
-
+call initracer(ngrid,nq,qsurf)
 ! In 1D model physics are called every time step
 ! ovverride iphysiq value that has been set by conf_phys
@@ -509 +509 @@
 ! CO2 ice on the surface
 ! ----------------------
-! get the index of co2 tracer (not known at this stage)
-igcm_co2 = 0
-do iq = 1,nq
-    if (trim(tname(iq)) == "co2") igcm_co2 = iq
-enddo
-if (igcm_co2 == 0) error stop "init_testphys1d error, missing co2 tracer!"
 
 if (.not. therestartfi) then