Changeset 3130 for trunk/LMDZ.MARS/libf/phymars

Timestamp:

Nov 20, 2023, 1:25:31 PM (20 months ago)

Author:

jbclement

Message:

PEM:
The perennial co2 ice is now taken into account with co2 frost (qsurf) to compute the tendency and to make the update + Rework of how co2 frost is converted to perennial co2 ice at the end of the PEM run + Correction of the value of 'threshold_co2_frost2perennial' to correspond to 10 m + Perennial co2 ice is now handled outside 'paleoclimate' in "phyetat0_mod.F90" of the Mars PCM + Some cleanings.

/!\ Commit for the PEM management of co2 ice before a rework of ice management in the PEM!
JBC

Location:

trunk/LMDZ.MARS/libf/phymars

Files:

• albedocaps.F90 (modified) (3 diffs)
• co2condens_mod.F (modified) (4 diffs)
• conf_phys.F (modified) (2 diffs)
• dyn1d/testphys1d.F90 (modified) (2 diffs)
• paleoclimate_mod.F90 (modified) (1 diff)
• phyetat0_mod.F90 (modified) (5 diffs)
• phyredem.F90 (modified) (3 diffs)
• physiq_mod.F (modified) (5 diffs)
• surfdat_h.F90 (modified) (5 diffs)

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

@@ -11 +11 @@
 USE mod_phys_lmdz_transfert_para, ONLY: bcast
 USE mod_phys_lmdz_para, ONLY: is_master
-USE paleoclimate_mod, ONLY: paleoclimate,albedo_perenialco2
+USE paleoclimate_mod, ONLY: paleoclimate,albedo_perennialco2
 
 implicit none
@@ -21 +21 @@
 integer,intent(in) :: ngrid
 real,intent(in) :: piceco2(ngrid) ! amount of CO2 ice on the surface (kg/m2)
-real,intent(inout) :: piceco2_peren(ngrid) ! amount of perenial co2 ice (kg/m^2)
+real,intent(inout) :: piceco2_peren(ngrid) ! amount of perennial co2 ice (kg/m^2)
 real,intent(out) :: psolaralb(ngrid,2) ! albedo of the surface
 real,intent(out) :: emisref(ngrid) ! emissivity of the surface
@@ -109 +109 @@
       if(paleoclimate) then
          if((piceco2_peren(ig).gt.0.).and.(piceco2(ig).lt.piceco2_peren(ig))) then
-             psolaralb(ig,1) = albedo_perenialco2
-             psolaralb(ig,2) = albedo_perenialco2
+             psolaralb(ig,1) = albedo_perennialco2
+             psolaralb(ig,2) = albedo_perennialco2
              piceco2_peren(ig) = piceco2(ig)
          endif

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,perenial_co2ice,
+     $                  piceco2,perennial_co2ice,
      $                  psolaralb,pemisurf,rdust,
      $                  pdtc,pdtsrfc,pdpsrf,pduc,pdvc,pdqc,
@@ -90 +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) :: perennial_co2ice(ngrid,nslope) ! perennial 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
@@ -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   :: perenial_co2ice_tmp(ngrid)    ! temporary amount of perenial 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]
@@ -581 +581 @@
         alb_tmp(:,:) = psolaralb(:,:,islope)
         emisref_tmp(:) = 0.
