Ignore:
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:
1 edited

Legend:

Unmodified
Added
Removed
  • LMDZ6/trunk/libf/dynphy_lonlat/phylmd/etat0phys_netcdf.F90

    r3630 r3771  
    107107  REAL, DIMENSION(SIZE(masque,1),SIZE(masque,2)) :: masque_tmp,phiso
    108108  REAL, DIMENSION(klon)               :: sn, rugmer, run_off_lic_0, fder
    109   REAL, DIMENSION(klon,nbsrf)         :: qsolsrf, snsrf
     109  REAL, DIMENSION(klon,nbsrf)         :: qsurf, snsrf
    110110  REAL, DIMENSION(klon,nsoilmx,nbsrf) :: tsoil
    111111
     
    192192  WRITE(lunout,*)'phystep =', phystep, radpas
    193193
    194 ! Init: ftsol, snsrf, qsolsrf, tsoil, rain_fall, snow_fall, solsw, sollw, z0
     194! Init: ftsol, snsrf, qsurf, tsoil, rain_fall, snow_fall, solsw, sollw, z0
    195195!*******************************************************************************
    196196  DO i=1,nbsrf; ftsol(:,i) = tsol; END DO
     
    209209
    210210  fevap(:,:) = 0.
    211   DO i=1,nbsrf; qsolsrf(:,i)=150.; END DO
     211  qsurf = 0.
    212212  DO i=1,nbsrf; DO j=1,nsoilmx; tsoil(:,j,i) = tsol; END DO; END DO
    213213  rain_fall  = 0.
     
    271271
    272272  CALL fonte_neige_init(run_off_lic_0)
    273   CALL pbl_surface_init( fder, snsrf, qsolsrf, tsoil )
     273  CALL pbl_surface_init( fder, snsrf, qsurf, tsoil )
    274274  CALL phyredem( "startphy.nc" )
    275275
Note: See TracChangeset for help on using the changeset viewer.