Changeset 2999 for trunk/LMDZ.MARS/libf/dynphy_lonlat/phymars

Timestamp:

Jul 19, 2023, 11:40:38 AM (18 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

Location:

trunk/LMDZ.MARS/libf/dynphy_lonlat/phymars

Files:

• lect_start_archive.F (modified) (9 diffs)
• newstart.F (modified) (10 diffs)
• start2archive.F (modified) (8 diffs)
• start2archive_SSO.F (modified) (7 diffs)

trunk/LMDZ.MARS/libf/dynphy_lonlat/phymars/lect_start_archive.F

@@ -3 +3 @@
      &     t,ucov,vcov,ps,h,phisold_newgrid,
      &     q,qsurf,tauscaling,totcloudfrac,surfith,nid,
-     &     watercap)
+     &     watercap,peren_co2ice)
 c=======================================================================
 c
@@ -57 +57 @@
       REAL,INTENT(OUT) :: totcloudfrac(ngrid) ! sub grid cloud fraction
       REAL,INTENT(OUT) :: watercap(ngrid,nslope) ! infinite polar cap
+      REAL,INTENT(OUT) :: peren_co2ice(ngrid,nslope) ! infinite co2 polar cap
       REAL,INTENT(OUT) :: phisold_newgrid(iip1,jjp1)
       REAL,INTENT(OUT) :: t(iip1,jjp1,llm)
@@ -140 +141 @@
       real totcloudfracS(iip1,jjp1)
       real watercapS(iip1,jjp1,nslope)
+      real peren_co2iceS(iip1,jjp1,nslope)
       real watercaptagS(iip1,jjp1)
       real albedoS(iip1,jjp1,nslope)
@@ -178 +180 @@
       real, dimension(:,:), allocatable :: totcloudfracold
       real, dimension(:,:,:), allocatable :: watercapold
+      real, dimension(:,:,:), allocatable :: peren_co2iceold
       real, dimension(:,:), allocatable :: watercapold_noslope
+      real, dimension(:,:), allocatable :: peren_co2iceold_noslope
       real, dimension(:,:), allocatable :: watercaptagold
       real, dimension(:), allocatable :: watercaptag_tmp
@@ -430 +434 @@
       allocate(totcloudfracold(imold+1,jmold+1))
       allocate(watercapold(imold+1,jmold+1,nslope))
+      allocate(peren_co2iceold(imold+1,jmold+1,nslope))
       allocate(watercaptagold(imold+1,jmold+1))
       allocate(watercaptag_tmp(ngrid))
@@ -439 +444 @@
       allocate(emisold_noslope(imold+1,jmold+1))
       allocate(watercapold_noslope(imold+1,jmold+1))
+      allocate(peren_co2iceold_noslope(imold+1,jmold+1))
       allocate(albedoold_noslope(imold+1,jmold+1))
       allocate(qsurfold_noslope(imold+1,jmold+1,nqtot))
@@ -914 +920 @@
         IF (ierr .NE. NF_NOERR) THEN
            PRINT*, "lect_start_archive: Failed loading <watercap>"
+           PRINT*, NF_STRERROR(ierr)
+           CALL abort
+        ENDIF
+      ENDIF
+
+      ierr = NF_INQ_VARID (nid, "perenial_co2ice", nvarid)
+      IF (ierr .NE. NF_NOERR) THEN
+         PRINT*, "lect_start_archive: <perenial_co2ice> not in file"
+         peren_co2iceold(:,:,:) = 0.
+      ELSE
+      if(no_slope) then
+       ierr = nf90_get_var(nid, nvarid,peren_co2iceold_noslope)
+       peren_co2iceold(:,:,1)=peren_co2iceold_noslope(:,:)
+      else
+       ierr = nf90_get_var(nid, nvarid,peren_co2iceold)
+      endif
+        IF (ierr .NE. NF_NOERR) THEN
+           PRINT*, "lect_start_archive:"
+           PRINT*, "Failed loading<peren_co2iceold>"
            PRINT*, NF_STRERROR(ierr)
            CALL abort
