Changeset 1201 for LMDZ4/branches/LMDZ4-dev/libf/dyn3dpar/sortvarc.F

Timestamp:

Jul 7, 2009, 4:01:00 PM (15 years ago)

Author:

Laurent Fairhead

Message:

Modifications nécessaires a l'inclusion d'un calendrier réaliste.
La date courante est calculée dans leapfrog.F et exprimée en Jour Julien
(modifié). On en a profité pour faire un peu de ménage dans la gestion des dates
du modèle.
Dans la physique, on utilise les routines de passages entre calendrier Julien et
Gregorien incluses dans IOIPSL pour calculer le nombre de jours écoulés depuis le
1er janvier (pour les conditions aux limites) ou l'equinoxe (pour le calcul de
la longitude solaire). Le calcul de l'orbite reprend celui du gcm planétaire
(codé par FH)
On décide du calendrier à utiliser à l'aide du paramètre calend du run.def. Par
défaut celui-ci est à earth_360d
LF

File:

• LMDZ4/branches/LMDZ4-dev/libf/dyn3dpar/sortvarc.F (modified) (2 diffs)

LMDZ4/branches/LMDZ4-dev/libf/dyn3dpar/sortvarc.F

@@ -129 +129 @@
       ang   = SSUM(     llm,  angl, 1 )
 
-      rday = FLOAT(INT ( day_ini + time ))
+c      rday = FLOAT(INT ( day_ini + time ))
 c
+       rday = FLOAT(INT(time-jD_ref-jH_ref))
       IF(ptot0.eq.0.)  THEN
          PRINT 3500, itau, rday, heure,time
@@ -156 +157 @@
       RETURN
 
-3500   FORMAT('0'10(1h*),4x,'pas'i7,5x,'jour'f5.0,'heure'f5.1,4x 
-     *   ,'date',f10.5,4x,10(1h*))
+3500   FORMAT('0'10(1h*),4x,'pas'i7,5x,'jour'f9.0,'heure'f5.1,4x 
+     *   ,'date',f14.4,4x,10(1h*))
 4000   FORMAT(10x,'masse',4x,'rmsdpdt',7x,'energie',2x,'enstrophie'
      * ,2x,'entropie',3x,'rmsv',4x,'mt.ang',/,'GLOB  '