Changeset 3201 for trunk/LMDZ.MARS

Timestamp:

Feb 5, 2024, 11:31:02 AM (10 months ago)

Author:

jbclement

Message:

Mars PCM:
Small modification of the change introduced in r3200 to make the code simpler and faster.
JBC

Location:

trunk/LMDZ.MARS

Files:

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

trunk/LMDZ.MARS/changelog.txt

@@ -4463 +4463 @@
 == 02/02/2024 == JBC
 Small update following r3188 to keep the CO2 mixing ration constant if "CO2cond_ps = 0" + Some cleanings.
+
+== 05/02/2024 == JBC
+Small modification of the change introduced in r3200 to make the code simpler and faster.

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

@@ -6 +6 @@
 use phyredem,            only: physdem0, physdem1
 use watersat_mod,        only: watersat
-use tracer_mod,          only: igcm_h2o_vap, igcm_h2o_ice, noms
+use tracer_mod,          only: igcm_h2o_vap, igcm_h2o_ice, igcm_co2, noms
 use comcstfi_h,          only: pi, rad, omeg, g, mugaz, rcp, r, cpp
 use time_phylmdz_mod,    only: daysec, day_step
@@ -270 +270 @@
 
     ! Increment tracers
-    if (abs(CO2cond_ps) < 1.e-10) then
-        do iq = 1,nq
-            if (iq == igcm_co2) cycle
-            q(1,:,iq) = q(1,:,iq) + dttestphys*dq(1,:,iq)
-        enddo
-    else
-        q = q + dttestphys*dq
-    endif
+    if (abs(CO2cond_ps) < 1.e-10) dq(1,:,igcm_co2) = 0.
+    q = q + dttestphys*dq
 enddo ! End of time stepping loop (idt=1,ndt)