Changeset 3030


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

Legend:

Unmodified
Added
Removed
  • trunk/LMDZ.COMMON/libf/evolution/evol_h2o_ice_s_mod.F90

    r2980 r3030  
    88                             iim_input,jjm_input,ngrid,cell_area,STOPPING,nslope)
    99
    10   USE temps_mod_evol, ONLY: dt_pem
    11   use comslope_mod, ONLY: subslope_dist
     10  use temps_mod_evol, only: dt_pem
     11  use comslope_mod, only: subslope_dist,def_slope_mean
    1212  use criterion_pem_stop_mod, only: criterion_waterice_stop
     13  use comconst_mod,only: pi
    1314
    1415  IMPLICIT NONE
     
    3940  REAL :: pos_tend, neg_tend, real_coefficient,negative_part    ! Variable to conserve water
    4041  REAL ::  new_tendencies(ngrid,nslope)                         ! Tendencies computed in order to conserve water ice on the surface, only exchange between surface are done
    41 
     42 
    4243!=======================================================================
    4344
     
    5253     if (qsurf(i,islope).GT.0) then
    5354         if (tendencies_h2o_ice_phys(i,islope).GT.0) then
    54             pos_tend=pos_tend+tendencies_h2o_ice_phys(i,islope)*cell_area(i)*subslope_dist(i,islope)
     55            pos_tend=pos_tend+tendencies_h2o_ice_phys(i,islope)*cell_area(i)*subslope_dist(i,islope)/cos(def_slope_mean(islope)*pi/180.)
    5556         else
    56             neg_tend=neg_tend-tendencies_h2o_ice_phys(i,islope)*cell_area(i)*subslope_dist(i,islope)
     57            neg_tend=neg_tend-tendencies_h2o_ice_phys(i,islope)*cell_area(i)*subslope_dist(i,islope)/cos(def_slope_mean(islope)*pi/180.)
    5758         endif
    5859     endif
  • trunk/LMDZ.COMMON/libf/evolution/pem.F90

    r3028 r3030  
    8686    use control_mod,              only: iphysiq, day_step, nsplit_phys
    8787#else
    88     use time_phylmdz_mod,         only: daysec, iphysiq, day_step, dtphys
     88    use time_phylmdz_mod,         only:  iphysiq, day_step
    8989#endif
    90 
     90use time_phylmdz_mod, only: daysec,dtphys
    9191#ifdef CPP_1D
    9292    use regular_lonlat_mod,       only: init_regular_lonlat
Note: See TracChangeset for help on using the changeset viewer.