Changeset 348 for LMDZ.3.3/trunk


Ignore:
Timestamp:
Mar 19, 2002, 5:12:44 PM (22 years ago)
Author:
lmdz
Message:

Regle le probleme de decalage de 1 jour au debut de chaque simulation
LF

Location:
LMDZ.3.3/trunk/libf
Files:
6 edited

Legend:

Unmodified
Added
Removed

  • LMDZ.3.3/trunk/libf/bibio/initdynav.F

    r2 r348  
    6262      real rlong(iip1,jjp1), rlat(iip1,jjp1)
    6363      integer ii,jj
    64       integer zan
     64      integer zan, idayref
    6565C
    6666C  Initialisations
    … …  
    7474
    7575      zan = anne0
    76       CALL ymds2ju(zan, 1, 1, 0.0, zjulian)
    77       zjulian = zjulian + day0
     76      idayref = day0
     77      CALL ymds2ju(zan, 1, idayref, 0.0, zjulian)
    7878      tau0 = 0
    7979     

  • LMDZ.3.3/trunk/libf/bibio/initfluxsto.F

    r54 r348  
    6565      integer uhoriid, vhoriid, thoriid, zvertiid,dhoriid,dvertiid
    6666      integer ii,jj
    67       integer zan
     67      integer zan, idayref
    6868      logical ok_sync
    6969C
    … …  
    7979
    8080      zan = anne_ini
    81       CALL ymds2ju(zan, 1, 1, 0.0, zjulian)
    82       zjulian = zjulian + day_ini
     81      idayref = day_ini
     82      CALL ymds2ju(zan, 1, idayref, 0.0, zjulian)
    8383      tau0 = 0
    8484       

  • LMDZ.3.3/trunk/libf/bibio/inithist.F

    r2 r348  
    6363      integer uhoriid, vhoriid, thoriid, zvertiid
    6464      integer ii,jj
    65       integer zan
     65      integer zan, idayref
    6666C
    6767C  Initialisations
    … …  
    7575
    7676      zan = anne0
    77       CALL ymds2ju(zan, 1, 1, 0.0, zjulian)
    78       zjulian = zjulian + day0
     77      idayref = day0
     78      CALL ymds2ju(zan, 1, idayref, 0.0, zjulian)
    7979      tau0 = 0
    8080     

  • LMDZ.3.3/trunk/libf/phylmd/initphysto.F

    r214 r348  
    6868      integer uhoriid, vhoriid, thoriid, zvertiid
    6969      integer ii,jj
    70       integer zan
     70      integer zan, idayref
    7171      logical ok_sync
    7272      REAL zx_lon(iim,jjm+1), zx_lat(iim,jjm+1)
    … …  
    8585
    8686      zan = anneeref
    87       CALL ymds2ju(zan, 1, 1, 0.0, zjulian)
    88       zjulian = zjulian + dayref
     87      idayref = dayref
     88      CALL ymds2ju(zan, 1, idayref, 0.0, zjulian)
    8989      tau0 = 0
    9090       

  • LMDZ.3.3/trunk/libf/phylmd/physiq.F

    r344 r348  
    632632      INTEGER nhori, nvert
    633633      REAL zsto, zout, zjulian
     634      integer idayref
    634635
    635636      character*20 modname
    … …  
    794795      IF (ok_journe) THEN
    795796c
    796          CALL ymds2ju(anne_ini, 1, 1, 0.0, zjulian)
    797          zjulian = zjulian + day_ini
     797         idayref = day_ini
     798         CALL ymds2ju(anne_ini, 1, idayref, 0.0, zjulian)
    798799c
    799800         CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlon,zx_lon)
    … …  
    950951      IF (ok_mensuel) THEN
    951952c
    952          CALL ymds2ju(anne_ini, 1, 1, 0.0, zjulian)
    953          zjulian = zjulian + day_ini
     953         idayref = day_ini
     954         CALL ymds2ju(anne_ini, 1, idayref, 0.0, zjulian)
    954955c
    955956         CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlon,zx_lon)
    … …  
    13431344      IF (ok_instan) THEN
    13441345c
    1345          CALL ymds2ju(anne_ini, 1, 1, 0.0, zjulian)
    1346          zjulian = zjulian + day_ini
     1346         idayref = day_ini
     1347         CALL ymds2ju(anne_ini, 1, idayref, 0.0, zjulian)
    13471348c
    13481349         CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlon,zx_lon)

  • LMDZ.3.3/trunk/libf/phylmd/phytrac.F

    r344 r348  
    210210      INTEGER nhori, nvert, nverta, nvertb, nverts
    211211      REAL zsto, zout, zjulian
     212      integer idayref
    212213      INTEGER nid_tra
    213214      SAVE nid_tra
    … …  
    331332
    332333C         
    333          CALL ymds2ju(anne_ini, 1, 1, 0.0, zjulian)
     334         idayref = 1
     335         CALL ymds2ju(anne_ini, 1, idayref, 0.0, zjulian)
    334336         itra = NINT(FLOAT(day_ini)*86400./pdtphys)
    335337         itap = NINT(FLOAT(day_ini)*86400./pdtphys)
Note: See TracChangeset for help on using the changeset viewer.