Changeset 3878 for trunk/LMDZ.PLUTO/libf/dynphy_lonlat/phypluto/lect_start_archive.F

Timestamp:

Aug 5, 2025, 1:29:32 PM (5 weeks ago)

Author:

tbertrand

Message:

LMDZ.PLUTO:
For newstart : reading in the intial albedo map and the new topography map
TB

File:

• trunk/LMDZ.PLUTO/libf/dynphy_lonlat/phypluto/lect_start_archive.F (modified) (8 diffs)

trunk/LMDZ.PLUTO/libf/dynphy_lonlat/phypluto/lect_start_archive.F

@@ -1 +1 @@
       SUBROUTINE lect_start_archive(ngrid,nlayer,
-     &     date,tsurf,tsoil,emis,q2,
+     &     date,tsurf,tsoil,emis,alb,q2,
      &     t,ucov,vcov,ps,h,phisold_newgrid,
      &     q,qsurf,surfith,nid)
@@ -68 +68 @@
       REAL n2ice(ngrid) ! N2 ice layer
       REAL,INTENT(OUT) :: emis(ngrid)
+      REAL,INTENT(OUT) :: alb(ngrid)
       REAL,INTENT(OUT) :: q2(ngrid,llm+1),qsurf(ngrid,nqtot)
       ! REAL,INTENT(OUT) :: tslab(ngrid,nslay)
@@ -100 +101 @@
       real n2iceS(iip1,jjp1)
       real emisS(iip1,jjp1)
+      real albS(iip1,jjp1)
       REAL q2S(iip1,jjp1,llm+1),qsurfS(iip1,jjp1,nqtot)
       ! real tslabS(iip1,jjp1,nslay)
@@ -135 +137 @@
       real, dimension(:,:), allocatable :: tsurfold
       real, dimension(:,:), allocatable :: emisold
+      real, dimension(:,:), allocatable :: albold
       real, dimension(:,:,:,:), allocatable :: qold
       ! real, dimension(:,:,:), allocatable :: tslabold
@@ -284 +287 @@
       allocate(tsurfold(imold+1,jmold+1))
       allocate(emisold(imold+1,jmold+1))
+      allocate(albold(imold+1,jmold+1))
       allocate(q2old(imold+1,jmold+1,lmold+1))
 !      allocate(tsoilold(imold+1,jmold+1,nsoilmx))
@@ -645 +649 @@
       IF (ierr .NE. NF_NOERR) THEN
          PRINT*, "lect_start_archive: Failed loading <emis>"
+         CALL abort
+      ENDIF
+c
+      ierr = NF_INQ_VARID (nid, "albedodat", nvarid)
+      IF (ierr .NE. NF_NOERR) THEN
+         PRINT*, "lect_start_archive: Field <albedodat> not found"
+         CALL abort
+      ENDIF
+#ifdef NC_DOUBLE
+      ierr = NF_GET_VARA_DOUBLE(nid, nvarid,start,count,albold)
+#else
+      ierr = NF_GET_VARA_REAL(nid, nvarid,start,count,albold)
+#endif
+      IF (ierr .NE. NF_NOERR) THEN
+         PRINT*, "lect_start_archive: Failed loading <albedodat>"
          CALL abort
       ENDIF
@@ -1086 +1105 @@
       call gr_dyn_fi (1,iim+1,jjm+1,ngrid,emiss,emis)
 c     write(46,*) 'emis',emis
+c Albedo de la surface
+      call interp_horiz (albold,albs,imold,jmold,iim,jjm,1,
+     &                   rlonuold,rlatvold,rlonu,rlatv)
+      call gr_dyn_fi (1,iim+1,jjm+1,ngrid,albs,alb)
+c     write(46,*) 'alb',alb
 
 
@@ -1503 +1527 @@
       deallocate(tsurfold)
       deallocate(emisold)
+      deallocate(albold)
       deallocate(q2old)
       deallocate(tsoilold)