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

Timestamp:

Mar 21, 2002, 4:09:32 PM (22 years ago)

Author:

lmdzadmin

Message:

2 changements pour les fichiers histoire:

• utilisation de l'entree "rectilineaire" de IOIPSL pour ne plus avoir

a

lancer ncregular a chaque fois

• le calendrier des fichiers histoire est maintenant base sur la date d'initialisation de la simulation plutot que sur la date de depart du

job

en cours

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

@@ +1 @@
+c
+c $Header$
+c
       SUBROUTINE phyetat0 (fichnom,dtime,co2_ppm,solaire,
      .            rlat,rlon, pctsrf, tsol,tsoil,deltat,qsol,snow,
@@ -16 +19 @@
 #include "dimsoil.h"
 #include "clesphys.h"
+#include "temps.h"
 c======================================================================
       CHARACTER*(*) fichnom
@@ -37 +41 @@
       real fder(klon)
       REAL frugs(klon,nbsrf)
-      REAL agesno(klon)
+      REAL agesno(klon,nbsrf)
       REAL zmea(klon)
       REAL zstd(klon)
@@ -161 +165 @@
 
       ENDIF
+
+      itau_phy = tab_cntrl(15)
+
 c
 c Lecture des latitudes (coordonnees):
@@ -874 +881 @@
             xmax = MAX(frugs(i,1),xmax)
          ENDDO
-         PRINT*,'Neige du sol <RUG>', xmin, xmax
+         PRINT*,'rugosite <RUG>', xmin, xmax
          DO nsrf = 2, nbsrf
          DO i = 1, klon
@@ -888 +895 @@
       IF (ierr.NE.NF_NOERR) THEN
          PRINT*, 'phyetat0: Le champ <AGESNO> est absent'
-         PRINT*, "          Valeur par default: 50"
-         DO i = 1, klon
-            agesno(i) = 50.0
-         ENDDO
-      ELSE
-#ifdef NC_DOUBLE
-         ierr = NF_GET_VAR_DOUBLE(nid, nvarid, agesno)
-#else
-         ierr = NF_GET_VAR_REAL(nid, nvarid, agesno)
+         PRINT*, '          Mais je vais essayer de lire AGESNO**'
+         DO nsrf = 1, nbsrf
+           IF (nsrf.GT.99) THEN
+             PRINT*, "Trop de sous-mailles"
+             CALL abort
+           ENDIF
+           WRITE(str2,'(i2.2)') nsrf
+           ierr = NF_INQ_VARID (nid, "AGESNO"//str2, nvarid)
+           IF (ierr.NE.NF_NOERR) THEN
+              PRINT*, "phyetat0: Le champ <AGESNO"//str2//"> est absent"
+              agesno = 50.0
+           ENDIF
+#ifdef NC_DOUBLE
+           ierr = NF_GET_VAR_DOUBLE(nid, nvarid, agesno(1,nsrf))
+#else
+           ierr = NF_GET_VAR_REAL(nid, nvarid, agesno(1,nsrf))
+#endif
+           IF (ierr.NE.NF_NOERR) THEN
+             PRINT*, "phyetat0: Lecture echouee pour <AGESNO"//str2//">"
+             CALL abort
+           ENDIF
+           xmin = 1.0E+20
+           xmax = -1.0E+20
+           DO i = 1, klon
+              xmin = MIN(agesno(i,nsrf),xmin)
+              xmax = MAX(agesno(i,nsrf),xmax)
+           ENDDO
+           PRINT*,'Age de la neige AGESNO**:', nsrf, xmin, xmax
+         ENDDO
+      ELSE
+         PRINT*, 'phyetat0: Le champ <AGESNO> est present'
+         PRINT*, '          J ignore donc les autres AGESNO**'
+#ifdef NC_DOUBLE
+         ierr = NF_GET_VAR_DOUBLE(nid, nvarid, agesno(1,1))
+#else
+         ierr = NF_GET_VAR_REAL(nid, nvarid, agesno(1,1))
 #endif
          IF (ierr.NE.NF_NOERR) THEN
-            PRINT*, 'phyetat0: Lecture echouee pour <AGESNO>'
+            PRINT*, "phyetat0: Lecture echouee pour <AGESNO>"
             CALL abort
          ENDIF
@@ -905 +939 @@
          xmax = -1.0E+20
          DO i = 1, klon
-            xmin = MIN(agesno(i),xmin)
-            xmax = MAX(agesno(i),xmax)
-         ENDDO
-         PRINT*,'Age de la neige agesno:', xmin, xmax
-      ENDIF
-c
+            xmin = MIN(agesno(i,1),xmin)
+            xmax = MAX(agesno(i,1),xmax)
+         ENDDO
+         PRINT*,'Age de la neige <AGESNO>', xmin, xmax
+         DO nsrf = 2, nbsrf
+         DO i = 1, klon
+            agesno(i,nsrf) = agesno(i,1)
+         ENDDO
+         ENDDO
+      ENDIF
+
 c
       ierr = NF_INQ_VARID (nid, "ZMEA", nvarid)