Changeset 4035 for trunk/LMDZ.MARS/libf/aeronomars/molvis.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/aeronomars/molvis.F (modified) (4 diffs)

trunk/LMDZ.MARS/libf/aeronomars/molvis.F

@@ -6 +6 @@
 
       SUBROUTINE molvis(ngrid,nlayer,ptimestep,
-     &            pplay,pplev,pt,pdteuv,pdtconduc
-     $           ,pvel,tsurf,zzlev,zzlay,zdvelmolvis)
+     &            pplay,pplev,pt,pdt,
+     $            pvel,tsurf,zzlev,zzlay,zdvelmolvis)
       
       use conc_mod, only: cpnew, Akknew, rnew
@@ -31 +31 @@
       integer,intent(in) :: ngrid ! number of atmospheric columns
       integer,intent(in) :: nlayer ! number of atmospheric layers
-      REAL ptimestep
-      REAL pplay(ngrid,nlayer)
-      REAL pplev(ngrid,nlayer+1)
-      REAL zzlay(ngrid,nlayer)
-      REAL zzlev(ngrid,nlayer+1)
-      real pt(ngrid,nlayer)
-      real tsurf(ngrid)
-      REAL pvel(ngrid,nlayer)
-      REAL pdvel(ngrid,nlayer)
-      real pdteuv(ngrid,nlayer)
-      real pdtconduc(ngrid,nlayer)
+      REAL,INTENT(IN) :: ptimestep ! physics time step (s)
+      REAL,INTENT(IN) :: pplay(ngrid,nlayer)
+      REAL,INTENT(IN) :: pplev(ngrid,nlayer+1)
+      REAL,INTENT(IN) :: zzlay(ngrid,nlayer)
+      REAL,INTENT(IN) :: zzlev(ngrid,nlayer+1)
+      real,intent(in) :: pt(ngrid,nlayer) ! input temperature (K)
+      real,intent(in) :: pdt(ngrid,nlayer) ! input tendency on temperature (K/s)
+      real,intent(in) :: tsurf(ngrid) ! input ssurface temperature (K)
+      REAL,INTENT(IN) :: pvel(ngrid,nlayer) ! wind velocity (m/s)
 
-      real zdvelmolvis(ngrid,nlayer)
+      real,intent(out) :: zdvelmolvis(ngrid,nlayer) ! tendency on velocity (m/s/s)
 
 c   local:
@@ -99 +97 @@
       do ig=1,ngrid
 
-        zt(1)=pt(ig,1)+(pdteuv(ig,1)+pdtconduc(ig,1))*ptimestep
+        zt(1)=pt(ig,1)+(pdt(ig,1))*ptimestep
         zvel(1)=pvel(ig,1)
         zlay(1)=zzlay(ig,1)
@@ -105 +103 @@
 
         do l=2,nz
-          zt(l)=pt(ig,l)+(pdteuv(ig,l)+pdtconduc(ig,l))*ptimestep
+          zt(l)=pt(ig,l)+(pdt(ig,l))*ptimestep
           zvel(l)=pvel(ig,l)
           zlay(l)=zzlay(ig,l)