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

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

@@ -10 +10 @@
       subroutine thermosphere(ngrid,nlayer,nq,
      &     pplev,pplay,dist_sol,
-     $     mu0,ptimestep,ptime,zday,tsurf,zzlev,zzlay,
-     &     pt,pq,pu,pv,pdt,pdq,
+     $     mu0,ptimestep,zday,tsurf,zzlev,zzlay,
+     &     pt,pq,pu,pv,pdt,pdq,pdu,pdv,
      $     zdteuv,zdtconduc,zdumolvis,zdvmolvis,zdqmoldiff,
      $     PhiEscH,PhiEscH2,PhiEscD)
@@ -31 +31 @@
       integer,intent(in) :: nlayer ! number of atmospheric layers
       integer,intent(in) :: nq ! number of advected tracers
-      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)
-      REAL,INTENT(in) :: zday
-      REAL,INTENT(in) :: dist_sol
-      REAL,INTENT(in) :: mu0(ngrid)
-      REAL,INTENT(in) :: pq(ngrid,nlayer,nq)
-      REAL,INTENT(in) :: ptimestep
-      REAL,INTENT(in) :: ptime
-      REAL,INTENT(in) :: tsurf(ngrid)
-      REAL,INTENT(in) :: pu(ngrid,nlayer),pv(ngrid,nlayer)
-      REAL,INTENT(in) :: pdt(ngrid,nlayer),pdq(ngrid,nlayer,nq)
+      REAL,INTENT(in) :: pplay(ngrid,nlayer) ! mid-layer pressure (Pa)
+      REAL,INTENT(in) :: pplev(ngrid,nlayer+1) ! inter-layer pressure (Pa)
+      REAL,INTENT(in) :: zzlay(ngrid,nlayer) ! altitude of mid-layer (m)
+      REAL,INTENT(in) :: zzlev(ngrid,nlayer+1) ! altitude of layer boundaries (m)
+      REAL,INTENT(in) :: pt(ngrid,nlayer) ! input temperature (K)
+      REAL,INTENT(in) :: zday ! Mars date (sols since Ls=0)
+      REAL,INTENT(in) :: dist_sol ! Sun-Mars distance (AU)
+      REAL,INTENT(in) :: mu0(ngrid) ! cosine of solar zenith angle
+      REAL,INTENT(in) :: pq(ngrid,nlayer,nq) ! input tracer mixing ratio (kg/kg_air)
+      REAL,INTENT(in) :: ptimestep ! physics time step (s)
+      REAL,INTENT(in) :: tsurf(ngrid) ! surface temperature (K)
+      REAL,INTENT(in) :: pu(ngrid,nlayer) ! input zonal wind (m/s)
+      REAL,INTENT(in) :: pv(ngrid,nlayer) ! input meridional wind (m/s)
 
-      REAL,INTENT(out) :: zdteuv(ngrid,nlayer)
-      REAL,INTENT(out) :: zdtconduc(ngrid,nlayer)
-      REAL,INTENT(out) :: zdumolvis(ngrid,nlayer)
-      REAL,INTENT(out) :: zdvmolvis(ngrid,nlayer)
-      REAL,INTENT(out) :: zdqmoldiff(ngrid,nlayer,nq)
-      REAL*8,INTENT(out) :: PhiEscH,PhiEscH2,PhiEscD
+      REAL,INTENT(out) :: zdteuv(ngrid,nlayer) ! tendency due to EUV (K/s)
+      REAL,INTENT(out) :: zdtconduc(ngrid,nlayer) ! tendency due to conduction (K/s)
+      REAL,INTENT(out) :: zdumolvis(ngrid,nlayer) ! tendency due to viscosity (m/s/s)
+      REAL,INTENT(out) :: zdvmolvis(ngrid,nlayer) ! tendency due to viscosity (m/s/s)
+      REAL,INTENT(out) :: zdqmoldiff(ngrid,nlayer,nq) ! tendency due to diffusion (kg/kg_air/s)
+      REAL*8,INTENT(out) :: PhiEscH ! total escape flux of H (s-1)
+      REAL*8,INTENT(out) :: PhiEscH2 ! total escape flux of H2 (s-1)
+      REAL*8,INTENT(out) :: PhiEscD ! total escape flux of D (s-1)
+
+      REAL,INTENT(inout) :: pdt(ngrid,nlayer) ! tendency on temperature (K/s)
+      REAL,INTENT(inout) :: pdq(ngrid,nlayer,nq) ! tendency on tracers (kg/kg_air/s)
+      REAL,INTENT(inout) :: pdu(ngrid,nlayer) ! tendency on zonal wind (m/s/s)
+      REAL,INTENT(inout) :: pdv(ngrid,nlayer) ! tendency on meridional wind (m/s/s)
 
       INTEGER :: l,ig
