Changeset 795 for trunk/LMDZ.MARS/libf

Timestamp:

Sep 25, 2012, 10:33:38 AM (12 years ago)

Author:

emillour

Message:

Mars GCM:

Minor improvement for running fractions of sols; enforce some rounding of

'Time' (ie: fraction of day) read from the start.nc file in order to limit
the accumultation of roundoffs from multiple runs.

EM

File:

• trunk/LMDZ.MARS/libf/dyn3d/gcm.F (modified) (1 diff)

trunk/LMDZ.MARS/libf/dyn3d/gcm.F

@@ -149 +149 @@
 
       CALL defrun_new( 99, .TRUE. )
+      ! in case time_0 (because of roundoffs) is close to zero,
+      ! set it to zero to avoid roundoff propagation issues
+      if ((time_0.gt.0.).and.(time_0.lt.(1./day_step))) then
+        write(*,*)"GCM: In start.nc, time=",time_0
+        write(*,*)"     but day_step=",day_step
+        write(*,*)"     and 1./day_step=",1./day_step
+        write(*,*)"     fix this drift by setting time=0"
+        time_0=0.
+      endif