Changeset 651 for LMDZ4/trunk/libf

Timestamp:

Jun 3, 2005, 12:03:07 PM (20 years ago)

Author:

Laurent Fairhead

Message:

On fait dépendre la lecture de tslab et seaice dans phyetat0 du flag ocean
(on n'a pas besoin de les lire si on n'est pas en ocean slab)
LF

Location:

LMDZ4/trunk/libf/phylmd

Files:

• phyetat0.F (modified) (3 diffs)
• physiq.F (modified) (1 diff)

LMDZ4/trunk/libf/phylmd/phyetat0.F

@@ -7 +7 @@
      .            rlat,rlon, pctsrf, tsol,tsoil,
 cIM "slab" ocean
-     .           tslab,seaice,
+     .           ocean, tslab,seaice,
      .           qsurf,qsol,snow,
      .           albe, alblw, evap, rain_fall, snow_fall, solsw, sollw,
@@ -67 +67 @@
       LOGICAL ancien_ok
 
+      CHARACTER*6 ocean
+
       INTEGER        longcles
       PARAMETER    ( longcles = 20 )
@@ -426 +428 @@
 c Lecture de tslab (pour slab ocean seulement):      
 c
-      ierr = NF_INQ_VARID (nid, "TSLAB", nvarid)
-      IF (ierr.NE.NF_NOERR) THEN
-         PRINT*, "phyetat0: Le champ <TSLAB> est absent"
-         CALL abort
-      ENDIF
-#ifdef NC_DOUBLE
-      ierr = NF_GET_VAR_DOUBLE(nid, nvarid, tslab)
-#else
-      ierr = NF_GET_VAR_REAL(nid, nvarid, tslab)
-#endif
-      IF (ierr.NE.NF_NOERR) THEN
-         PRINT*, "phyetat0: Lecture echouee pour <TSLAB>"
-         CALL abort
-      ENDIF
-      xmin = 1.0E+20
-      xmax = -1.0E+20
-      DO i = 1, klon
-         xmin = MIN(tslab(i),xmin)
-         xmax = MAX(tslab(i),xmax)
-      ENDDO
-      PRINT*,'Ecart de la SST tslab:', xmin, xmax
+      IF (ocean .eq. 'slab  ') then
+        ierr = NF_INQ_VARID (nid, "TSLAB", nvarid)
+        IF (ierr.NE.NF_NOERR) THEN
+          PRINT*, "phyetat0: Le champ <TSLAB> est absent"
+          CALL abort
+        ENDIF
+#ifdef NC_DOUBLE
+        ierr = NF_GET_VAR_DOUBLE(nid, nvarid, tslab)
+#else
+        ierr = NF_GET_VAR_REAL(nid, nvarid, tslab)
+#endif
+        IF (ierr.NE.NF_NOERR) THEN
+          PRINT*, "phyetat0: Lecture echouee pour <TSLAB>"
+          CALL abort
+        ENDIF
+        xmin = 1.0E+20
+        xmax = -1.0E+20
+        DO i = 1, klon
+          xmin = MIN(tslab(i),xmin)
+          xmax = MAX(tslab(i),xmax)
+        ENDDO
+        PRINT*,'Ecart de la SST tslab:', xmin, xmax
 c
 c Lecture de seaice (pour slab ocean seulement):
 c
-      ierr = NF_INQ_VARID (nid, "SEAICE", nvarid)
-      IF (ierr.NE.NF_NOERR) THEN
-         PRINT*, "phyetat0: Le champ <SEAICE> est absent"
-         CALL abort
-      ENDIF
-#ifdef NC_DOUBLE
-      ierr = NF_GET_VAR_DOUBLE(nid, nvarid, seaice)
-#else
-      ierr = NF_GET_VAR_REAL(nid, nvarid, seaice)
-#endif
-      IF (ierr.NE.NF_NOERR) THEN
-         PRINT*, "phyetat0: Lecture echouee pour <SEAICE>"
-         CALL abort
-      ENDIF
-      xmin = 1.0E+20
-      xmax = -1.0E+20
-      DO i = 1, klon
-         xmin = MIN(seaice(i),xmin)
-         xmax = MAX(seaice(i),xmax)
-      ENDDO
-      PRINT*,'Masse de la glace de mer seaice:', xmin, xmax
+        ierr = NF_INQ_VARID (nid, "SEAICE", nvarid)
+        IF (ierr.NE.NF_NOERR) THEN
+          PRINT*, "phyetat0: Le champ <SEAICE> est absent"
+          CALL abort
+        ENDIF
+#ifdef NC_DOUBLE
+        ierr = NF_GET_VAR_DOUBLE(nid, nvarid, seaice)
+#else
+        ierr = NF_GET_VAR_REAL(nid, nvarid, seaice)
+#endif
+        IF (ierr.NE.NF_NOERR) THEN
+          PRINT*, "phyetat0: Lecture echouee pour <SEAICE>"
+          CALL abort
+        ENDIF
+        xmin = 1.0E+20
+        xmax = -1.0E+20
+        DO i = 1, klon
+          xmin = MIN(seaice(i),xmin)
+          xmax = MAX(seaice(i),xmax)
+        ENDDO
+        PRINT*,'Masse de la glace de mer seaice:', xmin, xmax
+      ELSE
+        tslab = 0.
+        seaice = 0.
+      ENDIF
 c
 c Lecture de l'humidite de l'air juste au dessus du sol:

LMDZ4/trunk/libf/phylmd/physiq.F

@@ -1310 +1310 @@
      .       rlat,rlon,pctsrf, ftsol,ftsoil,
 cIM "slab" ocean
-     .       tslab,seaice,
+     .       ocean, tslab,seaice,
      .       fqsurf,qsol,fsnow,
      .       falbe, falblw, fevap, rain_fall,snow_fall,solsw, sollwdown,