-        perenial_co2ice_tmp(:) =  perenial_co2ice(:,islope) 
-        CALL albedocaps(zls,ngrid,piceco2_tmp,perenial_co2ice_tmp,
+        perennial_co2ice_tmp(:) =  perennial_co2ice(:,islope) 
+        CALL albedocaps(zls,ngrid,piceco2_tmp,perennial_co2ice_tmp,
      &                  alb_tmp,emisref_tmp)
-        perenial_co2ice(:,islope) = perenial_co2ice_tmp(:)
+        perennial_co2ice(:,islope) = perennial_co2ice_tmp(:)
         psolaralb(:,1,islope) =  alb_tmp(:,1)
         psolaralb(:,2,islope) =  alb_tmp(:,2)

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

@@ -37 +37 @@
       use aeropacity_mod, only: iddist, topdustref
       USE mod_phys_lmdz_transfert_para, ONLY: bcast
-      USE paleoclimate_mod,ONLY: paleoclimate,albedo_perenialco2,
+      USE paleoclimate_mod,ONLY: paleoclimate,albedo_perennialco2,
      &                           lag_layer 
       use microphys_h, only: mteta
@@ -351 +351 @@
 
 
-         write(*,*)"Albedo for perenial CO2 ice?"
-         albedo_perenialco2 = 0.85 ! default value
-         call getin_p("albedo_perenialco2",albedo_perenialco2)
-         write(*,*)"albedo_perenialco2 = ",albedo_perenialco2
+         write(*,*)"Albedo for perennial CO2 ice?"
+         albedo_perennialco2 = 0.85 ! default value
+         call getin_p("albedo_perennialco2",albedo_perennialco2)
+         write(*,*)"albedo_perennialco2 = ",albedo_perennialco2
 
 ! TRACERS:

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

@@ -2 +2 @@
 
 use comsoil_h,           only: inertiedat, inertiesoil, nsoilmx, tsoil, nqsoil, qsoil
-use surfdat_h,           only: albedodat, perenial_co2ice, watercap, tsurf, emis, qsurf
+use surfdat_h,           only: albedodat, perennial_co2ice, watercap, tsurf, emis, qsurf
 use comslope_mod,        only: def_slope, subslope_dist
 use phyredem,            only: physdem0, physdem1
@@ -161 +161 @@
     call physdem1("startfi.nc",nsoilmx,ngrid,llm,nq,nqsoil,dttestphys,time,      &
                   tsurf,tsoil,inertiesoil,albedo,emis,q2,qsurf,qsoil,tauscaling, &
-                  totcloudfrac,wstar,watercap,perenial_co2ice)
+                  totcloudfrac,wstar,watercap,perennial_co2ice)
 endif !(.not. therestartfi)
 

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

@@ -17 +17 @@
     real, save, allocatable :: h2o_ice_depth(:,:)  ! Thickness of the lag before H2O ice [m]
     real, save, allocatable :: lag_co2_ice(:,:)  ! Thickness of the lag before CO2 ice [m]
-    real, save :: albedo_perenialco2             ! Albedo for perenial co2 ice [1]
+    real, save :: albedo_perennialco2             ! Albedo for perennial co2 ice [1]
     real, save, allocatable :: d_coef(:,:)  ! Diffusion coeficent
     LOGICAL,SAVE :: lag_layer ! does lag layer is present?
-!$OMP THREADPRIVATE(h2o_ice_depth,d_coef,lag_co2_ice,albedo_perenialco2)
+!$OMP THREADPRIVATE(h2o_ice_depth,d_coef,lag_co2_ice,albedo_perennialco2)
 
     CONTAINS

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

@@ -10 +10 @@
 subroutine phyetat0 (fichnom,tab0,Lmodif,nsoil,ngrid,nlay,nq,nqsoil, &
                      day_ini,time0,tsurf,tsoil,albedo,emis,q2,qsurf,qsoil, &
-                     tauscaling,totcloudfrac,wstar,watercap,perenial_co2ice, &
+                     tauscaling,totcloudfrac,wstar,watercap,perennial_co2ice, &
                      def_slope,def_slope_mean,subslope_dist)
 
@@ -75 +75 @@
   real,intent(out) :: wstar(ngrid) ! Max vertical velocity in thermals (m/s)
   real,intent(out) :: watercap(ngrid,nslope) ! h2o_ice_cover
-  real,intent(out) :: perenial_co2ice(ngrid,nslope) ! perenial co2 ice(kg/m^2)
+  real,intent(out) :: perennial_co2ice(ngrid,nslope) ! perennial co2 ice(kg/m^2)
   real,intent(out) :: def_slope(nslope+1) !boundaries for bining of the slopes
   real,intent(out) :: def_slope_mean(nslope)
@@ -753 +753 @@
    endif
 
-! Diffuesion coeficent 
+! Diffusion coeficent 
    call get_field("d_coef",d_coef,found,indextime)
    if (.not.found) then
@@ -761 +761 @@
    endif
 
-
-!  Depth of CO2 lag
+! Depth of CO2 lag
    call get_field("lag_co2_ice",lag_co2_ice,found,indextime)
    if (.not.found) then
@@ -769 +768 @@
      lag_co2_ice(:,:) = -1.
    endif
