Changeset 3513 for LMDZ6/trunk/libf

Timestamp:

May 20, 2019, 12:37:25 PM (6 years ago)

Author:

fhourdin

Message:

Correction pour le 1D

File:

• LMDZ6/trunk/libf/phylmd/dyn1d/1DUTILS.h (modified) (2 diffs)

LMDZ6/trunk/libf/phylmd/dyn1d/1DUTILS.h

@@ -814 +814 @@
       character*80 abort_message
 !
-      INTEGER nb
-      SAVE nb
-      DATA nb / 0 /
+      INTEGER pass
 
       CALL open_restartphy(fichnom)
@@ -885 +883 @@
        tab_cntrl(31) = FLOAT(itau_dyn + itaufin)
 !
-      CALL put_var("controle","Param. de controle Dyn1D",tab_cntrl)
+      DO pass=1,2
+      CALL put_var(pass,"controle","Param. de controle Dyn1D",tab_cntrl)
 !
 
 !  Ecriture/extension de la coordonnee temps
 
-      nb = nb + 1
 
 !  Ecriture des champs
 !
-      CALL put_field("plev","p interfaces sauf la nulle",plev)
-      CALL put_field("play","",play)
-      CALL put_field("phi","geopotentielle",phi)
-      CALL put_field("phis","geopotentiell de surface",phis)
-      CALL put_field("presnivs","",presnivs)
-      CALL put_field("ucov","",ucov)
-      CALL put_field("vcov","",vcov)
-      CALL put_field("temp","",temp)
-      CALL put_field("omega2","",omega2)
+      CALL put_field(pass,"plev","p interfaces sauf la nulle",plev)
+      CALL put_field(pass,"play","",play)
+      CALL put_field(pass,"phi","geopotentielle",phi)
+      CALL put_field(pass,"phis","geopotentiell de surface",phis)
+      CALL put_field(pass,"presnivs","",presnivs)
+      CALL put_field(pass,"ucov","",ucov)
+      CALL put_field(pass,"vcov","",vcov)
+      CALL put_field(pass,"temp","",temp)
+      CALL put_field(pass,"omega2","",omega2)
 
       Do iq=1,nqtot
-        CALL put_field("q"//nmq(iq),"eau vap ou condens et traceurs",           &
+        CALL put_field(pass,"q"//nmq(iq),"eau vap ou condens et traceurs",           &
      &                                                      q(:,:,iq))
       EndDo
-      CALL close_restartphy
+    IF (pass==1) CALL enddef_restartphy
+    IF (pass==2) CALL close_restartphy
+
+
+      ENDDO
 
 !