Changeset 3341 for trunk/LMDZ.GENERIC

Timestamp:

May 23, 2024, 8:32:15 AM (6 months ago)

Author:

emillour

Message:

Generic PCM:
Follow-up of including surface albedo in restartfi.nc
Ensure initialization of this albedo to bare ground albedo in newstart
EM

Location:

trunk/LMDZ.GENERIC

Files:

• changelog.txt (modified) (1 diff)
• libf/dynphy_lonlat/phystd/newstart.F (modified) (3 diffs)

trunk/LMDZ.GENERIC/changelog.txt

@@ -1937 +1937 @@
 Turned "surfini" in a module in the process.
 Unrelated: added missing delarations in kcm1d so it compiles one again.
+
+== 23/05/2024 == EM
+Follow-up of including surface albedo in restartfi.nc
+Ensure initialization of this albedo to bare ground albedo in newstart

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

@@ -520 +520 @@
           CALL datareadnc(relief,surfacefile,phis,alb,surfith,
      &          zmeaS,zstdS,zsigS,zgamS,ztheS)
+          ! set surface albedo to bare ground albedo
+          albedodyn(:,:)=alb(:,:)
         else
         ! specific case when not using a "surface.nc" file
@@ -532 +534 @@
           read(*,*) alb(1,1)
           alb(:,:)=alb(1,1)
+          ! set surface albedo to bare ground albedo
+          albedodyn(:,:)=alb(:,:)
           
           write(*,*) "Enter value of thermal inertia of soil:"
@@ -542 +546 @@
         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,albedodyn,albedofi)
+        do i=1,L_NSPECTV
+          spectral_albedofi(:,i)=albedofi(:)
+        enddo
         CALL gr_dyn_fi(1,iip1,jjp1,ngridmx,zmeaS,zmea)
         CALL gr_dyn_fi(1,iip1,jjp1,ngridmx,zstdS,zstd)