Changeset 2999 for trunk


Ignore:
Timestamp:
Jul 19, 2023, 11:40:38 AM (17 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
Files:
12 edited

Legend:

Unmodified
Added
Removed

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

    r2959 r2999  
    33     &     t,ucov,vcov,ps,h,phisold_newgrid,
    44     &     q,qsurf,tauscaling,totcloudfrac,surfith,nid,
    5      &     watercap)
     5     &     watercap,peren_co2ice)
    66c=======================================================================
    77c
    … …  
    5757      REAL,INTENT(OUT) :: totcloudfrac(ngrid) ! sub grid cloud fraction
    5858      REAL,INTENT(OUT) :: watercap(ngrid,nslope) ! infinite polar cap
     59      REAL,INTENT(OUT) :: peren_co2ice(ngrid,nslope) ! infinite co2 polar cap
    5960      REAL,INTENT(OUT) :: phisold_newgrid(iip1,jjp1)
    6061      REAL,INTENT(OUT) :: t(iip1,jjp1,llm)
    … …  
    140141      real totcloudfracS(iip1,jjp1)
    141142      real watercapS(iip1,jjp1,nslope)
     143      real peren_co2iceS(iip1,jjp1,nslope)
    142144      real watercaptagS(iip1,jjp1)
    143145      real albedoS(iip1,jjp1,nslope)
    … …  
    178180      real, dimension(:,:), allocatable :: totcloudfracold
    179181      real, dimension(:,:,:), allocatable :: watercapold
     182      real, dimension(:,:,:), allocatable :: peren_co2iceold
    180183      real, dimension(:,:), allocatable :: watercapold_noslope
     184      real, dimension(:,:), allocatable :: peren_co2iceold_noslope
    181185      real, dimension(:,:), allocatable :: watercaptagold
    182186      real, dimension(:), allocatable :: watercaptag_tmp
    … …  
    430434      allocate(totcloudfracold(imold+1,jmold+1))
    431435      allocate(watercapold(imold+1,jmold+1,nslope))
     436      allocate(peren_co2iceold(imold+1,jmold+1,nslope))
    432437      allocate(watercaptagold(imold+1,jmold+1))
    433438      allocate(watercaptag_tmp(ngrid))
    … …  
    439444      allocate(emisold_noslope(imold+1,jmold+1))
    440445      allocate(watercapold_noslope(imold+1,jmold+1))
     446      allocate(peren_co2iceold_noslope(imold+1,jmold+1))
    441447      allocate(albedoold_noslope(imold+1,jmold+1))
    442448      allocate(qsurfold_noslope(imold+1,jmold+1,nqtot))
    … …  
    914920        IF (ierr .NE. NF_NOERR) THEN
    915921           PRINT*, "lect_start_archive: Failed loading <watercap>"
     922           PRINT*, NF_STRERROR(ierr)
     923           CALL abort
     924        ENDIF
     925      ENDIF
     926
     927      ierr = NF_INQ_VARID (nid, "perenial_co2ice", nvarid)
     928      IF (ierr .NE. NF_NOERR) THEN
     929         PRINT*, "lect_start_archive: <perenial_co2ice> not in file"
     930         peren_co2iceold(:,:,:) = 0.
     931      ELSE
     932      if(no_slope) then
     933       ierr = nf90_get_var(nid, nvarid,peren_co2iceold_noslope)
     934       peren_co2iceold(:,:,1)=peren_co2iceold_noslope(:,:)
     935      else
     936       ierr = nf90_get_var(nid, nvarid,peren_co2iceold)
     937      endif
     938        IF (ierr .NE. NF_NOERR) THEN
     939           PRINT*, "lect_start_archive:"
     940           PRINT*, "Failed loading<peren_co2iceold>"
    916941           PRINT*, NF_STRERROR(ierr)
    917942           CALL abort
    … …  
    12521277     &       watercap(:,1))
    12531278
     1279c Watercap
     1280      call interp_horiz (peren_co2iceold(:,:,1),peren_co2iceS(:,:,1),
     1281     &       imold,jmold,iim,jjm,1,rlonuold,rlatvold,rlonu,rlatv)
     1282      call gr_dyn_fi (1,iim+1,jjm+1,ngrid,peren_co2iceS(:,:,1),
     1283     &       peren_co2ice(:,1))
     1284
    12541285c Watercaptag
    12551286      if(imold.eq.iim .and. jmold.eq.jjm) then
    … …  
    16481679      deallocate(tauscalingold)
    16491680      deallocate(totcloudfracold)
     1681      deallocate(peren_co2iceold)
    16501682      deallocate(watercapold)
    16511683      deallocate(watercaptagold)

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

    r2952 r2999  
    3030     &                     albedodat, z0_default, qsurf, tsurf,
    3131     &                     emis, hmons, summit, base, watercap,
    32      &               ini_surfdat_h_slope_var,end_surfdat_h_slope_var
     32     &               ini_surfdat_h_slope_var,end_surfdat_h_slope_var,
     33     &               perenial_co2ice
    3334      use comsoil_h, only: inertiedat, inertiesoil,layer, mlayer,
    3435     & nsoilmx,tsoil,ini_comsoil_h_slope_var, end_comsoil_h_slope_var,
    … …  
    196197      REAL,ALLOCATABLE :: qsurf_old_slope(:,:,:) ! tracer on surface (e.g. kg.m-2)
    197198      REAL,ALLOCATABLE :: watercap_old_slope(:,:) ! Surface water ice (kg.m-2)
     199      REAL,ALLOCATABLE :: perenial_co2_old_slope(:,:) ! Surface water ice (kg.m-2)
    198200      REAL,ALLOCATABLE :: tsoil_old_slope(:,:,:)
    199201      REAL,ALLOCATABLE :: inertiesoil_old_slope(:,:,:)
    … …  
    464466     &        day_ini,time,tsurf,tsoil,albedo,emis,
    465467     &        q2,qsurf,tauscaling,totcloudfrac,
    466      &        wstar,watercap,def_slope,def_slope_mean,subslope_dist)
     468     &        wstar,watercap,perenial_co2ice,
     469     &        def_slope,def_slope_mean,subslope_dist)
    467470       
    468471        ! copy albedo and soil thermal inertia
    … …  
    18001803        allocate(emis_old_slope(ngridmx,nslope_old))   
    18011804        allocate(watercap_old_slope(ngridmx,nslope_old))
     1805        allocate(perenial_co2_old_slope(ngridmx,nslope_old))
    18021806
    18031807
    … …  
    18061810        emis_old_slope(:,:)=emis(:,:)
    18071811        watercap_old_slope(:,:)=watercap(:,:)
    1808 
     1812        perenial_co2_old_slope(:,:) = perenial_co2ice(:,:)
    18091813        call end_surfdat_h_slope_var
    18101814        call ini_surfdat_h_slope_var(ngridmx,nqtot,nslope_new)
    … …  
    18371841             emis(:,islope)=emis_old_slope(:,1)
    18381842             watercap(:,islope)=watercap_old_slope(:,1)
     1843             perenial_co2ice(:,islope)= perenial_co2_old_slope(:,1)
    18391844             tsoil(:,:,islope)=tsoil_old_slope(:,:,1)
    18401845             albedo(:,:,islope)=albedo_old_slope(:,:,1)
    … …  
    18471852             emis(:,1)=emis_old_slope(:,iflat)
    18481853             watercap(:,1)=watercap_old_slope(:,iflat)
     1854             perenial_co2ice(:,islope)= perenial_co2_old_slope(:,iflat)
    18491855             tsoil(:,:,1)=tsoil_old_slope(:,:,iflat)
    18501856             albedo(:,:,1)=albedo_old_slope(:,:,iflat)
    … …  
    18571863             emis(:,islope)=emis_old_slope(:,iflat)
    18581864             watercap(:,islope)=watercap_old_slope(:,iflat)
     1865             perenial_co2ice(:,islope)= perenial_co2_old_slope(:,iflat)
    18591866             tsoil(:,:,islope)=tsoil_old_slope(:,:,iflat)
    18601867             albedo(:,:,islope)=albedo_old_slope(:,:,iflat)
    … …  
    18681875             emis(:,islope)=emis_old_slope(:,iflat)
    18691876             watercap(:,islope)=watercap_old_slope(:,iflat)
     1877             perenial_co2ice(:,islope)= perenial_co2_old_slope(:,iflat)
    18701878             tsoil(:,:,islope)=tsoil_old_slope(:,:,iflat)
    18711879             albedo(:,:,islope)=albedo_old_slope(:,:,iflat)
    … …  
    20012009     &              dtphys,hour_ini,
    20022010     &              tsurf,tsoil,inertiesoil,albedo,emis,q2,qsurf,
    2003      &              tauscaling,totcloudfrac,wstar,watercap)
     2011     &              tauscaling,totcloudfrac,wstar,watercap,
     2012     &              perenial_co2ice)
    20042013
    20052014c=======================================================================

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

    r2959 r2999  
    6767      REAL,ALLOCATABLE :: tsoil(:,:,:) ! Soil temperature
    6868      REAL,ALLOCATABLE :: watercap(:,:)        ! h2o ice layer
     69      REAL,ALLOCATABLE :: perenial_co2ice(:,:) ! co2 ice layer
    6970      REAL :: tauscaling(ngridmx) ! dust conversion factor
    7071      REAL:: totcloudfrac(ngridmx) ! sub-grid cloud fraction
    … …  
    8687      REAL ithS(ip1jmp1,nsoilmx) ! Soil Thermal Inertia for inertie dat (present day climate)
    8788      REAL,ALLOCATABLE :: watercapS(:,:)
     89      REAL,ALLOCATABLE :: perenial_co2iceS(:,:)
    8890      REAL,ALLOCATABLE :: watercaptag_tmp(:)
    8991      REAL,ALLOCATABLE :: watercaptagS(:)
    … …  
    169171      allocate(tsoil(ngridmx,nsoilmx,nslope))
    170172      allocate(watercap(ngridmx,nslope))
     173      allocate(perenial_co2ice(ngridmx,nslope))
    171174      allocate(emis(ngridmx,nslope))
    172175      allocate(albedo(ngridmx,2,nslope))
    … …  
    177180      allocate(inertiesoilS(ip1jmp1,nsoilmx,nslope))
    178181      allocate(watercapS(ip1jmp1,nslope))
     182      allocate(perenial_co2iceS(ip1jmp1,nslope))
    179183      allocate(watercaptagS(ip1jmp1))
    180184      allocate(emisS(ip1jmp1,nslope))
    … …  
    184188      CALL phyetat0 (fichnom,0,Lmodif,nsoilmx,ngridmx,llm,nqtot,
    185189     &      day_ini_fi,timefi,tsurf,tsoil,albedo,emis,q2,qsurf,
    186      &      tauscaling,totcloudfrac,wstar,watercap,def_slope,
    187      &      def_slope_mean,subslope_dist)
     190     &      tauscaling,totcloudfrac,wstar,watercap,perenial_co2ice,
     191     &      def_slope, def_slope_mean,subslope_dist)
    188192
    189193       ierr = NF_OPEN (fichnom, NF_NOWRITE,nid1)
    … …  
    270274c tsurf --> tsurfS
    271275c watercap --> watercapS
     276c perenial_co2ice --> perenial_co2iceS
    272277c tsoil --> tsoilS
    273278c inertiesoil --> inertiesoilS
    … …  
    287292      call gr_fi_dyn(1,ngridmx,iip1,jjp1,watercap(:,islope),
    288293     &    watercapS(:,islope))
     294      call gr_fi_dyn(1,ngridmx,iip1,jjp1,perenial_co2ice(:,islope),
     295     &    perenial_co2iceS(:,islope))
    289296      call gr_fi_dyn(nsoilmx,ngridmx,iip1,jjp1,tsoil(:,:,islope),
    290297     &    tsoilS(:,:,islope))
    … …  
    417424      call write_archive(nid,ntime,'watercap','couche de glace h2o',
    418425     &  'kg/m2',2,watercapS)
     426      call write_archive(nid,ntime,'perenial_co2ice',
     427     &'couche de glace co2','kg/m2',2,perenial_co2iceS)
    419428      call write_archive(nid,ntime,'watercaptag','couche de glace h2o',
    420429     &  'kg/m2',2,watercaptagS)

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

    r2959 r2999  
    6969      REAL,ALLOCATABLE :: tsoil(:,:,:) ! Soil temperature
    7070      REAL,ALLOCATABLE :: watercap(:,:)        ! h2o ice layer
     71      REAL,ALLOCATABLE :: perenial_co2ice(:,:)        ! co2 ice layer
    7172      REAL tauscaling(ngridmx) ! dust conversion factor
    7273      REAL totcloudfrac(ngridmx) ! sub-grid cloud fraction
    … …  
    9192      REAL ithS(ip1jmp1,nsoilmx) ! Soil Thermal Inertia for inertie dat (present day climate)
    9293      REAL,ALLOCATABLE :: watercapS(:,:)
     94      REAL,ALLOCATABLE :: perenial_co2iceS(:,:)        ! co2 ice layer
    9395      REAL tauscalingS(ip1jmp1)
    9496      REAL totcloudfracS(ip1jmp1)
    … …  
    184186      allocate(tsoil(ngridmx,nsoilmx,nslope))
    185187      allocate(watercap(ngridmx,nslope))
     188      allocate(perenial_co2ice(ngridmx,nslope))
    186189      allocate(emis(ngridmx,nslope))
    187190      allocate(albedo(ngridmx,2,nslope))
    … …  
    192195      allocate(inertiesoilS(ip1jmp1,nsoilmx,nslope))
    193196      allocate(watercapS(ip1jmp1,nslope))
     197      allocate(perenial_co2iceS(ngridmx,nslope))
    194198      allocate(emisS(ip1jmp1,nslope))
    195199      allocate(albedoS(ip1jmp1,nslope))
    … …  
    198202      CALL phyetat0 (fichnom,0,Lmodif,nsoilmx,ngridmx,llm,nqtot,
    199203     &      day_ini_fi,timefi,tsurf,tsoil,albedo,emis,q2,qsurf,
    200      &      tauscaling,totcloudfrac,wstar,watercap,def_slope,
    201      &      def_slope_mean,subslope_dist)
     204     &      tauscaling,totcloudfrac,wstar,watercap,perenial_co2ice,
     205     &      def_slope,def_slope_mean,subslope_dist)
    202206
    203207       ierr = NF_OPEN (fichnom, NF_NOWRITE,nid1)
    … …  
    296300      call gr_fi_dyn(1,ngridmx,iip1,jjp1,watercap(:,islope),
    297301     &    watercapS(:,islope))
     302      call gr_fi_dyn(1,ngridmx,iip1,jjp1,perenial_co2ice(:,islope),
     303     &    perenial_co2iceS(:,islope))
    298304      call gr_fi_dyn(nsoilmx,ngridmx,iip1,jjp1,tsoil(:,:,islope),
    299305     &    tsoilS(:,:,islope))
    … …  
    427433      call write_archive(nid,ntime,'watercap','couche de glace h2o',
    428434     &  'kg/m2',2,watercapS)
     435      call write_archive(nid,ntime,'perenial_co2ice','couche de glace co2',
     436     &  'kg/m2',2,perenial_co2iceS)
    429437      call write_archive(nid,ntime,'tauscaling',
    430438     &  'dust conversion factor',' ',2,tauscalingS)

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

    r2609 r2999  
    1 subroutine albedocaps(zls,ngrid,piceco2,psolaralb,emisref)
     1subroutine albedocaps(zls,ngrid,piceco2,piceco2_peren,psolaralb,emisref)
    22
    33! routine which changes the albedo (and emissivity) of the surface
    … …  
    1111USE mod_phys_lmdz_transfert_para, ONLY: bcast
    1212USE mod_phys_lmdz_para, ONLY: is_master
     13USE paleoclimate_mod, ONLY: paleoclimate,albedo_perenialco2
     14
    1315implicit none
    1416
    … …  
    1921integer,intent(in) :: ngrid
    2022real,intent(in) :: piceco2(ngrid) ! amount of CO2 ice on the surface (kg/m2)
     23real,intent(inout) :: piceco2_peren(ngrid) ! amount of perenial co2 ice (kg/m^2)
    2124real,intent(out) :: psolaralb(ngrid,2) ! albedo of the surface
    2225real,intent(out) :: emisref(ngrid) ! emissivity of the surface
    … …  
    104107      psolaralb(ig,1)=albedice(icap)
    105108      psolaralb(ig,2)=albedice(icap)
     109      if(paleoclimate) then
     110         if((piceco2_peren(ig).gt.0.).and.(piceco2(ig).lt.piceco2_peren(ig))) then
     111             psolaralb(ig,1) = albedo_perenialco2
     112             psolaralb(ig,2) = albedo_perenialco2
     113             piceco2_peren(ig) = piceco2(ig)
     114         endif
     115      endif
    106116    endif
    107117  else if (watercaptag(ig) .and. water) then

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

    r2977 r2999  
    1111     $                  pcapcal,pplay,pplev,ptsrf,pt,
    1212     $                  pphi,pdt,pdu,pdv,pdtsrf,pu,pv,pq,pdq,
    13      $                  piceco2,psolaralb,pemisurf,rdust,
     13     $                  piceco2,perenial_co2ice,
     14     $                  psolaralb,pemisurf,rdust,
    1415     $                  pdtc,pdtsrfc,pdpsrf,pduc,pdvc,pdqc,
    1516     $                  fluxsurf_sw,zls,
    … …  
    3536#endif
    3637      use comslope_mod, ONLY: subslope_dist,def_slope_mean
     38      USE paleoclimate_mod, ONLY: paleoclimate
     39
    3740       IMPLICIT NONE
    3841c=======================================================================
    … …  
    8790
    8891      REAL,INTENT(INOUT) :: piceco2(ngrid,nslope) ! CO2 ice on the surface (kg.m-2)
     92      REAL,INTENT(INOUT) :: perenial_co2ice(ngrid,nslope) ! Perenial CO2 ice on the surface (kg.m-2)
    8993      REAL,INTENT(INOUT) :: psolaralb(ngrid,2,nslope) ! albedo of the surface
    9094      REAL,INTENT(INOUT) :: pemisurf(ngrid,nslope) ! emissivity of the surface
    … …  
    180184      REAL   :: alb_tmp(ngrid,2) ! local
    181185      REAL   :: zcondices_tmp(ngrid)    ! local 
    182       REAL   :: piceco2_tmp(ngrid)    ! local 
     186      REAL   :: piceco2_tmp(ngrid)    ! local
     187      REAL   :: perenial_co2ice_tmp(ngrid) ! perenial ice on one subslope (kg/m^2)
    183188      REAL   :: pemisurf_tmp(ngrid)! local
    184189      LOGICAL :: condsub_tmp(ngrid) !local
    … …  
    568573        piceco2_tmp(:) = piceco2(:,islope)
    569574        alb_tmp(:,:) = psolaralb(:,:,islope)
    570         emisref_tmp(:) = 0.
    571         CALL albedocaps(zls,ngrid,piceco2_tmp,alb_tmp,emisref_tmp)
     575        emisref_tmp(:) = 0.
     576        perenial_co2ice_tmp(:) =  perenial_co2ice(:,islope)
     577        CALL albedocaps(zls,ngrid,piceco2_tmp,perenial_co2ice_tmp,
     578     &                  alb_tmp,emisref_tmp)
     579        perenial_co2ice(:,islope) = perenial_co2ice_tmp(:)
    572580        psolaralb(:,1,islope) =  alb_tmp(:,1)
    573581        psolaralb(:,2,islope) =  alb_tmp(:,2)
    … …  
    823831! Extra special case for surface temperature tendency pdtsrfc:
    824832! we want to fix the south pole temperature to CO2 condensation temperature
    825          if (caps.and.(obliquit.lt.27.)) then
     833         if (caps.and.(obliquit.lt.27.).and.(.not.(paleoclimate))) then
    826834           ! check if last grid point is the south pole
    827835           if (abs(latitude(ngrid)-(-pi/2.)).lt.1.e-5) then

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

    r2994 r2999  
    5454      use aeropacity_mod, only: iddist, topdustref
    5555      USE mod_phys_lmdz_transfert_para, ONLY: bcast
    56       USE paleoclimate_mod,ONLY: paleoclimate
     56      USE paleoclimate_mod,ONLY: paleoclimate,albedo_perenialco2
    5757      IMPLICIT NONE
    5858      include "callkeys.h"
    … …  
    355355         write(*,*)" paleoclimate = ",paleoclimate
    356356
     357         write(*,*)"Albedo for perenial CO2 ice?"
     358         albedo_perenialco2 = 0.85 ! default value
     359         call getin_p("albedo_perenialco2",albedo_perenialco2)
     360         write(*,*)"albedo_perenialco2 = ",albedo_perenialco2
    357361
    358362! TRACERS:

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

    r2994 r2999  
    1111    IMPLICIT NONE
    1212
    13     LOGICAL :: paleoclimate                        ! False by default, is activate for paleoclimates specific processes (e.g., lag layer)
     13    LOGICAL :: paleoclimate                        ! False by default, is activate  for paleoclimates specific processes (e.g., lag layer)
    1414    real, save, allocatable :: lag_h2o_ice(:,:)  ! Thickness of the lag before H2O ice [m]
    1515    real, save, allocatable :: lag_co2_ice(:,:)  ! Thickness of the lag before CO2 ice [m]
    16 
     16    real, save :: albedo_perenialco2             ! Albedo for perenial co2 ice [1]
    1717    CONTAINS
    1818

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

    r2994 r2999  
    1010subroutine phyetat0 (fichnom,tab0,Lmodif,nsoil,ngrid,nlay,nq, &
    1111                     day_ini,time0,tsurf,tsoil,albedo,emis,q2,qsurf, &
    12                      tauscaling,totcloudfrac,wstar,watercap,def_slope, &
    13                      def_slope_mean,subslope_dist)
     12                     tauscaling,totcloudfrac,wstar,watercap,perenial_co2ice, &
     13                     def_slope,def_slope_mean,subslope_dist)
    1414
    1515  use tracer_mod, only: noms ! tracer names
    … …  
    2828  USE comslope_mod, ONLY: nslope, major_slope
    2929  USE paleoclimate_mod, ONLY: paleoclimate, lag_h2o_ice, lag_co2_ice
     30  USE comcstfi_h, only: pi
     31  use geometry_mod, only: latitude
    3032  implicit none
    3133 
    … …  
    7173  real,intent(out) :: wstar(ngrid) ! Max vertical velocity in thermals (m/s)
    7274  real,intent(out) :: watercap(ngrid,nslope) ! h2o_ice_cover
    73   real, intent(out) :: def_slope(nslope+1) !boundaries for bining of the slopes
    74   real, intent(out) :: def_slope_mean(nslope)
    75   real, intent(out) :: subslope_dist(ngrid,nslope) !undermesh statistics
     75  real,intent(out) :: perenial_co2ice(ngrid,nslope) ! perenial co2 ice(kg/m^2)
     76  real,intent(out) :: def_slope(nslope+1) !boundaries for bining of the slopes
     77  real,intent(out) :: def_slope_mean(nslope)
     78  real,intent(out) :: subslope_dist(ngrid,nslope) !undermesh statistics
    7679!======================================================================
    7780!  Local variables:
    … …  
    117120      REAL :: sum_dist
    118121      REAL :: current_max !var to find max distrib slope
     122! Variables for CO2 index
     123      INTEGER :: igcm_co2_tmp
    119124
    120125      CHARACTER(len=5) :: modname="phyetat0"
    … …  
    755760
    756761if (paleoclimate) then
     762  do iq=1,nq
     763   txt=noms(iq)
     764   if (txt.eq."co2") igcm_co2_tmp = iq
     765  enddo
    757766  if (startphy_file) then
     767! Depth of H2O lag
    758768   call get_field("lag_h2o_ice",lag_h2o_ice,found,indextime)
    759    write(*,*) 'paleo found start?',found
    760    
    761769   if (.not.found) then
    762770     write(*,*) "phyetat0: Failed loading <lag_h2o_ice> : ", &
    … …  
    765773   endif
    766774
     775!  Depth of CO2 lag
    767776   call get_field("lag_co2_ice",lag_co2_ice,found,indextime)
    768777   if (.not.found) then
    … …  
    771780     lag_co2_ice(:,:) = -1.
    772781   endif
    773   else
     782
     783! Perenial CO2 ice   
     784   call get_field("perenial_co2ice",perenial_co2ice,found,indextime)
     785   if (.not.found) then
     786     write(*,*) "phyetat0: Failed loading <perenial_co2ice> : ", &
     787                "<perenial_co2ice> is set as 10m at the South Pole"
     788     perenial_co2ice(:,:) = 0.
     789     if (abs(latitude(ngrid)-(-pi/2.)).lt.1.e-5) then
     790       DO islope = 1,nslope
     791        perenial_co2ice(ngrid,islope) = 10*1.6e3 ! 10m which is convert to kg/m^2
     792        qsurf(ngrid,igcm_co2_tmp,islope) = qsurf(ngrid-1,igcm_co2_tmp,islope) + & ! perenial ice + frost
     793                                        perenial_co2ice(ngrid,islope)
     794       ENDDO
     795     endif
     796   endif ! notfound
     797  else ! no startfiphyle
    774798     lag_h2o_ice(:,:) = -1.
    775799     lag_co2_ice(:,:) = -1.
     800     perenial_co2ice(:,:) = 0.
     801     if (abs(latitude(ngrid)-(-pi/2.)).lt.1.e-5) then
     802       DO islope = 1,nslope
     803        perenial_co2ice(ngrid,islope) = 10*1.6e3 ! 10m which is convert to kg/m^2
     804        qsurf(ngrid,igcm_co2_tmp,islope) = qsurf(ngrid-1,igcm_co2_tmp,islope) + & ! perenial ice + frost
     805                                        perenial_co2ice(ngrid,islope)
     806       ENDDO
     807     endif
    776808  endif !startphy_file
    777809else
    778    write(*,*) 'paleo found? nostart',found
    779    
    780810  lag_h2o_ice(:,:) = -1.
    781811  lag_co2_ice(:,:) = -1.
     812  perenial_co2ice(:,:) = 0.
    782813endif !paleoclimate
    783814

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

    r2994 r2999  
    168168                    albedo,emis,q2,qsurf,&
    169169                    tauscaling,totcloudfrac,wstar, &
    170                     watercap)
     170                    watercap,perenial_co2ice)
    171171  ! write time-dependent variable to restart file
    172172  use iostart, only : open_restartphy, close_restartphy, &
    … …  
    178178  use dust_param_mod, only: dustscaling_mode
    179179  use comsoil_h,only: flux_geo
    180   use comslope_mod, ONLY: nslope
     180  use comslope_mod, only: nslope
     181  use paleoclimate_mod, only: paleoclimate
    181182  implicit none
    182183 
    … …  
    201202  real,intent(in) :: wstar(ngrid)
    202203  real,intent(in) :: watercap(ngrid,nslope)
    203  
     204  real,intent(in) :: perenial_co2ice(ngrid,nslope)
     205
    204206  integer :: iq
    205207  character(len=30) :: txt ! to store some text
    … …  
    220222  ! Water ice layer
    221223  call put_field("watercap","H2O ice cover",watercap,time)
    222  
     224 
     225 ! Perenial CO2 ice layer
     226  if(paleoclimate) then
     227    call put_field("perenial_co2ice","CO2 ice cover",perenial_co2ice,time)
     228  endif
    223229  ! Surface temperature
    224230  call put_field("tsurf","Surface temperature",tsurf,time)

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

    r2965 r2999  
    5252     &                     tsurf, emis,
    5353     &                     capcal, fluxgrd, qsurf,
    54      &                     hmons,summit,base,watercap,watercaptag
     54     &                     hmons,summit,base,watercap,watercaptag,
     55     &                     perenial_co2ice
    5556      use comsaison_h, only: dist_sol, declin, zls,
    5657     &                       mu0, fract, local_time
    … …  
    595596     &         tsurf,tsoil,albedo,emis,
    596597     &         q2,qsurf,tauscaling,totcloudfrac,wstar,
    597      &         watercap,def_slope,def_slope_mean,subslope_dist)
     598     &         watercap,perenial_co2ice,
     599     &         def_slope,def_slope_mean,subslope_dist)
    598600
    599601!     Sky view:
    … …  
    21992201     $              capcal,zplay,zplev,tsurf,pt,
    22002202     $              pphi,pdt,pdu,pdv,zdtsurf,pu,pv,pq,pdq,
    2201      $              qsurf(:,igcm_co2,:),albedo,emis,rdust,
     2203     $              qsurf(:,igcm_co2,:),perenial_co2ice,
     2204     $              albedo,emis,rdust,
    22022205     $              zdtc,zdtsurfc,pdpsrf,zduc,zdvc,zdqc,
    22032206     $              fluxsurf_dn_sw,zls,
    … …  
    26142617     .                tsurf,tsoil,inertiesoil,albedo,
    26152618     .                emis,q2,qsurf,tauscaling,totcloudfrac,wstar,
    2616      .                watercap)
     2619     .                watercap,perenial_co2ice)
    26172620          ENDIF ! of IF (write_restart)
    26182621

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

    r2909 r2999  
    5050  REAL,ALLOCATABLE,SAVE :: qsurf(:,:,:) ! tracer on surface (e.g. kg.m-2)
    5151  REAL,SAVE,ALLOCATABLE :: watercap(:,:) ! Surface water ice (kg.m-2)
    52 
    53 !$OMP THREADPRIVATE(tsurf,emis,capcal,fluxgrd,qsurf,watercap)
     52  REAL,SAVE,ALLOCATABLE :: perenial_co2ice(:,:) ! Perenial CO2 ice (kg.m-2)
     53!$OMP THREADPRIVATE(tsurf,emis,capcal,fluxgrd,qsurf,watercap,perenial_co2ice)
    5454
    5555contains
    … …  
    7474    allocate(tsurf(ngrid,nslope))
    7575    allocate(watercap(ngrid,nslope))
     76    allocate(perenial_co2ice(ngrid,nslope))
    7677    allocate(emis(ngrid,nslope))
    7778    allocate(capcal(ngrid,nslope))
    … …  
    9192  implicit none
    9293
    93     if (allocated(albedodat))     deallocate(albedodat)
    94     if (allocated(phisfi))        deallocate(phisfi)
    95     if (allocated(watercaptag))   deallocate(watercaptag)
    96     if (allocated(dryness))       deallocate(dryness)
    97     if (allocated(zmea))          deallocate(zmea)
    98     if (allocated(zstd))          deallocate(zstd)
    99     if (allocated(zsig))          deallocate(zsig)
    100     if (allocated(zgam))          deallocate(zgam)
    101     if (allocated(zthe))          deallocate(zthe)
    102     if (allocated(z0))            deallocate(z0)
    103     if (allocated(qsurf))         deallocate(qsurf)
    104     if (allocated(tsurf))         deallocate(tsurf)
    105     if (allocated(watercap))      deallocate(watercap)
    106     if (allocated(emis))          deallocate(emis)
    107     if (allocated(capcal))        deallocate(capcal)
    108     if (allocated(fluxgrd))       deallocate(fluxgrd)
    109     if (allocated(hmons))         deallocate(hmons)
    110     if (allocated(summit))        deallocate(summit)
    111     if (allocated(base))          deallocate(base)
    112     if (allocated(alpha_hmons))   deallocate(alpha_hmons)
    113     if (allocated(hsummit))       deallocate(hsummit)
    114     if (allocated(contains_mons)) deallocate(contains_mons)
     94    if (allocated(albedodat))        deallocate(albedodat)
     95    if (allocated(phisfi))           deallocate(phisfi)
     96    if (allocated(watercaptag))      deallocate(watercaptag)
     97    if (allocated(dryness))          deallocate(dryness)
     98    if (allocated(zmea))             deallocate(zmea)
     99    if (allocated(zstd))             deallocate(zstd)
     100    if (allocated(zsig))             deallocate(zsig)
     101    if (allocated(zgam))             deallocate(zgam)
     102    if (allocated(zthe))             deallocate(zthe)
     103    if (allocated(z0))               deallocate(z0)
     104    if (allocated(qsurf))            deallocate(qsurf)
     105    if (allocated(tsurf))            deallocate(tsurf)
     106    if (allocated(watercap))         deallocate(watercap)
     107    if (allocated(perenial_co2ice))  deallocate(perenial_co2ice)
     108    if (allocated(emis))             deallocate(emis)
     109    if (allocated(capcal))           deallocate(capcal)
     110    if (allocated(fluxgrd))          deallocate(fluxgrd)
     111    if (allocated(hmons))            deallocate(hmons)
     112    if (allocated(summit))           deallocate(summit)
     113    if (allocated(base))             deallocate(base)
     114    if (allocated(alpha_hmons))      deallocate(alpha_hmons)
     115    if (allocated(hsummit))          deallocate(hsummit)
     116    if (allocated(contains_mons))    deallocate(contains_mons)
    115117   
    116118  end subroutine end_surfdat_h
    … …  
    125127    allocate(tsurf(ngrid,nslope))
    126128    allocate(watercap(ngrid,nslope))
     129    allocate(perenial_co2ice(ngrid,nslope))
    127130    allocate(emis(ngrid,nslope))
    128131    allocate(capcal(ngrid,nslope))
    … …  
    135138  implicit none
    136139
    137     if (allocated(qsurf))         deallocate(qsurf)
    138     if (allocated(tsurf))         deallocate(tsurf)
    139     if (allocated(watercap))      deallocate(watercap)
    140     if (allocated(emis))          deallocate(emis)
    141     if (allocated(capcal))        deallocate(capcal)
    142     if (allocated(fluxgrd))       deallocate(fluxgrd)
     140    if (allocated(qsurf))            deallocate(qsurf)
     141    if (allocated(tsurf))            deallocate(tsurf)
     142    if (allocated(watercap))         deallocate(watercap)
     143    if (allocated(perenial_co2ice))  deallocate(perenial_co2ice)
     144    if (allocated(emis))             deallocate(emis)
     145    if (allocated(capcal))           deallocate(capcal)
     146    if (allocated(fluxgrd))          deallocate(fluxgrd)
    143147   
    144148  end subroutine end_surfdat_h_slope_var
Note: See TracChangeset for help on using the changeset viewer.