Changeset 348 for LMDZ.3.3/trunk/libf/phylmd

Timestamp:

Mar 19, 2002, 5:12:44 PM (23 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/phylmd

Files:

• initphysto.F (modified) (2 diffs)
• physiq.F (modified) (4 diffs)
• phytrac.F (modified) (2 diffs)

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

@@ -68 +68 @@
       integer uhoriid, vhoriid, thoriid, zvertiid
       integer ii,jj
-      integer zan
+      integer zan, idayref
       logical ok_sync
       REAL zx_lon(iim,jjm+1), zx_lat(iim,jjm+1)
@@ -85 +85 @@
 
       zan = anneeref
-      CALL ymds2ju(zan, 1, 1, 0.0, zjulian)
-      zjulian = zjulian + dayref
+      idayref = dayref
+      CALL ymds2ju(zan, 1, idayref, 0.0, zjulian)
       tau0 = 0
         

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

@@ -632 +632 @@
       INTEGER nhori, nvert
       REAL zsto, zout, zjulian
+      integer idayref
 
       character*20 modname
@@ -794 +795 @@
       IF (ok_journe) THEN
 c
-         CALL ymds2ju(anne_ini, 1, 1, 0.0, zjulian)
-         zjulian = zjulian + day_ini
+         idayref = day_ini
+         CALL ymds2ju(anne_ini, 1, idayref, 0.0, zjulian)
 c
          CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlon,zx_lon)
@@ -950 +951 @@
       IF (ok_mensuel) THEN
 c
-         CALL ymds2ju(anne_ini, 1, 1, 0.0, zjulian)
-         zjulian = zjulian + day_ini
+         idayref = day_ini
+         CALL ymds2ju(anne_ini, 1, idayref, 0.0, zjulian)
 c
          CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlon,zx_lon)
@@ -1343 +1344 @@
       IF (ok_instan) THEN
 c
-         CALL ymds2ju(anne_ini, 1, 1, 0.0, zjulian)
-         zjulian = zjulian + day_ini
+         idayref = day_ini
+         CALL ymds2ju(anne_ini, 1, idayref, 0.0, zjulian)
 c
          CALL gr_fi_ecrit(1,klon,iim,jjmp1,rlon,zx_lon)

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

@@ -210 +210 @@
       INTEGER nhori, nvert, nverta, nvertb, nverts
       REAL zsto, zout, zjulian
+      integer idayref
       INTEGER nid_tra
       SAVE nid_tra
@@ -331 +332 @@
 
 C         
-         CALL ymds2ju(anne_ini, 1, 1, 0.0, zjulian)
+         idayref = 1
+         CALL ymds2ju(anne_ini, 1, idayref, 0.0, zjulian)
          itra = NINT(FLOAT(day_ini)*86400./pdtphys)
          itap = NINT(FLOAT(day_ini)*86400./pdtphys)