Changeset 1201 for LMDZ4/branches/LMDZ4-dev/libf/dyn3d/leapfrog.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/dyn3d/leapfrog.F (modified) (11 diffs)

LMDZ4/branches/LMDZ4-dev/libf/dyn3d/leapfrog.F

@@ +1 @@
+!
+! $Id$
 !
 c
@@ -112 +114 @@
 c
       INTEGER itau,itaufinp1,iav
-      INTEGER*4  iday ! jour julien
-      REAL       time ! Heure de la journee en fraction d'1 jour
+!      INTEGER  iday ! jour julien
+      REAL       time 
 
       REAL  SSUM
@@ -125 +127 @@
       real time_step, t_wrt, t_ops
 
-      REAL rdayvrai,rdaym_ini
+!      REAL rdayvrai,rdaym_ini
+! jD_cur: jour julien courant
+! jH_cur: heure julienne courante
+      REAL :: jD_cur, jH_cur
+      INTEGER :: an, mois, jour
+      REAL :: secondes
+
       LOGICAL first,callinigrads
 cIM : pour sortir les param. du modele dans un fis. netcdf 110106
       save first
       data first/.true./
-      real dt_cum, zjulian
+      real dt_cum
       character*10 infile
       integer zan, tau0, thoriid
@@ -167 +175 @@
       character*80 abort_message
 
-C Calendrier
-      LOGICAL true_calendar
-      PARAMETER (true_calendar = .false.)
-
       logical dissip_conservative
       save dissip_conservative
@@ -193 +197 @@
 
       itau = 0
-      iday = day_ini+itau/day_step
-      time = FLOAT(itau-(iday-day_ini)*day_step)/day_step+time_0
-         IF(time.GT.1.) THEN
-          time = time-1.
-          iday = iday+1
-         ENDIF
+c$$$      iday = day_ini+itau/day_step
+c$$$      time = FLOAT(itau-(iday-day_ini)*day_step)/day_step+time_0
+c$$$         IF(time.GT.1.) THEN
+c$$$          time = time-1.
+c$$$          iday = iday+1
+c$$$         ENDIF
 
 
@@ -214 +218 @@
 
    1  CONTINUE
+
+      jD_cur = jD_ref + (day_ini - 1) + int (itau * dtvr / daysec) 
+      jH_cur = jH_ref +                                                 &
+     &          (itau * dtvr / daysec - int(itau * dtvr / daysec)) 
 
 
@@ -284 +292 @@
       CALL geopot  ( ip1jmp1, teta  , pk , pks,  phis  , phi   )
 
+      time = jD_cur + jH_cur
       CALL caldyn 
      $  ( itau,ucov,vcov,teta,ps,masse,pk,pkf,phis ,
-     $    phi,conser,du,dv,dteta,dp,w, pbaru,pbarv, time+iday-day_ini )
+     $    phi,conser,du,dv,dteta,dp,w, pbaru,pbarv, time )
 
 
@@ -345 +354 @@
          CALL exner_hyb(  ip1jmp1, ps, p,alpha,beta,pks, pk, pkf )
 
-           rdaym_ini  = itau * dtvr / daysec
-           rdayvrai   = rdaym_ini  + day_ini
-
+!           rdaym_ini  = itau * dtvr / daysec
+!           rdayvrai   = rdaym_ini  + day_ini
+           jD_cur = jD_ref + (day_ini - 1) + int (itau * dtvr / daysec) 
+           jH_cur = jH_ref +                                            &
+     &              (itau * dtvr / daysec - int(itau * dtvr / daysec)) 
+!         write(lunout,*)'itau, jD_cur = ', itau, jD_cur, jH_cur
+!         call ju2ymds(jD_cur+jH_cur, an, mois, jour, secondes)
+!         write(lunout,*)'current date = ',an, mois, jour, secondes 
 
 c rajout debug
@@ -379 +393 @@
 #endif
 ! #endif of #ifdef CPP_IOIPSL
-         CALL calfis( lafin ,rdayvrai,time  ,
+         CALL calfis( lafin , jD_cur, jH_cur,
      $               ucov,vcov,teta,q,masse,ps,p,pk,phis,phi ,
      $               du,dv,dteta,dq,
@@ -507 +521 @@
             IF(forward. OR. leapf) THEN
               itau= itau + 1
-              iday= day_ini+itau/day_step
-              time= FLOAT(itau-(iday-day_ini)*day_step)/day_step+time_0
-                IF(time.GT.1.) THEN
-                  time = time-1.
-                  iday = iday+1
-                ENDIF
+c$$$              iday= day_ini+itau/day_step
+c$$$              time= FLOAT(itau-(iday-day_ini)*day_step)/day_step+time_0
+c$$$                IF(time.GT.1.) THEN
+c$$$                  time = time-1.
+c$$$                  iday = iday+1
+c$$$                ENDIF
             ENDIF
 
@@ -632 +646 @@
 
              itau =  itau + 1
-             iday = day_ini+itau/day_step
-             time = FLOAT(itau-(iday-day_ini)*day_step)/day_step+time_0
-
-                  IF(time.GT.1.) THEN
-                   time = time-1.
-                   iday = iday+1
-                  ENDIF
+c$$$             iday = day_ini+itau/day_step
+c$$$             time = FLOAT(itau-(iday-day_ini)*day_step)/day_step+time_0
+c$$$
+c$$$                  IF(time.GT.1.) THEN
+c$$$                   time = time-1.
+c$$$                   iday = iday+1
+c$$$                  ENDIF
 
                forward =  .FALSE.