Changeset 3159 for trunk/LMDZ.MARS

Timestamp:

Dec 14, 2023, 6:56:40 PM (11 months ago)

Author:

jbclement

Message:

PEM:

• A number of corrections related to r3149 for CO2 ice management:

  'co2_ice' was not transferred to 'perennial_co2ice' at the end of the PEM run;
  In the Mars PCM, 'perennial_co2ice' is now correctly handled regarding the frost in case of CO2 sublimation/condensation.
  

• Addition of (CO2 and H2O) ice metamorphism: if the frost is above a given value, then the excess frost is transferred to the perennial ice.
• Thresholds value for ice management can now be set in the "run_PEM.def".

JBC

Location:

trunk/LMDZ.MARS/libf/phymars

Files:

• albedocaps.F90 (modified) (1 diff)
• co2condens_mod.F (modified) (4 diffs)
• physiq_mod.F (modified) (1 diff)

trunk/LMDZ.MARS/libf/phymars/albedocaps.F90

@@ -105 +105 @@
       psolaralb(ig,2)=psolaralb(ig,1)
     else
-      psolaralb(ig,1)=albedice(icap)
-      psolaralb(ig,2)=albedice(icap)
-      if(paleoclimate) then
-         if((piceco2_peren(ig).gt.0.).and.(piceco2(ig).lt.piceco2_peren(ig))) then
-             psolaralb(ig,1) = albedo_perennialco2
-             psolaralb(ig,2) = albedo_perennialco2
-             piceco2_peren(ig) = piceco2(ig)
-         endif
-      endif
+      psolaralb(ig,:)=albedice(icap)
+      if (paleoclimate .and. piceco2_peren(ig) > 0. .and. abs(piceco2(ig)) < 1.e-10) psolaralb(ig,:) = albedo_perennialco2
     endif
   else if (watercaptag(ig) .and. water) then

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

@@ -5 +5 @@
       logical, save :: scavco2cond = .false. ! flag for using scavenging_by_co2
 !$OMP THREADPRIVATE(scavco2cond)
-      
+
       CONTAINS
 
@@ -186 +186 @@
       REAL   :: zcondices_tmp(ngrid)          ! temporary condensation rate [kg/m^2/s]
       REAL   :: piceco2_tmp(ngrid)            ! temporary amount of CO2 frost on the surface [kg/m^2]
-      REAL   :: perennial_co2ice_tmp(ngrid)    ! temporary amount of perennial CO2 frost on the surface [kg/m^2]
+      REAL   :: perennial_co2ice_tmp(ngrid)   ! temporary amount of perennial CO2 frost on the surface [kg/m^2]
       LOGICAL :: condsub_tmp(ngrid)           ! Boolean to check if CO2 ice is condensing / sublimating on the sub grid surface [1]
       REAL :: zfallice_tmp(ngrid,nlayer+1)    ! temporary amount of ice falling from layer l for a specific sub-grid surface [kg/m^2/s]
       REAL :: condens_layer_tmp(ngrid,nlayer) ! temporary condensation rate in layer l (co2 cloud) for a specific sub-grid surface [kg/m^2/s]
       INTEGER :: islope                       ! index for loop variables
-      REAL :: pdqsc_tmp(ngrid,nq)         ! tendency on surface tracers (grid-mesh average) after scavenging
+      REAL :: pdqsc_tmp(ngrid,nq)             ! tendency on surface tracers (grid-mesh average) after scavenging
 c----------------------------------------------------------------------
 
@@ -227 +227 @@
       ENDIF ! of IF (firstcall)
       zcpi=1./cpp
+      if (paleoclimate) piceco2 = piceco2 + perennial_co2ice
 
 c
@@ -564 +565 @@
 c  ********************************************************************
 
+!     Redistribute piceco2 into piceco2 (frost) and perennial_co2ice
+!     --------------------------------------------------------------
+      if (paleoclimate) then
+          where (piceco2 > perennial_co2ice) ! Perennial co2 ice has not been affected
+              ! It means:
+              !     - In case of sublimation, only frost is lost
+              !     - In case of condensation, only frost accumulates new ice
+              piceco2 = piceco2 - perennial_co2ice
+          else where ! Perennial co2 ice has been affected
+              ! It means that frost disappeared with sublimation and perennial ice is being lost
+              perennial_co2ice = piceco2
+              piceco2 = 0.
+          end where
+      endif
+
 !     Check that amont of CO2 ice is not problematic