-
-! Perenial CO2 ice   
-   call get_field("perenial_co2ice",perenial_co2ice,found,indextime)
-   if (.not.found) then
-     write(*,*) "phyetat0: Failed loading <perenial_co2ice> : ", &
-                "<perenial_co2ice> is set as 10m at the South Pole"
-     perenial_co2ice(:,:) = 0.
-     if (abs(latitude(ngrid)-(-pi/2.)).lt.1.e-5) then
-       DO islope = 1,nslope
-        perenial_co2ice(ngrid,islope) = 10*1.6e3 ! 10m which is convert to kg/m^2
-        qsurf(ngrid,igcm_co2_tmp,islope) = qsurf(ngrid-1,igcm_co2_tmp,islope) + & ! perenial ice + frost
-                                        perenial_co2ice(ngrid,islope)
-       ENDDO
-     endif
-   endif ! notfound
   else ! no startfiphyle
      h2o_ice_depth(:,:) = -1.
      lag_co2_ice(:,:) = -1.
-     perenial_co2ice(:,:) = 0.
      d_coef(:,:)= 4.e-4
-     if (abs(latitude(ngrid)-(-pi/2.)).lt.1.e-5) then
-       DO islope = 1,nslope
-        perenial_co2ice(ngrid,islope) = 10*1.6e3 ! 10m which is convert to kg/m^2
-        qsurf(ngrid,igcm_co2_tmp,islope) = qsurf(ngrid-1,igcm_co2_tmp,islope) + & ! perenial ice + frost
-                                        perenial_co2ice(ngrid,islope)
-       ENDDO
-     endif
   endif !startphy_file
 else
   h2o_ice_depth(:,:) = -1.
   lag_co2_ice(:,:) = -1.
-  perenial_co2ice(:,:) = 0.
   d_coef(:,:)= 4.e-4
 endif !paleoclimate
+
+! Perennial CO2 ice
+perennial_co2ice(:,:) = 0.
+if (startphy_file) then
+    call get_field("perennial_co2ice",perennial_co2ice,found,indextime)
+    if (.not. found) then
+        write(*,*) "phyetat0: Failed loading <perennial_co2ice> : ", &
+                   "<perennial_co2ice> is set as 10m at the South Pole"
+        if (abs(latitude(ngrid) - (-pi/2.)) < 1.e-5) then
+            do islope = 1,nslope
+                perennial_co2ice(ngrid,islope) = 10*1.6e3 ! 10m which is convert to kg/m^2
+               qsurf(ngrid,igcm_co2_tmp,islope) = qsurf(ngrid - 1,igcm_co2_tmp,islope) + perennial_co2ice(ngrid,islope) ! perennial ice + frost
+            enddo
+        endif
+    endif ! not found
+else
+    if (abs(latitude(ngrid) - (-pi/2.)) < 1.e-5) then
+        do islope = 1,nslope
+            perennial_co2ice(ngrid,islope) = 10*1.6e3 ! 10m which is convert to kg/m^2
+            qsurf(ngrid,igcm_co2_tmp,islope) = qsurf(ngrid - 1,igcm_co2_tmp,islope) + perennial_co2ice(ngrid,islope) ! perennial ice + frost
+        enddo
+    endif
+endif !startphy_file
 
 ! close file:

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

@@ -168 +168 @@
                     albedo,emis,q2,qsurf,qsoil,&
                     tauscaling,totcloudfrac,wstar, &
-                    watercap,perenial_co2ice)
+                    watercap,perennial_co2ice)
   ! write time-dependent variable to restart file
   use iostart, only : open_restartphy, close_restartphy, & 
@@ -205 +205 @@
   real,intent(in) :: wstar(ngrid)
   real,intent(in) :: watercap(ngrid,nslope)
-  real,intent(in) :: perenial_co2ice(ngrid,nslope)
+  real,intent(in) :: perennial_co2ice(ngrid,nslope)
 
   integer :: iq
@@ -226 +226 @@
   call put_field("watercap","H2O ice cover",watercap,time)
  
- ! Perenial CO2 ice layer
-  if (paleoclimate) call put_field("perenial_co2ice","CO2 ice cover",perenial_co2ice,time)
+  ! Perennial CO2 ice layer
+  call put_field("perennial_co2ice","CO2 ice cover",perennial_co2ice,time)
 
   ! Surface temperature

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

@@ -59 +59 @@
      &                     capcal, fluxgrd, qsurf,
      &                     hmons,summit,base,watercap,watercaptag,
