1 | ! |
---|
2 | ! $Header$ |
---|
3 | ! |
---|
4 | c |
---|
5 | c |
---|
6 | SUBROUTINE phyetat0 (fichnom, |
---|
7 | . ocean_in, ok_veget_in, |
---|
8 | . clesphy0, |
---|
9 | . tabcntr0) |
---|
10 | |
---|
11 | USE dimphy |
---|
12 | USE mod_grid_phy_lmdz |
---|
13 | USE mod_phys_lmdz_para |
---|
14 | USE iophy |
---|
15 | USE ocean_slab_mod, ONLY : ocean_slab_init |
---|
16 | USE ocean_cpl_mod, ONLY : ocean_cpl_init |
---|
17 | USE ocean_forced_mod, ONLY : ocean_forced_init |
---|
18 | USE fonte_neige_mod, ONLY : fonte_neige_init |
---|
19 | USE pbl_surface_mod, ONLY : pbl_surface_init |
---|
20 | USE surface_data, ONLY : ocean, ok_veget |
---|
21 | USE phys_state_var_mod |
---|
22 | |
---|
23 | IMPLICIT none |
---|
24 | c====================================================================== |
---|
25 | c Auteur(s) Z.X. Li (LMD/CNRS) date: 19930818 |
---|
26 | c Objet: Lecture de l'etat initial pour la physique |
---|
27 | c====================================================================== |
---|
28 | #include "dimensions.h" |
---|
29 | #include "netcdf.inc" |
---|
30 | #include "indicesol.h" |
---|
31 | #include "dimsoil.h" |
---|
32 | #include "clesphys.h" |
---|
33 | #include "temps.h" |
---|
34 | #include "thermcell.h" |
---|
35 | #include "compbl.h" |
---|
36 | c====================================================================== |
---|
37 | CHARACTER*(*) fichnom |
---|
38 | |
---|
39 | c les variables globales lues dans le fichier restart |
---|
40 | REAL rlat_glo(klon_glo), rlon_glo(klon_glo) |
---|
41 | REAL pctsrf_glo(klon_glo, nbsrf) |
---|
42 | REAL tsol_glo(klon_glo,nbsrf) |
---|
43 | REAL alb1_glo(klon_glo,nbsrf) |
---|
44 | REAL alb2_glo(klon_glo,nbsrf) |
---|
45 | REAL rain_fall_glo(klon_glo) |
---|
46 | REAL snow_fall_glo(klon_glo) |
---|
47 | real solsw_glo(klon_glo) |
---|
48 | real sollw_glo(klon_glo) |
---|
49 | REAL radsol_glo(klon_glo) |
---|
50 | REAL zmea_glo(klon_glo) |
---|
51 | REAL zstd_glo(klon_glo) |
---|
52 | REAL zsig_glo(klon_glo) |
---|
53 | REAL zgam_glo(klon_glo) |
---|
54 | REAL zthe_glo(klon_glo) |
---|
55 | REAL zpic_glo(klon_glo) |
---|
56 | REAL zval_glo(klon_glo) |
---|
57 | REAL rugsrel_glo(klon_glo) |
---|
58 | REAL t_ancien_glo(klon_glo,klev), q_ancien_glo(klon_glo,klev) |
---|
59 | REAL clwcon_glo(klon_glo,klev) |
---|
60 | REAL rnebcon_glo(klon_glo,klev) |
---|
61 | REAL ratqs_glo(klon_glo,klev) |
---|
62 | REAL pbl_tke_glo(klon_glo,klev+1,nbsrf) |
---|
63 | REAL zmax0_glo(klon_glo), f0_glo(klon) |
---|
64 | REAL ema_work1_glo(klon_glo, klev), ema_work2_glo(klon_glo, klev) |
---|
65 | REAL wake_deltat_glo(klon_glo,klev) |
---|
66 | REAL wake_deltaq_glo(klon_glo,klev) |
---|
67 | REAL wake_s_glo(klon_glo), wake_cstar_glo(klon_glo) |
---|
68 | REAL wake_fip_glo(klon_glo) |
---|
69 | REAL tsoil_p(klon,nsoilmx,nbsrf) |
---|
70 | REAL tslab_p(klon), seaice_p(klon) |
---|
71 | REAL qsurf_p(klon,nbsrf) |
---|
72 | REAL qsol_p(klon) |
---|
73 | REAL snow_p(klon,nbsrf) |
---|
74 | REAL evap_p(klon,nbsrf) |
---|
75 | real fder_p(klon) |
---|
76 | REAL frugs_p(klon,nbsrf) |
---|
77 | REAL agesno_p(klon,nbsrf) |
---|
78 | REAL run_off_lic_0_p(klon) |
---|
79 | |
---|
80 | LOGICAL,SAVE :: ancien_ok_glo |
---|
81 | !$OMP THREADPRIVATE(ancien_ok_glo) |
---|
82 | |
---|
83 | REAL zmasq_glo(klon_glo) |
---|
84 | REAL tsoil(klon_glo,nsoilmx,nbsrf) |
---|
85 | cIM "slab" ocean |
---|
86 | REAL tslab(klon_glo), seaice(klon_glo) |
---|
87 | REAL qsurf(klon_glo,nbsrf) |
---|
88 | REAL qsol(klon_glo) |
---|
89 | REAL snow(klon_glo,nbsrf) |
---|
90 | REAL evap(klon_glo,nbsrf) |
---|
91 | real fder(klon_glo) |
---|
92 | REAL frugs(klon_glo,nbsrf) |
---|
93 | REAL agesno(klon_glo,nbsrf) |
---|
94 | REAL fractint(klon_glo) |
---|
95 | REAL run_off_lic_0(klon_glo) |
---|
96 | |
---|
97 | CHARACTER*6 ocean_in |
---|
98 | LOGICAL ok_veget_in |
---|
99 | |
---|
100 | INTEGER longcles |
---|
101 | PARAMETER ( longcles = 20 ) |
---|
102 | REAL clesphy0( longcles ) |
---|
103 | c |
---|
104 | REAL xmin, xmax |
---|
105 | c |
---|
106 | INTEGER nid, nvarid |
---|
107 | INTEGER ierr, i, nsrf, isoil ,k |
---|
108 | INTEGER length |
---|
109 | PARAMETER (length=100) |
---|
110 | REAL tab_cntrl(length), tabcntr0(length) |
---|
111 | REAL,SAVE :: tab_cntrl_omp(length) |
---|
112 | CHARACTER*7 str7 |
---|
113 | CHARACTER*2 str2 |
---|
114 | |
---|
115 | c FH1D |
---|
116 | c real iolat(jjm+1) |
---|
117 | real iolat(jjm+1-1/iim) |
---|
118 | c |
---|
119 | c Ouvrir le fichier contenant l'etat initial: |
---|
120 | c |
---|
121 | |
---|
122 | c$OMP MASTER |
---|
123 | print *,'MASTER -x , omp_rank=',omp_rank |
---|
124 | c$OMP END MASTER |
---|
125 | |
---|
126 | c$OMP MASTER |
---|
127 | IF (is_mpi_root) THEN |
---|
128 | print*,'fichnom ',fichnom |
---|
129 | ierr = NF_OPEN (fichnom, NF_NOWRITE,nid) |
---|
130 | IF (ierr.NE.NF_NOERR) THEN |
---|
131 | write(6,*)' Pb d''ouverture du fichier '//fichnom |
---|
132 | write(6,*)' ierr = ', ierr |
---|
133 | CALL ABORT |
---|
134 | ENDIF |
---|
135 | ENDIF |
---|
136 | c$OMP END MASTER |
---|
137 | c |
---|
138 | c Lecture des parametres de controle: |
---|
139 | c |
---|
140 | c$OMP MASTER |
---|
141 | IF (is_mpi_root) THEN |
---|
142 | |
---|
143 | ierr = NF_INQ_VARID (nid, "controle", nvarid) |
---|
144 | IF (ierr.NE.NF_NOERR) THEN |
---|
145 | PRINT*, 'phyetat0: Le champ <controle> est absent' |
---|
146 | CALL abort |
---|
147 | ENDIF |
---|
148 | #ifdef NC_DOUBLE |
---|
149 | ierr = NF_GET_VAR_DOUBLE(nid, nvarid, tab_cntrl) |
---|
150 | #else |
---|
151 | ierr = NF_GET_VAR_REAL(nid, nvarid, tab_cntrl) |
---|
152 | #endif |
---|
153 | IF (ierr.NE.NF_NOERR) THEN |
---|
154 | PRINT*, 'phyetat0: Lecture echouee pour <controle>' |
---|
155 | CALL abort |
---|
156 | ENDIF |
---|
157 | ENDIF |
---|
158 | |
---|
159 | c$OMP END MASTER |
---|
160 | |
---|
161 | CALL bcast(tab_cntrl) |
---|
162 | |
---|
163 | c |
---|
164 | !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! |
---|
165 | ! FH 2008/05/09 On elimine toutes les clefs physiques dans la dynamique |
---|
166 | ! Les constantes de la physiques sont lues dans la physique seulement. |
---|
167 | ! Les egalites du type |
---|
168 | ! tab_cntrl( 5 )=clesphy0(1) |
---|
169 | ! sont remplacees par |
---|
170 | ! clesphy0(1)=tab_cntrl( 5 ) |
---|
171 | ! On inverse aussi la logique. |
---|
172 | ! On remplit les tab_cntrl avec les parametres lus dans les .def |
---|
173 | !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! |
---|
174 | |
---|
175 | DO i = 1, length |
---|
176 | tabcntr0( i ) = tab_cntrl( i ) |
---|
177 | ENDDO |
---|
178 | c |
---|
179 | tab_cntrl(1)=dtime |
---|
180 | tab_cntrl(2)=radpas |
---|
181 | co2_ppm_etat0 = tab_cntrl(3) |
---|
182 | solaire_etat0 = tab_cntrl(4) |
---|
183 | tab_cntrl(5)=iflag_con |
---|
184 | tab_cntrl(6)=nbapp_rad |
---|
185 | |
---|
186 | if (cycle_diurne) tab_cntrl( 7) =1. |
---|
187 | if (soil_model) tab_cntrl( 8) =1. |
---|
188 | if (new_oliq) tab_cntrl( 9) =1. |
---|
189 | if (ok_orodr) tab_cntrl(10) =1. |
---|
190 | if (ok_orolf) tab_cntrl(11) =1. |
---|
191 | if (ok_limitvrai) tab_cntrl(12) =1. |
---|
192 | |
---|
193 | |
---|
194 | itau_phy = tab_cntrl(15) |
---|
195 | |
---|
196 | |
---|
197 | |
---|
198 | IF( clesphy0(1).NE.tab_cntrl( 5 ) ) THEN |
---|
199 | clesphy0(1)=tab_cntrl( 5 ) |
---|
200 | ENDIF |
---|
201 | |
---|
202 | IF( clesphy0(2).NE.tab_cntrl( 6 ) ) THEN |
---|
203 | clesphy0(2)=tab_cntrl( 6 ) |
---|
204 | ENDIF |
---|
205 | |
---|
206 | IF( clesphy0(3).NE.tab_cntrl( 7 ) ) THEN |
---|
207 | clesphy0(3)=tab_cntrl( 7 ) |
---|
208 | ENDIF |
---|
209 | |
---|
210 | IF( clesphy0(4).NE.tab_cntrl( 8 ) ) THEN |
---|
211 | clesphy0(4)=tab_cntrl( 8 ) |
---|
212 | ENDIF |
---|
213 | |
---|
214 | IF( clesphy0(5).NE.tab_cntrl( 9 ) ) THEN |
---|
215 | clesphy0(5)=tab_cntrl( 9 ) |
---|
216 | ENDIF |
---|
217 | |
---|
218 | IF( clesphy0(6).NE.tab_cntrl( 10 ) ) THEN |
---|
219 | clesphy0(6)=tab_cntrl( 10 ) |
---|
220 | ENDIF |
---|
221 | |
---|
222 | IF( clesphy0(7).NE.tab_cntrl( 11 ) ) THEN |
---|
223 | clesphy0(7)=tab_cntrl( 11 ) |
---|
224 | ENDIF |
---|
225 | |
---|
226 | IF( clesphy0(8).NE.tab_cntrl( 12 ) ) THEN |
---|
227 | clesphy0(8)=tab_cntrl( 12 ) |
---|
228 | ENDIF |
---|
229 | |
---|
230 | |
---|
231 | c |
---|
232 | c Lecture des latitudes (coordonnees): |
---|
233 | c |
---|
234 | c$OMP MASTER |
---|
235 | IF (is_mpi_root) THEN |
---|
236 | |
---|
237 | ierr = NF_INQ_VARID (nid, "latitude", nvarid) |
---|
238 | IF (ierr.NE.NF_NOERR) THEN |
---|
239 | PRINT*, 'phyetat0: Le champ <latitude> est absent' |
---|
240 | CALL abort |
---|
241 | ENDIF |
---|
242 | #ifdef NC_DOUBLE |
---|
243 | ierr = NF_GET_VAR_DOUBLE(nid, nvarid, rlat_glo) |
---|
244 | #else |
---|
245 | ierr = NF_GET_VAR_REAL(nid, nvarid, rlat_glo) |
---|
246 | #endif |
---|
247 | IF (ierr.NE.NF_NOERR) THEN |
---|
248 | PRINT*, 'phyetat0: Lecture echouee pour <latitude>' |
---|
249 | CALL abort |
---|
250 | ENDIF |
---|
251 | |
---|
252 | c |
---|
253 | c Lecture des longitudes (coordonnees): |
---|
254 | c |
---|
255 | ierr = NF_INQ_VARID (nid, "longitude", nvarid) |
---|
256 | IF (ierr.NE.NF_NOERR) THEN |
---|
257 | PRINT*, 'phyetat0: Le champ <longitude> est absent' |
---|
258 | CALL abort |
---|
259 | ENDIF |
---|
260 | #ifdef NC_DOUBLE |
---|
261 | ierr = NF_GET_VAR_DOUBLE(nid, nvarid, rlon_glo) |
---|
262 | #else |
---|
263 | ierr = NF_GET_VAR_REAL(nid, nvarid, rlon_glo) |
---|
264 | #endif |
---|
265 | IF (ierr.NE.NF_NOERR) THEN |
---|
266 | PRINT*, 'phyetat0: Lecture echouee pour <latitude>' |
---|
267 | CALL abort |
---|
268 | ENDIF |
---|
269 | C |
---|
270 | C |
---|
271 | C Lecture du masque terre mer |
---|
272 | C |
---|
273 | |
---|
274 | ierr = NF_INQ_VARID (nid, "masque", nvarid) |
---|
275 | IF (ierr .EQ. NF_NOERR) THEN |
---|
276 | #ifdef NC_DOUBLE |
---|
277 | ierr = NF_GET_VAR_DOUBLE(nid, nvarid, zmasq_glo) |
---|
278 | #else |
---|
279 | ierr = NF_GET_VAR_REAL(nid, nvarid, zmasq_glo) |
---|
280 | #endif |
---|
281 | IF (ierr.NE.NF_NOERR) THEN |
---|
282 | PRINT*, 'phyetat0: Lecture echouee pour <masque>' |
---|
283 | CALL abort |
---|
284 | ENDIF |
---|
285 | else |
---|
286 | PRINT*, 'phyetat0: Le champ <masque> est absent' |
---|
287 | PRINT*, 'fichier startphy non compatible avec phyetat0' |
---|
288 | C CALL abort |
---|
289 | ENDIF |
---|
290 | |
---|
291 | |
---|
292 | C Lecture des fractions pour chaque sous-surface |
---|
293 | C |
---|
294 | C initialisation des sous-surfaces |
---|
295 | C |
---|
296 | pctsrf = 0. |
---|
297 | C |
---|
298 | C fraction de terre |
---|
299 | C |
---|
300 | |
---|
301 | ierr = NF_INQ_VARID (nid, "FTER", nvarid) |
---|
302 | IF (ierr .EQ. NF_NOERR) THEN |
---|
303 | #ifdef NC_DOUBLE |
---|
304 | ierr = NF_GET_VAR_DOUBLE(nid, nvarid, |
---|
305 | . pctsrf_glo(1 : klon_glo,is_ter)) |
---|
306 | #else |
---|
307 | ierr = NF_GET_VAR_REAL(nid, nvarid, |
---|
308 | . pctsrf_glo(1 : klon_glo,is_ter)) |
---|
309 | #endif |
---|
310 | IF (ierr.NE.NF_NOERR) THEN |
---|
311 | PRINT*, 'phyetat0: Lecture echouee pour <FTER>' |
---|
312 | CALL abort |
---|
313 | ENDIF |
---|
314 | else |
---|
315 | PRINT*, 'phyetat0: Le champ <FTER> est absent' |
---|
316 | c@$$ CALL abort |
---|
317 | ENDIF |
---|
318 | |
---|
319 | C |
---|
320 | C fraction de glace de terre |
---|
321 | C |
---|
322 | ierr = NF_INQ_VARID (nid, "FLIC", nvarid) |
---|
323 | IF (ierr .EQ. NF_NOERR) THEN |
---|
324 | #ifdef NC_DOUBLE |
---|
325 | ierr = NF_GET_VAR_DOUBLE(nid, nvarid, |
---|
326 | . pctsrf_glo(1 : klon_glo,is_lic)) |
---|
327 | #else |
---|
328 | ierr = NF_GET_VAR_REAL(nid, nvarid, |
---|
329 | . pctsrf_glo(1 : klon_glo,is_lic)) |
---|
330 | #endif |
---|
331 | IF (ierr.NE.NF_NOERR) THEN |
---|
332 | PRINT*, 'phyetat0: Lecture echouee pour <FLIC>' |
---|
333 | CALL abort |
---|
334 | ENDIF |
---|
335 | else |
---|
336 | PRINT*, 'phyetat0: Le champ <FLIC> est absent' |
---|
337 | c@$$ CALL abort |
---|
338 | ENDIF |
---|
339 | C |
---|
340 | C fraction d'ocean |
---|
341 | C |
---|
342 | ierr = NF_INQ_VARID (nid, "FOCE", nvarid) |
---|
343 | IF (ierr .EQ. NF_NOERR) THEN |
---|
344 | #ifdef NC_DOUBLE |
---|
345 | ierr = NF_GET_VAR_DOUBLE(nid, nvarid, |
---|
346 | . pctsrf_glo(1 : klon_glo,is_oce)) |
---|
347 | #else |
---|
348 | ierr = NF_GET_VAR_REAL(nid, nvarid, |
---|
349 | . pctsrf_glo(1 : klon_glo,is_oce)) |
---|
350 | #endif |
---|
351 | IF (ierr.NE.NF_NOERR) THEN |
---|
352 | PRINT*, 'phyetat0: Lecture echouee pour <FOCE>' |
---|
353 | CALL abort |
---|
354 | ENDIF |
---|
355 | else |
---|
356 | PRINT*, 'phyetat0: Le champ <FOCE> est absent' |
---|
357 | c@$$ CALL abort |
---|
358 | ENDIF |
---|
359 | |
---|
360 | C |
---|
361 | C fraction glace de mer |
---|
362 | C |
---|
363 | ierr = NF_INQ_VARID (nid, "FSIC", nvarid) |
---|
364 | IF (ierr .EQ. NF_NOERR) THEN |
---|
365 | #ifdef NC_DOUBLE |
---|
366 | ierr = NF_GET_VAR_DOUBLE(nid, nvarid, |
---|
367 | . pctsrf_glo(1 : klon_glo,is_sic)) |
---|
368 | #else |
---|
369 | ierr = NF_GET_VAR_REAL(nid, nvarid, |
---|
370 | . pctsrf_glo(1 : klon_glo, is_sic)) |
---|
371 | #endif |
---|
372 | IF (ierr.NE.NF_NOERR) THEN |
---|
373 | PRINT*, 'phyetat0: Lecture echouee pour <FSIC>' |
---|
374 | CALL abort |
---|
375 | ENDIF |
---|
376 | else |
---|
377 | PRINT*, 'phyetat0: Le champ <FSIC> est absent' |
---|
378 | c@$$ CALL abort |
---|
379 | ENDIF |
---|
380 | |
---|
381 | C |
---|
382 | C Verification de l'adequation entre le masque et les sous-surfaces |
---|
383 | C |
---|
384 | fractint( 1 : klon_glo) = pctsrf_glo(1 : klon_glo, is_ter) |
---|
385 | $ + pctsrf_glo(1 : klon_glo, is_lic) |
---|
386 | DO i = 1 , klon_glo |
---|
387 | IF ( abs(fractint(i) - zmasq_glo(i) ) .GT. EPSFRA ) THEN |
---|
388 | WRITE(*,*) 'phyetat0: attention fraction terre pas ', |
---|
389 | $ 'coherente ', i, zmasq_glo(i), pctsrf_glo(i, is_ter) |
---|
390 | $ ,pctsrf_glo(i, is_lic) |
---|
391 | ENDIF |
---|
392 | END DO |
---|
393 | fractint (1 : klon_glo) = pctsrf_glo(1 : klon_glo, is_oce) |
---|
394 | $ + pctsrf_glo(1 : klon_glo, is_sic) |
---|
395 | DO i = 1 , klon_glo |
---|
396 | IF ( abs( fractint(i) - (1. - zmasq_glo(i))) .GT. EPSFRA ) THEN |
---|
397 | WRITE(*,*) 'phyetat0 attention fraction ocean pas ', |
---|
398 | $ 'coherente ', i, zmasq_glo(i) , pctsrf_glo(i, is_oce) |
---|
399 | $ ,pctsrf_glo(i, is_sic) |
---|
400 | ENDIF |
---|
401 | END DO |
---|
402 | |
---|
403 | C |
---|
404 | c Lecture des temperatures du sol: |
---|
405 | c |
---|
406 | |
---|
407 | ierr = NF_INQ_VARID (nid, "TS", nvarid) |
---|
408 | IF (ierr.NE.NF_NOERR) THEN |
---|
409 | PRINT*, 'phyetat0: Le champ <TS> est absent' |
---|
410 | PRINT*, ' Mais je vais essayer de lire TS**' |
---|
411 | DO nsrf = 1, nbsrf |
---|
412 | IF (nsrf.GT.99) THEN |
---|
413 | PRINT*, "Trop de sous-mailles" |
---|
414 | CALL abort |
---|
415 | ENDIF |
---|
416 | WRITE(str2,'(i2.2)') nsrf |
---|
417 | ierr = NF_INQ_VARID (nid, "TS"//str2, nvarid) |
---|
418 | IF (ierr.NE.NF_NOERR) THEN |
---|
419 | PRINT*, "phyetat0: Le champ <TS"//str2//"> est absent" |
---|
420 | CALL abort |
---|
421 | ENDIF |
---|
422 | #ifdef NC_DOUBLE |
---|
423 | ierr = NF_GET_VAR_DOUBLE(nid, nvarid, tsol_glo(1,nsrf)) |
---|
424 | #else |
---|
425 | ierr = NF_GET_VAR_REAL(nid, nvarid, tsol_glo(1,nsrf)) |
---|
426 | #endif |
---|
427 | IF (ierr.NE.NF_NOERR) THEN |
---|
428 | PRINT*, "phyetat0: Lecture echouee pour <TS"//str2//">" |
---|
429 | CALL abort |
---|
430 | ENDIF |
---|
431 | |
---|
432 | xmin = 1.0E+20 |
---|
433 | xmax = -1.0E+20 |
---|
434 | DO i = 1, klon_glo |
---|
435 | xmin = MIN(tsol_glo(i,nsrf),xmin) |
---|
436 | xmax = MAX(tsol_glo(i,nsrf),xmax) |
---|
437 | ENDDO |
---|
438 | PRINT*,'Temperature du sol TS**:', nsrf, xmin, xmax |
---|
439 | ENDDO |
---|
440 | ELSE |
---|
441 | PRINT*, 'phyetat0: Le champ <TS> est present' |
---|
442 | PRINT*, ' J ignore donc les autres temperatures TS**' |
---|
443 | #ifdef NC_DOUBLE |
---|
444 | ierr = NF_GET_VAR_DOUBLE(nid, nvarid, tsol_glo(1,1)) |
---|
445 | #else |
---|
446 | ierr = NF_GET_VAR_REAL(nid, nvarid, tsol_glo(1,1)) |
---|
447 | #endif |
---|
448 | IF (ierr.NE.NF_NOERR) THEN |
---|
449 | PRINT*, "phyetat0: Lecture echouee pour <TS>" |
---|
450 | CALL abort |
---|
451 | ENDIF |
---|
452 | xmin = 1.0E+20 |
---|
453 | xmax = -1.0E+20 |
---|
454 | DO i = 1, klon_glo |
---|
455 | xmin = MIN(tsol_glo(i,1),xmin) |
---|
456 | xmax = MAX(tsol_glo(i,1),xmax) |
---|
457 | ENDDO |
---|
458 | PRINT*,'Temperature du sol <TS>', xmin, xmax |
---|
459 | DO nsrf = 2, nbsrf |
---|
460 | DO i = 1, klon_glo |
---|
461 | tsol_glo(i,nsrf) = tsol_glo(i,1) |
---|
462 | ENDDO |
---|
463 | ENDDO |
---|
464 | ENDIF |
---|
465 | |
---|
466 | c |
---|
467 | c Lecture des temperatures du sol profond: |
---|
468 | c |
---|
469 | DO nsrf = 1, nbsrf |
---|
470 | DO isoil=1, nsoilmx |
---|
471 | IF (isoil.GT.99 .AND. nsrf.GT.99) THEN |
---|
472 | PRINT*, "Trop de couches ou sous-mailles" |
---|
473 | CALL abort |
---|
474 | ENDIF |
---|
475 | WRITE(str7,'(i2.2,"srf",i2.2)') isoil, nsrf |
---|
476 | ierr = NF_INQ_VARID (nid, 'Tsoil'//str7, nvarid) |
---|
477 | IF (ierr.NE.NF_NOERR) THEN |
---|
478 | PRINT*, "phyetat0: Le champ <Tsoil"//str7//"> est absent" |
---|
479 | PRINT*, " Il prend donc la valeur de surface" |
---|
480 | DO i=1, klon_glo |
---|
481 | tsoil(i,isoil,nsrf)=tsol_glo(i,nsrf) |
---|
482 | ENDDO |
---|
483 | ELSE |
---|
484 | #ifdef NC_DOUBLE |
---|
485 | ierr = NF_GET_VAR_DOUBLE(nid, nvarid, tsoil(1,isoil,nsrf)) |
---|
486 | #else |
---|
487 | ierr = NF_GET_VAR_REAL(nid, nvarid, tsoil(1,isoil,nsrf)) |
---|
488 | #endif |
---|
489 | IF (ierr.NE.NF_NOERR) THEN |
---|
490 | PRINT*, "Lecture echouee pour <Tsoil"//str7//">" |
---|
491 | CALL abort |
---|
492 | ENDIF |
---|
493 | ENDIF |
---|
494 | ENDDO |
---|
495 | ENDDO |
---|
496 | c |
---|
497 | cIM "slab" ocean |
---|
498 | c |
---|
499 | c Lecture de tslab (pour slab ocean seulement): |
---|
500 | c |
---|
501 | IF (ocean_in .eq. 'slab ') then |
---|
502 | ierr = NF_INQ_VARID (nid, "TSLAB", nvarid) |
---|
503 | IF (ierr.NE.NF_NOERR) THEN |
---|
504 | PRINT*, "phyetat0: Le champ <TSLAB> est absent" |
---|
505 | CALL abort |
---|
506 | ENDIF |
---|
507 | #ifdef NC_DOUBLE |
---|
508 | ierr = NF_GET_VAR_DOUBLE(nid, nvarid, tslab) |
---|
509 | #else |
---|
510 | ierr = NF_GET_VAR_REAL(nid, nvarid, tslab) |
---|
511 | #endif |
---|
512 | IF (ierr.NE.NF_NOERR) THEN |
---|
513 | PRINT*, "phyetat0: Lecture echouee pour <TSLAB>" |
---|
514 | CALL abort |
---|
515 | ENDIF |
---|
516 | xmin = 1.0E+20 |
---|
517 | xmax = -1.0E+20 |
---|
518 | DO i = 1, klon_glo |
---|
519 | xmin = MIN(tslab(i),xmin) |
---|
520 | xmax = MAX(tslab(i),xmax) |
---|
521 | ENDDO |
---|
522 | PRINT*,'Min, Max tslab (utilise si OCEAN=slab )', xmin, xmax |
---|
523 | c |
---|
524 | c Lecture de seaice (pour slab ocean seulement): |
---|
525 | c |
---|
526 | ierr = NF_INQ_VARID (nid, "SEAICE", nvarid) |
---|
527 | IF (ierr.NE.NF_NOERR) THEN |
---|
528 | PRINT*, "phyetat0: Le champ <SEAICE> est absent" |
---|
529 | CALL abort |
---|
530 | ENDIF |
---|
531 | #ifdef NC_DOUBLE |
---|
532 | ierr = NF_GET_VAR_DOUBLE(nid, nvarid, seaice) |
---|
533 | #else |
---|
534 | ierr = NF_GET_VAR_REAL(nid, nvarid, seaice) |
---|
535 | #endif |
---|
536 | IF (ierr.NE.NF_NOERR) THEN |
---|
537 | PRINT*, "phyetat0: Lecture echouee pour <SEAICE>" |
---|
538 | CALL abort |
---|
539 | ENDIF |
---|
540 | xmin = 1.0E+20 |
---|
541 | xmax = -1.0E+20 |
---|
542 | DO i = 1, klon_glo |
---|
543 | xmin = MIN(seaice(i),xmin) |
---|
544 | xmax = MAX(seaice(i),xmax) |
---|
545 | ENDDO |
---|
546 | PRINT*,'Masse de la glace de mer (utilise si OCEAN=slab)', |
---|
547 | $ xmin, xmax |
---|
548 | ELSE |
---|
549 | tslab = 0. |
---|
550 | seaice = 0. |
---|
551 | ENDIF |
---|
552 | c |
---|
553 | c Lecture de l'humidite de l'air juste au dessus du sol: |
---|
554 | c |
---|
555 | ierr = NF_INQ_VARID (nid, "QS", nvarid) |
---|
556 | IF (ierr.NE.NF_NOERR) THEN |
---|
557 | PRINT*, 'phyetat0: Le champ <QS> est absent' |
---|
558 | PRINT*, ' Mais je vais essayer de lire QS**' |
---|
559 | DO nsrf = 1, nbsrf |
---|
560 | IF (nsrf.GT.99) THEN |
---|
561 | PRINT*, "Trop de sous-mailles" |
---|
562 | CALL abort |
---|
563 | ENDIF |
---|
564 | WRITE(str2,'(i2.2)') nsrf |
---|
565 | ierr = NF_INQ_VARID (nid, "QS"//str2, nvarid) |
---|
566 | IF (ierr.NE.NF_NOERR) THEN |
---|
567 | PRINT*, "phyetat0: Le champ <QS"//str2//"> est absent" |
---|
568 | CALL abort |
---|
569 | ENDIF |
---|
570 | #ifdef NC_DOUBLE |
---|
571 | ierr = NF_GET_VAR_DOUBLE(nid, nvarid, qsurf(1,nsrf)) |
---|
572 | #else |
---|
573 | ierr = NF_GET_VAR_REAL(nid, nvarid, qsurf(1,nsrf)) |
---|
574 | #endif |
---|
575 | IF (ierr.NE.NF_NOERR) THEN |
---|
576 | PRINT*, "phyetat0: Lecture echouee pour <QS"//str2//">" |
---|
577 | CALL abort |
---|
578 | ENDIF |
---|
579 | xmin = 1.0E+20 |
---|
580 | xmax = -1.0E+20 |
---|
581 | DO i = 1, klon_glo |
---|
582 | xmin = MIN(qsurf(i,nsrf),xmin) |
---|
583 | xmax = MAX(qsurf(i,nsrf),xmax) |
---|
584 | ENDDO |
---|
585 | PRINT*,'Humidite pres du sol QS**:', nsrf, xmin, xmax |
---|
586 | ENDDO |
---|
587 | ELSE |
---|
588 | PRINT*, 'phyetat0: Le champ <QS> est present' |
---|
589 | PRINT*, ' J ignore donc les autres humidites QS**' |
---|
590 | #ifdef NC_DOUBLE |
---|
591 | ierr = NF_GET_VAR_DOUBLE(nid, nvarid, qsurf(1,1)) |
---|
592 | #else |
---|
593 | ierr = NF_GET_VAR_REAL(nid, nvarid, qsurf(1,1)) |
---|
594 | #endif |
---|
595 | IF (ierr.NE.NF_NOERR) THEN |
---|
596 | PRINT*, "phyetat0: Lecture echouee pour <QS>" |
---|
597 | CALL abort |
---|
598 | ENDIF |
---|
599 | xmin = 1.0E+20 |
---|
600 | xmax = -1.0E+20 |
---|
601 | DO i = 1, klon_glo |
---|
602 | xmin = MIN(qsurf(i,1),xmin) |
---|
603 | xmax = MAX(qsurf(i,1),xmax) |
---|
604 | ENDDO |
---|
605 | PRINT*,'Humidite pres du sol <QS>', xmin, xmax |
---|
606 | DO nsrf = 2, nbsrf |
---|
607 | DO i = 1, klon_glo |
---|
608 | qsurf(i,nsrf) = qsurf(i,1) |
---|
609 | ENDDO |
---|
610 | ENDDO |
---|
611 | ENDIF |
---|
612 | C |
---|
613 | C Eau dans le sol (pour le modele de sol "bucket") |
---|
614 | C |
---|
615 | ierr = NF_INQ_VARID (nid, "QSOL", nvarid) |
---|
616 | IF (ierr .EQ. NF_NOERR) THEN |
---|
617 | #ifdef NC_DOUBLE |
---|
618 | ierr = NF_GET_VAR_DOUBLE(nid, nvarid, qsol) |
---|
619 | #else |
---|
620 | ierr = NF_GET_VAR_REAL(nid, nvarid, qsol) |
---|
621 | #endif |
---|
622 | IF (ierr.NE.NF_NOERR) THEN |
---|
623 | PRINT*, 'phyetat0: Lecture echouee pour <QSOL>' |
---|
624 | CALL abort |
---|
625 | ENDIF |
---|
626 | else |
---|
627 | PRINT*, 'phyetat0: Le champ <QSOL> est absent' |
---|
628 | PRINT*, ' Valeur par defaut nulle' |
---|
629 | qsol(:)=0. |
---|
630 | c@$$ CALL abort |
---|
631 | ENDIF |
---|
632 | xmin = 1.0E+20 |
---|
633 | xmax = -1.0E+20 |
---|
634 | DO i = 1, klon_glo |
---|
635 | xmin = MIN(qsol(i),xmin) |
---|
636 | xmax = MAX(qsol(i),xmax) |
---|
637 | ENDDO |
---|
638 | PRINT*,'Eau dans le sol (mm) <QSOL>', xmin, xmax |
---|
639 | c |
---|
640 | c Lecture de neige au sol: |
---|
641 | c |
---|
642 | ierr = NF_INQ_VARID (nid, "SNOW", nvarid) |
---|
643 | IF (ierr.NE.NF_NOERR) THEN |
---|
644 | PRINT*, 'phyetat0: Le champ <SNOW> est absent' |
---|
645 | PRINT*, ' Mais je vais essayer de lire SNOW**' |
---|
646 | DO nsrf = 1, nbsrf |
---|
647 | IF (nsrf.GT.99) THEN |
---|
648 | PRINT*, "Trop de sous-mailles" |
---|
649 | CALL abort |
---|
650 | ENDIF |
---|
651 | WRITE(str2,'(i2.2)') nsrf |
---|
652 | ierr = NF_INQ_VARID (nid, "SNOW"//str2, nvarid) |
---|
653 | IF (ierr.NE.NF_NOERR) THEN |
---|
654 | PRINT*, "phyetat0: Le champ <SNOW"//str2//"> est absent" |
---|
655 | CALL abort |
---|
656 | ENDIF |
---|
657 | #ifdef NC_DOUBLE |
---|
658 | ierr = NF_GET_VAR_DOUBLE(nid, nvarid, snow(1,nsrf)) |
---|
659 | #else |
---|
660 | ierr = NF_GET_VAR_REAL(nid, nvarid, snow(1,nsrf)) |
---|
661 | #endif |
---|
662 | IF (ierr.NE.NF_NOERR) THEN |
---|
663 | PRINT*, "phyetat0: Lecture echouee pour <SNOW"//str2//">" |
---|
664 | CALL abort |
---|
665 | ENDIF |
---|
666 | xmin = 1.0E+20 |
---|
667 | xmax = -1.0E+20 |
---|
668 | DO i = 1, klon_glo |
---|
669 | xmin = MIN(snow(i,nsrf),xmin) |
---|
670 | xmax = MAX(snow(i,nsrf),xmax) |
---|
671 | ENDDO |
---|
672 | PRINT*,'Neige du sol SNOW**:', nsrf, xmin, xmax |
---|
673 | ENDDO |
---|
674 | ELSE |
---|
675 | PRINT*, 'phyetat0: Le champ <SNOW> est present' |
---|
676 | PRINT*, ' J ignore donc les autres neiges SNOW**' |
---|
677 | #ifdef NC_DOUBLE |
---|
678 | ierr = NF_GET_VAR_DOUBLE(nid, nvarid, snow(1,1)) |
---|
679 | #else |
---|
680 | ierr = NF_GET_VAR_REAL(nid, nvarid, snow(1,1)) |
---|
681 | #endif |
---|
682 | IF (ierr.NE.NF_NOERR) THEN |
---|
683 | PRINT*, "phyetat0: Lecture echouee pour <SNOW>" |
---|
684 | CALL abort |
---|
685 | ENDIF |
---|
686 | xmin = 1.0E+20 |
---|
687 | xmax = -1.0E+20 |
---|
688 | DO i = 1, klon_glo |
---|
689 | xmin = MIN(snow(i,1),xmin) |
---|
690 | xmax = MAX(snow(i,1),xmax) |
---|
691 | ENDDO |
---|
692 | PRINT*,'Neige du sol <SNOW>', xmin, xmax |
---|
693 | DO nsrf = 2, nbsrf |
---|
694 | DO i = 1, klon_glo |
---|
695 | snow(i,nsrf) = snow(i,1) |
---|
696 | ENDDO |
---|
697 | ENDDO |
---|
698 | ENDIF |
---|
699 | c |
---|
700 | c Lecture de albedo de l'interval visible au sol: |
---|
701 | c |
---|
702 | ierr = NF_INQ_VARID (nid, "ALBE", nvarid) |
---|
703 | IF (ierr.NE.NF_NOERR) THEN |
---|
704 | PRINT*, 'phyetat0: Le champ <ALBE> est absent' |
---|
705 | PRINT*, ' Mais je vais essayer de lire ALBE**' |
---|
706 | DO nsrf = 1, nbsrf |
---|
707 | IF (nsrf.GT.99) THEN |
---|
708 | PRINT*, "Trop de sous-mailles" |
---|
709 | CALL abort |
---|
710 | ENDIF |
---|
711 | WRITE(str2,'(i2.2)') nsrf |
---|
712 | ierr = NF_INQ_VARID (nid, "ALBE"//str2, nvarid) |
---|
713 | IF (ierr.NE.NF_NOERR) THEN |
---|
714 | PRINT*, "phyetat0: Le champ <ALBE"//str2//"> est absent" |
---|
715 | CALL abort |
---|
716 | ENDIF |
---|
717 | #ifdef NC_DOUBLE |
---|
718 | ierr = NF_GET_VAR_DOUBLE(nid, nvarid, alb1_glo(1,nsrf)) |
---|
719 | #else |
---|
720 | ierr = NF_GET_VAR_REAL(nid, nvarid, alb1_glo(1,nsrf)) |
---|
721 | #endif |
---|
722 | IF (ierr.NE.NF_NOERR) THEN |
---|
723 | PRINT*, "phyetat0: Lecture echouee pour <ALBE"//str2//">" |
---|
724 | CALL abort |
---|
725 | ENDIF |
---|
726 | xmin = 1.0E+20 |
---|
727 | xmax = -1.0E+20 |
---|
728 | DO i = 1, klon_glo |
---|
729 | xmin = MIN(alb1_glo(i,nsrf),xmin) |
---|
730 | xmax = MAX(alb1_glo(i,nsrf),xmax) |
---|
731 | ENDDO |
---|
732 | PRINT*,'Albedo du sol ALBE**:', nsrf, xmin, xmax |
---|
733 | ENDDO |
---|
734 | ELSE |
---|
735 | PRINT*, 'phyetat0: Le champ <ALBE> est present' |
---|
736 | PRINT*, ' J ignore donc les autres ALBE**' |
---|
737 | #ifdef NC_DOUBLE |
---|
738 | ierr = NF_GET_VAR_DOUBLE(nid, nvarid, alb1_glo(1,1)) |
---|
739 | #else |
---|
740 | ierr = NF_GET_VAR_REAL(nid, nvarid, alb1_glo(1,1)) |
---|
741 | #endif |
---|
742 | IF (ierr.NE.NF_NOERR) THEN |
---|
743 | PRINT*, "phyetat0: Lecture echouee pour <ALBE>" |
---|
744 | CALL abort |
---|
745 | ENDIF |
---|
746 | xmin = 1.0E+20 |
---|
747 | xmax = -1.0E+20 |
---|
748 | DO i = 1, klon_glo |
---|
749 | xmin = MIN(alb1_glo(i,1),xmin) |
---|
750 | xmax = MAX(alb1_glo(i,1),xmax) |
---|
751 | ENDDO |
---|
752 | PRINT*,'Neige du sol <ALBE>', xmin, xmax |
---|
753 | DO nsrf = 2, nbsrf |
---|
754 | DO i = 1, klon_glo |
---|
755 | alb1_glo(i,nsrf) = alb1_glo(i,1) |
---|
756 | ENDDO |
---|
757 | ENDDO |
---|
758 | ENDIF |
---|
759 | |
---|
760 | c |
---|
761 | c Lecture de albedo au sol dans l'interval proche infra-rouge: |
---|
762 | c |
---|
763 | ierr = NF_INQ_VARID (nid, "ALBLW", nvarid) |
---|
764 | IF (ierr.NE.NF_NOERR) THEN |
---|
765 | PRINT*, 'phyetat0: Le champ <ALBLW> est absent' |
---|
766 | c PRINT*, ' Mais je vais essayer de lire ALBLW**' |
---|
767 | PRINT*, ' Mais je vais prendre ALBE**' |
---|
768 | DO nsrf = 1, nbsrf |
---|
769 | DO i = 1, klon_glo |
---|
770 | alb2_glo(i,nsrf) = alb1_glo(i,nsrf) |
---|
771 | ENDDO |
---|
772 | ENDDO |
---|
773 | ELSE |
---|
774 | PRINT*, 'phyetat0: Le champ <ALBLW> est present' |
---|
775 | PRINT*, ' J ignore donc les autres ALBLW**' |
---|
776 | #ifdef NC_DOUBLE |
---|
777 | ierr = NF_GET_VAR_DOUBLE(nid, nvarid, alb2_glo(1,1)) |
---|
778 | #else |
---|
779 | ierr = NF_GET_VAR_REAL(nid, nvarid, alb2_glo(1,1)) |
---|
780 | #endif |
---|
781 | IF (ierr.NE.NF_NOERR) THEN |
---|
782 | PRINT*, "phyetat0: Lecture echouee pour <ALBLW>" |
---|
783 | CALL abort |
---|
784 | ENDIF |
---|
785 | xmin = 1.0E+20 |
---|
786 | xmax = -1.0E+20 |
---|
787 | DO i = 1, klon_glo |
---|
788 | xmin = MIN(alb2_glo(i,1),xmin) |
---|
789 | xmax = MAX(alb2_glo(i,1),xmax) |
---|
790 | ENDDO |
---|
791 | PRINT*,'Neige du sol <ALBLW>', xmin, xmax |
---|
792 | DO nsrf = 2, nbsrf |
---|
793 | DO i = 1, klon_glo |
---|
794 | alb2_glo(i,nsrf) = alb2_glo(i,1) |
---|
795 | ENDDO |
---|
796 | ENDDO |
---|
797 | ENDIF |
---|
798 | c |
---|
799 | c Lecture de evaporation: |
---|
800 | c |
---|
801 | ierr = NF_INQ_VARID (nid, "EVAP", nvarid) |
---|
802 | IF (ierr.NE.NF_NOERR) THEN |
---|
803 | PRINT*, 'phyetat0: Le champ <EVAP> est absent' |
---|
804 | PRINT*, ' Mais je vais essayer de lire EVAP**' |
---|
805 | DO nsrf = 1, nbsrf |
---|
806 | IF (nsrf.GT.99) THEN |
---|
807 | PRINT*, "Trop de sous-mailles" |
---|
808 | CALL abort |
---|
809 | ENDIF |
---|
810 | WRITE(str2,'(i2.2)') nsrf |
---|
811 | ierr = NF_INQ_VARID (nid, "EVAP"//str2, nvarid) |
---|
812 | IF (ierr.NE.NF_NOERR) THEN |
---|
813 | PRINT*, "phyetat0: Le champ <EVAP"//str2//"> est absent" |
---|
814 | CALL abort |
---|
815 | ENDIF |
---|
816 | #ifdef NC_DOUBLE |
---|
817 | ierr = NF_GET_VAR_DOUBLE(nid, nvarid, evap(1,nsrf)) |
---|
818 | #else |
---|
819 | ierr = NF_GET_VAR_REAL(nid, nvarid, evap(1,nsrf)) |
---|
820 | #endif |
---|
821 | IF (ierr.NE.NF_NOERR) THEN |
---|
822 | PRINT*, "phyetat0: Lecture echouee pour <EVAP"//str2//">" |
---|
823 | CALL abort |
---|
824 | ENDIF |
---|
825 | xmin = 1.0E+20 |
---|
826 | xmax = -1.0E+20 |
---|
827 | DO i = 1, klon_glo |
---|
828 | xmin = MIN(evap(i,nsrf),xmin) |
---|
829 | xmax = MAX(evap(i,nsrf),xmax) |
---|
830 | ENDDO |
---|
831 | PRINT*,'evap du sol EVAP**:', nsrf, xmin, xmax |
---|
832 | ENDDO |
---|
833 | ELSE |
---|
834 | PRINT*, 'phyetat0: Le champ <EVAP> est present' |
---|
835 | PRINT*, ' J ignore donc les autres EVAP**' |
---|
836 | #ifdef NC_DOUBLE |
---|
837 | ierr = NF_GET_VAR_DOUBLE(nid, nvarid, evap(1,1)) |
---|
838 | #else |
---|
839 | ierr = NF_GET_VAR_REAL(nid, nvarid, evap(1,1)) |
---|
840 | #endif |
---|
841 | IF (ierr.NE.NF_NOERR) THEN |
---|
842 | PRINT*, "phyetat0: Lecture echouee pour <EVAP>" |
---|
843 | CALL abort |
---|
844 | ENDIF |
---|
845 | xmin = 1.0E+20 |
---|
846 | xmax = -1.0E+20 |
---|
847 | DO i = 1, klon_glo |
---|
848 | xmin = MIN(evap(i,1),xmin) |
---|
849 | xmax = MAX(evap(i,1),xmax) |
---|
850 | ENDDO |
---|
851 | PRINT*,'Evap du sol <EVAP>', xmin, xmax |
---|
852 | DO nsrf = 2, nbsrf |
---|
853 | DO i = 1, klon_glo |
---|
854 | evap(i,nsrf) = evap(i,1) |
---|
855 | ENDDO |
---|
856 | ENDDO |
---|
857 | ENDIF |
---|
858 | c |
---|
859 | c Lecture precipitation liquide: |
---|
860 | c |
---|
861 | ierr = NF_INQ_VARID (nid, "rain_f", nvarid) |
---|
862 | IF (ierr.NE.NF_NOERR) THEN |
---|
863 | PRINT*, 'phyetat0: Le champ <rain_f> est absent' |
---|
864 | CALL abort |
---|
865 | ENDIF |
---|
866 | #ifdef NC_DOUBLE |
---|
867 | ierr = NF_GET_VAR_DOUBLE(nid, nvarid, rain_fall_glo) |
---|
868 | #else |
---|
869 | ierr = NF_GET_VAR_REAL(nid, nvarid, rain_fall_glo) |
---|
870 | #endif |
---|
871 | IF (ierr.NE.NF_NOERR) THEN |
---|
872 | PRINT*, 'phyetat0: Lecture echouee pour <rain_f>' |
---|
873 | CALL abort |
---|
874 | ENDIF |
---|
875 | xmin = 1.0E+20 |
---|
876 | xmax = -1.0E+20 |
---|
877 | DO i = 1, klon_glo |
---|
878 | xmin = MIN(rain_fall_glo(i),xmin) |
---|
879 | xmax = MAX(rain_fall_glo(i),xmax) |
---|
880 | ENDDO |
---|
881 | PRINT*,'Precipitation liquide rain_f:', xmin, xmax |
---|
882 | c |
---|
883 | c Lecture precipitation solide: |
---|
884 | c |
---|
885 | ierr = NF_INQ_VARID (nid, "snow_f", nvarid) |
---|
886 | IF (ierr.NE.NF_NOERR) THEN |
---|
887 | PRINT*, 'phyetat0: Le champ <snow_f> est absent' |
---|
888 | CALL abort |
---|
889 | ENDIF |
---|
890 | #ifdef NC_DOUBLE |
---|
891 | ierr = NF_GET_VAR_DOUBLE(nid, nvarid, snow_fall_glo) |
---|
892 | #else |
---|
893 | ierr = NF_GET_VAR_REAL(nid, nvarid, snow_fall_glo) |
---|
894 | #endif |
---|
895 | IF (ierr.NE.NF_NOERR) THEN |
---|
896 | PRINT*, 'phyetat0: Lecture echouee pour <snow_f>' |
---|
897 | CALL abort |
---|
898 | ENDIF |
---|
899 | xmin = 1.0E+20 |
---|
900 | xmax = -1.0E+20 |
---|
901 | DO i = 1, klon_glo |
---|
902 | xmin = MIN(snow_fall_glo(i),xmin) |
---|
903 | xmax = MAX(snow_fall_glo(i),xmax) |
---|
904 | ENDDO |
---|
905 | PRINT*,'Precipitation solide snow_f:', xmin, xmax |
---|
906 | c |
---|
907 | c Lecture rayonnement solaire au sol: |
---|
908 | c |
---|
909 | ierr = NF_INQ_VARID (nid, "solsw", nvarid) |
---|
910 | IF (ierr.NE.NF_NOERR) THEN |
---|
911 | PRINT*, 'phyetat0: Le champ <solsw> est absent' |
---|
912 | PRINT*, 'mis a zero' |
---|
913 | solsw_glo = 0. |
---|
914 | ELSE |
---|
915 | #ifdef NC_DOUBLE |
---|
916 | ierr = NF_GET_VAR_DOUBLE(nid, nvarid, solsw_glo) |
---|
917 | #else |
---|
918 | ierr = NF_GET_VAR_REAL(nid, nvarid, solsw_glo) |
---|
919 | #endif |
---|
920 | IF (ierr.NE.NF_NOERR) THEN |
---|
921 | PRINT*, 'phyetat0: Lecture echouee pour <solsw>' |
---|
922 | CALL abort |
---|
923 | ENDIF |
---|
924 | ENDIF |
---|
925 | xmin = 1.0E+20 |
---|
926 | xmax = -1.0E+20 |
---|
927 | DO i = 1, klon_glo |
---|
928 | xmin = MIN(solsw_glo(i),xmin) |
---|
929 | xmax = MAX(solsw_glo(i),xmax) |
---|
930 | ENDDO |
---|
931 | PRINT*,'Rayonnement solaire au sol solsw:', xmin, xmax |
---|
932 | c |
---|
933 | c Lecture rayonnement IF au sol: |
---|
934 | c |
---|
935 | ierr = NF_INQ_VARID (nid, "sollw", nvarid) |
---|
936 | IF (ierr.NE.NF_NOERR) THEN |
---|
937 | PRINT*, 'phyetat0: Le champ <sollw> est absent' |
---|
938 | PRINT*, 'mis a zero' |
---|
939 | sollw_glo = 0. |
---|
940 | ELSE |
---|
941 | #ifdef NC_DOUBLE |
---|
942 | ierr = NF_GET_VAR_DOUBLE(nid, nvarid, sollw_glo) |
---|
943 | #else |
---|
944 | ierr = NF_GET_VAR_REAL(nid, nvarid, sollw_glo) |
---|
945 | #endif |
---|
946 | IF (ierr.NE.NF_NOERR) THEN |
---|
947 | PRINT*, 'phyetat0: Lecture echouee pour <sollw>' |
---|
948 | CALL abort |
---|
949 | ENDIF |
---|
950 | ENDIF |
---|
951 | xmin = 1.0E+20 |
---|
952 | xmax = -1.0E+20 |
---|
953 | DO i = 1, klon_glo |
---|
954 | xmin = MIN(sollw_glo(i),xmin) |
---|
955 | xmax = MAX(sollw_glo(i),xmax) |
---|
956 | ENDDO |
---|
957 | PRINT*,'Rayonnement IF au sol sollw:', xmin, xmax |
---|
958 | |
---|
959 | ENDIF ! is_mpi_root |
---|
960 | c$OMP END MASTER |
---|
961 | |
---|
962 | |
---|
963 | c$OMP MASTER |
---|
964 | IF (is_mpi_root) THEN |
---|
965 | c |
---|
966 | c Lecture derive des flux: |
---|
967 | c |
---|
968 | ierr = NF_INQ_VARID (nid, "fder", nvarid) |
---|
969 | IF (ierr.NE.NF_NOERR) THEN |
---|
970 | PRINT*, 'phyetat0: Le champ <fder> est absent' |
---|
971 | PRINT*, 'mis a zero' |
---|
972 | fder = 0. |
---|
973 | ELSE |
---|
974 | #ifdef NC_DOUBLE |
---|
975 | ierr = NF_GET_VAR_DOUBLE(nid, nvarid, fder) |
---|
976 | #else |
---|
977 | ierr = NF_GET_VAR_REAL(nid, nvarid, fder) |
---|
978 | #endif |
---|
979 | IF (ierr.NE.NF_NOERR) THEN |
---|
980 | PRINT*, 'phyetat0: Lecture echouee pour <fder>' |
---|
981 | CALL abort |
---|
982 | ENDIF |
---|
983 | ENDIF |
---|
984 | xmin = 1.0E+20 |
---|
985 | xmax = -1.0E+20 |
---|
986 | DO i = 1, klon_glo |
---|
987 | xmin = MIN(fder(i),xmin) |
---|
988 | xmax = MAX(fder(i),xmax) |
---|
989 | ENDDO |
---|
990 | PRINT*,'Derive des flux fder:', xmin, xmax |
---|
991 | |
---|
992 | c |
---|
993 | c Lecture du rayonnement net au sol: |
---|
994 | c |
---|
995 | ierr = NF_INQ_VARID (nid, "RADS", nvarid) |
---|
996 | IF (ierr.NE.NF_NOERR) THEN |
---|
997 | PRINT*, 'phyetat0: Le champ <RADS> est absent' |
---|
998 | CALL abort |
---|
999 | ENDIF |
---|
1000 | #ifdef NC_DOUBLE |
---|
1001 | ierr = NF_GET_VAR_DOUBLE(nid, nvarid, radsol_glo) |
---|
1002 | #else |
---|
1003 | ierr = NF_GET_VAR_REAL(nid, nvarid, radsol_glo) |
---|
1004 | #endif |
---|
1005 | IF (ierr.NE.NF_NOERR) THEN |
---|
1006 | PRINT*, 'phyetat0: Lecture echouee pour <RADS>' |
---|
1007 | CALL abort |
---|
1008 | ENDIF |
---|
1009 | xmin = 1.0E+20 |
---|
1010 | xmax = -1.0E+20 |
---|
1011 | DO i = 1, klon_glo |
---|
1012 | xmin = MIN(radsol_glo(i),xmin) |
---|
1013 | xmax = MAX(radsol_glo(i),xmax) |
---|
1014 | ENDDO |
---|
1015 | PRINT*,'Rayonnement net au sol radsol:', xmin, xmax |
---|
1016 | c |
---|
1017 | c Lecture de la longueur de rugosite |
---|
1018 | c |
---|
1019 | c |
---|
1020 | ierr = NF_INQ_VARID (nid, "RUG", nvarid) |
---|
1021 | IF (ierr.NE.NF_NOERR) THEN |
---|
1022 | PRINT*, 'phyetat0: Le champ <RUG> est absent' |
---|
1023 | PRINT*, ' Mais je vais essayer de lire RUG**' |
---|
1024 | DO nsrf = 1, nbsrf |
---|
1025 | IF (nsrf.GT.99) THEN |
---|
1026 | PRINT*, "Trop de sous-mailles" |
---|
1027 | CALL abort |
---|
1028 | ENDIF |
---|
1029 | WRITE(str2,'(i2.2)') nsrf |
---|
1030 | ierr = NF_INQ_VARID (nid, "RUG"//str2, nvarid) |
---|
1031 | IF (ierr.NE.NF_NOERR) THEN |
---|
1032 | PRINT*, "phyetat0: Le champ <RUG"//str2//"> est absent" |
---|
1033 | CALL abort |
---|
1034 | ENDIF |
---|
1035 | #ifdef NC_DOUBLE |
---|
1036 | ierr = NF_GET_VAR_DOUBLE(nid, nvarid, frugs(1,nsrf)) |
---|
1037 | #else |
---|
1038 | ierr = NF_GET_VAR_REAL(nid, nvarid, frugs(1,nsrf)) |
---|
1039 | #endif |
---|
1040 | IF (ierr.NE.NF_NOERR) THEN |
---|
1041 | PRINT*, "phyetat0: Lecture echouee pour <RUG"//str2//">" |
---|
1042 | CALL abort |
---|
1043 | ENDIF |
---|
1044 | xmin = 1.0E+20 |
---|
1045 | xmax = -1.0E+20 |
---|
1046 | DO i = 1, klon_glo |
---|
1047 | xmin = MIN(frugs(i,nsrf),xmin) |
---|
1048 | xmax = MAX(frugs(i,nsrf),xmax) |
---|
1049 | ENDDO |
---|
1050 | PRINT*,'rugosite du sol RUG**:', nsrf, xmin, xmax |
---|
1051 | ENDDO |
---|
1052 | ELSE |
---|
1053 | PRINT*, 'phyetat0: Le champ <RUG> est present' |
---|
1054 | PRINT*, ' J ignore donc les autres RUG**' |
---|
1055 | #ifdef NC_DOUBLE |
---|
1056 | ierr = NF_GET_VAR_DOUBLE(nid, nvarid, frugs(1,1)) |
---|
1057 | #else |
---|
1058 | ierr = NF_GET_VAR_REAL(nid, nvarid, frugs(1,1)) |
---|
1059 | #endif |
---|
1060 | IF (ierr.NE.NF_NOERR) THEN |
---|
1061 | PRINT*, "phyetat0: Lecture echouee pour <RUG>" |
---|
1062 | CALL abort |
---|
1063 | ENDIF |
---|
1064 | xmin = 1.0E+20 |
---|
1065 | xmax = -1.0E+20 |
---|
1066 | DO i = 1, klon_glo |
---|
1067 | xmin = MIN(frugs(i,1),xmin) |
---|
1068 | xmax = MAX(frugs(i,1),xmax) |
---|
1069 | ENDDO |
---|
1070 | PRINT*,'rugosite <RUG>', xmin, xmax |
---|
1071 | DO nsrf = 2, nbsrf |
---|
1072 | DO i = 1, klon_glo |
---|
1073 | frugs(i,nsrf) = frugs(i,1) |
---|
1074 | ENDDO |
---|
1075 | ENDDO |
---|
1076 | ENDIF |
---|
1077 | |
---|
1078 | c |
---|
1079 | c Lecture de l'age de la neige: |
---|
1080 | c |
---|
1081 | ierr = NF_INQ_VARID (nid, "AGESNO", nvarid) |
---|
1082 | IF (ierr.NE.NF_NOERR) THEN |
---|
1083 | PRINT*, 'phyetat0: Le champ <AGESNO> est absent' |
---|
1084 | PRINT*, ' Mais je vais essayer de lire AGESNO**' |
---|
1085 | DO nsrf = 1, nbsrf |
---|
1086 | IF (nsrf.GT.99) THEN |
---|
1087 | PRINT*, "Trop de sous-mailles" |
---|
1088 | CALL abort |
---|
1089 | ENDIF |
---|
1090 | WRITE(str2,'(i2.2)') nsrf |
---|
1091 | ierr = NF_INQ_VARID (nid, "AGESNO"//str2, nvarid) |
---|
1092 | IF (ierr.NE.NF_NOERR) THEN |
---|
1093 | PRINT*, "phyetat0: Le champ <AGESNO"//str2//"> est absent" |
---|
1094 | agesno = 50.0 |
---|
1095 | ENDIF |
---|
1096 | #ifdef NC_DOUBLE |
---|
1097 | ierr = NF_GET_VAR_DOUBLE(nid, nvarid, agesno(1,nsrf)) |
---|
1098 | #else |
---|
1099 | ierr = NF_GET_VAR_REAL(nid, nvarid, agesno(1,nsrf)) |
---|
1100 | #endif |
---|
1101 | IF (ierr.NE.NF_NOERR) THEN |
---|
1102 | PRINT*, "phyetat0: Lecture echouee pour <AGESNO"//str2//">" |
---|
1103 | CALL abort |
---|
1104 | ENDIF |
---|
1105 | xmin = 1.0E+20 |
---|
1106 | xmax = -1.0E+20 |
---|
1107 | DO i = 1, klon_glo |
---|
1108 | xmin = MIN(agesno(i,nsrf),xmin) |
---|
1109 | xmax = MAX(agesno(i,nsrf),xmax) |
---|
1110 | ENDDO |
---|
1111 | PRINT*,'Age de la neige AGESNO**:', nsrf, xmin, xmax |
---|
1112 | ENDDO |
---|
1113 | ELSE |
---|
1114 | PRINT*, 'phyetat0: Le champ <AGESNO> est present' |
---|
1115 | PRINT*, ' J ignore donc les autres AGESNO**' |
---|
1116 | #ifdef NC_DOUBLE |
---|
1117 | ierr = NF_GET_VAR_DOUBLE(nid, nvarid, agesno(1,1)) |
---|
1118 | #else |
---|
1119 | ierr = NF_GET_VAR_REAL(nid, nvarid, agesno(1,1)) |
---|
1120 | #endif |
---|
1121 | IF (ierr.NE.NF_NOERR) THEN |
---|
1122 | PRINT*, "phyetat0: Lecture echouee pour <AGESNO>" |
---|
1123 | CALL abort |
---|
1124 | ENDIF |
---|
1125 | xmin = 1.0E+20 |
---|
1126 | xmax = -1.0E+20 |
---|
1127 | DO i = 1, klon_glo |
---|
1128 | xmin = MIN(agesno(i,1),xmin) |
---|
1129 | xmax = MAX(agesno(i,1),xmax) |
---|
1130 | ENDDO |
---|
1131 | PRINT*,'Age de la neige <AGESNO>', xmin, xmax |
---|
1132 | DO nsrf = 2, nbsrf |
---|
1133 | DO i = 1, klon_glo |
---|
1134 | agesno(i,nsrf) = agesno(i,1) |
---|
1135 | ENDDO |
---|
1136 | ENDDO |
---|
1137 | ENDIF |
---|
1138 | |
---|
1139 | c |
---|
1140 | ierr = NF_INQ_VARID (nid, "ZMEA", nvarid) |
---|
1141 | IF (ierr.NE.NF_NOERR) THEN |
---|
1142 | PRINT*, 'phyetat0: Le champ <ZMEA> est absent' |
---|
1143 | CALL abort |
---|
1144 | ENDIF |
---|
1145 | #ifdef NC_DOUBLE |
---|
1146 | ierr = NF_GET_VAR_DOUBLE(nid, nvarid, zmea_glo) |
---|
1147 | #else |
---|
1148 | ierr = NF_GET_VAR_REAL(nid, nvarid, zmea_glo) |
---|
1149 | #endif |
---|
1150 | IF (ierr.NE.NF_NOERR) THEN |
---|
1151 | PRINT*, 'phyetat0: Lecture echouee pour <ZMEA>' |
---|
1152 | CALL abort |
---|
1153 | ENDIF |
---|
1154 | xmin = 1.0E+20 |
---|
1155 | xmax = -1.0E+20 |
---|
1156 | DO i = 1, klon_glo |
---|
1157 | xmin = MIN(zmea_glo(i),xmin) |
---|
1158 | xmax = MAX(zmea_glo(i),xmax) |
---|
1159 | ENDDO |
---|
1160 | PRINT*,'OROGRAPHIE SOUS-MAILLE zmea:', xmin, xmax |
---|
1161 | c |
---|
1162 | c |
---|
1163 | ierr = NF_INQ_VARID (nid, "ZSTD", nvarid) |
---|
1164 | IF (ierr.NE.NF_NOERR) THEN |
---|
1165 | PRINT*, 'phyetat0: Le champ <ZSTD> est absent' |
---|
1166 | CALL abort |
---|
1167 | ENDIF |
---|
1168 | #ifdef NC_DOUBLE |
---|
1169 | ierr = NF_GET_VAR_DOUBLE(nid, nvarid, zstd_glo) |
---|
1170 | #else |
---|
1171 | ierr = NF_GET_VAR_REAL(nid, nvarid, zstd_glo) |
---|
1172 | #endif |
---|
1173 | IF (ierr.NE.NF_NOERR) THEN |
---|
1174 | PRINT*, 'phyetat0: Lecture echouee pour <ZSTD>' |
---|
1175 | CALL abort |
---|
1176 | ENDIF |
---|
1177 | xmin = 1.0E+20 |
---|
1178 | xmax = -1.0E+20 |
---|
1179 | DO i = 1, klon_glo |
---|
1180 | xmin = MIN(zstd_glo(i),xmin) |
---|
1181 | xmax = MAX(zstd_glo(i),xmax) |
---|
1182 | ENDDO |
---|
1183 | PRINT*,'OROGRAPHIE SOUS-MAILLE zstd:', xmin, xmax |
---|
1184 | c |
---|
1185 | c |
---|
1186 | ierr = NF_INQ_VARID (nid, "ZSIG", nvarid) |
---|
1187 | IF (ierr.NE.NF_NOERR) THEN |
---|
1188 | PRINT*, 'phyetat0: Le champ <ZSIG> est absent' |
---|
1189 | CALL abort |
---|
1190 | ENDIF |
---|
1191 | #ifdef NC_DOUBLE |
---|
1192 | ierr = NF_GET_VAR_DOUBLE(nid, nvarid, zsig_glo) |
---|
1193 | #else |
---|
1194 | ierr = NF_GET_VAR_REAL(nid, nvarid, zsig_glo) |
---|
1195 | #endif |
---|
1196 | IF (ierr.NE.NF_NOERR) THEN |
---|
1197 | PRINT*, 'phyetat0: Lecture echouee pour <ZSIG>' |
---|
1198 | CALL abort |
---|
1199 | ENDIF |
---|
1200 | xmin = 1.0E+20 |
---|
1201 | xmax = -1.0E+20 |
---|
1202 | DO i = 1, klon_glo |
---|
1203 | xmin = MIN(zsig_glo(i),xmin) |
---|
1204 | xmax = MAX(zsig_glo(i),xmax) |
---|
1205 | ENDDO |
---|
1206 | PRINT*,'OROGRAPHIE SOUS-MAILLE zsig:', xmin, xmax |
---|
1207 | c |
---|
1208 | c |
---|
1209 | ierr = NF_INQ_VARID (nid, "ZGAM", nvarid) |
---|
1210 | IF (ierr.NE.NF_NOERR) THEN |
---|
1211 | PRINT*, 'phyetat0: Le champ <ZGAM> est absent' |
---|
1212 | CALL abort |
---|
1213 | ENDIF |
---|
1214 | #ifdef NC_DOUBLE |
---|
1215 | ierr = NF_GET_VAR_DOUBLE(nid, nvarid, zgam_glo) |
---|
1216 | #else |
---|
1217 | ierr = NF_GET_VAR_REAL(nid, nvarid, zgam_glo) |
---|
1218 | #endif |
---|
1219 | IF (ierr.NE.NF_NOERR) THEN |
---|
1220 | PRINT*, 'phyetat0: Lecture echouee pour <ZGAM>' |
---|
1221 | CALL abort |
---|
1222 | ENDIF |
---|
1223 | xmin = 1.0E+20 |
---|
1224 | xmax = -1.0E+20 |
---|
1225 | DO i = 1, klon_glo |
---|
1226 | xmin = MIN(zgam_glo(i),xmin) |
---|
1227 | xmax = MAX(zgam_glo(i),xmax) |
---|
1228 | ENDDO |
---|
1229 | PRINT*,'OROGRAPHIE SOUS-MAILLE zgam:', xmin, xmax |
---|
1230 | c |
---|
1231 | c |
---|
1232 | ierr = NF_INQ_VARID (nid, "ZTHE", nvarid) |
---|
1233 | IF (ierr.NE.NF_NOERR) THEN |
---|
1234 | PRINT*, 'phyetat0: Le champ <ZTHE> est absent' |
---|
1235 | CALL abort |
---|
1236 | ENDIF |
---|
1237 | #ifdef NC_DOUBLE |
---|
1238 | ierr = NF_GET_VAR_DOUBLE(nid, nvarid, zthe_glo) |
---|
1239 | #else |
---|
1240 | ierr = NF_GET_VAR_REAL(nid, nvarid, zthe_glo) |
---|
1241 | #endif |
---|
1242 | IF (ierr.NE.NF_NOERR) THEN |
---|
1243 | PRINT*, 'phyetat0: Lecture echouee pour <ZTHE>' |
---|
1244 | CALL abort |
---|
1245 | ENDIF |
---|
1246 | xmin = 1.0E+20 |
---|
1247 | xmax = -1.0E+20 |
---|
1248 | DO i = 1, klon_glo |
---|
1249 | xmin = MIN(zthe_glo(i),xmin) |
---|
1250 | xmax = MAX(zthe_glo(i),xmax) |
---|
1251 | ENDDO |
---|
1252 | PRINT*,'OROGRAPHIE SOUS-MAILLE zthe:', xmin, xmax |
---|
1253 | c |
---|
1254 | c |
---|
1255 | ierr = NF_INQ_VARID (nid, "ZPIC", nvarid) |
---|
1256 | IF (ierr.NE.NF_NOERR) THEN |
---|
1257 | PRINT*, 'phyetat0: Le champ <ZPIC> est absent' |
---|
1258 | CALL abort |
---|
1259 | ENDIF |
---|
1260 | #ifdef NC_DOUBLE |
---|
1261 | ierr = NF_GET_VAR_DOUBLE(nid, nvarid, zpic_glo) |
---|
1262 | #else |
---|
1263 | ierr = NF_GET_VAR_REAL(nid, nvarid, zpic_glo) |
---|
1264 | #endif |
---|
1265 | IF (ierr.NE.NF_NOERR) THEN |
---|
1266 | PRINT*, 'phyetat0: Lecture echouee pour <ZPIC>' |
---|
1267 | CALL abort |
---|
1268 | ENDIF |
---|
1269 | xmin = 1.0E+20 |
---|
1270 | xmax = -1.0E+20 |
---|
1271 | DO i = 1, klon_glo |
---|
1272 | xmin = MIN(zpic_glo(i),xmin) |
---|
1273 | xmax = MAX(zpic_glo(i),xmax) |
---|
1274 | ENDDO |
---|
1275 | PRINT*,'OROGRAPHIE SOUS-MAILLE zpic:', xmin, xmax |
---|
1276 | c |
---|
1277 | ierr = NF_INQ_VARID (nid, "ZVAL", nvarid) |
---|
1278 | IF (ierr.NE.NF_NOERR) THEN |
---|
1279 | PRINT*, 'phyetat0: Le champ <ZVAL> est absent' |
---|
1280 | CALL abort |
---|
1281 | ENDIF |
---|
1282 | #ifdef NC_DOUBLE |
---|
1283 | ierr = NF_GET_VAR_DOUBLE(nid, nvarid, zval_glo) |
---|
1284 | #else |
---|
1285 | ierr = NF_GET_VAR_REAL(nid, nvarid, zval_glo) |
---|
1286 | #endif |
---|
1287 | IF (ierr.NE.NF_NOERR) THEN |
---|
1288 | PRINT*, 'phyetat0: Lecture echouee pour <ZVAL>' |
---|
1289 | CALL abort |
---|
1290 | ENDIF |
---|
1291 | xmin = 1.0E+20 |
---|
1292 | xmax = -1.0E+20 |
---|
1293 | DO i = 1, klon_glo |
---|
1294 | xmin = MIN(zval_glo(i),xmin) |
---|
1295 | xmax = MAX(zval_glo(i),xmax) |
---|
1296 | ENDDO |
---|
1297 | PRINT*,'OROGRAPHIE SOUS-MAILLE zval:', xmin, xmax |
---|
1298 | c |
---|
1299 | c |
---|
1300 | ierr = NF_INQ_VARID (nid, "RUGSREL", nvarid) |
---|
1301 | IF (ierr.NE.NF_NOERR) THEN |
---|
1302 | PRINT*, 'phyetat0: Le champ <RUGSREL> est absent' |
---|
1303 | CALL abort |
---|
1304 | ENDIF |
---|
1305 | #ifdef NC_DOUBLE |
---|
1306 | ierr = NF_GET_VAR_DOUBLE(nid, nvarid, rugsrel_glo) |
---|
1307 | #else |
---|
1308 | ierr = NF_GET_VAR_REAL(nid, nvarid, rugsrel_glo) |
---|
1309 | #endif |
---|
1310 | IF (ierr.NE.NF_NOERR) THEN |
---|
1311 | PRINT*, 'phyetat0: Lecture echouee pour <RUGSREL>' |
---|
1312 | CALL abort |
---|
1313 | ENDIF |
---|
1314 | xmin = 1.0E+20 |
---|
1315 | xmax = -1.0E+20 |
---|
1316 | DO i = 1, klon_glo |
---|
1317 | xmin = MIN(rugsrel_glo(i),xmin) |
---|
1318 | xmax = MAX(rugsrel_glo(i),xmax) |
---|
1319 | ENDDO |
---|
1320 | PRINT*,'Rugosite relief (ecart-type) rugsrel:', xmin, xmax |
---|
1321 | c |
---|
1322 | c |
---|
1323 | ancien_ok_glo = .TRUE. |
---|
1324 | c |
---|
1325 | ierr = NF_INQ_VARID (nid, "TANCIEN", nvarid) |
---|
1326 | IF (ierr.NE.NF_NOERR) THEN |
---|
1327 | PRINT*, "phyetat0: Le champ <TANCIEN> est absent" |
---|
1328 | PRINT*, "Depart legerement fausse. Mais je continue" |
---|
1329 | ancien_ok_glo = .FALSE. |
---|
1330 | ELSE |
---|
1331 | #ifdef NC_DOUBLE |
---|
1332 | ierr = NF_GET_VAR_DOUBLE(nid, nvarid, t_ancien_glo) |
---|
1333 | #else |
---|
1334 | ierr = NF_GET_VAR_REAL(nid, nvarid, t_ancien_glo) |
---|
1335 | #endif |
---|
1336 | IF (ierr.NE.NF_NOERR) THEN |
---|
1337 | PRINT*, "phyetat0: Lecture echouee pour <TANCIEN>" |
---|
1338 | CALL abort |
---|
1339 | ENDIF |
---|
1340 | ENDIF |
---|
1341 | c |
---|
1342 | ierr = NF_INQ_VARID (nid, "QANCIEN", nvarid) |
---|
1343 | IF (ierr.NE.NF_NOERR) THEN |
---|
1344 | PRINT*, "phyetat0: Le champ <QANCIEN> est absent" |
---|
1345 | PRINT*, "Depart legerement fausse. Mais je continue" |
---|
1346 | ancien_ok_glo = .FALSE. |
---|
1347 | ELSE |
---|
1348 | #ifdef NC_DOUBLE |
---|
1349 | ierr = NF_GET_VAR_DOUBLE(nid, nvarid, q_ancien_glo) |
---|
1350 | #else |
---|
1351 | ierr = NF_GET_VAR_REAL(nid, nvarid, q_ancien_glo) |
---|
1352 | #endif |
---|
1353 | IF (ierr.NE.NF_NOERR) THEN |
---|
1354 | PRINT*, "phyetat0: Lecture echouee pour <QANCIEN>" |
---|
1355 | CALL abort |
---|
1356 | ENDIF |
---|
1357 | ENDIF |
---|
1358 | c |
---|
1359 | clwcon_glo=0. |
---|
1360 | ierr = NF_INQ_VARID (nid, "CLWCON", nvarid) |
---|
1361 | IF (ierr.NE.NF_NOERR) THEN |
---|
1362 | PRINT*, "phyetat0: Le champ CLWCON est absent" |
---|
1363 | PRINT*, "Depart legerement fausse. Mais je continue" |
---|
1364 | c clwcon_glo = 0. |
---|
1365 | ELSE |
---|
1366 | #ifdef NC_DOUBLE |
---|
1367 | ierr = NF_GET_VAR_DOUBLE(nid, nvarid, clwcon_glo) |
---|
1368 | #else |
---|
1369 | ierr = NF_GET_VAR_REAL(nid, nvarid, clwcon_glo) |
---|
1370 | #endif |
---|
1371 | IF (ierr.NE.NF_NOERR) THEN |
---|
1372 | PRINT*, "phyetat0: Lecture echouee pour <CLWCON>" |
---|
1373 | CALL abort |
---|
1374 | ENDIF |
---|
1375 | ENDIF |
---|
1376 | xmin = 1.0E+20 |
---|
1377 | xmax = -1.0E+20 |
---|
1378 | xmin = MINval(clwcon_glo) |
---|
1379 | xmax = MAXval(clwcon_glo) |
---|
1380 | PRINT*,'Eau liquide convective (ecart-type) clwcon:', xmin, xmax |
---|
1381 | c |
---|
1382 | rnebcon_glo = 0. |
---|
1383 | ierr = NF_INQ_VARID (nid, "RNEBCON", nvarid) |
---|
1384 | IF (ierr.NE.NF_NOERR) THEN |
---|
1385 | PRINT*, "phyetat0: Le champ RNEBCON est absent" |
---|
1386 | PRINT*, "Depart legerement fausse. Mais je continue" |
---|
1387 | c rnebcon_glo = 0. |
---|
1388 | ELSE |
---|
1389 | #ifdef NC_DOUBLE |
---|
1390 | ierr = NF_GET_VAR_DOUBLE(nid, nvarid, rnebcon_glo) |
---|
1391 | #else |
---|
1392 | ierr = NF_GET_VAR_REAL(nid, nvarid, rnebcon_glo) |
---|
1393 | #endif |
---|
1394 | IF (ierr.NE.NF_NOERR) THEN |
---|
1395 | PRINT*, "phyetat0: Lecture echouee pour <RNEBCON>" |
---|
1396 | CALL abort |
---|
1397 | ENDIF |
---|
1398 | ENDIF |
---|
1399 | xmin = 1.0E+20 |
---|
1400 | xmax = -1.0E+20 |
---|
1401 | xmin = MINval(rnebcon_glo) |
---|
1402 | xmax = MAXval(rnebcon_glo) |
---|
1403 | PRINT*,'Nebulosite convective (ecart-type) rnebcon:', xmin, xmax |
---|
1404 | |
---|
1405 | c |
---|
1406 | c Lecture ratqs |
---|
1407 | c |
---|
1408 | ratqs_glo=0. |
---|
1409 | ierr = NF_INQ_VARID (nid, "RATQS", nvarid) |
---|
1410 | IF (ierr.NE.NF_NOERR) THEN |
---|
1411 | PRINT*, "phyetat0: Le champ <RATQS> est absent" |
---|
1412 | PRINT*, "Depart legerement fausse. Mais je continue" |
---|
1413 | ratqs_glo = 0. |
---|
1414 | ELSE |
---|
1415 | #ifdef NC_DOUBLE |
---|
1416 | ierr = NF_GET_VAR_DOUBLE(nid, nvarid, ratqs_glo) |
---|
1417 | #else |
---|
1418 | ierr = NF_GET_VAR_REAL(nid, nvarid, ratqs_glo) |
---|
1419 | #endif |
---|
1420 | IF (ierr.NE.NF_NOERR) THEN |
---|
1421 | PRINT*, "phyetat0: Lecture echouee pour <RATQS>" |
---|
1422 | CALL abort |
---|
1423 | ENDIF |
---|
1424 | ENDIF |
---|
1425 | xmin = 1.0E+20 |
---|
1426 | xmax = -1.0E+20 |
---|
1427 | xmin = MINval(ratqs_glo) |
---|
1428 | xmax = MAXval(ratqs_glo) |
---|
1429 | PRINT*,'(ecart-type) ratqs:', xmin, xmax |
---|
1430 | c |
---|
1431 | c Lecture run_off_lic_0 |
---|
1432 | c |
---|
1433 | ierr = NF_INQ_VARID (nid, "RUNOFFLIC0", nvarid) |
---|
1434 | IF (ierr.NE.NF_NOERR) THEN |
---|
1435 | PRINT*, "phyetat0: Le champ <RUNOFFLIC0> est absent" |
---|
1436 | PRINT*, "Depart legerement fausse. Mais je continue" |
---|
1437 | run_off_lic_0 = 0. |
---|
1438 | ELSE |
---|
1439 | #ifdef NC_DOUBLE |
---|
1440 | ierr = NF_GET_VAR_DOUBLE(nid, nvarid, run_off_lic_0) |
---|
1441 | #else |
---|
1442 | ierr = NF_GET_VAR_REAL(nid, nvarid, run_off_lic_0) |
---|
1443 | #endif |
---|
1444 | IF (ierr.NE.NF_NOERR) THEN |
---|
1445 | PRINT*, "phyetat0: Lecture echouee pour <RUNOFFLIC0>" |
---|
1446 | CALL abort |
---|
1447 | ENDIF |
---|
1448 | ENDIF |
---|
1449 | xmin = 1.0E+20 |
---|
1450 | xmax = -1.0E+20 |
---|
1451 | xmin = MINval(run_off_lic_0) |
---|
1452 | xmax = MAXval(run_off_lic_0) |
---|
1453 | PRINT*,'(ecart-type) run_off_lic_0:', xmin, xmax |
---|
1454 | |
---|
1455 | |
---|
1456 | c Lecture de l'energie cinetique turbulente |
---|
1457 | c |
---|
1458 | |
---|
1459 | IF (iflag_pbl>1) then |
---|
1460 | PRINT*, 'phyetat0: Le champ <TKE> est absent' |
---|
1461 | PRINT*, ' Mais je vais essayer de lire TKE**' |
---|
1462 | DO nsrf = 1, nbsrf |
---|
1463 | IF (nsrf.GT.99) THEN |
---|
1464 | PRINT*, "Trop de sous-mailles" |
---|
1465 | CALL abort |
---|
1466 | ENDIF |
---|
1467 | WRITE(str2,'(i2.2)') nsrf |
---|
1468 | ierr = NF_INQ_VARID (nid, "TKE"//str2, nvarid) |
---|
1469 | IF (ierr.NE.NF_NOERR) THEN |
---|
1470 | PRINT*, "WARNING phyetat0: <TKE"//str2//"> est absent" |
---|
1471 | pbl_tke_glo(:,:,nsrf)=1.e-8 |
---|
1472 | ELSE |
---|
1473 | #ifdef NC_DOUBLE |
---|
1474 | ierr = NF_GET_VAR_DOUBLE(nid, nvarid, |
---|
1475 | & pbl_tke_glo(1,1,nsrf)) |
---|
1476 | #else |
---|
1477 | ierr = NF_GET_VAR_REAL(nid, nvarid, pbl_tke_glo(1,1,nsrf)) |
---|
1478 | #endif |
---|
1479 | IF (ierr.NE.NF_NOERR) THEN |
---|
1480 | PRINT*, "WARNING phyetat0: echec lect <TKE"//str2//">" |
---|
1481 | CALL abort |
---|
1482 | ENDIF |
---|
1483 | ENDIF |
---|
1484 | |
---|
1485 | xmin = 1.0E+20 |
---|
1486 | xmax = -1.0E+20 |
---|
1487 | DO k = 1, klev |
---|
1488 | DO i = 1, klon_glo |
---|
1489 | xmin = MIN(pbl_tke_glo(i,k,nsrf),xmin) |
---|
1490 | xmax = MAX(pbl_tke_glo(i,k,nsrf),xmax) |
---|
1491 | ENDDO |
---|
1492 | ENDDO |
---|
1493 | PRINT*,'Temperature du sol TKE**:', nsrf, xmin, xmax |
---|
1494 | ENDDO |
---|
1495 | ENDIF |
---|
1496 | c |
---|
1497 | c zmax0 |
---|
1498 | ierr = NF_INQ_VARID (nid, "ZMAX0", nvarid) |
---|
1499 | IF (ierr.NE.NF_NOERR) THEN |
---|
1500 | PRINT*, "phyetat0: Le champ <ZMAX0> est absent" |
---|
1501 | PRINT*, "Depart legerement fausse. Mais je continue" |
---|
1502 | zmax0_glo=40. |
---|
1503 | ELSE |
---|
1504 | #ifdef NC_DOUBLE |
---|
1505 | ierr = NF_GET_VAR_DOUBLE(nid, nvarid, zmax0_glo) |
---|
1506 | #else |
---|
1507 | ierr = NF_GET_VAR_REAL(nid, nvarid, zmax0_glo) |
---|
1508 | #endif |
---|
1509 | IF (ierr.NE.NF_NOERR) THEN |
---|
1510 | PRINT*, "phyetat0: Lecture echouee pour <ZMAX0>" |
---|
1511 | CALL abort |
---|
1512 | ENDIF |
---|
1513 | ENDIF |
---|
1514 | xmin = 1.0E+20 |
---|
1515 | xmax = -1.0E+20 |
---|
1516 | xmin = MINval(zmax0_glo) |
---|
1517 | xmax = MAXval(zmax0_glo) |
---|
1518 | PRINT*,'(ecart-type) zmax0:', xmin, xmax |
---|
1519 | c |
---|
1520 | c f0(ig)=1.e-5 |
---|
1521 | c f0 |
---|
1522 | ierr = NF_INQ_VARID (nid, "f0", nvarid) |
---|
1523 | IF (ierr.NE.NF_NOERR) THEN |
---|
1524 | PRINT*, "phyetat0: Le champ <f0> est absent" |
---|
1525 | PRINT*, "Depart legerement fausse. Mais je continue" |
---|
1526 | f0_glo=1.e-5 |
---|
1527 | ELSE |
---|
1528 | #ifdef NC_DOUBLE |
---|
1529 | ierr = NF_GET_VAR_DOUBLE(nid, nvarid, f0_glo) |
---|
1530 | #else |
---|
1531 | ierr = NF_GET_VAR_REAL(nid, nvarid, f0_glo) |
---|
1532 | #endif |
---|
1533 | IF (ierr.NE.NF_NOERR) THEN |
---|
1534 | PRINT*, "phyetat0: Lecture echouee pour <f0>" |
---|
1535 | CALL abort |
---|
1536 | ENDIF |
---|
1537 | ENDIF |
---|
1538 | xmin = 1.0E+20 |
---|
1539 | xmax = -1.0E+20 |
---|
1540 | xmin = MINval(f0_glo) |
---|
1541 | xmax = MAXval(f0_glo) |
---|
1542 | PRINT*,'(ecart-type) f0:', xmin, xmax |
---|
1543 | c |
---|
1544 | c ema_work1 |
---|
1545 | c |
---|
1546 | ierr = NF_INQ_VARID (nid, "EMA_WORK1", nvarid) |
---|
1547 | IF (ierr.NE.NF_NOERR) THEN |
---|
1548 | PRINT*, "phyetat0: Le champ <EMA_WORK1> est absent" |
---|
1549 | PRINT*, "Depart legerement fausse. Mais je continue" |
---|
1550 | ema_work1_glo=0. |
---|
1551 | ELSE |
---|
1552 | #ifdef NC_DOUBLE |
---|
1553 | ierr = NF_GET_VAR_DOUBLE(nid, nvarid, ema_work1_glo) |
---|
1554 | #else |
---|
1555 | ierr = NF_GET_VAR_REAL(nid, nvarid, ema_work1_glo) |
---|
1556 | #endif |
---|
1557 | IF (ierr.NE.NF_NOERR) THEN |
---|
1558 | PRINT*, "phyetat0: Lecture echouee pour <EMA_WORK1>" |
---|
1559 | CALL abort |
---|
1560 | ENDIF |
---|
1561 | xmin = 1.0E+20 |
---|
1562 | xmax = -1.0E+20 |
---|
1563 | DO k = 1, klev |
---|
1564 | DO i = 1, klon |
---|
1565 | xmin = MIN(ema_work1_glo(i,k),xmin) |
---|
1566 | xmax = MAX(ema_work1_glo(i,k),xmax) |
---|
1567 | ENDDO |
---|
1568 | ENDDO |
---|
1569 | PRINT*,'ema_work1:', xmin, xmax |
---|
1570 | ENDIF |
---|
1571 | c |
---|
1572 | c ema_work2 |
---|
1573 | c |
---|
1574 | ierr = NF_INQ_VARID (nid, "EMA_WORK2", nvarid) |
---|
1575 | IF (ierr.NE.NF_NOERR) THEN |
---|
1576 | PRINT*, "phyetat0: Le champ <EMA_WORK2> est absent" |
---|
1577 | PRINT*, "Depart legerement fausse. Mais je continue" |
---|
1578 | ema_work2_glo=0. |
---|
1579 | ELSE |
---|
1580 | #ifdef NC_DOUBLE |
---|
1581 | ierr = NF_GET_VAR_DOUBLE(nid, nvarid, ema_work2_glo) |
---|
1582 | #else |
---|
1583 | ierr = NF_GET_VAR_REAL(nid, nvarid, ema_work2_glo) |
---|
1584 | #endif |
---|
1585 | IF (ierr.NE.NF_NOERR) THEN |
---|
1586 | PRINT*, "phyetat0: Lecture echouee pour <EMA_WORK2>" |
---|
1587 | CALL abort |
---|
1588 | ENDIF |
---|
1589 | xmin = 1.0E+20 |
---|
1590 | xmax = -1.0E+20 |
---|
1591 | DO k = 1, klev |
---|
1592 | DO i = 1, klon |
---|
1593 | xmin = MIN(ema_work2_glo(i,k),xmin) |
---|
1594 | xmax = MAX(ema_work2_glo(i,k),xmax) |
---|
1595 | ENDDO |
---|
1596 | ENDDO |
---|
1597 | PRINT*,'ema_work2:', xmin, xmax |
---|
1598 | ENDIF |
---|
1599 | c |
---|
1600 | c wake_deltat |
---|
1601 | c |
---|
1602 | ierr = NF_INQ_VARID (nid, "WAKE_DELTAT", nvarid) |
---|
1603 | IF (ierr.NE.NF_NOERR) THEN |
---|
1604 | PRINT*, "phyetat0: Le champ <WAKE_DELTAT> est absent" |
---|
1605 | PRINT*, "Depart legerement fausse. Mais je continue" |
---|
1606 | wake_deltat_glo=0. |
---|
1607 | ELSE |
---|
1608 | #ifdef NC_DOUBLE |
---|
1609 | ierr = NF_GET_VAR_DOUBLE(nid, nvarid, wake_deltat_glo) |
---|
1610 | #else |
---|
1611 | ierr = NF_GET_VAR_REAL(nid, nvarid, wake_deltat_glo) |
---|
1612 | #endif |
---|
1613 | IF (ierr.NE.NF_NOERR) THEN |
---|
1614 | PRINT*, "phyetat0: Lecture echouee pour <WAKE_DELTAT>" |
---|
1615 | CALL abort |
---|
1616 | ENDIF |
---|
1617 | xmin = 1.0E+20 |
---|
1618 | xmax = -1.0E+20 |
---|
1619 | DO k = 1, klev |
---|
1620 | DO i = 1, klon_glo |
---|
1621 | xmin = MIN(wake_deltat_glo(i,k),xmin) |
---|
1622 | xmax = MAX(wake_deltat_glo(i,k),xmax) |
---|
1623 | ENDDO |
---|
1624 | ENDDO |
---|
1625 | PRINT*,'wake_deltat:', xmin, xmax |
---|
1626 | ENDIF |
---|
1627 | c |
---|
1628 | c wake_deltaq |
---|
1629 | c |
---|
1630 | ierr = NF_INQ_VARID (nid, "WAKE_DELTAQ", nvarid) |
---|
1631 | IF (ierr.NE.NF_NOERR) THEN |
---|
1632 | PRINT*, "phyetat0: Le champ <WAKE_DELTAQ> est absent" |
---|
1633 | PRINT*, "Depart legerement fausse. Mais je continue" |
---|
1634 | wake_deltaq_glo=0. |
---|
1635 | ELSE |
---|
1636 | #ifdef NC_DOUBLE |
---|
1637 | ierr = NF_GET_VAR_DOUBLE(nid, nvarid, wake_deltaq_glo) |
---|
1638 | #else |
---|
1639 | ierr = NF_GET_VAR_REAL(nid, nvarid, wake_deltaq_glo) |
---|
1640 | #endif |
---|
1641 | IF (ierr.NE.NF_NOERR) THEN |
---|
1642 | PRINT*, "phyetat0: Lecture echouee pour <WAKE_DELTAQ>" |
---|
1643 | CALL abort |
---|
1644 | ENDIF |
---|
1645 | xmin = 1.0E+20 |
---|
1646 | xmax = -1.0E+20 |
---|
1647 | DO k = 1, klev |
---|
1648 | DO i = 1, klon_glo |
---|
1649 | xmin = MIN(wake_deltaq_glo(i,k),xmin) |
---|
1650 | xmax = MAX(wake_deltaq_glo(i,k),xmax) |
---|
1651 | ENDDO |
---|
1652 | ENDDO |
---|
1653 | PRINT*,'wake_deltaq:', xmin, xmax |
---|
1654 | ENDIF |
---|
1655 | c |
---|
1656 | c wake_s |
---|
1657 | c |
---|
1658 | ierr = NF_INQ_VARID (nid, "WAKE_S", nvarid) |
---|
1659 | IF (ierr.NE.NF_NOERR) THEN |
---|
1660 | PRINT*, "phyetat0: Le champ <WAKE_S> est absent" |
---|
1661 | PRINT*, "Depart legerement fausse. Mais je continue" |
---|
1662 | wake_s_glo=0. |
---|
1663 | ELSE |
---|
1664 | #ifdef NC_DOUBLE |
---|
1665 | ierr = NF_GET_VAR_DOUBLE(nid, nvarid, wake_s_glo) |
---|
1666 | #else |
---|
1667 | ierr = NF_GET_VAR_REAL(nid, nvarid, wake_s_glo) |
---|
1668 | #endif |
---|
1669 | IF (ierr.NE.NF_NOERR) THEN |
---|
1670 | PRINT*, "phyetat0: Lecture echouee pour <WAKE_S>" |
---|
1671 | CALL abort |
---|
1672 | ENDIF |
---|
1673 | ENDIF |
---|
1674 | xmin = 1.0E+20 |
---|
1675 | xmax = -1.0E+20 |
---|
1676 | xmin = MINval(wake_s_glo) |
---|
1677 | xmax = MAXval(wake_s_glo) |
---|
1678 | PRINT*,'(ecart-type) wake_s:', xmin, xmax |
---|
1679 | c |
---|
1680 | c wake_cstar |
---|
1681 | c |
---|
1682 | ierr = NF_INQ_VARID (nid, "WAKE_CSTAR", nvarid) |
---|
1683 | IF (ierr.NE.NF_NOERR) THEN |
---|
1684 | PRINT*, "phyetat0: Le champ <WAKE_CSTAR> est absent" |
---|
1685 | PRINT*, "Depart legerement fausse. Mais je continue" |
---|
1686 | wake_cstar_glo=0. |
---|
1687 | ELSE |
---|
1688 | #ifdef NC_DOUBLE |
---|
1689 | ierr = NF_GET_VAR_DOUBLE(nid, nvarid, wake_cstar_glo) |
---|
1690 | #else |
---|
1691 | ierr = NF_GET_VAR_REAL(nid, nvarid, wake_cstar_glo) |
---|
1692 | #endif |
---|
1693 | IF (ierr.NE.NF_NOERR) THEN |
---|
1694 | PRINT*, "phyetat0: Lecture echouee pour <WAKE_CSTAR>" |
---|
1695 | CALL abort |
---|
1696 | ENDIF |
---|
1697 | ENDIF |
---|
1698 | xmin = 1.0E+20 |
---|
1699 | xmax = -1.0E+20 |
---|
1700 | xmin = MINval(wake_cstar_glo) |
---|
1701 | xmax = MAXval(wake_cstar_glo) |
---|
1702 | PRINT*,'(ecart-type) wake_cstar:', xmin, xmax |
---|
1703 | c |
---|
1704 | c wake_fip |
---|
1705 | c |
---|
1706 | ierr = NF_INQ_VARID (nid, "WAKE_FIP", nvarid) |
---|
1707 | IF (ierr.NE.NF_NOERR) THEN |
---|
1708 | PRINT*, "phyetat0: Le champ <WAKE_FIP> est absent" |
---|
1709 | PRINT*, "Depart legerement fausse. Mais je continue" |
---|
1710 | wake_fip_glo=0. |
---|
1711 | ELSE |
---|
1712 | #ifdef NC_DOUBLE |
---|
1713 | ierr = NF_GET_VAR_DOUBLE(nid, nvarid, wake_fip_glo) |
---|
1714 | #else |
---|
1715 | ierr = NF_GET_VAR_REAL(nid, nvarid, wake_fip_glo) |
---|
1716 | #endif |
---|
1717 | IF (ierr.NE.NF_NOERR) THEN |
---|
1718 | PRINT*, "phyetat0: Lecture echouee pour <WAKE_FIP>" |
---|
1719 | CALL abort |
---|
1720 | ENDIF |
---|
1721 | ENDIF |
---|
1722 | xmin = 1.0E+20 |
---|
1723 | xmax = -1.0E+20 |
---|
1724 | xmin = MINval(wake_fip_glo) |
---|
1725 | xmax = MAXval(wake_fip_glo) |
---|
1726 | PRINT*,'(ecart-type) wake_fip:', xmin, xmax |
---|
1727 | c |
---|
1728 | c Fermer le fichier: |
---|
1729 | c |
---|
1730 | ierr = NF_CLOSE(nid) |
---|
1731 | ENDIF ! is_mpi_root |
---|
1732 | c |
---|
1733 | c$OMP END MASTER |
---|
1734 | |
---|
1735 | c$OMP MASTER |
---|
1736 | cym en attendant mieux |
---|
1737 | iolat(1)=rlat_glo(1) |
---|
1738 | |
---|
1739 | !FH1D |
---|
1740 | !iolat(jjm+1)=rlat(klon_glo) |
---|
1741 | iolat(jjm+1-1/iim)=rlat_glo(klon_glo) |
---|
1742 | if (iim.gt.1) then |
---|
1743 | do i=2,jjm |
---|
1744 | iolat(i)=rlat_glo(2+(i-2)*iim) |
---|
1745 | enddo |
---|
1746 | endif |
---|
1747 | |
---|
1748 | CALL bcast_mpi(iolat) |
---|
1749 | CALL bcast_mpi(rlon_glo) |
---|
1750 | |
---|
1751 | !FH1D |
---|
1752 | ! call init_iophy(iolat,rlon(2:iim+1)) |
---|
1753 | call init_iophy(iolat,rlon_glo(2-1/iim:iim+1-1/iim)) |
---|
1754 | |
---|
1755 | c$OMP END MASTER |
---|
1756 | |
---|
1757 | call Scatter( rlat_glo,rlat) |
---|
1758 | call Scatter( rlon_glo,rlon) |
---|
1759 | call Scatter( tsol_glo,ftsol) |
---|
1760 | IF (iflag_pbl>1) then |
---|
1761 | call Scatter( pbl_tke_glo,pbl_tke) |
---|
1762 | endif |
---|
1763 | call Scatter( zmax0_glo,zmax0) |
---|
1764 | call Scatter( f0_glo,f0) |
---|
1765 | call Scatter( ema_work1_glo, ema_work1) |
---|
1766 | call Scatter( ema_work2_glo, ema_work2) |
---|
1767 | call Scatter( wake_deltat_glo, wake_deltat) |
---|
1768 | call Scatter( wake_deltaq_glo, wake_deltaq) |
---|
1769 | call Scatter( wake_s_glo, wake_s) |
---|
1770 | call Scatter( wake_cstar_glo, wake_cstar) |
---|
1771 | call Scatter( wake_fip_glo, wake_fip) |
---|
1772 | call Scatter( tsoil,tsoil_p) |
---|
1773 | call Scatter( tslab,tslab_p) |
---|
1774 | call Scatter( seaice,seaice_p) |
---|
1775 | call Scatter( qsurf,qsurf_p) |
---|
1776 | call Scatter( qsol,qsol_p) |
---|
1777 | call Scatter( snow,snow_p) |
---|
1778 | call Scatter( alb1_glo,falb1) |
---|
1779 | call Scatter( alb2_glo,falb2) |
---|
1780 | call Scatter( evap,evap_p) |
---|
1781 | call Scatter( radsol_glo,radsol) |
---|
1782 | call Scatter( rain_fall_glo,rain_fall) |
---|
1783 | call Scatter( snow_fall_glo,snow_fall) |
---|
1784 | call Scatter( sollw_glo,sollw) |
---|
1785 | call Scatter( solsw_glo,solsw) |
---|
1786 | call Scatter( fder,fder_p) |
---|
1787 | call Scatter( frugs,frugs_p) |
---|
1788 | call Scatter( agesno,agesno_p) |
---|
1789 | call Scatter( zmea_glo,zmea) |
---|
1790 | call Scatter( zstd_glo,zstd) |
---|
1791 | call Scatter( zsig_glo,zsig) |
---|
1792 | call Scatter( zgam_glo,zgam) |
---|
1793 | call Scatter( zthe_glo,zthe) |
---|
1794 | call Scatter( zpic_glo,zpic) |
---|
1795 | call Scatter( zval_glo,zval) |
---|
1796 | call Scatter( rugsrel_glo,rugoro) |
---|
1797 | call Scatter( pctsrf_glo,pctsrf) |
---|
1798 | call Scatter( run_off_lic_0,run_off_lic_0_p) |
---|
1799 | call Scatter( t_ancien_glo,t_ancien) |
---|
1800 | call Scatter( q_ancien_glo,q_ancien) |
---|
1801 | call Scatter( rnebcon_glo,rnebcon) |
---|
1802 | call Scatter( clwcon_glo,clwcon) |
---|
1803 | call Scatter( ratqs_glo,ratqs) |
---|
1804 | call Scatter( zmasq_glo,zmasq) |
---|
1805 | |
---|
1806 | c |
---|
1807 | c Initilalize variables in module surface_data |
---|
1808 | c |
---|
1809 | ok_veget = ok_veget_in |
---|
1810 | ocean = ocean_in |
---|
1811 | c |
---|
1812 | c Initialize module pbl_surface_mod |
---|
1813 | c |
---|
1814 | CALL pbl_surface_init(qsol_p, fder_p, snow_p, qsurf_p, |
---|
1815 | $ evap_p, frugs_p, agesno_p, tsoil_p) |
---|
1816 | |
---|
1817 | c Initialize ocean module according to ocean type |
---|
1818 | IF ( ocean == 'slab' ) THEN |
---|
1819 | c initilalize module ocean_slab_init |
---|
1820 | CALL ocean_slab_init(dtime, tslab_p, seaice_p, pctsrf) |
---|
1821 | ELSEIF ( ocean == 'couple' ) THEN |
---|
1822 | c initilalize module ocean_cpl_init |
---|
1823 | CALL ocean_cpl_init(dtime, rlon, rlat) |
---|
1824 | ELSE |
---|
1825 | c initilalize module ocean_forced_init |
---|
1826 | CALL ocean_forced_init |
---|
1827 | ENDIF |
---|
1828 | c |
---|
1829 | c Initilialize module fonte_neige_mod |
---|
1830 | c |
---|
1831 | CALL fonte_neige_init(run_off_lic_0_p) |
---|
1832 | |
---|
1833 | |
---|
1834 | RETURN |
---|
1835 | END |
---|