-      DO ig=1,ngrid
+!     ----------------------------------------------
+      DO ig = 1,ngrid
          DO islope = 1,nslope
-           if(.not.piceco2(ig,islope).ge.0.) THEN
-              if(piceco2(ig,islope).le.-5.e-8) print*,
-     $         'WARNING co2condens piceco2(',ig,')=', piceco2(ig,islope)
-               piceco2(ig,islope)=0.
+           if (piceco2(ig,islope) < 0.) then
+              if (piceco2(ig,islope) <= -5.e-8) print*,
+     $         'WARNING co2condens piceco2(',ig,',',islope,
+     $         ') =',piceco2(ig,islope)
+               piceco2(ig,islope) = 0.
            endif
         ENDDO
       ENDDO
+      if (paleoclimate) then
+      DO ig = 1,ngrid
+         DO islope = 1,nslope
+           if (perennial_co2ice(ig,islope) < 0.) then
+              if (perennial_co2ice(ig,islope) <= -5.e-8) print*,
+     $         'WARNING co2condens perennial_co2ice(',ig,',',islope,
+     $         ') =',perennial_co2ice(ig,islope)
+               perennial_co2ice(ig,islope) = 0.
+           endif
+        ENDDO
+      ENDDO
+      endif
       
 !     Set albedo and emissivity of the surface
 !     ----------------------------------------
       DO islope = 1,nslope
-        piceco2_tmp(:) = piceco2(:,islope)
-        alb_tmp(:,:) = psolaralb(:,:,islope)
-        emisref_tmp(:) = 0.
-        perennial_co2ice_tmp(:) =  perennial_co2ice(:,islope) 
+        piceco2_tmp = piceco2(:,islope)
+        alb_tmp = psolaralb(:,:,islope)
+        emisref_tmp = 0.
+        perennial_co2ice_tmp =  perennial_co2ice(:,islope) 
         CALL albedocaps(zls,ngrid,piceco2_tmp,perennial_co2ice_tmp,
      &                  alb_tmp,emisref_tmp)
-        perennial_co2ice(:,islope) = perennial_co2ice_tmp(:)
+        perennial_co2ice(:,islope) = perennial_co2ice_tmp
         psolaralb(:,1,islope) =  alb_tmp(:,1)
         psolaralb(:,2,islope) =  alb_tmp(:,2)
-        emisref(:,islope) = emisref_tmp(:)
+        emisref(:,islope) = emisref_tmp
       ENDDO
 
-! set pemisurf() to emissiv when there is bare surface (needed for co2snow)
-      DO ig=1,ngrid
-        DO islope = 1,nslope
-          if (piceco2(ig,islope).eq.0) then
-            pemisurf(ig,islope)=emissiv
-          endif
-        ENDDO
-      ENDDO
+!     Set pemisurf() to emissiv when there is bare surface (needed for co2snow)
+!     -------------------------------------------------------------------------
+      where (piceco2 == 0.) pemisurf = emissiv
+
 
 !         firstcall2=.false.

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

@@ -2473 +2473 @@
       IF (refill_watercap) THEN
 
-        DO ig=1,ngrid
+        DO ig = 1,ngrid
          DO islope = 1,nslope
-           if (watercaptag(ig).and. (qsurf(ig,igcm_h2o_ice,islope)
-     &        .gt.frost_metam_threshold)) then
-
-                watercap(ig,islope)=watercap(ig,islope)
-     &          +qsurf(ig,igcm_h2o_ice,islope)
-     &         - frost_metam_threshold
+           if (watercaptag(ig) .and. (qsurf(ig,igcm_h2o_ice,islope)
+     &        > frost_metam_threshold)) then
+
+                watercap(ig,islope) = watercap(ig,islope)
+     &          + qsurf(ig,igcm_h2o_ice,islope)
+     &          - frost_metam_threshold
                 qsurf(ig,igcm_h2o_ice,islope) = frost_metam_threshold
-           endif ! (watercaptag(ig).and.
+           endif ! watercaptag
           ENDDO
         ENDDO
 
-      ENDIF ! (refill_watercap) THEN
+      ENDIF ! refill_watercap