Changeset 2999 for trunk/LMDZ.MARS/libf/phymars/co2condens_mod.F

Timestamp:

Jul 19, 2023, 11:40:38 AM (19 months ago)

Author:

llange

Message:

Mars PCM
Include perenial_co2ice (equivalent of watercap) to distinguich between CO2 frost and perenial CO2 ice for paleoclimate studies.
When no frost is present and we dig into perenial ice, the surface albedo is changed. The albedo for seasonal ice is set to 0.65, and the perenial ice albedo can be fixed in the callphys.def. I recommand values between 0.8 and 0.9.
To use this, paleoclimate must be set to True and TESalbedo to false in the callphys.def. Else, it runs as usual with TES albedo
LL

File:

• trunk/LMDZ.MARS/libf/phymars/co2condens_mod.F (modified) (6 diffs)

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

@@ -11 +11 @@
      $                  pcapcal,pplay,pplev,ptsrf,pt,
      $                  pphi,pdt,pdu,pdv,pdtsrf,pu,pv,pq,pdq,
-     $                  piceco2,psolaralb,pemisurf,rdust,
+     $                  piceco2,perenial_co2ice,
+     $                  psolaralb,pemisurf,rdust,
      $                  pdtc,pdtsrfc,pdpsrf,pduc,pdvc,pdqc,
      $                  fluxsurf_sw,zls,
@@ -35 +36 @@
 #endif
       use comslope_mod, ONLY: subslope_dist,def_slope_mean
+      USE paleoclimate_mod, ONLY: paleoclimate
+
        IMPLICIT NONE
 c=======================================================================
@@ -87 +90 @@
 
       REAL,INTENT(INOUT) :: piceco2(ngrid,nslope) ! CO2 ice on the surface (kg.m-2)
+      REAL,INTENT(INOUT) :: perenial_co2ice(ngrid,nslope) ! Perenial CO2 ice on the surface (kg.m-2)
       REAL,INTENT(INOUT) :: psolaralb(ngrid,2,nslope) ! albedo of the surface
       REAL,INTENT(INOUT) :: pemisurf(ngrid,nslope) ! emissivity of the surface
@@ -180 +184 @@
       REAL   :: alb_tmp(ngrid,2) ! local
       REAL   :: zcondices_tmp(ngrid)    ! local  
-      REAL   :: piceco2_tmp(ngrid)    ! local  
+      REAL   :: piceco2_tmp(ngrid)    ! local
+      REAL   :: perenial_co2ice_tmp(ngrid) ! perenial ice on one subslope (kg/m^2)
       REAL   :: pemisurf_tmp(ngrid)! local
       LOGICAL :: condsub_tmp(ngrid) !local
@@ -568 +573 @@
         piceco2_tmp(:) = piceco2(:,islope)
         alb_tmp(:,:) = psolaralb(:,:,islope)
-        emisref_tmp(:) = 0. 
-        CALL albedocaps(zls,ngrid,piceco2_tmp,alb_tmp,emisref_tmp)
+        emisref_tmp(:) = 0.
+        perenial_co2ice_tmp(:) =  perenial_co2ice(:,islope) 
+        CALL albedocaps(zls,ngrid,piceco2_tmp,perenial_co2ice_tmp,
+     &                  alb_tmp,emisref_tmp)
+        perenial_co2ice(:,islope) = perenial_co2ice_tmp(:)
         psolaralb(:,1,islope) =  alb_tmp(:,1)
         psolaralb(:,2,islope) =  alb_tmp(:,2)
@@ -823 +831 @@
 ! Extra special case for surface temperature tendency pdtsrfc:
 ! we want to fix the south pole temperature to CO2 condensation temperature
-         if (caps.and.(obliquit.lt.27.)) then
+         if (caps.and.(obliquit.lt.27.).and.(.not.(paleoclimate))) then
            ! check if last grid point is the south pole
            if (abs(latitude(ngrid)-(-pi/2.)).lt.1.e-5) then