Changeset 3878 for trunk

Timestamp:

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

Author:

tbertrand

Message:

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

Location:

trunk/LMDZ.PLUTO/libf/dynphy_lonlat/phypluto

Files:

• datareadnc.F (modified) (3 diffs)
• grid_noro1.F (modified) (2 diffs)
• lect_start_archive.F (modified) (8 diffs)
• newstart.F (modified) (2 diffs)

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

@@ -58 +58 @@
 
       INTEGER    imd,jmd,imdp1,jmdp1
-      parameter    (imd=1440,jmd=719,imdp1=1441,jmdp1=720)
+      !parameter    (imd=1440,jmd=719,imdp1=1441,jmdp1=720)
+      parameter    (imd=360,jmd=179,imdp1=361,jmdp1=180)
 
       INTEGER    iimp1
@@ -85 +86 @@
 
       INTEGER klatdat,ngridmixgdat
-      PARAMETER (klatdat=720,ngridmixgdat=1440)
+      PARAMETER (klatdat=jmd+1,ngridmixgdat=imd)
 
 c    on passe une grille en rlonu rlatv et im+1 jm a interp_horiz)
@@ -221 +222 @@
           latitude(:)=(pi/180.)*latitude(:)
 
-          call grid_noro1(1440,720, longitude, latitude, zdata,
+          call grid_noro1(imd,jmd+1, longitude, latitude, zdata,
      .         iim, jjp1, rlonv, rlatu, zmea,zstd,zsig,zgam,zthe)
 

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

@@ -56 +56 @@
        implicit none
        integer iusn,jusn,iext
-       parameter(iusn=1440,jusn=720,iext=40)
+       !parameter(iusn=1440,jusn=720,iext=40)
+       parameter(iusn=360,jusn=180,iext=40)
 c!-*-      include 'param1'
 c!-*-      include 'comcstfi.h'
@@ -113 +114 @@
 
       IF(imdep.ne.iusn.or.jmdep.ne.jusn)then
-         print *,' imdep ou jmdep mal dimensionnes:',imdep,jmdep
+         print *,' imdep ou jmdep mal dimensionnes:',imdep,jmdep,
+     &                    iusn,jusn
          call abort
       ENDIF

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)

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

@@ -567 +567 @@
         CALL gr_dyn_fi(1,iip1,jjp1,ngridmx,phis,phisfi)
         CALL gr_dyn_fi(1,iip1,jjp1,ngridmx,surfith,surfithfi)
-        CALL gr_dyn_fi(1,iip1,jjp1,ngridmx,alb,albfi)
+        CALL gr_dyn_fi(1,iip1,jjp1,ngridmx,alb,albfi) ! erased at next step
         CALL gr_dyn_fi(1,iip1,jjp1,ngridmx,zmeaS,zmea)
         CALL gr_dyn_fi(1,iip1,jjp1,ngridmx,zstdS,zstd)
@@ -584 +584 @@
         write(*,*) 'Reading file START_ARCHIVE'
         CALL lect_start_archive(ngridmx,llm,
-     &   date,tsurf,tsoil,emis,q2,
+     &   date,tsurf,tsoil,emis,albfi,q2,
      &   t,ucov,vcov,ps,teta,phisold_newgrid,q,qsurf,
      &   surfith,nid)