Changeset 775

Timestamp:

Sep 6, 2012, 4:14:51 PM (12 years ago)

Author:

jleconte

Message:

05/09/2012 == JL

• Corrects typo in previous commit to account for possible eccentricity in

tlocked case

File:

• trunk/LMDZ.GENERIC/libf/phystd/physiq.F90 (modified) (1 diff)

trunk/LMDZ.GENERIC/libf/phystd/physiq.F90

@@ -645 +645 @@
 !              ztim3=-COS(declin)*SIN(2.*pi*(zday/year_day) - zls*nres)
 ! JL12 corrects tidally resonant cases. nres=omega_rot/omega_orb
-              ztim2=COS(declin)*COS(2.*pi*(zday/year_day)*(nres-1.))
-              ztim3=-COS(declin)*SIN(2.*pi*(zday/year_day)*(nres-1.))
+              ztim2=COS(declin)*COS(2.*pi*(zday/year_day)*nres - zls)
+              ztim3=-COS(declin)*SIN(2.*pi*(zday/year_day)*nres - zls)
 
               call stelang(ngrid,sinlon,coslon,sinlat,coslat,    &