@@ -1252 +1277 @@
      &       watercap(:,1))
 
+c Watercap
+      call interp_horiz (peren_co2iceold(:,:,1),peren_co2iceS(:,:,1),
+     &       imold,jmold,iim,jjm,1,rlonuold,rlatvold,rlonu,rlatv)
+      call gr_dyn_fi (1,iim+1,jjm+1,ngrid,peren_co2iceS(:,:,1),
+     &       peren_co2ice(:,1))
+
 c Watercaptag
       if(imold.eq.iim .and. jmold.eq.jjm) then
@@ -1648 +1679 @@
       deallocate(tauscalingold)
       deallocate(totcloudfracold)
+      deallocate(peren_co2iceold)
       deallocate(watercapold)
       deallocate(watercaptagold)

trunk/LMDZ.MARS/libf/dynphy_lonlat/phymars/newstart.F

@@ -30 +30 @@
      &                     albedodat, z0_default, qsurf, tsurf,
      &                     emis, hmons, summit, base, watercap,
-     &               ini_surfdat_h_slope_var,end_surfdat_h_slope_var
+     &               ini_surfdat_h_slope_var,end_surfdat_h_slope_var,
+     &               perenial_co2ice
       use comsoil_h, only: inertiedat, inertiesoil,layer, mlayer,
      & nsoilmx,tsoil,ini_comsoil_h_slope_var, end_comsoil_h_slope_var,
@@ -196 +197 @@
       REAL,ALLOCATABLE :: qsurf_old_slope(:,:,:) ! tracer on surface (e.g. kg.m-2)
       REAL,ALLOCATABLE :: watercap_old_slope(:,:) ! Surface water ice (kg.m-2)
+      REAL,ALLOCATABLE :: perenial_co2_old_slope(:,:) ! Surface water ice (kg.m-2)
       REAL,ALLOCATABLE :: tsoil_old_slope(:,:,:) 
       REAL,ALLOCATABLE :: inertiesoil_old_slope(:,:,:) 
@@ -464 +466 @@
      &        day_ini,time,tsurf,tsoil,albedo,emis,
      &        q2,qsurf,tauscaling,totcloudfrac,
-     &        wstar,watercap,def_slope,def_slope_mean,subslope_dist)
+     &        wstar,watercap,perenial_co2ice,
+     &        def_slope,def_slope_mean,subslope_dist)
         
         ! copy albedo and soil thermal inertia
@@ -1800 +1803 @@
         allocate(emis_old_slope(ngridmx,nslope_old))    
         allocate(watercap_old_slope(ngridmx,nslope_old)) 
+        allocate(perenial_co2_old_slope(ngridmx,nslope_old)) 
 
 
@@ -1806 +1810 @@
         emis_old_slope(:,:)=emis(:,:)
         watercap_old_slope(:,:)=watercap(:,:)
-
+        perenial_co2_old_slope(:,:) = perenial_co2ice(:,:)
         call end_surfdat_h_slope_var
         call ini_surfdat_h_slope_var(ngridmx,nqtot,nslope_new)
@@ -1837 +1841 @@
              emis(:,islope)=emis_old_slope(:,1)
              watercap(:,islope)=watercap_old_slope(:,1)
+             perenial_co2ice(:,islope)= perenial_co2_old_slope(:,1)
              tsoil(:,:,islope)=tsoil_old_slope(:,:,1)
              albedo(:,:,islope)=albedo_old_slope(:,:,1)
@@ -1847 +1852 @@
              emis(:,1)=emis_old_slope(:,iflat)
              watercap(:,1)=watercap_old_slope(:,iflat)
+             perenial_co2ice(:,islope)= perenial_co2_old_slope(:,iflat)
              tsoil(:,:,1)=tsoil_old_slope(:,:,iflat)
              albedo(:,:,1)=albedo_old_slope(:,:,iflat)
