Changeset 442 for LMDZ.3.3

Timestamp:

Jan 27, 2003, 11:45:56 AM (21 years ago)

Author:

lmdzadmin

Message:

Rajout qsol JLD
IM/LF

File:

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

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

@@ -3 +3 @@
 c
       SUBROUTINE phyetat0 (fichnom,dtime,co2_ppm_etat0,solaire_etat0,
-     .            rlat,rlon, pctsrf, tsol,tsoil,deltat,qsol,snow,
+     .            rlat,rlon, pctsrf, tsol,tsoil,deltat,qsurf,qsol,snow,
      .           albe, evap, rain_fall, snow_fall, solsw, sollw,
      .           fder,radsol,frugs,agesno,clesphy0,
@@ -30 +30 @@
       REAL tsoil(klon,nsoilmx,nbsrf)
       REAL deltat(klon)
-      REAL qsol(klon,nbsrf)
+      REAL qsurf(klon,nbsrf)
+      REAL qsol(klon)
       REAL snow(klon,nbsrf)
       REAL albe(klon,nbsrf)
@@ -436 +437 @@
       PRINT*,'Ecart de la SST deltat:', xmin, xmax
 c
-c Lecture de l'humidite du sol:
+c Lecture de l'humidite de l'air juste au dessus du sol:
 c
       ierr = NF_INQ_VARID (nid, "QS", nvarid)
@@ -454 +455 @@
            ENDIF
 #ifdef NC_DOUBLE
-           ierr = NF_GET_VAR_DOUBLE(nid, nvarid, qsol(1,nsrf))
-#else
-           ierr = NF_GET_VAR_REAL(nid, nvarid, qsol(1,nsrf))
+           ierr = NF_GET_VAR_DOUBLE(nid, nvarid, qsurf(1,nsrf))
+#else
+           ierr = NF_GET_VAR_REAL(nid, nvarid, qsurf(1,nsrf))
 #endif
            IF (ierr.NE.NF_NOERR) THEN
@@ -465 +466 @@
            xmax = -1.0E+20
            DO i = 1, klon
-              xmin = MIN(qsol(i,nsrf),xmin)
-              xmax = MAX(qsol(i,nsrf),xmax)
+              xmin = MIN(qsurf(i,nsrf),xmin)
+              xmax = MAX(qsurf(i,nsrf),xmax)
            ENDDO
-           PRINT*,'Humidite du sol QS**:', nsrf, xmin, xmax
+           PRINT*,'Humidite pres du sol QS**:', nsrf, xmin, xmax
          ENDDO
       ELSE
@@ -474 +475 @@
          PRINT*, '          J ignore donc les autres humidites QS**'
 #ifdef NC_DOUBLE
-         ierr = NF_GET_VAR_DOUBLE(nid, nvarid, qsol(1,1))
-#else
-         ierr = NF_GET_VAR_REAL(nid, nvarid, qsol(1,1))
+         ierr = NF_GET_VAR_DOUBLE(nid, nvarid, qsurf(1,1))
+#else
+         ierr = NF_GET_VAR_REAL(nid, nvarid, qsurf(1,1))
 #endif
          IF (ierr.NE.NF_NOERR) THEN
@@ -485 +486 @@
          xmax = -1.0E+20
          DO i = 1, klon
-            xmin = MIN(qsol(i,1),xmin)
-            xmax = MAX(qsol(i,1),xmax)
-         ENDDO
-         PRINT*,'Humidite du sol <QS>', xmin, xmax
+            xmin = MIN(qsurf(i,1),xmin)
+            xmax = MAX(qsurf(i,1),xmax)
+         ENDDO
+         PRINT*,'Humidite pres du sol <QS>', xmin, xmax
          DO nsrf = 2, nbsrf
          DO i = 1, klon
-            qsol(i,nsrf) = qsol(i,1)
-         ENDDO
-         ENDDO
-      ENDIF
+            qsurf(i,nsrf) = qsurf(i,1)
+         ENDDO
+         ENDDO
+      ENDIF
+C
+C Eau dans le sol (pour le modele de sol "bucket")
+C
+      ierr = NF_INQ_VARID (nid, "QSOL", nvarid)
+      IF (ierr .EQ.  NF_NOERR) THEN
+#ifdef NC_DOUBLE
+          ierr = NF_GET_VAR_DOUBLE(nid, nvarid, qsol)
+#else
+          ierr = NF_GET_VAR_REAL(nid, nvarid, qsol)
+#endif
+          IF (ierr.NE.NF_NOERR) THEN
+              PRINT*, 'phyetat0: Lecture echouee pour <QSOL>'
+              CALL abort
+          ENDIF
+      else
+          PRINT*, 'phyetat0: Le champ <QSOL> est absent'
+          PRINT*, '          Valeur par defaut nulle'
+          qsol(:)=0.
+c$$$         CALL abort
+      ENDIF
+      xmin = 1.0E+20
+      xmax = -1.0E+20
+      DO i = 1, klon
+        xmin = MIN(qsol(i),xmin)
+        xmax = MAX(qsol(i),xmax)
+      ENDDO
+      PRINT*,'Eau dans le sol (mm) <QSOL>', xmin, xmax
 c
 c Lecture de neige au sol: