Changeset 1912 for LMDZ5/trunk/libf/phylmd/ini_undefSTD.F

Timestamp:

Dec 5, 2013, 6:32:35 PM (10 years ago)

Author:

musat

Message:

1) Modifications pour faire des simulations par an avec un calendrier realiste (365 jours ou autre).

Il faut mettre une frequence de sortie de -1 (variable phys_out_filetimesteps dans config.def) pour
que IOIPSL calcule les moyennes mensuels en prenant en compte des longuers variables de chaque
mois. Par exemple, pour le fichier histmth (1er fichier) et histmthNMC (7eme).
phys_out_filetimesteps= -1 1day 6hr 6hr 6hr 1d -1 1day 6hr

2) Corrections titres variables niveaux de pression des fichiers histmth, histday, etc

File:

• LMDZ5/trunk/libf/phylmd/ini_undefSTD.F (modified) (5 diffs)

LMDZ5/trunk/libf/phylmd/ini_undefSTD.F

@@ -3 +3 @@
 !
 
-      SUBROUTINE ini_undefSTD(itap,
-     $           freq_outNMC)
+      SUBROUTINE ini_undefSTD(itap,itapm1)
       USE dimphy
       USE phys_state_var_mod ! Variables sauvegardees de la physique
+      USE phys_cal_mod, only : mth_len
       IMPLICIT none
+      include "clesphys.h"
 c
 c====================================================================
@@ -36 +37 @@
 c variables Input/Output
 c     INTEGER nlevSTD, klevSTD, itap
-      INTEGER itap
+      INTEGER itap, itapm1
 c     PARAMETER(klevSTD=17)
 c     REAL dtime
@@ -44 +45 @@
       INTEGER i, k, n
 c     PARAMETER(nout=3) !nout=1 day/nout=2 mth/nout=3 NMC
-      REAL freq_outNMC(nout)
+c     REAL freq_outNMC(nout)
+      REAL un_jour
+      PARAMETER(un_jour=86400.)
 c
 c variables Output
@@ -52 +55 @@
       DO n=1, nout
 c
-c initialisation variables en debut de la journee ou du mois
+c initialisation variables en debut du mois, de la journee ou des 6h
 c
-       IF(MOD(itap,NINT(freq_outNMC(n)/dtime)).EQ.1.) THEN
+       IF(n.EQ.1.AND.itap-itapm1.EQ.1.OR.
+     &n.GT.1.AND.MOD(itap,NINT(freq_outNMC(n)/dtime)).EQ.1) THEN
+c       print*,'ini_undefSTD n itap',n,itap
         DO k=1, nlevSTD
          DO i=1, klon
@@ -78 +83 @@
          ENDDO !i
         ENDDO !k
-c
-       ENDIF !MOD(itap,NINT(freq_outNMC(n)/dtime)).EQ.1.
-c
+       ENDIF !
       ENDDO !n
-c
       RETURN
       END