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

Timestamp:

May 6, 2004, 10:57:36 AM (20 years ago)

Author:

lmdzadmin

Message:

Modifications pour la fermeture en eau (fonte des glaciers) JLD, OM
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

@@ -7 +7 @@
      .           fder,radsol,frugs,agesno,clesphy0,
      .           zmea,zstd,zsig,zgam,zthe,zpic,zval,rugsrel,tabcntr0,
-     .           t_ancien,q_ancien,ancien_ok, rnebcon, ratqs,clwcon)
+     .           t_ancien,q_ancien,ancien_ok, rnebcon, ratqs,clwcon,
+     .           run_off_lic_0)
       IMPLICIT none
 c======================================================================
@@ -56 +57 @@
       REAL pctsrf(klon, nbsrf)
       REAL fractint(klon)
+      REAL run_off_lic_0(klon)
 
       REAL t_ancien(klon,klev), q_ancien(klon,klev)
@@ -1333 +1335 @@
       ENDIF
 c
+c Lecture ratqs
+c
       ierr = NF_INQ_VARID (nid, "RATQS", nvarid)
       IF (ierr.NE.NF_NOERR) THEN
@@ -1355 +1359 @@
       PRINT*,'(ecart-type) ratqs:', xmin, xmax
 c
+c Lecture run_off_lic_0
+c
+      ierr = NF_INQ_VARID (nid, "RUNOFFLIC0", nvarid)
+      IF (ierr.NE.NF_NOERR) THEN
+         PRINT*, "phyetat0: Le champ <RUNOFFLIC0> est absent"
+         PRINT*, "Depart legerement fausse. Mais je continue"
+         run_off_lic_0 = 0.
+      ELSE
+#ifdef NC_DOUBLE
+         ierr = NF_GET_VAR_DOUBLE(nid, nvarid, run_off_lic_0)
+#else
+         ierr = NF_GET_VAR_REAL(nid, nvarid, run_off_lic_0)
+#endif
+         IF (ierr.NE.NF_NOERR) THEN
+            PRINT*, "phyetat0: Lecture echouee pour <RUNOFFLIC0>"
+            CALL abort
+         ENDIF
+      ENDIF
+      xmin = 1.0E+20
+      xmax = -1.0E+20
+      xmin = MINval(run_off_lic_0)
+      xmax = MAXval(run_off_lic_0)
+      PRINT*,'(ecart-type) run_off_lic_0:', xmin, xmax
+c
 c Fermer le fichier:
 c