Changeset 3030 for trunk/LMDZ.COMMON/libf/evolution

Timestamp:

Aug 24, 2023, 11:40:30 AM (16 months ago)

Author:

llange

Message:

MARS PEM

• Following -r3028, added the declaration of daysec and dtphys that has been removed during the cleaning
• Add the corrective factor /cos(slope) when compute H2O ice tendencies to ensure mass conservation

LL

Location:

trunk/LMDZ.COMMON/libf/evolution

Files:

• evol_h2o_ice_s_mod.F90 (modified) (3 diffs)
• pem.F90 (modified) (1 diff)

trunk/LMDZ.COMMON/libf/evolution/evol_h2o_ice_s_mod.F90

@@ -8 +8 @@
                              iim_input,jjm_input,ngrid,cell_area,STOPPING,nslope)
 
-  USE temps_mod_evol, ONLY: dt_pem
-  use comslope_mod, ONLY: subslope_dist
+  use temps_mod_evol, only: dt_pem
+  use comslope_mod, only: subslope_dist,def_slope_mean
   use criterion_pem_stop_mod, only: criterion_waterice_stop
+  use comconst_mod,only: pi
 
   IMPLICIT NONE
@@ -39 +40 @@
   REAL :: pos_tend, neg_tend, real_coefficient,negative_part    ! Variable to conserve water
   REAL ::  new_tendencies(ngrid,nslope)                         ! Tendencies computed in order to conserve water ice on the surface, only exchange between surface are done
-
+  
 !=======================================================================
 
@@ -52 +53 @@
      if (qsurf(i,islope).GT.0) then
          if (tendencies_h2o_ice_phys(i,islope).GT.0) then
-            pos_tend=pos_tend+tendencies_h2o_ice_phys(i,islope)*cell_area(i)*subslope_dist(i,islope)
+            pos_tend=pos_tend+tendencies_h2o_ice_phys(i,islope)*cell_area(i)*subslope_dist(i,islope)/cos(def_slope_mean(islope)*pi/180.)
          else
-            neg_tend=neg_tend-tendencies_h2o_ice_phys(i,islope)*cell_area(i)*subslope_dist(i,islope)
+            neg_tend=neg_tend-tendencies_h2o_ice_phys(i,islope)*cell_area(i)*subslope_dist(i,islope)/cos(def_slope_mean(islope)*pi/180.)
          endif
      endif

trunk/LMDZ.COMMON/libf/evolution/pem.F90

@@ -86 +86 @@
     use control_mod,              only: iphysiq, day_step, nsplit_phys
 #else
-    use time_phylmdz_mod,         only: daysec, iphysiq, day_step, dtphys
+    use time_phylmdz_mod,         only:  iphysiq, day_step
 #endif
-
+use time_phylmdz_mod, only: daysec,dtphys
 #ifdef CPP_1D
     use regular_lonlat_mod,       only: init_regular_lonlat