Changeset 468 for LMDZ.3.3

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:

• etat0_netcdf.F (modified) (4 diffs)
• gcm.F (modified) (1 diff)

LMDZ.3.3/branches/rel-LF/libf/dyn3d/etat0_netcdf.F

@@ -36 +36 @@
       REAL :: tsol(klon), qsol(klon), sn(klon)
       REAL :: tsolsrf(klon,nbsrf), qsurf(klon,nbsrf),snsrf(klon,nbsrf) 
-      REAL :: albe(klon,nbsrf), evap(klon,nbsrf)
+      REAL :: albe(klon,nbsrf), evap(klon,nbsrf), alblw(klon,nbsrf)
       REAL :: tsoil(klon,nsoilmx,nbsrf) 
       REAL :: radsol(klon),rain_fall(klon), snow_fall(klon)
@@ -634 +634 @@
       albe(:,is_oce) = 0.5
       albe(:,is_sic) = 0.6
+      alblw(:,is_ter) = albe(:,is_ter)
+      alblw(:,is_lic) = albe(:,is_lic)
+      alblw(:,is_oce) = albe(:,is_oce)
+      alblw(:,is_sic) = albe(:,is_sic)
+
       evap(:,:) = 0.
 C     vapeur d'eau a la surface =- vapeur d'eau premier niveau
@@ -661 +666 @@
       frugs(1:klon,is_ter) = MAX(1.0e-05, zstd(1:klon)*zsig(1:klon)/2.0)
       frugs(1:klon,is_lic) = MAX(1.0e-05, zstd(1:klon)*zsig(1:klon)/2.0)
-      frugs(1:klon,is_sic) = 0.001
+      frugs(1:klon,is_sic) = 0.002
       fder = 0.0
       clwcon = 0.0
@@ -670 +675 @@
       call phyredem("startphy.nc",phystep,radpas,
      $    latfi, lonfi, pctsrf, tsolsrf, tsoil, deltat, qsurf, qsol, 
-     $    snsrf, albe, evap, rain_fall, snow_fall, solsw, sollw, fder,
+     $    snsrf, albe, alblw,
+     $    evap, rain_fall, snow_fall, solsw, sollw, fder,
      $    radsol, frugs,  agesno, 
      $    zmea, zstd, zsig, zgam, zthe, zpic, zval, rugsrel, 

LMDZ.3.3/branches/rel-LF/libf/dyn3d/gcm.F

@@ -185 +185 @@
 c         et eventuellement      , iq = 3, nqmx pour les autres traceurs
 c
-c     iadv(1): choix pour l'eau vap. et  iadv(2) : choix pour l'eau liq.
-c
       DO iq = 1, nqmx
        iadv( iq ) = 3
       ENDDO
+c
+c     iadv(1): choix pour l'eau vap. et  iadv(2) : choix pour l'eau liq.
+      iadv( 1 ) = 4
 c
        IF( iadv(1).EQ.1 ) PRINT *,' Choix du shema humidite specifique'