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

Timestamp:

May 18, 2024, 8:21:24 PM (13 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/phystd/dyn1d

Files:

• kcm1d.F90 (modified) (1 diff)
• rcm1d.F (modified) (4 diffs)

trunk/LMDZ.GENERIC/libf/phystd/dyn1d/kcm1d.F90

@@ -106 +106 @@
   LOGICAL :: moderntracdef=.false. ! JVO, YJ : modern traceur.def
 
+  character(len=100) :: dt_file
+  integer :: ios
+  integer :: k
+  
   ! --------------
   ! Initialisation

trunk/LMDZ.GENERIC/libf/phystd/dyn1d/rcm1d.F

@@ -7 +7 @@
       use infotrac, only: nqtot, tname
       use tracer_h, only: noms, is_condensable
+      use radinc_h, only : L_NSPECTV
       use surfdat_h, only: albedodat, phisfi, dryness,
      &                     zmea, zstd, zsig, zgam, zthe,
@@ -94 +95 @@
       integer :: i_h2o_ice=0     ! tracer index of h2o ice
       integer :: i_h2o_vap=0     ! tracer index of h2o vapor
-      REAL emis(1)               ! surface layer
+      REAL emis(1)               ! emissivity of surface
+      real :: albedo(1,L_NSPECTV) ! surface albedo in various spectral bands
       REAL q2(llm+1)             ! Turbulent Kinetic Energy
       REAL zlay(llm)             ! altitude estimee dans les couches (km)
@@ -883 +885 @@
       call getin("albedo",albedodat(1))
       write(*,*) " albedo = ",albedodat(1)
+      ! Initialize surface albedo to that of bare ground
+      albedo(1,:)=albedodat(1)
 
       inertiedat(1,1)=400 ! default value for inertiedat
@@ -954 +958 @@
       call physdem1("startfi.nc",nsoilmx,1,llm,nq,
      &                dtphys,time,
-     &                tsurf,tsoil,emis,q2,qsurf,
+     &                tsurf,tsoil,emis,albedo,q2,qsurf,
      &                cloudfrac,totcloudfrac,hice, 
      &                rnat,pctsrf_sic,tslab,tsea_ice,sea_ice)