@@ -68 +74 @@
 
       ! initialize tendencies to zero in all cases
-      ! (tendencies are added later on, even if parametrization is not called)
       zdteuv(1:ngrid,1:nlayer)=0
       zdtconduc(1:ngrid,1:nlayer)=0
@@ -77 +82 @@
       if (calleuv) then
         call euvheat(ngrid,nlayer,nq,pt,pdt,pplev,pplay,zzlay,
-     $               mu0,ptimestep,ptime,zday,pq,pdq,zdteuv)
+     $               mu0,ptimestep,zday,pq,pdq,zdteuv)
+        ! update tendency on temperature
+        pdt(:,:)=pdt(:,:)+zdteuv(:,:)
       endif
 
       if (callconduct) THEN
-        call conduction(ngrid,nlayer,ptimestep,pplay,pplev,pt,zdteuv,
+        call conduction(ngrid,nlayer,ptimestep,pplay,pplev,pt,pdt,
      $                   tsurf,zzlev,zzlay,zdtconduc)
+        ! update tendency on temperature
+        pdt(:,:)=pdt(:,:)+zdtconduc(:,:)
       endif
 
       if (callmolvis) THEN
         call molvis(ngrid,nlayer,ptimestep,pplay,pplev,pt,
-     &                zdteuv,zdtconduc,pu,
-     $                   tsurf,zzlev,zzlay,zdumolvis)
+     &                pdt,pu,tsurf,zzlev,zzlay,zdumolvis)
         call molvis(ngrid,nlayer,ptimestep,pplay,pplev,pt,
-     &                zdteuv,zdtconduc,pv,
-     $                   tsurf,zzlev,zzlay,zdvmolvis)
+     &                pdt,pv,tsurf,zzlev,zzlay,zdvmolvis)
+        ! update tendencies on winds
+        pdu(:,:)=pdu(:,:)+zdumolvis(:,:)
+        pdv(:,:)=pdv(:,:)+zdvmolvis(:,:)
       endif
 
@@ -99 +109 @@
         call moldiff_red(ngrid,nlayer,nq,
      &                   pplay,pplev,pt,pdt,pq,pdq,ptimestep,
-     &                   zzlay,zdteuv,zdtconduc,zdqmoldiff,
-     &                   PhiEscH,PhiEscH2,PhiEscD)
+     &                   zdqmoldiff,PhiEscH,PhiEscH2,PhiEscD)
        else
         ! new MPF (modified pass flow) scheme
         call moldiff_MPF(ngrid,nlayer,nq,
      &                   pplay,pplev,pt,pdt,pq,pdq,ptimestep,
-     &                   zzlay,zdteuv,zdtconduc,zdqmoldiff,
-     &                   PhiEscH,PhiEscH2,PhiEscD)
+     &                   zdqmoldiff,PhiEscH,PhiEscH2,PhiEscD)
        endif ! of if (moldiff_scheme==1)
-      endif
+
+       ! update tendencies on tracers
+       pdq(:,:,:)=pdq(:,:,:)+zdqmoldiff(:,:,:)
+       
+      endif ! of if (callmoldiff)
 
       end subroutine thermosphere