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

Timestamp:

May 25, 2005, 3:10:09 PM (20 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/phyredem.F (modified) (4 diffs)

LMDZ4/trunk/libf/phylmd/phyredem.F

@@ -4 +4 @@
 c
       SUBROUTINE phyredem (fichnom,dtime,radpas,
-     .           rlat,rlon, pctsrf,tsol,tsoil,deltat,qsurf,qsol,snow,
+     .           rlat,rlon, pctsrf,tsol,tsoil,
+cIM "slab" ocean
+     .           tslab,seaice,
+     .           qsurf,qsol,snow,
      .           albedo, alblw, evap, rain_fall, snow_fall,
      .           solsw, sollw,fder,
@@ -31 +34 @@
       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)
@@ -288 +292 @@
       ENDDO
 c
-c
-      ierr = NF_REDEF (nid)
-#ifdef NC_DOUBLE
-      ierr = NF_DEF_VAR (nid, "DELTAT", NF_DOUBLE, 1, idim2,nvarid)
-#else
-      ierr = NF_DEF_VAR (nid, "DELTAT", NF_FLOAT, 1, idim2,nvarid)
+cIM "slab" ocean
+      ierr = NF_REDEF (nid)
+#ifdef NC_DOUBLE
+      ierr = NF_DEF_VAR (nid, "TSLAB", NF_DOUBLE, 1, idim2,nvarid)
+#else
+      ierr = NF_DEF_VAR (nid, "TSLAB", NF_FLOAT, 1, idim2,nvarid)
 #endif
       ierr = NF_PUT_ATT_TEXT (nid,nvarid,"title", 33,
@@ -299 +303 @@
       ierr = NF_ENDDEF(nid)
 #ifdef NC_DOUBLE
-      ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,deltat)
-#else
-      ierr = NF_PUT_VAR_REAL (nid,nvarid,deltat)
+      ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,tslab)
+#else
+      ierr = NF_PUT_VAR_REAL (nid,nvarid,tslab)
+#endif
+c
+      ierr = NF_REDEF (nid)
+#ifdef NC_DOUBLE
+      ierr = NF_DEF_VAR (nid, "SEAICE", NF_DOUBLE, 1, idim2,nvarid)
+#else
+      ierr = NF_DEF_VAR (nid, "SEAICE", NF_FLOAT, 1, idim2,nvarid)
+#endif
+      ierr = NF_PUT_ATT_TEXT (nid,nvarid,"title", 33,
+     .                        "Glace de mer kg/m2 (pour slab-ocean)")
+      ierr = NF_ENDDEF(nid)
+#ifdef NC_DOUBLE
+      ierr = NF_PUT_VAR_DOUBLE (nid,nvarid,seaice)
+#else
+      ierr = NF_PUT_VAR_REAL (nid,nvarid,seaice)
 #endif
 c