Ignore:
Timestamp:
Aug 6, 2003, 4:51:49 PM (21 years ago)
Author:
lmdzadmin
Message:

Choix sur le schéma d'advection d'eau et sur les seuils
LF

Location:
LMDZ.3.3/branches/rel-LF/libf/dyn3d
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • LMDZ.3.3/branches/rel-LF/libf/dyn3d/etat0_netcdf.F

    r437 r468  
    3636      REAL :: tsol(klon), qsol(klon), sn(klon)
    3737      REAL :: tsolsrf(klon,nbsrf), qsurf(klon,nbsrf),snsrf(klon,nbsrf)
    38       REAL :: albe(klon,nbsrf), evap(klon,nbsrf)
     38      REAL :: albe(klon,nbsrf), evap(klon,nbsrf), alblw(klon,nbsrf)
    3939      REAL :: tsoil(klon,nsoilmx,nbsrf)
    4040      REAL :: radsol(klon),rain_fall(klon), snow_fall(klon)
     
    634634      albe(:,is_oce) = 0.5
    635635      albe(:,is_sic) = 0.6
     636      alblw(:,is_ter) = albe(:,is_ter)
     637      alblw(:,is_lic) = albe(:,is_lic)
     638      alblw(:,is_oce) = albe(:,is_oce)
     639      alblw(:,is_sic) = albe(:,is_sic)
     640
    636641      evap(:,:) = 0.
    637642C     vapeur d'eau a la surface =- vapeur d'eau premier niveau
     
    661666      frugs(1:klon,is_ter) = MAX(1.0e-05, zstd(1:klon)*zsig(1:klon)/2.0)
    662667      frugs(1:klon,is_lic) = MAX(1.0e-05, zstd(1:klon)*zsig(1:klon)/2.0)
    663       frugs(1:klon,is_sic) = 0.001
     668      frugs(1:klon,is_sic) = 0.002
    664669      fder = 0.0
    665670      clwcon = 0.0
     
    670675      call phyredem("startphy.nc",phystep,radpas,
    671676     $    latfi, lonfi, pctsrf, tsolsrf, tsoil, deltat, qsurf, qsol,
    672      $    snsrf, albe, evap, rain_fall, snow_fall, solsw, sollw, fder,
     677     $    snsrf, albe, alblw,
     678     $    evap, rain_fall, snow_fall, solsw, sollw, fder,
    673679     $    radsol, frugs,  agesno,
    674680     $    zmea, zstd, zsig, zgam, zthe, zpic, zval, rugsrel,
  • LMDZ.3.3/branches/rel-LF/libf/dyn3d/gcm.F

    r462 r468  
    185185c         et eventuellement      , iq = 3, nqmx pour les autres traceurs
    186186c
    187 c     iadv(1): choix pour l'eau vap. et  iadv(2) : choix pour l'eau liq.
    188 c
    189187      DO iq = 1, nqmx
    190188       iadv( iq ) = 3
    191189      ENDDO
     190c
     191c     iadv(1): choix pour l'eau vap. et  iadv(2) : choix pour l'eau liq.
     192      iadv( 1 ) = 4
    192193c
    193194       IF( iadv(1).EQ.1 ) PRINT *,' Choix du shema humidite specifique'
Note: See TracChangeset for help on using the changeset viewer.