1 | ! $Id: initdynav_p.F 1279 2009-12-10 09:02:56Z fairhead $ |
---|
2 | |
---|
3 | SUBROUTINE 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) |
---|
277 | END SUBROUTINE initdynav_loc |
---|