Changeset 1912 for LMDZ5/trunk/libf/phylmd/calcul_divers.h

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/calcul_divers.h (modified) (3 diffs)

LMDZ5/trunk/libf/phylmd/calcul_divers.h

@@ -2 +2 @@
 ! $Header$
 !
-
-!     Initialisations diverses au "debut" du mois
-      IF(debut) THEN
-         nday_rain(:)=0.
-
+! Initialisations diverses au tout debut
+      IF(itap.EQ.1) THEN
+         itapm1=0
 !        surface terre
-         paire_ter(:)=0.
          DO i=1, klon
             IF(pctsrf(i,is_ter).GT.0.) THEN
@@ -16 +13 @@
       ENDIF
 
-!IM   Calcul une fois par jour : total_rain, nday_rain
-      IF(MOD(itap,INT(un_jour/dtime)).EQ.0) THEN
+! Initialisation debut de mois
+      IF(itap.EQ.itapm1+1) THEN
+        nday_rain(:)=0.
+!       print*,'initialisation mois suivants day_rain itap',itap
+      ENDIF
+
+! Calcul fin de journee : total_rain, nday_rain
+      IF(MOD(itap,NINT(un_jour/dtime)).EQ.0) THEN
+!        print*,'calcul nday_rain itap ',itap
          DO i = 1, klon
             total_rain(i)=rain_fall(i)+snow_fall(i)  
@@ -23 +27 @@
          ENDDO
       ENDIF
+
+! Initialisation fin de mois
+      IF(MOD(itap-itapm1,NINT(mth_len*un_jour/dtime)).EQ.0) THEN
+        itapm1=itapm1+NINT(mth_len*un_jour/dtime)
+!       print*,'initialisation itapm1 ',itapm1
+      ENDIF