Changeset 3423 for trunk


Ignore:
Timestamp:
Sep 2, 2024, 5:20:13 PM (4 months ago)
Author:
bhatnags
Message:

Generic-PCM
Including variable "tice" in startarchive and newstart
SB

Location:
trunk/LMDZ.GENERIC
Files:
4 edited

Legend:

Unmodified
Added
Removed
  • trunk/LMDZ.GENERIC/changelog.txt

    r3397 r3423  
    19701970Fix/clarification of the usage of tice (sea ice surface temperature) vs.
    19711971tsea_ice (temperature of the top layer (ice or snow)).
     1972
     1973== 02/09/2024 == SB
     1974Including variable "tice" in startarchive and newstart
  • trunk/LMDZ.GENERIC/libf/dynphy_lonlat/phystd/lect_start_archive.F

    r3100 r3423  
    33     &     t,ucov,vcov,ps,h,phisold_newgrid,
    44     &     q,qsurf,surfith,nid,
    5      &     rnat,pctsrf_sic,tslab,tsea_ice,sea_ice,
     5     &     rnat,pctsrf_sic,tslab,tsea_ice,tice,sea_ice,
    66     &     du_nonoro_gwd,dv_nonoro_gwd,east_gwstress,west_gwstress)
    77
     
    7171      REAL,INTENT(OUT) :: tslab(ngrid,nslay)
    7272      REAL ,INTENT(OUT) ::rnat(ngrid),pctsrf_sic(ngrid)
    73       REAL,INTENT(OUT) :: tsea_ice(ngrid),sea_ice(ngrid)
     73      REAL,INTENT(OUT) :: tsea_ice(ngrid),tice(ngrid),sea_ice(ngrid)
    7474c     REAL phisfi(ngrid)
    7575      REAL,INTENT(OUT):: du_nonoro_gwd(ngrid,llm)
     
    100100      REAL q2S(iip1,jjp1,llm+1),qsurfS(iip1,jjp1,nqtot)
    101101      real tslabS(iip1,jjp1,nslay)
    102       real pctsrf_sicS(iip1,jjp1),tsea_iceS(iip1,jjp1)
     102      real pctsrf_sicS(iip1,jjp1),tsea_iceS(iip1,jjp1),ticeS(iip1,jjp1)
    103103      real rnatS(iip1,jjp1), sea_iceS(iip1,jjp1)
    104104      real du_nonoro_gwdS(iip1,jjp1,llm),dv_nonoro_gwdS(iip1,jjp1,llm)
     
    134134      real, dimension(:,:,:), allocatable :: tslabold
    135135      real, dimension(:,:), allocatable :: rnatold,pctsrf_sicold
    136       real, dimension(:,:), allocatable :: tsea_iceold,sea_iceold
     136      real, dimension(:,:), allocatable :: tsea_iceold,ticeold
     137      real, dimension(:,:), allocatable :: sea_iceold
    137138      real,allocatable :: du_nonoro_gwdold(:,:,:)
    138139      real,allocatable :: dv_nonoro_gwdold(:,:,:)
     
    292293      allocate(pctsrf_sicold(imold+1,jmold+1))
    293294      allocate(tsea_iceold(imold+1,jmold+1))
     295      allocate(ticeold(imold+1,jmold+1))
    294296      allocate(sea_iceold(imold+1,jmold+1))
    295297
     
    742744      IF (ierr .NE. NF_NOERR) THEN
    743745         PRINT*, "lect_start_archive: Failed loading <tsea_ice>"
     746      ENDIF
     747c     
     748      ierr = NF_INQ_VARID (nid, "tice", nvarid)
     749      IF (ierr .NE. NF_NOERR) THEN
     750         PRINT*, "lect_start_archive: Field <tice> not found"
     751      ENDIF
     752#ifdef NC_DOUBLE
     753      ierr = NF_GET_VARA_DOUBLE(nid, nvarid,start,count,ticeold)
     754#else
     755      ierr = NF_GET_VARA_REAL(nid, nvarid,start,count,ticeold)
     756#endif
     757      IF (ierr .NE. NF_NOERR) THEN
     758         PRINT*, "lect_start_archive: Failed loading <tice>"
    744759      ENDIF
    745760c
     
    13141329     &                   rlonuold,rlatvold,rlonu,rlatv)
    13151330      call gr_dyn_fi (1,iim+1,jjm+1,ngrid,tsea_ices,tsea_ice)
     1331
     1332      call interp_horiz (ticeold,tices,imold,jmold,iim,jjm,1,
     1333     &                   rlonuold,rlatvold,rlonu,rlatv)
     1334      call gr_dyn_fi (1,iim+1,jjm+1,ngrid,tices,tice)
    13161335
    13171336      call interp_horiz (sea_iceold,sea_ices,imold,jmold,iim,jjm,1,
     
    15221541      deallocate(pctsrf_sicold)
    15231542      deallocate(tsea_iceold)
     1543      deallocate(ticeold)
    15241544      deallocate(sea_iceold)
    15251545
  • trunk/LMDZ.GENERIC/libf/dynphy_lonlat/phystd/newstart.F

    r3341 r3423  
    128128      REAL,ALLOCATABLE :: tslab(:,:) ! slab ocean temperature
    129129      REAL pctsrf_sic(ngridmx) ! sea ice cover
    130       REAL tsea_ice(ngridmx) ! temperature sea_ice
     130      REAL tsea_ice(ngridmx) ! sea ice temperature
     131      REAL tice(ngridmx) ! top layer temp, snow/seaice
    131132      REAL sea_ice(ngridmx) ! mass of sea ice (kg/m2)
    132133
     
    367368     &        tsurf,tsoil,emis,spectral_albedofi,q2,qsurf,
    368369     &        cloudfrac,totalfrac,hice,rnat,pctsrf_sic,tslab,tsea_ice,
    369      &        sea_ice)
     370     &        tice,sea_ice)
    370371
    371372        ! copy albedo and soil thermal inertia on (local) physics grid
     
    570571     &   t,ucov,vcov,ps,teta,phisold_newgrid,q,qsurf,
    571572     &   surfith,nid,
    572      &   rnat,pctsrf_sic,tslab,tsea_ice,sea_ice,
     573     &   rnat,pctsrf_sic,tslab,tsea_ice,tice,sea_ice,
    573574     &   du_nonoro_gwd,dv_nonoro_gwd,east_gwstress,west_gwstress)
    574575        write(*,*) "OK, read start_archive file"
     
    620621      write(*,*) 'noacglac  : H2O ice across Noachis Terra'
    621622      write(*,*) 'oborealis : H2O ice across Vastitas Borealis'
     623      write(*,*) 'oceanball : H2O liquid everywhere'
    622624      write(*,*) 'iceball   : Thick ice layer all over surface'
    623625      write(*,*) 'supercontinent: Create a continent of given Ab and TI'
     
    16771679     &                tsurf,tsoil,emis,spectral_albedofi,q2,qsurf,
    16781680     &                cloudfrac,totalfrac,hice,
    1679      &                rnat,pctsrf_sic,tslab,tsea_ice,sea_ice)
     1681     &                rnat,pctsrf_sic,tslab,tsea_ice,tice,sea_ice)
    16801682
    16811683c=======================================================================
  • trunk/LMDZ.GENERIC/libf/dynphy_lonlat/phystd/start2archive.F

    r3335 r3423  
    8787      REAL rnat(ngridmx),pctsrf_sic(ngridmx),sea_ice(ngridmx)
    8888      REAL, ALLOCATABLE :: tslab(:,:)
    89       REAL tsea_ice(ngridmx)
     89      REAL tsea_ice(ngridmx),tice(ngridmx)
    9090
    9191
     
    109109      REAL rnatS(ip1jmp1),pctsrf_sicS(ip1jmp1),sea_iceS(ip1jmp1)
    110110      REAL, ALLOCATABLE :: tslabS(:,:)
    111       REAL tsea_iceS(ip1jmp1)
     111      REAL tsea_iceS(ip1jmp1),ticeS(ip1jmp1)
    112112
    113113!     For non-orographic GW
     
    250250!       change FF 05/2011
    251251     .       cloudfrac,totalcloudfrac,hice,
    252 !       change BC 05/2014
    253      .       rnat,pctsrf_sic,tslab,tsea_ice,sea_ice)
     252!       change BC 05/2014, SB 08/2024
     253     .       rnat,pctsrf_sic,tslab,tsea_ice,tice,sea_ice)
    254254
    255255
     
    363363      call gr_fi_dyn(1,ngridmx,iip1,jjp1,pctsrf_sic,pctsrf_sicS)
    364364      call gr_fi_dyn(1,ngridmx,iip1,jjp1,tsea_ice,tsea_iceS)
     365      call gr_fi_dyn(1,ngridmx,iip1,jjp1,tice,ticeS)
    365366      call gr_fi_dyn(1,ngridmx,iip1,jjp1,sea_ice,sea_iceS)
    366367      call gr_fi_dyn(nslay,ngridmx,iip1,jjp1,tslab,tslabS)
     
    550551          call write_archive(nid,ntime,'tsea_ice'
    551552     &            ,'tsea_ice','',2,tsea_iceS)
     553          call write_archive(nid,ntime,'tice'
     554     &            ,'tice','',2,ticeS)
    552555        endif !ok_slab_ocean
    553556     
Note: See TracChangeset for help on using the changeset viewer.