Changeset 3285 for trunk

Timestamp:

Mar 28, 2024, 9:54:07 AM (8 months ago)

Author:

evos

Message:

small changes in the subsufrace lag model.
summing the small subtimestep instead of the ptimestep
and setting zqsurf to 0 if we remove the ice from the surface.
changing the interaction with frost threshold to 2 microns.

File:

• trunk/LMDZ.MARS/libf/phymars/vdifc_mod.F (modified) (5 diffs)

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

@@ -247 +247 @@
       REAL zdqsdif_ssi_tot(ngrid,nslope)       ! Total flux of the interactions with SSI kg/m^2/s^-1)
       
-      REAL :: tol_frost = 1e-4 ! tolerence for frost thicnkess (kg/m^2) to avoid numerical noise effect
+      REAL :: tol_frost = 2e-3 ! tolerence for frost thicnkess (kg/m^2) to avoid numerical noise effect
 c    ** un petit test de coherence
 c       --------------------------
@@ -1120 +1120 @@
      
                     zdqsdif_tot(ig)=-zqsurf(ig)/subtimestep
-!                    zqsurf(ig) = 0
                     if((h2o_ice_depth(ig,islope).gt.0)
      &                  .and.lag_layer) then
                   ! Atm <-> subsurface exchange, we need to update the exchange coefficient zb by a factor 1/1+R; R = zice*Cd,h*/D and add the flux from the subsurface
-                  
+                         zqsurf(ig) = 0
                          if(old_wsublimation_scheme) then
                             resist(ig,islope)=h2o_ice_depth(ig,islope)
@@ -1146 +1145 @@
                          qsat_ssi(ig,islope)=ztsrf(ig)/Tice(ig,islope)
      &                                       *qsat_ssi(ig,islope)
-                  ! Flux from the subsurface      
-                        if(old_wsublimation_scheme) then
-                           zdqsdif_ssi_atm(ig,islope)=rho(ig)
-     &                     *dryness(ig)*zcdv(ig,islope)
-     &                     *1/(1+resist(ig,islope))    
-     &                     *(zq1temp(ig)-qsat_ssi(ig,islope))
-                        else
-                           zdqsdif_ssi_atm(ig,islope)=rho(ig)*
-     &                     *dryness(ig) *zcdh(ig,islope) 
-     &                     *1/(1+resist(ig,islope))   
-     &                     *(zq1temp(ig)-qsat_ssi(ig,islope))                         
-                        endif
+     
                   ! And now we solve correctly the system      
                         z1(ig)=1./(za(ig,1)+zb(ig,1)+
@@ -1267 +1255 @@
               zdqsdif_ssi_atm_tot(ig,islope) = 
      &                              zdqsdif_ssi_atm_tot(ig,islope)
-     &                            + zdqsdif_ssi_atm(ig,islope)
+     &                     + zdqsdif_ssi_atm(ig,islope)*subtimestep
               zdqsdif_ssi_frost_tot(ig,islope) = 
      &                                zdqsdif_ssi_frost_tot(ig,islope) 
-     &                              + zdqsdif_ssi_frost(ig,islope)
+     &                     +zdqsdif_ssi_frost(ig,islope)*subtimestep
 c             Monitoring instantaneous latent heat flux in W.m-2 :
               zsurf_h2o_lh(ig,islope) = zsurf_h2o_lh(ig,islope)+
@@ -1299 +1287 @@
      &                       /ptimestep
      
-            zdqsdif_ssi_tot(ig,islope) = 
+            zdqsdif_ssi_atm_tot(ig,islope)=
+     &           zdqsdif_ssi_atm_tot(ig,islope)/ptimestep
+            zdqsdif_ssi_frost_tot(ig,islope)=
+     &           zdqsdif_ssi_frost_tot(ig,islope)/ptimestep
+                   zdqsdif_ssi_tot(ig,islope) = 
      &                        zdqsdif_ssi_atm_tot(ig,islope) 
      &                      + zdqsdif_ssi_frost_tot(ig,islope)