Changeset 3423 for trunk/LMDZ.GENERIC/libf/dynphy_lonlat/phystd/lect_start_archive.F

Timestamp:

Sep 2, 2024, 5:20:13 PM (6 months ago)

Author:

bhatnags

Message:

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

File:

• trunk/LMDZ.GENERIC/libf/dynphy_lonlat/phystd/lect_start_archive.F (modified) (8 diffs)

trunk/LMDZ.GENERIC/libf/dynphy_lonlat/phystd/lect_start_archive.F

@@ -3 +3 @@
      &     t,ucov,vcov,ps,h,phisold_newgrid,
      &     q,qsurf,surfith,nid,
-     &     rnat,pctsrf_sic,tslab,tsea_ice,sea_ice,
+     &     rnat,pctsrf_sic,tslab,tsea_ice,tice,sea_ice,
      &     du_nonoro_gwd,dv_nonoro_gwd,east_gwstress,west_gwstress)
 
@@ -71 +71 @@
       REAL,INTENT(OUT) :: tslab(ngrid,nslay)
       REAL ,INTENT(OUT) ::rnat(ngrid),pctsrf_sic(ngrid)
-      REAL,INTENT(OUT) :: tsea_ice(ngrid),sea_ice(ngrid)
+      REAL,INTENT(OUT) :: tsea_ice(ngrid),tice(ngrid),sea_ice(ngrid)
 c     REAL phisfi(ngrid)
       REAL,INTENT(OUT):: du_nonoro_gwd(ngrid,llm)
@@ -100 +100 @@
       REAL q2S(iip1,jjp1,llm+1),qsurfS(iip1,jjp1,nqtot)
       real tslabS(iip1,jjp1,nslay)
-      real pctsrf_sicS(iip1,jjp1),tsea_iceS(iip1,jjp1)
+      real pctsrf_sicS(iip1,jjp1),tsea_iceS(iip1,jjp1),ticeS(iip1,jjp1)
       real rnatS(iip1,jjp1), sea_iceS(iip1,jjp1)
       real du_nonoro_gwdS(iip1,jjp1,llm),dv_nonoro_gwdS(iip1,jjp1,llm)
@@ -134 +134 @@
       real, dimension(:,:,:), allocatable :: tslabold
       real, dimension(:,:), allocatable :: rnatold,pctsrf_sicold
-      real, dimension(:,:), allocatable :: tsea_iceold,sea_iceold
+      real, dimension(:,:), allocatable :: tsea_iceold,ticeold
+      real, dimension(:,:), allocatable :: sea_iceold
       real,allocatable :: du_nonoro_gwdold(:,:,:)
       real,allocatable :: dv_nonoro_gwdold(:,:,:)
@@ -292 +293 @@
       allocate(pctsrf_sicold(imold+1,jmold+1))
       allocate(tsea_iceold(imold+1,jmold+1))
+      allocate(ticeold(imold+1,jmold+1))
       allocate(sea_iceold(imold+1,jmold+1))
 
@@ -742 +744 @@
       IF (ierr .NE. NF_NOERR) THEN
          PRINT*, "lect_start_archive: Failed loading <tsea_ice>"
+      ENDIF
+c      
+      ierr = NF_INQ_VARID (nid, "tice", nvarid)
+      IF (ierr .NE. NF_NOERR) THEN
+         PRINT*, "lect_start_archive: Field <tice> not found"
+      ENDIF
+#ifdef NC_DOUBLE
+      ierr = NF_GET_VARA_DOUBLE(nid, nvarid,start,count,ticeold)
+#else
+      ierr = NF_GET_VARA_REAL(nid, nvarid,start,count,ticeold)
+#endif
+      IF (ierr .NE. NF_NOERR) THEN
+         PRINT*, "lect_start_archive: Failed loading <tice>"
       ENDIF
 c
@@ -1314 +1329 @@
      &                   rlonuold,rlatvold,rlonu,rlatv)
       call gr_dyn_fi (1,iim+1,jjm+1,ngrid,tsea_ices,tsea_ice)
+
+      call interp_horiz (ticeold,tices,imold,jmold,iim,jjm,1,
+     &                   rlonuold,rlatvold,rlonu,rlatv)
+      call gr_dyn_fi (1,iim+1,jjm+1,ngrid,tices,tice)
 
       call interp_horiz (sea_iceold,sea_ices,imold,jmold,iim,jjm,1,
@@ -1522 +1541 @@
       deallocate(pctsrf_sicold)
       deallocate(tsea_iceold)
+      deallocate(ticeold)
       deallocate(sea_iceold)