Changeset 644 for LMDZ4/trunk/libf/phylmd/phyetat0.F

Timestamp:

May 25, 2005, 3:10:09 PM (19 years ago)

Author:

Laurent Fairhead

Message:

Synchronisation avec tous les diagnostiques de Ionela IM
Inclusion du slab ocean IM
LF

File:

• LMDZ4/trunk/libf/phylmd/phyetat0.F (modified) (5 diffs)

LMDZ4/trunk/libf/phylmd/phyetat0.F

@@ -5 +5 @@
 c
       SUBROUTINE phyetat0 (fichnom,dtime,co2_ppm_etat0,solaire_etat0,
-     .            rlat,rlon, pctsrf, tsol,tsoil,deltat,qsurf,qsol,snow,
+     .            rlat,rlon, pctsrf, tsol,tsoil,
+cIM "slab" ocean
+     .           tslab,seaice,
+     .           qsurf,qsol,snow,
      .           albe, alblw, evap, rain_fall, snow_fall, solsw, sollw,
      .           fder,radsol,frugs,agesno,clesphy0,
@@ -32 +35 @@
       REAL tsol(klon,nbsrf)
       REAL tsoil(klon,nsoilmx,nbsrf)
-      REAL deltat(klon)
+cIM "slab" ocean
+      REAL tslab(klon), seaice(klon)
       REAL qsurf(klon,nbsrf)
       REAL qsol(klon)
       REAL snow(klon,nbsrf)
       REAL albe(klon,nbsrf)
-cIM BEG alblw
       REAL alblw(klon,nbsrf)
-cIM END alblw
       REAL evap(klon,nbsrf)
       REAL radsol(klon)
@@ -420 +422 @@
       ENDDO
 c
-c Lecture de deltat (pour slab ocean seulement):
-c
-      ierr = NF_INQ_VARID (nid, "DELTAT", nvarid)
-      IF (ierr.NE.NF_NOERR) THEN
-         PRINT*, "phyetat0: Le champ <DELTAT> est absent"
-         CALL abort
-      ENDIF
-#ifdef NC_DOUBLE
-      ierr = NF_GET_VAR_DOUBLE(nid, nvarid, deltat)
-#else
-      ierr = NF_GET_VAR_REAL(nid, nvarid, deltat)
-#endif
-      IF (ierr.NE.NF_NOERR) THEN
-         PRINT*, "phyetat0: Lecture echouee pour <DELTAT>"
-         CALL abort
-      ENDIF
-      xmin = 1.0E+20
-      xmax = -1.0E+20
-      DO i = 1, klon
-         xmin = MIN(deltat(i),xmin)
-         xmax = MAX(deltat(i),xmax)
-      ENDDO
-      PRINT*,'Ecart de la SST deltat:', xmin, xmax
+cIM "slab" ocean 
+c
+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
+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
 c
 c Lecture de l'humidite de l'air juste au dessus du sol:
@@ -652 +680 @@
 
 c
-cIM BEG alblw
 c Lecture de albedo au sol LW:
 c
@@ -690 +717 @@
          ENDDO
       ENDIF
-
-cIM END alblw
-
 c
 c Lecture de evaporation: