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

Timestamp:

Oct 30, 2000, 11:57:01 AM (24 years ago)

Author:

lmdzadmin

Message:

Rajout de la derivee des fluxs dans le fichier restart de la physique
LF

File:

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

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

@@ -2 +2 @@
      .            rlat,rlon, pctsrf, tsol,tsoil,deltat,qsol,snow,
      .           albe, evap, rain_fall, snow_fall, solsw, sollw,
-     .           radsol,frugs,agesno,clesphy0,
+     .           fder,radsol,frugs,agesno,clesphy0,
      .           zmea,zstd,zsig,zgam,zthe,zpic,zval,rugsrel,tabcntr0,
      .           t_ancien,q_ancien,ancien_ok)
@@ -35 +35 @@
       REAL sollw(klon)
       real solsw(klon)
+      real fder(klon)
       REAL frugs(klon,nbsrf)
       REAL agesno(klon)
@@ -769 +770 @@
 
 c
+c Lecture derive des flux:
+c
+      ierr = NF_INQ_VARID (nid, "fder", nvarid)
+      IF (ierr.NE.NF_NOERR) THEN
+         PRINT*, 'phyetat0: Le champ <fder> est absent'
+         PRINT*, 'mis a zero'
+         fder = 0.
+      ELSE
+#ifdef NC_DOUBLE
+        ierr = NF_GET_VAR_DOUBLE(nid, nvarid, fder)
+#else
+        ierr = NF_GET_VAR_REAL(nid, nvarid, fder)
+#endif
+        IF (ierr.NE.NF_NOERR) THEN
+          PRINT*, 'phyetat0: Lecture echouee pour <fder>'
+          CALL abort
+        ENDIF
+      ENDIF
+      xmin = 1.0E+20
+      xmax = -1.0E+20
+      DO i = 1, klon
+         xmin = MIN(fder(i),xmin)
+         xmax = MAX(fder(i),xmax)
+      ENDDO
+      PRINT*,'Derive des flux fder:', xmin, xmax
+
+c
 c Lecture du rayonnement net au sol:
 c
@@ -826 +854 @@
               xmax = MAX(frugs(i,nsrf),xmax)
            ENDDO
-           PRINT*,'evap du sol RUG**:', nsrf, xmin, xmax
+           PRINT*,'rugosite du sol RUG**:', nsrf, xmin, xmax
          ENDDO
       ELSE