@@ -1857 +1863 @@
              emis(:,islope)=emis_old_slope(:,iflat)
              watercap(:,islope)=watercap_old_slope(:,iflat)
+             perenial_co2ice(:,islope)= perenial_co2_old_slope(:,iflat)
              tsoil(:,:,islope)=tsoil_old_slope(:,:,iflat)
              albedo(:,:,islope)=albedo_old_slope(:,:,iflat)
@@ -1868 +1875 @@
              emis(:,islope)=emis_old_slope(:,iflat)
              watercap(:,islope)=watercap_old_slope(:,iflat)
+             perenial_co2ice(:,islope)= perenial_co2_old_slope(:,iflat)
              tsoil(:,:,islope)=tsoil_old_slope(:,:,iflat)
              albedo(:,:,islope)=albedo_old_slope(:,:,iflat)
@@ -2001 +2009 @@
      &              dtphys,hour_ini,
      &              tsurf,tsoil,inertiesoil,albedo,emis,q2,qsurf,
-     &              tauscaling,totcloudfrac,wstar,watercap)
+     &              tauscaling,totcloudfrac,wstar,watercap,
+     &              perenial_co2ice)
 
 c=======================================================================

trunk/LMDZ.MARS/libf/dynphy_lonlat/phymars/start2archive.F

@@ -67 +67 @@
       REAL,ALLOCATABLE :: tsoil(:,:,:) ! Soil temperature
       REAL,ALLOCATABLE :: watercap(:,:)        ! h2o ice layer
+      REAL,ALLOCATABLE :: perenial_co2ice(:,:) ! co2 ice layer
       REAL :: tauscaling(ngridmx) ! dust conversion factor
       REAL:: totcloudfrac(ngridmx) ! sub-grid cloud fraction
@@ -86 +87 @@
       REAL ithS(ip1jmp1,nsoilmx) ! Soil Thermal Inertia for inertie dat (present day climate)
       REAL,ALLOCATABLE :: watercapS(:,:)
+      REAL,ALLOCATABLE :: perenial_co2iceS(:,:)
       REAL,ALLOCATABLE :: watercaptag_tmp(:)
       REAL,ALLOCATABLE :: watercaptagS(:)
@@ -169 +171 @@
       allocate(tsoil(ngridmx,nsoilmx,nslope))
       allocate(watercap(ngridmx,nslope))
+      allocate(perenial_co2ice(ngridmx,nslope))
       allocate(emis(ngridmx,nslope))
       allocate(albedo(ngridmx,2,nslope))
@@ -177 +180 @@
       allocate(inertiesoilS(ip1jmp1,nsoilmx,nslope))
       allocate(watercapS(ip1jmp1,nslope))
+      allocate(perenial_co2iceS(ip1jmp1,nslope))
       allocate(watercaptagS(ip1jmp1))
       allocate(emisS(ip1jmp1,nslope))
@@ -184 +188 @@
       CALL phyetat0 (fichnom,0,Lmodif,nsoilmx,ngridmx,llm,nqtot,
      &      day_ini_fi,timefi,tsurf,tsoil,albedo,emis,q2,qsurf,
-     &      tauscaling,totcloudfrac,wstar,watercap,def_slope,
-     &      def_slope_mean,subslope_dist)
+     &      tauscaling,totcloudfrac,wstar,watercap,perenial_co2ice,
+     &      def_slope, def_slope_mean,subslope_dist)
 
        ierr = NF_OPEN (fichnom, NF_NOWRITE,nid1)
@@ -270 +274 @@
 c tsurf --> tsurfS
 c watercap --> watercapS
+c perenial_co2ice --> perenial_co2iceS
 c tsoil --> tsoilS
 c inertiesoil --> inertiesoilS
@@ -287 +292 @@
       call gr_fi_dyn(1,ngridmx,iip1,jjp1,watercap(:,islope),
      &    watercapS(:,islope))
