Changeset 4035 for trunk/LMDZ.MARS/libf/phymars/physiq_mod.F

Timestamp:

Jan 30, 2026, 12:40:55 PM (4 weeks ago)

Author:

emillour

Message:

Mars PCM:
Cleanup and correction around thermospheric processes: ensure that the
tendencies from previous processes are included before being sent to
the next process (was not the case for conduction).
EM

File:

• trunk/LMDZ.MARS/libf/phymars/physiq_mod.F (modified) (3 diffs)

trunk/LMDZ.MARS/libf/phymars/physiq_mod.F

@@ -2130 +2130 @@
 c     ------------------
 
-!#ifndef MESOSCALE
 c        --------------
 c        photochemistry :
@@ -2200 +2199 @@
 
          END IF  ! of IF (photochem)
-!#endif
-
-
-!#ifndef MESOSCALE
+
+
 c-----------------------------------------------------------------------
 c   10. THERMOSPHERE CALCULATION
@@ -2210 +2207 @@
       if (callthermos) then
         call thermosphere(ngrid,nlayer,nq,zplev,zplay,dist_sol,
-     $     mu0,ptimestep,ptime,zday,tsurf_meshavg,zzlev,zzlay,
-     &     pt,pq,pu,pv,pdt,pdq,
-     $     zdteuv,zdtconduc,zdumolvis,zdvmolvis,zdqmoldiff,
-     $     PhiEscH,PhiEscH2,PhiEscD)
-
-        DO l=1,nlayer
-          DO ig=1,ngrid
-            dtrad(ig,l)=dtrad(ig,l)+zdteuv(ig,l)
-            pdt(ig,l)=pdt(ig,l)+zdtconduc(ig,l)+zdteuv(ig,l)
-            pdv(ig,l)=pdv(ig,l)+zdvmolvis(ig,l)
-            pdu(ig,l)=pdu(ig,l)+zdumolvis(ig,l)
-            DO iq=1, nq
-              pdq(ig,l,iq)=pdq(ig,l,iq)+zdqmoldiff(ig,l,iq)
-            ENDDO
-          ENDDO
-        ENDDO
+     &     mu0,ptimestep,zday,tsurf_meshavg,zzlev,zzlay,
+     &     pt,pq,pu,pv,pdt,pdq,pdu,pdv,
+     &     zdteuv,zdtconduc,zdumolvis,zdvmolvis,zdqmoldiff,
+     &     PhiEscH,PhiEscH2,PhiEscD)
+        
+        ! tendencies zdteuv,zdtconduc,zdumolvis,zdvmolvis,zdqmoldiff
+        ! already included in pdt,pdu,pdv,pdq by thermosphere
+        
+        dtrad(1:ngrid,1:nlayer)=dtrad(1:ngrid,1:nlayer)
+     &                               +zdteuv(1:ngrid,1:nlayer)
 
       endif ! of if (callthermos)
-!#endif
 
 c-----------------------------------------------------------------------