Changeset 3091 for trunk/LMDZ.MARS/libf/phymars

Timestamp:

Oct 18, 2023, 3:20:35 PM (15 months ago)

Author:

llange

Message:

MARS PCM
Fixing a bug in co2condens: pdqsc was in the wrong dimension (nslope
missing)
LL & TP

File:

• trunk/LMDZ.MARS/libf/phymars/co2condens_mod.F (modified) (4 diffs)

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

@@ -102 +102 @@
       REAL,INTENT(OUT) :: pdvc(ngrid,nlayer) ! tendency on meridional wind (m.s-2)
       REAL,INTENT(OUT) :: pdqc(ngrid,nlayer,nq) ! tendency on tracers
-      REAL,INTENT(OUT) :: pdqsc(ngrid,nq) ! tendency on surface tracers
+      REAL,INTENT(OUT) :: pdqsc(ngrid,nq,nslope) ! tendency on surface tracers
       
       ! added to calculate flux dependent albedo:
@@ -181 +181 @@
 c Subslope:
 
-      REAL   :: emisref_tmp(ngrid) 
-      REAL   :: alb_tmp(ngrid,2) ! local
-      REAL   :: zcondices_tmp(ngrid)    ! local  
-      REAL   :: piceco2_tmp(ngrid)    ! local
-      REAL   :: perenial_co2ice_tmp(ngrid) ! perenial ice on one subslope (kg/m^2)
-      REAL   :: pemisurf_tmp(ngrid)! local
-      LOGICAL :: condsub_tmp(ngrid) !local
-      REAL :: zfallice_tmp(ngrid,nlayer+1)
-      REAL :: condens_layer_tmp(ngrid,nlayer) ! co2clouds: condensation rate in layer  l  (kg/m2/s)
-      INTEGER :: islope
+      REAL   :: emisref_tmp(ngrid)            ! temporary emissiviy of a sub-grid surface before computing emissivity with snow  [1]
+      REAL   :: pemisurf_tmp(ngrid)           ! temporary emissiviy of a sub-grid surface after computing emissivity with snow   [1]
+      REAL   :: alb_tmp(ngrid,2)              ! temporary albedo before computing albedo with TES/fixed value [1]
+      REAL   :: zcondices_tmp(ngrid)          ! temporary condensation rate [kg/m^2/s]
+      REAL   :: piceco2_tmp(ngrid)            ! temporary amount of CO2 frost on the surface [kg/m^2]
+      REAL   :: perenial_co2ice_tmp(ngrid)    ! temporary amount of perenial CO2 frost on the surface [kg/m^2]
+      LOGICAL :: condsub_tmp(ngrid)           ! Boolean to check if CO2 ice is condensing / sublimating on the sub grid surface [1]
+      REAL :: zfallice_tmp(ngrid,nlayer+1)    ! temporary amount of ice falling from layer l for a specific sub-grid surface [kg/m^2/s]
+      REAL :: condens_layer_tmp(ngrid,nlayer) ! temporary condensation rate in layer l (co2 cloud) for a specific sub-grid surface [kg/m^2/s]
+      INTEGER :: islope                       ! index for loop variables
+      REAL :: pdqsc_tmp(ngrid,nq)         ! tendency on surface tracers (grid-mesh average) after scavenging
 c----------------------------------------------------------------------
 
@@ -248 +249 @@
       pdvc(1:ngrid,1:nlayer)  = 0
       pdtc(1:ngrid,1:nlayer) = 0.
-      pdqsc(1:ngrid,1:nq) = 0
-         
+      pdqsc(1:ngrid,1:nq,1:nslope) = 0
+      pdqsc_tmp(1:ngrid,1:nq) = 0
       pdqc(1:ngrid,1:nlayer,1:nq)  = 0
 
@@ -353 +354 @@
       if (scavco2cond) then
         call scavenging_by_co2(ngrid,nlayer,nq,ptimestep,pplev,zq,
-     &                       rdust,zcondicea,zfallice,zdq_scav,pdqsc)
+     &                       rdust,zcondicea,zfallice,zdq_scav,
+     &                       pdqsc_tmp)
+
+        DO islope = 1,nslope
+              pdqsc(:,:,islope)= pdqsc_tmp(:,:)*
+     &       cos(pi*def_slope_mean(islope)/180.)
+        ENDDO
       endif
              call write_output("co2condens_zfallice",