Changeset 3400 for trunk/LMDZ.MARS/libf

Timestamp:

Aug 6, 2024, 6:39:29 PM (5 months ago)

Author:

llange

Message:

Mars PCM
Back to the commit
Correction of the computation of the water vapor flux from the subsurface ice when thin frost is at the surface: When the frost was too thin (thickness < tol_frost), the flux from the subsurface ice to the surface/atmosphere was not computed. It is now corrected, assuming that the frost is too thin to prevent exchange between subsurface ice and the atmosphere
cleaning testphys1d (unused comments)
LL, who should be writing his PhD dissertation

Location:

trunk/LMDZ.MARS/libf/phymars

Files:

• dyn1d/testphys1d.F90 (modified) (1 diff)
• vdifc_mod.F (modified) (4 diffs)

trunk/LMDZ.MARS/libf/phymars/dyn1d/testphys1d.F90

@@ -207 +207 @@
             endif
         endif
-    !EV profile
-!            IF(atm_wat_profile.eq.2) THEN
-!             DO ilayer=1,16
-!              q(1,ilayer,igcm_h2o_vap)=min(zqsat(ilayer),6e-5)
-!             ENDDO! ilayer=1,16
-!             DO ilayer=17,22
-!              q(1,ilayer,igcm_h2o_vap)=min(zqsat(ilayer),3e-5)
-!             ENDDO! ilayer=17,22
-!             DO ilayer=23,nlayer
-!              q(1,ilayer,igcm_h2o_vap)=0
-!             ENDDO! ilayer=23,nlayer
-!             endif
      endif
 

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

@@ -1113 +1113 @@
 
               if(.not.watercaptag(ig).and.(.not.adsorption_soil)) then
-                 if ((-zdqsdif_tot(ig)*subtimestep) 
-     &                .ge.(zqsurf(ig))) then
+                 if ((((-zdqsdif_tot(ig)*subtimestep) 
+     &                .ge.(zqsurf(ig)))) 
+     &              .or.
+     &                (((-zdqsdif_tot(ig)*subtimestep)
+     &                .lt.(zqsurf(ig)))
+     &                .and. ((zqsurf(ig).lt.tol_frost))
+     &                .and.lag_layer)) then
+
      
-                    zdqsdif_tot(ig)=-zqsurf(ig)/subtimestep
+!                     zdqsdif_tot(ig)=-zqsurf(ig)/subtimestep
                     if((h2o_ice_depth(ig,islope).gt.0)
-     &                  .and. (zqsurf(ig) .lt. 1e-9)
+     &                  .and. (zqsurf(ig) .lt. tol_frost)
      &                  .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
+!                         zqsurf(ig) = 0
                          if(old_wsublimation_scheme) then
                             resist(ig,islope)=h2o_ice_depth(ig,islope)
@@ -1214 +1220 @@
 ! 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
@@ -1248 +1252 @@
               zqsurf(ig)= zqsurf(ig)+(
      &                       zdqsdif_tot(ig))*subtimestep
+
               if (zqsurf(ig)<0 .and.(.not.watercaptag(ig))) then
                       zqsurf(ig)=0
@@ -1412 +1417 @@
      &                     "Ground ice latent heat flux",
      &                     "W.m-2",surf_h2o_lh(:,iflat))
-         call write_output("nsubtimestep",
-     &                     "Number of subtimesteps in vdifc",
-     &                     "int",nsubtimestep(:))
 
          call write_output('zdqsdif_ssi_frost_tot',