Changeset 3273 for trunk/LMDZ.MARS/libf

Timestamp:

Mar 19, 2024, 11:06:57 AM (10 months ago)

Author:

llange

Message:

Mars PCM
Small fixes in vdifc for the interaction between the subsurface and frost (the flux from the surface frost to the subsurface was oo large if one removes all the frost at the surface)
LL

File:

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

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

@@ -1169 +1169 @@
      &                              *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
 
                     else 
@@ -1214 +1225 @@
      &                                 *(qsat(ig)-qsat_ssi(ig,islope))
 
-! Line to comment for now since we don't have a mass of subsurface frost in our computation (otherwise, we would not conserve the H2O mass in the system)               
+! Line to comment for now since we don't have a mass of subsurface frost
+! in our computation (otherwise, we would not conserve the H2O mass in
+! the system)     
+
+
+                  if ((-zdqsdif_tot(ig)+zdqsdif_ssi_frost(ig,islope))
+     &            *subtimestep.ge.(zqsurf(ig))) then
+                     zdqsdif_ssi_frost(ig,islope) =
+     &                         zqsurf(ig)/subtimestep
+     &                         +   zdqsdif_tot(ig)
+                  endif
+
                     zdqsdif_tot(ig) = zdqsdif_tot(ig) - 
      &                        zdqsdif_ssi_frost(ig,islope)