My Project
 All Classes Files Functions Variables Macros
ini_undefSTD.F
Go to the documentation of this file.
1 !
2 ! $Id: ini_undefSTD.F 1403 2010-07-01 09:02:53Z fairhead $
3 !
4 
5  SUBROUTINE ini_undefstd(itap,
6  $ freq_outnmc)
7  USE dimphy
8  USE phys_state_var_mod ! Variables sauvegardees de la physique
9  IMPLICIT none
10 c
11 c====================================================================
12 c
13 c I. Musat : 09.2004
14 c
15 c Initialisation - a des frequences differentes :
16 c
17 c 1) des variables moyennees sur la journee "day" ou sur le mois "mth"
18 c calculees a partir des valeurs "instantannees" de la physique
19 c
20 c 2) des variables moyennes mensuelles "NMC" calculees a partir des val.
21 c toutes les 6 heures
22 c
23 c nout=1 !var. journaliere "day" moyenne sur tous les pas de temps
24 c ! de la physique
25 c nout=2 !var. mensuelle "mth" moyennee sur tous les pas de temps
26 c ! de la physique
27 c nout=3 !var. mensuelle "NMC" moyennee toutes les 6heures
28 c
29 c NB: mettre "inst(X)" dans le write_hist*NMC.h !
30 c====================================================================
31 c
32 cym #include "dimensions.h"
33 cym integer jjmp1
34 cym parameter (jjmp1=jjm+1-1/jjm)
35 cym #include "dimphy.h"
36 c variables Input/Output
37 c INTEGER nlevSTD, klevSTD, itap
38  INTEGER itap
39 c PARAMETER(klevSTD=17)
40 c REAL dtime
41 c
42 c variables locales
43 c INTEGER i, k, nout, n
44  INTEGER i, k, n
45 c PARAMETER(nout=3) !nout=1 day/nout=2 mth/nout=3 NMC
46  REAL freq_outnmc(nout)
47 c
48 c variables Output
49 c REAL tnondef(klon,klevSTD,nout)
50 c REAL tsumSTD(klon,klevSTD,nout)
51 c
52  DO n=1, nout
53 c
54 c initialisation variables en debut de la journee ou du mois
55 c
56  IF(mod(itap,nint(freq_outnmc(n)/dtime)).EQ.1.) THEN
57  DO k=1, nlevstd
58  DO i=1, klon
59  tnondef(i,k,n)=0.
60  tsumstd(i,k,n)=0.
61  usumstd(i,k,n)=0.
62  vsumstd(i,k,n)=0.
63  wsumstd(i,k,n)=0.
64  phisumstd(i,k,n)=0.
65  qsumstd(i,k,n)=0.
66  rhsumstd(i,k,n)=0.
67  uvsumstd(i,k,n)=0.
68  vqsumstd(i,k,n)=0.
69  vtsumstd(i,k,n)=0.
70  wqsumstd(i,k,n)=0.
71  vphisumstd(i,k,n)=0.
72  wtsumstd(i,k,n)=0.
73  u2sumstd(i,k,n)=0.
74  v2sumstd(i,k,n)=0.
75  t2sumstd(i,k,n)=0.
76  o3sumstd(i,k,n)=0.
77  o3daysumstd(i,k,n)=0.
78  ENDDO !i
79  ENDDO !k
80 c
81  ENDIF !MOD(itap,NINT(freq_outNMC(n)/dtime)).EQ.1.
82 c
83  ENDDO !n
84 c
85  RETURN
86  END