Changeset 3146 for trunk/LMDZ.MARS/libf/phymars

Timestamp:

Dec 2, 2023, 3:01:42 PM (19 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

File:

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

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