source: LMDZ5/trunk/libf/phylmd/moy_undefSTD.F @ 1907

Last change on this file since 1907 was 1907, checked in by lguez, 10 years ago

Added a copyright property to every file of the distribution, except
for the fcm files (which have their own copyright). Use svn propget on
a file to see the copyright. For instance:

$ svn propget copyright libf/phylmd/physiq.F90
Name of program: LMDZ
Creation date: 1984
Version: LMDZ5
License: CeCILL version 2
Holder: Laboratoire de m\'et\'eorologie dynamique, CNRS, UMR 8539
See the license file in the root directory

Also added the files defining the CeCILL version 2 license, in French
and English, at the top of the LMDZ tree.
  • Property copyright set to
    Name of program: LMDZ
    Creation date: 1984
    Version: LMDZ5
    License: CeCILL version 2
    Holder: Laboratoire de m\'et\'eorologie dynamique, CNRS, UMR 8539
    See the license file in the root directory
  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 4.0 KB
Line 
1!
2! $Id: moy_undefSTD.F 1907 2013-11-26 13:10:46Z lguez $
3!
4      SUBROUTINE moy_undefSTD(itap,freq_outNMC,freq_moyNMC)
5      USE netcdf
6      USE dimphy
7      USE phys_state_var_mod ! Variables sauvegardees de la physique
8      IMPLICIT none
9c
10c====================================================================
11c
12c I. Musat : 09.2004
13c
14c Moyenne - a des frequences differentes - des valeurs bien definies
15c         (.NE.missing_val) des variables interpolees a un niveau de
16c         pression.
17c 1) les variables de type "day" (nout=1) ou "mth" (nout=2) sont sommees
18c    tous les pas de temps de la physique
19c
20c 2) les variables de type "NMC" (nout=3) sont calculees a partir
21c    des valeurs instantannees toutes les 6 heures
22c
23c
24c NB: mettre "inst(X)" dans le write_hist*NMC.h !
25c====================================================================
26cym#include "dimensions.h"
27cym      integer jjmp1
28cym      parameter (jjmp1=jjm+1-1/jjm)
29cym#include "dimphy.h"
30c
31c
32c variables Input
33c     INTEGER nlevSTD, klevSTD, itap
34c     PARAMETER(klevSTD=17)
35      INTEGER itap
36c
37c variables locales
38c     INTEGER i, k, nout, n
39c     PARAMETER(nout=3) !nout=1 day/nout=2 mth/nout=3 NMC
40      INTEGER i, k, n
41c     REAL dtime, freq_outNMC(nout), freq_moyNMC(nout)
42      REAL freq_outNMC(nout), freq_moyNMC(nout)
43c
44c variables Output
45c     REAL tnondef(klon,klevSTD,nout)
46c     REAL tsumSTD(klon,klevSTD,nout)
47c
48      REAL missing_val
49c
50      missing_val=nf90_fill_real
51c
52      DO n=1, nout
53c
54c calcul 1 fois par jour
55c
56       IF(MOD(itap,NINT(freq_outNMC(n)/dtime)).EQ.0) THEN
57c
58        DO k=1, nlevSTD
59         DO i=1, klon
60          IF(tnondef(i,k,n).NE.(freq_moyNMC(n))) THEN
61           tsumSTD(i,k,n)=tsumSTD(i,k,n)/
62     $     (freq_moyNMC(n)-tnondef(i,k,n))
63cIM BEG
64          usumSTD(i,k,n)=usumSTD(i,k,n)/
65     $     (freq_moyNMC(n)-tnondef(i,k,n))
66          vsumSTD(i,k,n)=vsumSTD(i,k,n)/
67     $     (freq_moyNMC(n)-tnondef(i,k,n))
68          wsumSTD(i,k,n)=wsumSTD(i,k,n)/
69     $     (freq_moyNMC(n)-tnondef(i,k,n))
70          phisumSTD(i,k,n)=phisumSTD(i,k,n)/
71     $     (freq_moyNMC(n)-tnondef(i,k,n))
72          qsumSTD(i,k,n)=qsumSTD(i,k,n)/
73     $     (freq_moyNMC(n)-tnondef(i,k,n))
74          rhsumSTD(i,k,n)=rhsumSTD(i,k,n)/
75     $     (freq_moyNMC(n)-tnondef(i,k,n))
76          uvsumSTD(i,k,n)=uvsumSTD(i,k,n)/
77     $     (freq_moyNMC(n)-tnondef(i,k,n))
78          vqsumSTD(i,k,n)=vqsumSTD(i,k,n)/
79     $     (freq_moyNMC(n)-tnondef(i,k,n))
80          vTsumSTD(i,k,n)=vTsumSTD(i,k,n)/
81     $     (freq_moyNMC(n)-tnondef(i,k,n))
82          wqsumSTD(i,k,n)=wqsumSTD(i,k,n)/
83     $     (freq_moyNMC(n)-tnondef(i,k,n))
84          vphisumSTD(i,k,n)=vphisumSTD(i,k,n)/
85     $     (freq_moyNMC(n)-tnondef(i,k,n))
86          wTsumSTD(i,k,n)=wTsumSTD(i,k,n)/
87     $     (freq_moyNMC(n)-tnondef(i,k,n))
88          u2sumSTD(i,k,n)=u2sumSTD(i,k,n)/
89     $     (freq_moyNMC(n)-tnondef(i,k,n))
90          v2sumSTD(i,k,n)=v2sumSTD(i,k,n)/
91     $     (freq_moyNMC(n)-tnondef(i,k,n))
92          T2sumSTD(i,k,n)=T2sumSTD(i,k,n)/
93     $     (freq_moyNMC(n)-tnondef(i,k,n))
94          O3sumSTD(i,k,n)=O3sumSTD(i,k,n)/
95     $     (freq_moyNMC(n)-tnondef(i,k,n))
96          O3daysumSTD(i,k,n)=O3daysumSTD(i,k,n)/
97     $     (freq_moyNMC(n)-tnondef(i,k,n))
98cIM END
99          ELSE
100           tsumSTD(i,k,n)=missing_val
101           usumSTD(i,k,n)=missing_val
102           vsumSTD(i,k,n)=missing_val
103           wsumSTD(i,k,n)=missing_val
104           phisumSTD(i,k,n)=missing_val
105           qsumSTD(i,k,n)=missing_val
106           rhsumSTD(i,k,n)=missing_val
107           uvsumSTD(i,k,n)=missing_val
108           vqsumSTD(i,k,n)=missing_val
109           vTsumSTD(i,k,n)=missing_val
110           wqsumSTD(i,k,n)=missing_val
111           vphisumSTD(i,k,n)=missing_val
112           wTsumSTD(i,k,n)=missing_val
113           u2sumSTD(i,k,n)=missing_val
114           v2sumSTD(i,k,n)=missing_val
115           T2sumSTD(i,k,n)=missing_val
116           O3sumSTD(i,k,n)=missing_val
117           O3daysumSTD(i,k,n)=missing_val
118          ENDIF !tnondef(i,k,n).NE.(freq_moyNMC(n))
119         ENDDO !i
120        ENDDO !k
121       ENDIF !MOD(itap,NINT(freq_outNMC(n)/dtime)).EQ.0
122c
123      ENDDO !n
124c
125      RETURN
126      END 
Note: See TracBrowser for help on using the repository browser.