Changeset 1921 for LMDZ5/branches/testing/libf/phylmd/moy_undefSTD.F

Timestamp:

Dec 20, 2013, 10:04:40 AM (10 years ago)

Author:

Laurent Fairhead

Message:

Merged trunk changes r1909:1920 into testing branch

Location:

LMDZ5/branches/testing

Files:

• . (modified) (1 prop)
• libf/phylmd/moy_undefSTD.F (modified) (7 diffs)

LMDZ5/branches/testing/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)
+       ELSE
+         freq_moyNMC(n)=freq_outNMC(n)/freq_calNMC(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 +76 @@
            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 +110 @@
           O3daysumSTD(i,k,n)=O3daysumSTD(i,k,n)/
      $     (freq_moyNMC(n)-tnondef(i,k,n))
-cIM END
           ELSE
            tsumSTD(i,k,n)=missing_val