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

Timestamp:

Mar 31, 2008, 5:59:23 PM (17 years ago)

Author:

lmdzadmin

Message:

Ajout variables convection (ema_work1, ema_work2) dans startphy.nc
IM

File:

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

LMDZ4/trunk/libf/phylmd/phyetat0.F

@@ -13 +13 @@
      .           rugsrel_p,tabcntr0,
      .           t_ancien_p,q_ancien_p,ancien_ok_p, rnebcon_p, ratqs_p,
-     .           clwcon_p,pbl_tke_p, zmax0_p, f0_p)
+     .           clwcon_p,pbl_tke_p, zmax0_p, f0_p, 
+     .           ema_work1_p, ema_work2_p)
 
       USE dimphy
@@ -49 +50 @@
       REAL pbl_tke_p(klon,klev,nbsrf)
       REAL zmax0_p(klon), f0_p(klon)
+      REAL ema_work1_p(klon,klev), ema_work2_p(klon,klev)
       REAL tsoil_p(klon,nsoilmx,nbsrf)
       REAL tslab_p(klon), seaice_p(klon)
@@ -87 +89 @@
       REAL pbl_tke(klon_glo,klev,nbsrf)
       REAL zmax0(klon), f0(klon)
+      REAL ema_work1(klon,klev), ema_work2(klon,klev)
       REAL tsoil(klon_glo,nsoilmx,nbsrf)
 cIM "slab" ocean
@@ -1434 +1437 @@
 
 c
-      ierr = NF_INQ_VARID (nid, "QANCIEN", nvarid)
-      IF (ierr.NE.NF_NOERR) THEN
-         PRINT*, "phyetat0: Le champ <QANCIEN> est absent"
-         PRINT*, "Depart legerement fausse. Mais je continue"
-         ancien_ok = .FALSE.
-      ELSE
-#ifdef NC_DOUBLE
-         ierr = NF_GET_VAR_DOUBLE(nid, nvarid, q_ancien)
-#else
-         ierr = NF_GET_VAR_REAL(nid, nvarid, q_ancien)
-#endif
-         IF (ierr.NE.NF_NOERR) THEN
-            PRINT*, "phyetat0: Lecture echouee pour <QANCIEN>"
-            CALL abort
-         ENDIF
-      ENDIF
-c
 c Lecture ratqs
 c
@@ -1587 +1573 @@
       xmax = MAXval(f0)
       PRINT*,'(ecart-type) f0:', xmin, xmax
+c
+c ema_work1
+c
+      ierr = NF_INQ_VARID (nid, "EMA_WORK1", nvarid)
+      IF (ierr.NE.NF_NOERR) THEN
+         PRINT*, "phyetat0: Le champ <EMA_WORK1> est absent"
+         PRINT*, "Depart legerement fausse. Mais je continue"
+         ema_work1=0.
+      ELSE
+#ifdef NC_DOUBLE
+         ierr = NF_GET_VAR_DOUBLE(nid, nvarid, ema_work1)
+#else
+         ierr = NF_GET_VAR_REAL(nid, nvarid, ema_work1)
+#endif
+         IF (ierr.NE.NF_NOERR) THEN
+            PRINT*, "phyetat0: Lecture echouee pour <EMA_WORK1>"
+            CALL abort
+         ENDIF
+           xmin = 1.0E+20
+           xmax = -1.0E+20
+           DO k = 1, klev
+           DO i = 1, klon
+              xmin = MIN(ema_work1(i,k),xmin)
+              xmax = MAX(ema_work1(i,k),xmax)
+           ENDDO
+           ENDDO
+           PRINT*,'ema_work1:', xmin, xmax
+      ENDIF
+c
+c ema_work2
+c
+      ierr = NF_INQ_VARID (nid, "EMA_WORK2", nvarid)
+      IF (ierr.NE.NF_NOERR) THEN
+         PRINT*, "phyetat0: Le champ <EMA_WORK2> est absent"
+         PRINT*, "Depart legerement fausse. Mais je continue"
+         ema_work2=0.
+      ELSE
+#ifdef NC_DOUBLE
+         ierr = NF_GET_VAR_DOUBLE(nid, nvarid, ema_work2)
+#else
+         ierr = NF_GET_VAR_REAL(nid, nvarid, ema_work2)
+#endif
+         IF (ierr.NE.NF_NOERR) THEN
+            PRINT*, "phyetat0: Lecture echouee pour <EMA_WORK2>"
+            CALL abort
+         ENDIF
+           xmin = 1.0E+20
+           xmax = -1.0E+20
+           DO k = 1, klev
+           DO i = 1, klon
+              xmin = MIN(ema_work2(i,k),xmin)
+              xmax = MAX(ema_work2(i,k),xmax)
+           ENDDO
+           ENDDO
+           PRINT*,'ema_work2:', xmin, xmax
+      ENDIF
 c
 c Fermer le fichier:
@@ -1625 +1667 @@
       call Scatter( zmax0,zmax0_p)
       call Scatter( f0,f0_p)
+      call Scatter( ema_work1, ema_work1_p)
+      call Scatter( ema_work2, ema_work2_p)
       call Scatter( tsoil,tsoil_p)
       call Scatter( tslab,tslab_p)