source: LMDZ.3.3/branches/rel-LF/libf/bibio/initdynav.F @ 402

Last change on this file since 402 was 362, checked in by lmdzadmin, 22 years ago

Resynchronisation de la sortie de la temperature potentielle
LF

  • Property svn:eol-style set to native
  • Property svn:keywords set to Author Date Id Revision
File size: 4.4 KB
Line 
1c
2c $Header$
3c
4      subroutine initdynav(infile,day0,anne0,tstep,t_ops,t_wrt
5     .                     ,nq,fileid)
6
7       USE IOIPSL
8       USE histcom
9
10      implicit none
11
12C
13C   Routine d'initialisation des ecritures des fichiers histoires LMDZ
14C   au format IOIPSL. Initialisation du fichier histoire moyenne.
15C
16C   Appels succesifs des routines: histbeg
17C                                  histhori
18C                                  histver
19C                                  histdef
20C                                  histend
21C
22C   Entree:
23C
24C      infile: nom du fichier histoire a creer
25C      day0,anne0: date de reference
26C      tstep : frequence d'ecriture
27C      t_ops: frequence de l'operation pour IOIPSL
28C      t_wrt: frequence d'ecriture sur le fichier
29C      nq: nombre de traceurs
30C
31C   Sortie:
32C      fileid: ID du fichier netcdf cree
33C
34C   L. Fairhead, LMD, 03/99
35C
36C =====================================================================
37C
38C   Declarations
39#include "dimensions.h"
40#include "paramet.h"
41#include "comconst.h"
42#include "comvert.h"
43#include "comgeom.h"
44#include "temps.h"
45#include "ener.h"
46#include "logic.h"
47#include "description.h"
48#include "serre.h"
49
50C   Arguments
51C
52      character*(*) infile
53      integer*4 day0, anne0
54      real tstep, t_ops, t_wrt
55      integer fileid
56      integer nq
57      integer thoriid, zvertiid
58
59C   Variables locales
60C
61      integer tau0
62      real zjulian
63      character*3 str
64      character*10 ctrac
65      integer iq
66      real rlong(iip1,jjp1), rlat(iip1,jjp1)
67      integer ii,jj
68      integer zan, dayref
69C
70C  Initialisations
71C
72      pi = 4. * atan (1.)
73      str='q  '
74      ctrac = 'traceur   '
75C
76C  Appel a histbeg: creation du fichier netcdf et initialisations diverses
77C         
78
79      zan = anne0
80      dayref = day0
81      CALL ymds2ju(zan, 1, dayref, 0.0, zjulian)
82      tau0 = itau_dyn
83     
84      do jj = 1, jjp1
85        do ii = 1, iip1
86          rlong(ii,jj) = rlonv(ii) * 180. / pi
87          rlat(ii,jj)  = rlatu(jj) * 180. / pi
88        enddo
89      enddo
90       
91      call histbeg(infile, iip1, rlong(:,1), jjp1, rlat(1,:),
92     .             1, iip1, 1, jjp1,
93     .             tau0, zjulian, tstep, thoriid, fileid)
94
95C
96C  Appel a histvert pour la grille verticale
97C
98      call histvert(fileid, 'sigss', 'Niveaux sigma','Pa',
99     .              llm, nivsigs, zvertiid)
100C
101C  Appels a histdef pour la definition des variables a sauvegarder
102C
103C  Vents U
104C
105      write(6,*)'inithistave',tstep
106      call histdef(fileid, 'u', 'vents u scalaires moyennes',
107     .             'm/s', iip1, jjp1, thoriid, llm, 1, llm, zvertiid,
108     .             32, 'ave(X)', t_ops, t_wrt)
109
110C
111C  Vents V
112C
113      call histdef(fileid, 'v', 'vents v scalaires moyennes',
114     .             'm/s', iip1, jjp1, thoriid, llm, 1, llm, zvertiid,
115     .             32, 'ave(X)', t_ops, t_wrt)
116
117C
118C  Temperature
119C
120      call histdef(fileid, 'temp', 'temperature moyennee', 'K',
121     .             iip1, jjp1, thoriid, llm, 1, llm, zvertiid,
122     .             32, 'ave(X)', t_ops, t_wrt)
123C
124C  Temperature potentielle
125C
126      call histdef(fileid, 'theta', 'temperature potentielle', 'K',
127     .             iip1, jjp1, thoriid, llm, 1, llm, zvertiid,
128     .             32, 'ave(X)', t_ops, t_wrt)
129
130
131C
132C  Geopotentiel
133C
134      call histdef(fileid, 'phi', 'geopotentiel moyenne', '-',
135     .             iip1, jjp1, thoriid, llm, 1, llm, zvertiid,
136     .             32, 'ave(X)', t_ops, t_wrt)
137C
138C  Traceurs
139C
140      IF(nq.GE.1) THEN
141        DO iq=1,nq
142          IF ( iq.LE.9 ) THEN
143            WRITE(str(2:2),'(i1.1)') iq
144            write(ctrac(9:9),'(i1.1)') iq
145          ELSE
146            WRITE(str(2:3),'(i2.2)') iq
147            write(ctrac(9:10),'(i2.2)') iq
148          ENDIF
149          call histdef(fileid, str, ctrac, '-',
150     .             iip1, jjp1, thoriid, llm, 1, llm, zvertiid,
151     .             32, 'ave(X)', t_ops, t_wrt)
152        enddo
153      endif
154C
155C  Masse
156C
157      call histdef(fileid, 'masse', 'masse', 'kg',
158     .             iip1, jjp1, thoriid, 1, 1, 1, -99,
159     .             32, 'ave(X)', t_ops, t_wrt)
160C
161C  Pression au sol
162C
163      call histdef(fileid, 'ps', 'pression naturelle au sol', 'Pa',
164     .             iip1, jjp1, thoriid, 1, 1, 1, -99,
165     .             32, 'ave(X)', t_ops, t_wrt)
166C
167C  Pression au sol
168C
169      call histdef(fileid, 'phis', 'geopotentiel au sol', '-',
170     .             iip1, jjp1, thoriid, 1, 1, 1, -99,
171     .             32, 'ave(X)', t_ops, t_wrt)
172C
173C  Fin
174C
175      call histend(fileid)
176      return
177      end
Note: See TracBrowser for help on using the repository browser.