source: LMDZ6/branches/Amaury_dev/libf/dyn3dmem/initdynav_loc.f90

Last change on this file was 5195, checked in by abarral, 8 weeks ago

Correct r5192, some lmdz_description cases were missing

  • 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
File size: 7.0 KB
Line 
1! $Id: initdynav_p.F 1279 2009-12-10 09:02:56Z fairhead $
2
3SUBROUTINE initdynav_loc(day0, anne0, tstep, t_ops, t_wrt)
4
5  ! This routine needs IOIPSL
6  USE IOIPSL
7  USE parallel_lmdz
8  USE lmdz_write_field
9  USE misc_mod
10  ! USE lmdz_infotrac
11  USE com_io_dyn_mod, ONLY: histaveid, histvaveid, histuaveid, &
12          dynhistave_file, dynhistvave_file, dynhistuave_file
13  USE comconst_mod, ONLY: pi
14  USE comvert_mod, ONLY: presnivs
15  USE temps_mod, ONLY: itau_dyn
16  USE lmdz_iniprint, ONLY: lunout, prt_level
17  USE lmdz_comgeom
18
19  USE lmdz_dimensions, ONLY: iim, jjm, llm, ndm
20  USE lmdz_paramet
21  IMPLICIT NONE
22
23
24  !   Routine d'initialisation des ecritures des fichiers histoires LMDZ
25  !   au format IOIPSL. Initialisation du fichier histoire moyenne.
26
27  !   Appels succesifs des routines: histbeg
28  !                              histhori
29  !                              histver
30  !                              histdef
31  !                              histend
32
33  !   Entree:
34
35  !  day0,anne0: date de reference
36  !  tstep : frequence d'ecriture
37  !  t_ops: frequence de l'operation pour IOIPSL
38  !  t_wrt: frequence d'ecriture sur le fichier
39
40  !   Sortie:
41  !  fileid: ID du fichier netcdf cree
42
43  !   L. Fairhead, LMD, 03/99
44
45  ! =====================================================================
46
47  !   Declarations
48
49
50
51  !   Arguments
52
53  INTEGER(kind = 4) day0, anne0
54  REAL :: tstep, t_ops, t_wrt
55
56  ! This routine needs IOIPSL
57  !   Variables locales
58
59  INTEGER :: tau0
60  REAL :: zjulian
61  INTEGER :: iq
62  REAL :: rlong(iip1, jjp1), rlat(iip1, jjp1)
63  INTEGER :: uhoriid, vhoriid, thoriid
64  INTEGER :: zvertiid, zvertiidv, zvertiidu
65  INTEGER :: ii, jj
66  INTEGER :: zan, dayref
67  INTEGER :: jjb, jje, jjn
68
69  ! definition du domaine d'ecriture pour le rebuild
70
71  INTEGER, DIMENSION(2) :: ddid
72  INTEGER, DIMENSION(2) :: dsg
73  INTEGER, DIMENSION(2) :: dsl
74  INTEGER, DIMENSION(2) :: dpf
75  INTEGER, DIMENSION(2) :: dpl
76  INTEGER, DIMENSION(2) :: dhs
77  INTEGER, DIMENSION(2) :: dhe
78
79  INTEGER :: dynhistave_domain_id
80  INTEGER :: dynhistvave_domain_id
81  INTEGER :: dynhistuave_domain_id
82
83  IF (adjust) return
84
85
86  !  Initialisations
87
88  pi = 4. * atan (1.)
89
90  !  Appel a histbeg: creation du fichier netcdf et initialisations diverses
91  !
92
93  zan = anne0
94  dayref = day0
95  CALL ymds2ju(zan, 1, dayref, 0.0, zjulian)
96  tau0 = itau_dyn
97
98  DO jj = 1, jjp1
99    DO ii = 1, iip1
100      rlong(ii, jj) = rlonv(ii) * 180. / pi
101      rlat(ii, jj) = rlatu(jj) * 180. / pi
102    enddo
103  enddo
104
105
106  ! Creation de 3 fichiers pour les differentes grilles horizontales
107  ! Restriction de IOIPSL: seulement 2 coordonnees dans le meme fichier
108  ! Grille Scalaire
109
110  jjb = jj_begin
111  jje = jj_end
112  jjn = jj_nb
113
114  ddid = (/ 1, 2 /)
115  dsg = (/ iip1, jjp1 /)
116  dsl = (/ iip1, jjn /)
117  dpf = (/ 1, jjb /)
118  dpl = (/ iip1, jje /)
119  dhs = (/ 0, 0 /)
120  dhe = (/ 0, 0 /)
121
122  CALL flio_dom_set(mpi_size, mpi_rank, ddid, dsg, dsl, dpf, dpl, dhs, dhe, &
123          'box', dynhistave_domain_id)
124
125  CALL histbeg(dynhistave_file, iip1, rlong(:, 1), jjn, &
126          rlat(1, jjb:jje), 1, iip1, 1, jjn, tau0, &
127          zjulian, tstep, thoriid, &
128          histaveid, dynhistave_domain_id)
129
130
131  !  Creation du fichier histoire pour les grilles en V et U (oblige pour l'instant,
132  !  IOIPSL ne permet pas de grilles avec des nombres de point differents dans
133  !  un meme fichier)
134  ! Grille V
135
136  jjb = jj_begin
137  jje = jj_end
138  jjn = jj_nb
139  IF (pole_sud) jjn = jjn - 1
140  IF (pole_sud) jje = jje - 1
141
142  DO jj = jjb, jje
143    DO ii = 1, iip1
144      rlong(ii, jj) = rlonv(ii) * 180. / pi
145      rlat(ii, jj) = rlatv(jj) * 180. / pi
146    enddo
147  enddo
148
149  ddid = (/ 1, 2 /)
150  dsg = (/ iip1, jjm /)
151  dsl = (/ iip1, jjn /)
152  dpf = (/ 1, jjb /)
153  dpl = (/ iip1, jje /)
154  dhs = (/ 0, 0 /)
155  dhe = (/ 0, 0 /)
156
157  CALL flio_dom_set(mpi_size, mpi_rank, ddid, dsg, dsl, dpf, dpl, dhs, dhe, &
158          'box', dynhistvave_domain_id)
159
160  CALL histbeg(dynhistvave_file, iip1, rlong(:, 1), jjn, &
161          rlat(1, jjb:jje), 1, iip1, 1, jjn, tau0, &
162          zjulian, tstep, vhoriid, &
163          histvaveid, dynhistvave_domain_id)
164
165  ! Grille U
166
167  DO jj = 1, jjp1
168    DO ii = 1, iip1
169      rlong(ii, jj) = rlonu(ii) * 180. / pi
170      rlat(ii, jj) = rlatu(jj) * 180. / pi
171    enddo
172  enddo
173
174  jjb = jj_begin
175  jje = jj_end
176  jjn = jj_nb
177
178  ddid = (/ 1, 2 /)
179  dsg = (/ iip1, jjp1 /)
180  dsl = (/ iip1, jjn /)
181  dpf = (/ 1, jjb /)
182  dpl = (/ iip1, jje /)
183  dhs = (/ 0, 0 /)
184  dhe = (/ 0, 0 /)
185
186  CALL flio_dom_set(mpi_size, mpi_rank, ddid, dsg, dsl, dpf, dpl, dhs, dhe, &
187          'box', dynhistuave_domain_id)
188
189  CALL histbeg(dynhistuave_file, iip1, rlong(:, 1), jjn, &
190          rlat(1, jjb:jje), 1, iip1, 1, jjn, tau0, &
191          zjulian, tstep, uhoriid, &
192          histuaveid, dynhistuave_domain_id)
193
194
195
196  !  Appel a histvert pour la grille verticale
197
198  CALL histvert(histaveid, 'presnivs', 'Niveaux Pression&
199          &     approximatifs', 'mb', llm, presnivs / 100., zvertiid, 'down')
200  CALL histvert(histuaveid, 'presnivs', 'Niveaux Pression&
201          &     approximatifs', 'mb', llm, presnivs / 100., zvertiidv, 'down')
202  CALL histvert(histvaveid, 'presnivs', 'Niveaux Pression&
203          &     approximatifs', 'mb', llm, presnivs / 100., zvertiidu, 'down')
204
205
206  !  Appels a histdef pour la definition des variables a sauvegarder
207
208  !  Vents U
209
210  jjn = jj_nb
211  CALL histdef(histuaveid, 'u', 'vent u moyen ', &
212          'm/s', iip1, jjn, uhoriid, llm, 1, llm, zvertiidu, &
213          32, 'ave(X)', t_ops, t_wrt)
214
215
216  !  Vents V
217
218  IF (pole_sud) jjn = jj_nb - 1
219  CALL histdef(histvaveid, 'v', 'vent v moyen', &
220          'm/s', iip1, jjn, vhoriid, llm, 1, llm, zvertiidv, &
221          32, 'ave(X)', t_ops, t_wrt)
222
223
224  !  Temperature
225
226  jjn = jj_nb
227  CALL histdef(histaveid, 'temp', 'temperature moyenne', 'K', &
228          iip1, jjn, thoriid, llm, 1, llm, zvertiid, &
229          32, 'ave(X)', t_ops, t_wrt)
230
231  !  Temperature potentielle
232
233  CALL histdef(histaveid, 'theta', 'temperature potentielle', 'K', &
234          iip1, jjn, thoriid, llm, 1, llm, zvertiid, &
235          32, 'ave(X)', t_ops, t_wrt)
236
237
238
239  !  Geopotentiel
240
241  CALL histdef(histaveid, 'phi', 'geopotentiel moyen', '-', &
242          iip1, jjn, thoriid, llm, 1, llm, zvertiid, &
243          32, 'ave(X)', t_ops, t_wrt)
244
245  !  Traceurs
246
247  !    DO iq=1,nqtot
248  !      CALL histdef(histaveid, tracers(iq)%name,
249  ! .                            tracers(iq)%longName, '-',
250  ! .             iip1, jjn, thoriid, llm, 1, llm, zvertiid,
251  ! .             32, 'ave(X)', t_ops, t_wrt)
252  !    enddo
253
254  !  Masse
255
256  CALL histdef(histaveid, 'masse', 'masse moyenne', 'kg', &
257          iip1, jjn, thoriid, llm, 1, llm, zvertiid, &
258          32, 'ave(X)', t_ops, t_wrt)
259
260  !  Pression au sol
261
262  CALL histdef(histaveid, 'ps', 'pression naturelle au sol', 'Pa', &
263          iip1, jjn, thoriid, 1, 1, 1, -99, &
264          32, 'ave(X)', t_ops, t_wrt)
265
266  !  Geopotentiel au sol
267
268  !  CALL histdef(histaveid, 'phis', 'geopotentiel au sol', '-',
269  ! .             iip1, jjn, thoriid, 1, 1, 1, -99,
270  ! .             32, 'ave(X)', t_ops, t_wrt)
271
272  !  Fin
273
274  CALL histend(histaveid)
275  CALL histend(histuaveid)
276  CALL histend(histvaveid)
277END SUBROUTINE initdynav_loc
Note: See TracBrowser for help on using the repository browser.