Ignore:
Timestamp:
May 20, 2019, 12:37:25 PM (6 years ago)
Author:
fhourdin
Message:

Correction pour le 1D

File:
1 edited

Legend:

Unmodified
Added
Removed
  • LMDZ6/trunk/libf/phylmd/dyn1d/1DUTILS.h

    r3316 r3513  
    814814      character*80 abort_message
    815815!
    816       INTEGER nb
    817       SAVE nb
    818       DATA nb / 0 /
     816      INTEGER pass
    819817
    820818      CALL open_restartphy(fichnom)
     
    885883       tab_cntrl(31) = FLOAT(itau_dyn + itaufin)
    886884!
    887       CALL put_var("controle","Param. de controle Dyn1D",tab_cntrl)
     885      DO pass=1,2
     886      CALL put_var(pass,"controle","Param. de controle Dyn1D",tab_cntrl)
    888887!
    889888
    890889!  Ecriture/extension de la coordonnee temps
    891890
    892       nb = nb + 1
    893891
    894892!  Ecriture des champs
    895893!
    896       CALL put_field("plev","p interfaces sauf la nulle",plev)
    897       CALL put_field("play","",play)
    898       CALL put_field("phi","geopotentielle",phi)
    899       CALL put_field("phis","geopotentiell de surface",phis)
    900       CALL put_field("presnivs","",presnivs)
    901       CALL put_field("ucov","",ucov)
    902       CALL put_field("vcov","",vcov)
    903       CALL put_field("temp","",temp)
    904       CALL put_field("omega2","",omega2)
     894      CALL put_field(pass,"plev","p interfaces sauf la nulle",plev)
     895      CALL put_field(pass,"play","",play)
     896      CALL put_field(pass,"phi","geopotentielle",phi)
     897      CALL put_field(pass,"phis","geopotentiell de surface",phis)
     898      CALL put_field(pass,"presnivs","",presnivs)
     899      CALL put_field(pass,"ucov","",ucov)
     900      CALL put_field(pass,"vcov","",vcov)
     901      CALL put_field(pass,"temp","",temp)
     902      CALL put_field(pass,"omega2","",omega2)
    905903
    906904      Do iq=1,nqtot
    907         CALL put_field("q"//nmq(iq),"eau vap ou condens et traceurs",           &
     905        CALL put_field(pass,"q"//nmq(iq),"eau vap ou condens et traceurs",           &
    908906     &                                                      q(:,:,iq))
    909907      EndDo
    910       CALL close_restartphy
     908    IF (pass==1) CALL enddef_restartphy
     909    IF (pass==2) CALL close_restartphy
     910
     911
     912      ENDDO
    911913
    912914!
Note: See TracChangeset for help on using the changeset viewer.