Changeset 4035 for trunk/LMDZ.MARS/libf/aeronomars/moldiff_red.F90

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/moldiff_red.F90 (modified) (4 diffs)

trunk/LMDZ.MARS/libf/aeronomars/moldiff_red.F90

@@ -8 +8 @@
 
 subroutine moldiff_red(ngrid,nlayer,nq,pplay,pplev,pt,pdt,pq,pdq,&
-                       ptimestep,zzlay,pdteuv,pdtconduc,pdqdiff,&
+                       ptimestep,pdqdiff,&
                        PhiEscH,PhiEscH2,PhiEscD)
 
@@ -31 +31 @@
       integer,intent(in) :: nlayer ! number of atmospheric layers
       integer,intent(in) :: nq ! number of advected tracers
-      real ptimestep
-      real pplay(ngrid,nlayer)
-      real zzlay(ngrid,nlayer)
-      real pplev(ngrid,nlayer+1)
-      real pq(ngrid,nlayer,nq)
-      real pdq(ngrid,nlayer,nq)
-      real pt(ngrid,nlayer)
-      real pdt(ngrid,nlayer)
-      real pdteuv(ngrid,nlayer)
-      real pdtconduc(ngrid,nlayer)
-      real pdqdiff(ngrid,nlayer,nq)
-      real*8 PhiEscH,PhiEscH2,PhiEscD
+      real,intent(in) :: ptimestep ! physics time step (s)
+      real,intent(in) :: pplay(ngrid,nlayer) ! mid-layer pressure (Pa)
+      real,intent(in) :: pplev(ngrid,nlayer+1) ! pressure (Pa) at layer boundaries 
+      real,intent(in) :: pq(ngrid,nlayer,nq) ! input tracer mmr (kg/kg_air)
+      real,intent(in) :: pdq(ngrid,nlayer,nq) ! input tendency on tracers (kg/kg_air/s)
+      real,intent(in) :: pt(ngrid,nlayer) ! input temperature (K)
+      real,intent(in) :: pdt(ngrid,nlayer) ! input tendency on temperature (K/s)
+      real,intent(out) :: pdqdiff(ngrid,nlayer,nq) ! tendency on tracers (kg/kg_air/s)
+      real*8,intent(out) :: PhiEscH,PhiEscH2,PhiEscD
 !
 ! Local
@@ -255 +252 @@
 !     &  ,tt,ptimestep,nlayer,ig)
         do l=1,nlayer
-          tt(l)=pt(ig,l)*1D0+(pdt(ig,l)*dble(ptimestep)+                &
-                              pdtconduc(ig,l)*dble(ptimestep)+          &
-                              pdteuv(ig,l)*dble(ptimestep))
+          tt(l)=pt(ig,l)*1D0+(pdt(ig,l)*dble(ptimestep))
           ! to cach Nans...
           if (tt(l).ne.tt(l)) then
             print*,'Err TMNEW',ig,l,tt(l),pt(ig,l),                     &
-              pdt(ig,l),pdtconduc(ig,l),pdteuv(ig,l),dble(ptimestep)
+              pdt(ig,l),dble(ptimestep)
           endif
         enddo ! of do l=1,nlayer
@@ -322 +317 @@
         Print*,'Zmax too high',ig,zmax,zmin
         do l=1,nlayer
-        print*,'old',zz(l),pt(ig,l),pdteuv(ig,l),pdq(ig,l,:)
+        print*,'old',zz(l),pt(ig,l),pdt(ig,l),pdq(ig,l,:)
         print*,'l',l,rhot(l),tt(l),pp(l),massemoy(l),qq(l,:)
         enddo