Ignore:
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:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/LMDZ.MARS/libf/aeronomars/euvheat.F90

    r4023 r4035  
    66
    77      SUBROUTINE euvheat(ngrid,nlayer,nq,pt,pdt,pplev,pplay,zzlay, &
    8            mu0,ptimestep,ptime,zday,pq,pdq,pdteuv)
     8           mu0,ptimestep,zday,pq,pdq,pdteuv)
    99
    1010      use tracer_mod, only: igcm_co2, igcm_co, igcm_o, igcm_o1d,         &
     
    4747      integer,intent(in) :: nlayer ! number of atmospheric layers
    4848      integer,intent(in) :: nq ! number of advected tracers
    49       real :: pt(ngrid,nlayer)
    50       real :: pdt(ngrid,nlayer)
    51       real :: pplev(ngrid,nlayer+1)
    52       real :: pplay(ngrid,nlayer)
    53       real :: zzlay(ngrid,nlayer)
    54       real :: mu0(ngrid)
    55       real :: ptimestep,ptime
    56       real :: zday
    57       real :: pq(ngrid,nlayer,nq)
    58       real :: pdq(ngrid,nlayer,nq)
    59 
    60       real :: pdteuv(ngrid,nlayer)
     49      real,intent(in) :: pt(ngrid,nlayer) ! input temperature (K)
     50      real,intent(in) :: pdt(ngrid,nlayer) ! input tendency on temperatur (K/s)
     51      real,intent(in) :: pplev(ngrid,nlayer+1) ! pressure (Pa) at layer interfaces
     52      real,intent(in) :: pplay(ngrid,nlayer) ! pressure (Pa) at mid-layer
     53      real,intent(in) :: zzlay(ngrid,nlayer) ! mid-layer altitude (m)
     54      real,intent(in) :: mu0(ngrid) ! cosine of solar zenith angle
     55      real,intent(in) :: ptimestep ! physics time step (s)
     56      real,intent(in) :: zday ! date (sols since Ls=0)
     57      real,intent(in) :: pq(ngrid,nlayer,nq) ! input tracer mmr (kg/kg_air)
     58      real,intent(in) :: pdq(ngrid,nlayer,nq) ! input tendencies on tracers (kg/kg_air/s)
     59
     60      real,intent(out) :: pdteuv(ngrid,nlayer) ! temperature tendency (K/s)
    6161!
    6262!    Local variables :
Note: See TracChangeset for help on using the changeset viewer.