Changeset 3771

Timestamp:

Sep 30, 2020, 5:02:50 PM (4 years ago)

Author:

lguez

Message:

Bug fix: initial value of surface humidity

Rename qsolsrf to qsurf as it should contain surface humidity. Fix the
value, which should be between 0 and 1. A better initial value would
be saturation humidity but we do not have yet the surface pressure in
etat0phys_netcdf, it will be defined later in etat0dyn_netcdf. The
old value, 150, probably came from a confusion with soil water
content. This initial value of surface humidity is used in the first
time step of physics to compute a virtual temperature at the surface,
in cdrag.

File:

• LMDZ6/trunk/libf/dynphy_lonlat/phylmd/etat0phys_netcdf.F90 (modified) (4 diffs)

LMDZ6/trunk/libf/dynphy_lonlat/phylmd/etat0phys_netcdf.F90

@@ -107 +107 @@
   REAL, DIMENSION(SIZE(masque,1),SIZE(masque,2)) :: masque_tmp,phiso
   REAL, DIMENSION(klon)               :: sn, rugmer, run_off_lic_0, fder
-  REAL, DIMENSION(klon,nbsrf)         :: qsolsrf, snsrf
+  REAL, DIMENSION(klon,nbsrf)         :: qsurf, snsrf
   REAL, DIMENSION(klon,nsoilmx,nbsrf) :: tsoil
 
@@ -192 +192 @@
   WRITE(lunout,*)'phystep =', phystep, radpas
 
-! Init: ftsol, snsrf, qsolsrf, tsoil, rain_fall, snow_fall, solsw, sollw, z0
+! Init: ftsol, snsrf, qsurf, tsoil, rain_fall, snow_fall, solsw, sollw, z0
 !*******************************************************************************
   DO i=1,nbsrf; ftsol(:,i) = tsol; END DO
@@ -209 +209 @@
 
   fevap(:,:) = 0.
-  DO i=1,nbsrf; qsolsrf(:,i)=150.; END DO
+  qsurf = 0.
   DO i=1,nbsrf; DO j=1,nsoilmx; tsoil(:,j,i) = tsol; END DO; END DO
   rain_fall  = 0.
@@ -271 +271 @@
 
   CALL fonte_neige_init(run_off_lic_0)
-  CALL pbl_surface_init( fder, snsrf, qsolsrf, tsoil )
+  CALL pbl_surface_init( fder, snsrf, qsurf, tsoil )
   CALL phyredem( "startphy.nc" )