Changeset 444 for LMDZ.3.3/branches/rel-LF/libf/phylmd

Timestamp:

Jan 27, 2003, 11:57:08 AM (22 years ago)

Author:

lmdzadmin

Message:

Rajout qsol, slp JLD/FH
IM

File:

• LMDZ.3.3/branches/rel-LF/libf/phylmd/physiq.F (modified) (8 diffs)

LMDZ.3.3/branches/rel-LF/libf/phylmd/physiq.F

@@ -267 +267 @@
       SAVE deltat                 ! ecart avec la SST de reference
 c
-      REAL fqsol(klon,nbsrf)
-      SAVE fqsol                  ! humidite du sol
+      REAL fqsurf(klon,nbsrf)
+      SAVE fqsurf                 ! humidite de l'air au contact de la surface
+c
+      REAL qsol(klon)
+      SAVE qsol                  ! hauteur d'eau dans le sol
 c
       REAL fsnow(klon,nbsrf)
@@ -502 +505 @@
       REAL cldt(klon),cldq(klon) !nuage total, eau liquide integree
 c
-      REAL zxtsol(klon), zxqsol(klon), zxsnow(klon), zxfluxlat(klon)
+      REAL zxtsol(klon), zxqsurf(klon), zxsnow(klon), zxfluxlat(klon)
 c
       REAL dist, rmu0(klon), fract(klon)
@@ -746 +749 @@
 c
          CALL phyetat0 ("startphy.nc",dtime,co2_ppm_etat0,solaire_etat0,
-     .       rlat,rlon,pctsrf, ftsol,ftsoil,deltat,fqsol,fsnow,
+     .       rlat,rlon,pctsrf, ftsol,ftsoil,deltat,fqsurf,fqsol,fsnow,
      .       falbe, fevap, rain_fall,snow_fall,solsw, sollwdown,
      .       dlw,radsol,frugs,agesno,clesphy0,
@@ -1099 +1102 @@
      e            julien, rmu0, 
      e            ok_veget, ocean, npas, nexca, ftsol,
-     $            soil_model,ftsoil,
-     $            paprs,pplay,radsol, fsnow,fqsol,fevap,falbe,falblw,
+     $            soil_model,ftsoil, qsol, 
+     $            paprs,pplay,radsol, fsnow,fqsurf,fevap,falbe,falblw,
      $            fluxlat,
 cIM cf. JLD  e            rain_fall, snow_fall, solsw, sollw, sollwdown, fder,
@@ -1186 +1189 @@
 cIM cf. JLD
             wfbils(i,nsrf) = ( fsolsw(i,nsrf) + fsollw(i,nsrf)
-     $         + fluxt(i,nsrf) + fluxlat(i,nsrf) ) * pctsrf(i,nsrf)
+     $         + fluxt(i,1,nsrf) + fluxlat(i,nsrf) ) * pctsrf(i,nsrf)
             zxtsol(i) = zxtsol(i) + ftsol(i,nsrf)*pctsrf(i,nsrf)
             zxfluxlat(i) = zxfluxlat(i) + fluxlat(i,nsrf)*pctsrf(i,nsrf)
@@ -1718 +1721 @@
 c
 c      CALL hydrol(dtime,pctsrf,rain_fall, snow_fall, zxevap,
-c     .            agesno, ftsol,fqsol,fsnow, ruis)
-c
-      DO i = 1, klon
-         zxqsol(i) = 0.0
+c     .            agesno, ftsol,fqsurf,fsnow, ruis)
+c
+      DO i = 1, klon
+         zxqsurf(i) = 0.0
          zxsnow(i) = 0.0
       ENDDO
       DO nsrf = 1, nbsrf
       DO i = 1, klon
-         zxqsol(i) = zxqsol(i) + fqsol(i,nsrf)*pctsrf(i,nsrf)
+         zxqsurf(i) = zxqsurf(i) + fqsurf(i,nsrf)*pctsrf(i,nsrf)
          zxsnow(i) = zxsnow(i) + fsnow(i,nsrf)*pctsrf(i,nsrf)
       ENDDO
@@ -1736 +1739 @@
 cXXX      DO i = 1, klon
 cXXX         IF (pctsrf(i,nsrf).LT.epsfra) THEN
-cXXX            fqsol(i,nsrf) = zxqsol(i)
+cXXX            fqsurf(i,nsrf) = zxqsurf(i)
 cXXX            fsnow(i,nsrf) = zxsnow(i)
 cXXX         ENDIF
@@ -1993 +1996 @@
 ccc         IF (ok_oasis) CALL quitcpl
          CALL phyredem ("restartphy.nc",dtime,radpas,
-     .      rlat, rlon, pctsrf, ftsol, ftsoil, deltat, fqsol, fsnow,
-     .      falbe, fevap, rain_fall, snow_fall,
+     .      rlat, rlon, pctsrf, ftsol, ftsoil, deltat, fqsurf, qsol,
+     .      fsnow, falbe, fevap, rain_fall, snow_fall,
      .      solsw, sollwdown,dlw,
      .      radsol,frugs,agesno,