+      call gr_fi_dyn(1,ngridmx,iip1,jjp1,perenial_co2ice(:,islope),
+     &    perenial_co2iceS(:,islope))
       call gr_fi_dyn(nsoilmx,ngridmx,iip1,jjp1,tsoil(:,:,islope),
      &    tsoilS(:,:,islope))
@@ -417 +424 @@
       call write_archive(nid,ntime,'watercap','couche de glace h2o',
      &  'kg/m2',2,watercapS)
+      call write_archive(nid,ntime,'perenial_co2ice',
+     &'couche de glace co2','kg/m2',2,perenial_co2iceS)
       call write_archive(nid,ntime,'watercaptag','couche de glace h2o',
      &  'kg/m2',2,watercaptagS)

trunk/LMDZ.MARS/libf/dynphy_lonlat/phymars/start2archive_SSO.F

@@ -69 +69 @@
       REAL,ALLOCATABLE :: tsoil(:,:,:) ! Soil temperature
       REAL,ALLOCATABLE :: watercap(:,:)        ! h2o ice layer
+      REAL,ALLOCATABLE :: perenial_co2ice(:,:)        ! co2 ice layer
       REAL tauscaling(ngridmx) ! dust conversion factor
       REAL totcloudfrac(ngridmx) ! sub-grid cloud fraction
@@ -91 +92 @@
       REAL ithS(ip1jmp1,nsoilmx) ! Soil Thermal Inertia for inertie dat (present day climate)
       REAL,ALLOCATABLE :: watercapS(:,:)
+      REAL,ALLOCATABLE :: perenial_co2iceS(:,:)        ! co2 ice layer
       REAL tauscalingS(ip1jmp1)
       REAL totcloudfracS(ip1jmp1)
@@ -184 +186 @@
       allocate(tsoil(ngridmx,nsoilmx,nslope))
       allocate(watercap(ngridmx,nslope))
+      allocate(perenial_co2ice(ngridmx,nslope))
       allocate(emis(ngridmx,nslope))
       allocate(albedo(ngridmx,2,nslope))
@@ -192 +195 @@
       allocate(inertiesoilS(ip1jmp1,nsoilmx,nslope))
       allocate(watercapS(ip1jmp1,nslope))
+      allocate(perenial_co2iceS(ngridmx,nslope))
       allocate(emisS(ip1jmp1,nslope))
       allocate(albedoS(ip1jmp1,nslope))
@@ -198 +202 @@
       CALL phyetat0 (fichnom,0,Lmodif,nsoilmx,ngridmx,llm,nqtot,
      &      day_ini_fi,timefi,tsurf,tsoil,albedo,emis,q2,qsurf,
-     &      tauscaling,totcloudfrac,wstar,watercap,def_slope,
-     &      def_slope_mean,subslope_dist)
+     &      tauscaling,totcloudfrac,wstar,watercap,perenial_co2ice,
+     &      def_slope,def_slope_mean,subslope_dist)
 
        ierr = NF_OPEN (fichnom, NF_NOWRITE,nid1)
@@ -296 +300 @@
       call gr_fi_dyn(1,ngridmx,iip1,jjp1,watercap(:,islope),
      &    watercapS(:,islope))
+      call gr_fi_dyn(1,ngridmx,iip1,jjp1,perenial_co2ice(:,islope),
+     &    perenial_co2iceS(:,islope))
       call gr_fi_dyn(nsoilmx,ngridmx,iip1,jjp1,tsoil(:,:,islope),
      &    tsoilS(:,:,islope))
@@ -427 +433 @@
       call write_archive(nid,ntime,'watercap','couche de glace h2o',
      &  'kg/m2',2,watercapS)
+      call write_archive(nid,ntime,'perenial_co2ice','couche de glace co2',
+     &  'kg/m2',2,perenial_co2iceS)
       call write_archive(nid,ntime,'tauscaling',
      &  'dust conversion factor',' ',2,tauscalingS)