Changeset 795 for trunk

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

Location:

trunk/LMDZ.MARS

Files:

• README (modified) (1 diff)
• libf/dyn3d/gcm.F (modified) (1 diff)

trunk/LMDZ.MARS/README

@@ -1806 +1806 @@
   dynamics/physics cycles. The fraction of the sol that a (re)start.nc
   file corresponds to is (read)written as 'Time' in the file.
+
+== 25/09/2012 == EM
+>> 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.

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