-     &                     perenial_co2ice
+     &                     perennial_co2ice
       use comsaison_h, only: dist_sol, declin, zls, 
      &                       mu0, fract, local_time
@@ -604 +604 @@
      &         tsurf,tsoil,albedo,emis,
      &         q2,qsurf,qsoil,tauscaling,totcloudfrac,wstar,
-     &         watercap,perenial_co2ice,
+     &         watercap,perennial_co2ice,
      &         def_slope,def_slope_mean,subslope_dist)
 
@@ -2209 +2209 @@
      $              capcal,zplay,zplev,tsurf,pt,
      $              pphi,pdt,pdu,pdv,zdtsurf,pu,pv,pq,pdq,
-     $              qsurf(:,igcm_co2,:),perenial_co2ice,
+     $              qsurf(:,igcm_co2,:),perennial_co2ice,
      $              albedo,emis,rdust,
      $              zdtc,zdtsurfc,pdpsrf,zduc,zdvc,zdqc,
@@ -2625 +2625 @@
      .                tsurf,tsoil,inertiesoil,albedo,
      .                emis,q2,qsurf,qsoil,tauscaling,totcloudfrac,
-     .                wstar,watercap,perenial_co2ice)
+     .                wstar,watercap,perennial_co2ice)
           ENDIF ! of IF (write_restart)
 
@@ -3205 +3205 @@
      &         "Perennial water ice thickness"
      &         ,"kg.m-2",watercap(:,islope))
+         ENDDO
+         call write_output("perennial_co2ice",
+     &         "Perennial co2 ice thickness","kg.m-2",
+     &         perennial_co2ice(:,iflat))
+         do islope=1,nslope
+          write(str2(1:2),'(i2.2)') islope
+         call write_output("perennial_co2ice_slope"//str2,
+     &         "Perennial c ice thickness"
+     &         ,"kg.m-2",perennial_co2ice(:,islope))
          ENDDO
          call write_output("temp_layer1","temperature in layer 1",

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

@@ -53 +53 @@
   REAL,ALLOCATABLE,SAVE :: qsurf(:,:,:) ! tracer on surface (e.g. kg.m-2)
   REAL,SAVE,ALLOCATABLE :: watercap(:,:) ! Surface water ice (kg.m-2)
-  REAL,SAVE,ALLOCATABLE :: perenial_co2ice(:,:) ! Perenial CO2 ice (kg.m-2)
-!$OMP THREADPRIVATE(tsurf,emis,capcal,fluxgrd,qsurf,watercap,perenial_co2ice)
+  REAL,SAVE,ALLOCATABLE :: perennial_co2ice(:,:) ! Perennial CO2 ice (kg.m-2)
+!$OMP THREADPRIVATE(tsurf,emis,capcal,fluxgrd,qsurf,watercap,perennial_co2ice)
 
 contains
@@ -77 +77 @@
     allocate(tsurf(ngrid,nslope))
     allocate(watercap(ngrid,nslope))
-    allocate(perenial_co2ice(ngrid,nslope))
+    allocate(perennial_co2ice(ngrid,nslope))
     allocate(emis(ngrid,nslope))
     allocate(capcal(ngrid,nslope))
@@ -108 +108 @@
     if (allocated(tsurf))            deallocate(tsurf)
     if (allocated(watercap))         deallocate(watercap)
-    if (allocated(perenial_co2ice))  deallocate(perenial_co2ice)
+    if (allocated(perennial_co2ice))  deallocate(perennial_co2ice)
     if (allocated(emis))             deallocate(emis)
     if (allocated(capcal))           deallocate(capcal)
@@ -130 +130 @@
     allocate(tsurf(ngrid,nslope))
     allocate(watercap(ngrid,nslope))
-    allocate(perenial_co2ice(ngrid,nslope))
+    allocate(perennial_co2ice(ngrid,nslope))
     allocate(emis(ngrid,nslope))
     allocate(capcal(ngrid,nslope))
@@ -144 +144 @@
     if (allocated(tsurf))            deallocate(tsurf)
     if (allocated(watercap))         deallocate(watercap)
-    if (allocated(perenial_co2ice))  deallocate(perenial_co2ice)
+    if (allocated(perennial_co2ice)) deallocate(perennial_co2ice)
     if (allocated(emis))             deallocate(emis)
     if (allocated(capcal))           deallocate(capcal)