Changeset 152

Timestamp:

Jun 8, 2011, 4:49:35 PM (14 years ago)

Author:

slebonnois

Message:

SL: correction bug lie a itau_phy pour phytitan et phyvenus

Location:

trunk

Files:

• LMDZ.TITAN/libf/phytitan/physiq.F (modified) (2 diffs)
• LMDZ.VENUS/libf/phyvenus/physiq.F (modified) (3 diffs)

trunk/LMDZ.TITAN/libf/phytitan/physiq.F

@@ -474 +474 @@
      .       t_ancien, ancien_ok)
 
-c dtime est lu dans startphy, alors que 
+c dtime est lu dans startphy 
 c pdtphys est calcule a partir des nouvelles conditions:
 c Reinitialisation du pas de temps physique quand changement
@@ -482 +482 @@
 c           abort_message='Pas physique n est pas correct '
 c           call abort_gcm(modname,abort_message,1)
+c----------------
+c pour initialiser convenablement le time_counter, il faut tenir compte
+c du changement de dtime en changeant itau_phy (point de depart)
+            itau_phy = NINT(itau_phy*dtime/pdtphys)
+c----------------
             dtime=pdtphys
          ENDIF

trunk/LMDZ.VENUS/libf/phyvenus/physiq.F

@@ -340 +340 @@
 c
       INTEGER nhori, nvert, idayref
+      INTEGER itau_phy_modif
       REAL zsto, zout, zsto1, zsto2, zero
       parameter (zero=0.0e0)
@@ -446 +447 @@
      .       t_ancien, ancien_ok)
 
-c dtime est lu dans startphy, alors que 
+c dtime est lu dans startphy 
 c pdtphys est calcule a partir des nouvelles conditions:
 c Reinitialisation du pas de temps physique quand changement
@@ -454 +455 @@
 c           abort_message='Pas physique n est pas correct '
 c           call abort_gcm(modname,abort_message,1)
+c----------------
+c pour initialiser convenablement le time_counter, il faut tenir compte
+c du changement de dtime en changeant itau_phy (point de depart)
+            itau_phy = NINT(itau_phy*dtime/pdtphys)
+c----------------
             dtime=pdtphys
          ENDIF