Changeset 3335 for trunk/LMDZ.GENERIC/libf/dynphy_lonlat/phystd

Timestamp:

May 18, 2024, 8:21:24 PM (9 months ago)

Author:

emillour

Message:

Generic PCM:
Add reading/writing of surface albedo in (re)startfi.nc to
improve model restartability. For now only the simpler case
of non-spectral dependent surface albedo is handled.
Turned "surfini" in a module in the process.
Unrelated: added missing delarations in kcm1d so it compiles one again.
EM

Location:

trunk/LMDZ.GENERIC/libf/dynphy_lonlat/phystd

Files:

• newstart.F (modified) (5 diffs)
• start2archive.F (modified) (6 diffs)

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

@@ -21 +21 @@
       USE tracer_h, ONLY: igcm_co2_ice, igcm_h2o_vap, igcm_h2o_ice
       USE comsoil_h, ONLY: nsoilmx, layer, mlayer, inertiedat
+      USE radinc_h, only : L_NSPECTV ! number of spectral bands in the visible
       USE surfdat_h, ONLY: phisfi, albedodat,
      &                     zmea, zstd, zsig, zgam, zthe
@@ -112 +113 @@
       REAL q2(ngridmx,llm+1)
 !      REAL rnaturfi(ngridmx)
-      real alb(iip1,jjp1),albfi(ngridmx) ! albedos
+      real alb(iip1,jjp1),albfi(ngridmx) ! bare ground albedos
+      real albedodyn(iip1,jjp1),albedofi(ngridmx) ! surface albedos
+      real spectral_albedofi(ngridmx,L_NSPECTV) ! spectral surface albedo
       real,ALLOCATABLE :: ith(:,:,:),ithfi(:,:) ! thermal inertia (3D)
       real surfith(iip1,jjp1),surfithfi(ngridmx) ! surface thermal inertia (2D)
@@ -361 +364 @@
         fichnom = 'startfi.nc'
         CALL phyetat0(.true.,ngridmx,llm,fichnom,tab0,Lmodif,nsoilmx,
-     .        nqtot,day_ini,time,
-     .        tsurf,tsoil,emis,q2,qsurf,   !) ! temporary modif by RDW
-     .        cloudfrac,totalfrac,hice,rnat,pctsrf_sic,tslab,tsea_ice,
-     .        sea_ice)
+     &        nqtot,day_ini,time,
+     &        tsurf,tsoil,emis,spectral_albedofi,q2,qsurf,
+     &        cloudfrac,totalfrac,hice,rnat,pctsrf_sic,tslab,tsea_ice,
+     &        sea_ice)
 
         ! copy albedo and soil thermal inertia on (local) physics grid
         do i=1,ngridmx
           albfi(i) = albedodat(i)
+          albedofi(i)= spectral_albedofi(i,1) ! assume same albedo at all wavelenghts
           do j=1,nsoilmx
            ithfi(i,j) = inertiedat(i,j)
@@ -380 +384 @@
         ! to correctly recast things on physics grid)
         call gr_fi_dyn(1,ngridmx,iip1,jjp1,albfi,alb)
+        call gr_fi_dyn(1,ngridmx,iip1,jjp1,albedofi,albedodyn)
         call gr_fi_dyn(nsoilmx,ngridmx,iip1,jjp1,ithfi,ith)
         call gr_fi_dyn(1,ngridmx,iip1,jjp1,surfithfi,surfith)
@@ -1662 +1667 @@
       call physdem1("restartfi.nc",nsoilmx,ngridmx,llm,nqtot,
      &                dtphys,real(day_ini),
-     &                tsurf,tsoil,emis,q2,qsurf,
+     &                tsurf,tsoil,emis,spectral_albedofi,q2,qsurf,
      &                cloudfrac,totalfrac,hice,
      &                rnat,pctsrf_sic,tslab,tsea_ice,sea_ice)

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

@@ -21 +21 @@
       use infotrac, only: infotrac_init, nqtot, tname
       USE comsoil_h
-      
+      USE radinc_h, only : L_NSPECTV ! number of spectral bands in the visible
 !      use slab_ice_h, only: noceanmx
       USE ocean_slab_mod, ONLY: nslay
@@ -73 +73 @@
       REAL,ALLOCATABLE :: qsurf(:,:)
       REAL emis(ngridmx)
+      REAL :: albedo(ngridmx,L_NSPECTV) ! spectral surface albedo
       INTEGER start,length
       PARAMETER (length = 100)
@@ -99 +100 @@
       REAL,ALLOCATABLE :: qsurfS(:,:)
       REAL emisS(ip1jmp1)
+      REAL :: albedoS(ngridmx) ! surface albedo assumed same at all wavelengths
 
 !     added by FF for cloud fraction setup
@@ -245 +247 @@
       CALL phyetat0(.true.,ngridmx,llm,fichnom,0,Lmodif,nsoilmx,nqtot,
      .      day_ini_fi,timefi,
-     .      tsurf,tsoil,emis,q2,qsurf,
+     .      tsurf,tsoil,emis,albedo,q2,qsurf,
 !       change FF 05/2011
      .       cloudfrac,totalcloudfrac,hice,
@@ -351 +353 @@
       call gr_fi_dyn(nsoilmx,ngridmx,iip1,jjp1,inertiedat,ithS)
       call gr_fi_dyn(1,ngridmx,iip1,jjp1,emis,emisS)
+      call gr_fi_dyn(1,ngridmx,iip1,jjp1,albedo(1,1),albedoS)
       call gr_fi_dyn(llm+1,ngridmx,iip1,jjp1,q2,q2S)
       call gr_fi_dyn(nqtot,ngridmx,iip1,jjp1,qsurf,qsurfS)
@@ -455 +458 @@
 !     &  'kg/m2',2,co2iceS)
       call write_archive(nid,ntime,'emis','grd emis',' ',2,emisS)
+      call write_archive(nid,ntime,'albedo','surface albedo',' ',
+     &                   2,albedoS)
       call write_archive(nid,ntime,'ps','Psurf','Pa',2,ps)
       call write_archive(nid,ntime,'tsurf','surf T','K',2,tsurfS)