Changeset 1912 for LMDZ5/trunk/libf/phylmd/moy_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/moy_undefSTD.F (modified) (7 diffs)

LMDZ5/trunk/libf/phylmd/moy_undefSTD.F

@@ -2 +2 @@
 ! $Id$
 !
-      SUBROUTINE moy_undefSTD(itap,freq_outNMC,freq_moyNMC)
+      SUBROUTINE moy_undefSTD(itap,itapm1)
       USE netcdf
       USE dimphy
-      USE phys_state_var_mod ! Variables sauvegardees de la physique
+      USE phys_state_var_mod
+      USE phys_cal_mod, only : mth_len
       IMPLICIT none
+      include "clesphys.h"
 c
 c====================================================================
@@ -33 +35 @@
 c     INTEGER nlevSTD, klevSTD, itap
 c     PARAMETER(klevSTD=17)
-      INTEGER itap
+      INTEGER itap, itapm1
 c
 c variables locales
@@ -40 +42 @@
       INTEGER i, k, n
 c     REAL dtime, freq_outNMC(nout), freq_moyNMC(nout)
-      REAL freq_outNMC(nout), freq_moyNMC(nout)
+c     REAL freq_outNMC(nout), freq_calNMC(nout)
+      REAL freq_moyNMC(nout)
 c
 c variables Output
@@ -46 +49 @@
 c     REAL tsumSTD(klon,klevSTD,nout)
 c
+      REAL un_jour
+      PARAMETER(un_jour=86400.)
       REAL missing_val
 c
@@ -51 +56 @@
 c
       DO n=1, nout
+       IF(freq_outNMC(n).LT.0) THEN
+         freq_moyNMC(n)=(mth_len*un_jour)/freq_calNMC(n)
+c        print*,'moy_undefSTD n freq_out freq_moy =',
+c    $n,freq_moyNMC(n)
+       ENDIF         
 c 
-c calcul 1 fois par jour
+c calcul 1 fois pas mois, 1 fois par jour ou toutes les 6h
 c
-       IF(MOD(itap,NINT(freq_outNMC(n)/dtime)).EQ.0) THEN
+       IF(n.EQ.1.AND.itap.EQ.itapm1.OR.
+     $n.GT.1.AND.MOD(itap,NINT(freq_outNMC(n)/dtime)).EQ.0) THEN
+c
+c       print*,'moy_undefSTD n itap itapm1',n,itap,itapm1
 c
         DO k=1, nlevSTD
@@ -61 +74 @@
            tsumSTD(i,k,n)=tsumSTD(i,k,n)/
      $     (freq_moyNMC(n)-tnondef(i,k,n))
-cIM BEG
           usumSTD(i,k,n)=usumSTD(i,k,n)/
      $     (freq_moyNMC(n)-tnondef(i,k,n))
@@ -96 +108 @@
           O3daysumSTD(i,k,n)=O3daysumSTD(i,k,n)/
      $     (freq_moyNMC(n)-tnondef(i,k,n))
-cIM END
           ELSE
            tsumSTD(i,k,n)=missing_val