source: LMDZ5/trunk/libf/phylmd/ini_undefSTD.F @ 1912

Last change on this file since 1912 was 1912, checked in by musat, 10 years ago

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
  • 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: 2.4 KB
Line 
1!
2! $Id: ini_undefSTD.F 1912 2013-12-05 17:32:35Z musat $
3!
4
5      SUBROUTINE ini_undefSTD(itap,itapm1)
6      USE dimphy
7      USE phys_state_var_mod ! Variables sauvegardees de la physique
8      USE phys_cal_mod, only : mth_len
9      IMPLICIT none
10      include "clesphys.h"
11c
12c====================================================================
13c
14c I. Musat : 09.2004
15c
16c Initialisation - a des frequences differentes :
17c
18c 1) des variables moyennees sur la journee "day" ou sur le mois "mth"
19c    calculees a partir des valeurs "instantannees" de la physique
20c
21c 2) des variables moyennes mensuelles "NMC" calculees a partir des val.
22c    toutes les 6 heures
23c
24c nout=1 !var. journaliere "day" moyenne sur tous les pas de temps
25c              ! de la physique
26c nout=2 !var. mensuelle "mth" moyennee sur tous les pas de temps
27c              ! de la physique
28c nout=3 !var. mensuelle "NMC" moyennee toutes les 6heures
29c
30c NB: mettre "inst(X)" dans le write_hist*NMC.h !
31c====================================================================
32c
33cym #include "dimensions.h"
34cym      integer jjmp1
35cym      parameter (jjmp1=jjm+1-1/jjm)
36cym #include "dimphy.h"
37c variables Input/Output
38c     INTEGER nlevSTD, klevSTD, itap
39      INTEGER itap, itapm1
40c     PARAMETER(klevSTD=17)
41c     REAL dtime
42c
43c variables locales
44c     INTEGER i, k, nout, n
45      INTEGER i, k, n
46c     PARAMETER(nout=3) !nout=1 day/nout=2 mth/nout=3 NMC
47c     REAL freq_outNMC(nout)
48      REAL un_jour
49      PARAMETER(un_jour=86400.)
50c
51c variables Output
52c     REAL tnondef(klon,klevSTD,nout)
53c     REAL tsumSTD(klon,klevSTD,nout)
54c
55      DO n=1, nout
56c
57c initialisation variables en debut du mois, de la journee ou des 6h
58c
59       IF(n.EQ.1.AND.itap-itapm1.EQ.1.OR.
60     &n.GT.1.AND.MOD(itap,NINT(freq_outNMC(n)/dtime)).EQ.1) THEN
61c       print*,'ini_undefSTD n itap',n,itap
62        DO k=1, nlevSTD
63         DO i=1, klon
64          tnondef(i,k,n)=0.
65          tsumSTD(i,k,n)=0.
66          usumSTD(i,k,n)=0.
67          vsumSTD(i,k,n)=0.
68          wsumSTD(i,k,n)=0.
69          phisumSTD(i,k,n)=0.
70          qsumSTD(i,k,n)=0.
71          rhsumSTD(i,k,n)=0.
72          uvsumSTD(i,k,n)=0.
73          vqsumSTD(i,k,n)=0.
74          vTsumSTD(i,k,n)=0.
75          wqsumSTD(i,k,n)=0.
76          vphisumSTD(i,k,n)=0.
77          wTsumSTD(i,k,n)=0.
78          u2sumSTD(i,k,n)=0.
79          v2sumSTD(i,k,n)=0.
80          T2sumSTD(i,k,n)=0.
81          O3sumSTD(i,k,n)=0.
82          O3daysumSTD(i,k,n)=0.
83         ENDDO !i
84        ENDDO !k
85       ENDIF !
86      ENDDO !n
87      RETURN
88      END 
Note: See TracBrowser for help on using the repository browser.