/[lmdze]/trunk/libf/phylmd/phyetat0.f90
ViewVC logotype

Contents of /trunk/libf/phylmd/phyetat0.f90

Parent Directory Parent Directory | Revision Log Revision Log

Revision 49 - (show annotations)
Wed Aug 24 11:43:14 2011 UTC (12 years, 9 months ago) by guez
File size: 29432 byte(s) 
LMDZE now uses library Jumble.

Removed all calls to "flinget". Replaced calls to "flinget",
"flininfo", "flinopen_nozoom" by calls to NetCDF95 and Jumble.

Split file "cv_driver.f" into "cv_driver.f90", "cv_flag.f90" and
"cv_thermo.f90".

Bug fix: "QANCIEN" was read twice in "phyeytat0".

In "physiq", initialization of "d_t", "d_u", "d_v" was useless.
1 module phyetat0_m
2
3 use dimphy, only: klon
4
5 IMPLICIT none
6
7 REAL, save:: rlat(klon), rlon(klon) ! latitude and longitude, in degrees
8
9 private klon
10
11 contains
12
13 SUBROUTINE phyetat0(fichnom, pctsrf, tsol, tsoil, ocean, tslab, seaice, &
14 qsurf, qsol, snow, albe, alblw, evap, rain_fall, snow_fall, solsw, &
15 sollw, fder, radsol, frugs, agesno, zmea, zstd, zsig, zgam, zthe, &
16 zpic, zval, t_ancien, q_ancien, ancien_ok, rnebcon, ratqs, clwcon, &
17 run_off_lic_0)
18
19 ! From phylmd/phyetat0.F, version 1.4 2005/06/03 10:03:07
20 ! Author: Z.X. Li (LMD/CNRS)
21 ! Date: 19930818
22 ! Objet : Lecture de l'état initial pour la physique
23
24 USE indicesol, ONLY : epsfra, is_lic, is_oce, is_sic, is_ter, nbsrf
25 USE dimsoil, ONLY : nsoilmx
26 USE temps, ONLY : itau_phy
27 use netcdf, only: nf90_get_att, nf90_global, nf90_inq_varid, NF90_NOERR, &
28 nf90_get_var, NF90_NOWRITE
29 use netcdf95, only: handle_err, nf95_get_var, nf95_close, NF95_OPEN
30 use dimphy, only: zmasq, klev
31
32 CHARACTER(len=*), intent(in):: fichnom
33 REAL tsol(klon, nbsrf)
34 REAL tsoil(klon, nsoilmx, nbsrf)
35 REAL tslab(klon), seaice(klon)
36 REAL qsurf(klon, nbsrf)
37 REAL qsol(klon)
38 REAL snow(klon, nbsrf)
39 REAL albe(klon, nbsrf)
40 REAL alblw(klon, nbsrf)
41 REAL evap(klon, nbsrf)
42 REAL radsol(klon)
43 REAL rain_fall(klon)
44 REAL snow_fall(klon)
45 REAL sollw(klon)
46 real solsw(klon)
47 real fder(klon)
48 REAL frugs(klon, nbsrf)
49 REAL agesno(klon, nbsrf)
50 REAL zmea(klon)
51 REAL, intent(out):: zstd(klon)
52 REAL, intent(out):: zsig(klon)
53 REAL zgam(klon)
54 REAL zthe(klon)
55 REAL zpic(klon)
56 REAL zval(klon)
57 REAL pctsrf(klon, nbsrf)
58 REAL fractint(klon)
59 REAL run_off_lic_0(klon)
60
61 REAL t_ancien(klon, klev), q_ancien(klon, klev)
62 real rnebcon(klon, klev), clwcon(klon, klev), ratqs(klon, klev)
63 LOGICAL, intent(out):: ancien_ok
64
65 CHARACTER(len=*), intent(in):: ocean
66
67 REAL xmin, xmax
68
69 INTEGER ncid, varid
70 INTEGER ierr, i, nsrf, isoil
71 CHARACTER*7 str7
72 CHARACTER*2 str2
73
74 !---------------------------------------------------------------
75
76 print *, "Call sequence information: phyetat0"
77
78 ! Ouvrir le fichier contenant l'etat initial:
79 print *, 'fichnom = ', fichnom
80 call NF95_OPEN(fichnom, NF90_NOWRITE, ncid)
81
82 ierr = nf90_get_att(ncid, nf90_global, "itau_phy", itau_phy)
83 call handle_err("phyetat0 itau_phy", ierr, ncid, nf90_global)
84
85 ! Lecture des latitudes (coordonnees):
86
87 ierr = NF90_INQ_VARID(ncid, "latitude", varid)
88 IF (ierr /= NF90_NOERR) THEN
89 PRINT *, 'phyetat0: Le champ <latitude> est absent'
90 stop 1
91 ENDIF
92 ierr = NF90_GET_VAR(ncid, varid, rlat)
93 IF (ierr /= NF90_NOERR) THEN
94 PRINT *, 'phyetat0: Lecture echouee pour <latitude>'
95 stop 1
96 ENDIF
97
98 ! Lecture des longitudes (coordonnees):
99
100 ierr = NF90_INQ_VARID(ncid, "longitude", varid)
101 IF (ierr /= NF90_NOERR) THEN
102 PRINT *, 'phyetat0: Le champ <longitude> est absent'
103 stop 1
104 ENDIF
105 ierr = NF90_GET_VAR(ncid, varid, rlon)
106 IF (ierr /= NF90_NOERR) THEN
107 PRINT *, 'phyetat0: Lecture echouee pour <latitude>'
108 stop 1
109 ENDIF
110
111 ! Lecture du masque terre mer
112
113 ierr = NF90_INQ_VARID(ncid, "masque", varid)
114 IF (ierr == NF90_NOERR) THEN
115 call nf95_get_var(ncid, varid, zmasq)
116 else
117 PRINT *, 'phyetat0: Le champ <masque> est absent'
118 PRINT *, 'fichier startphy non compatible avec phyetat0'
119 ! stop 1
120 ENDIF
121 ! Lecture des fractions pour chaque sous-surface
122
123 ! initialisation des sous-surfaces
124
125 pctsrf = 0.
126
127 ! fraction de terre
128
129 ierr = NF90_INQ_VARID(ncid, "FTER", varid)
130 IF (ierr == NF90_NOERR) THEN
131 call nf95_get_var(ncid, varid, pctsrf(1 : klon, is_ter))
132 else
133 PRINT *, 'phyetat0: Le champ <FTER> est absent'
134 !$$$ stop 1
135 ENDIF
136
137 ! fraction de glace de terre
138
139 ierr = NF90_INQ_VARID(ncid, "FLIC", varid)
140 IF (ierr == NF90_NOERR) THEN
141 call nf95_get_var(ncid, varid, pctsrf(1 : klon, is_lic))
142 else
143 PRINT *, 'phyetat0: Le champ <FLIC> est absent'
144 !$$$ stop 1
145 ENDIF
146
147 ! fraction d'ocean
148
149 ierr = NF90_INQ_VARID(ncid, "FOCE", varid)
150 IF (ierr == NF90_NOERR) THEN
151 call nf95_get_var(ncid, varid, pctsrf(1 : klon, is_oce))
152 else
153 PRINT *, 'phyetat0: Le champ <FOCE> est absent'
154 !$$$ stop 1
155 ENDIF
156
157 ! fraction glace de mer
158
159 ierr = NF90_INQ_VARID(ncid, "FSIC", varid)
160 IF (ierr == NF90_NOERR) THEN
161 call nf95_get_var(ncid, varid, pctsrf(1 : klon, is_sic))
162 else
163 PRINT *, 'phyetat0: Le champ <FSIC> est absent'
164 !$$$ stop 1
165 ENDIF
166
167 ! Verification de l'adequation entre le masque et les sous-surfaces
168
169 fractint( 1 : klon) = pctsrf(1 : klon, is_ter) &
170 + pctsrf(1 : klon, is_lic)
171 DO i = 1 , klon
172 IF ( abs(fractint(i) - zmasq(i) ) .GT. EPSFRA ) THEN
173 WRITE(*, *) 'phyetat0: attention fraction terre pas ', &
174 'coherente ', i, zmasq(i), pctsrf(i, is_ter) &
175 , pctsrf(i, is_lic)
176 ENDIF
177 END DO
178 fractint(1 : klon) = pctsrf(1 : klon, is_oce) &
179 + pctsrf(1 : klon, is_sic)
180 DO i = 1 , klon
181 IF ( abs( fractint(i) - (1. - zmasq(i))) .GT. EPSFRA ) THEN
182 WRITE(*, *) 'phyetat0 attention fraction ocean pas ', &
183 'coherente ', i, zmasq(i) , pctsrf(i, is_oce) &
184 , pctsrf(i, is_sic)
185 ENDIF
186 END DO
187
188 ! Lecture des temperatures du sol:
189
190 ierr = NF90_INQ_VARID(ncid, "TS", varid)
191 IF (ierr /= NF90_NOERR) THEN
192 PRINT *, 'phyetat0 : Le champ <TS> est absent'
193 PRINT *, ' Mais je vais essayer de lire TS**'
194 DO nsrf = 1, nbsrf
195 IF (nsrf.GT.99) THEN
196 PRINT *, "Trop de sous-mailles"
197 stop 1
198 ENDIF
199 WRITE(str2, '(i2.2)') nsrf
200 ierr = NF90_INQ_VARID(ncid, "TS"//str2, varid)
201 IF (ierr /= NF90_NOERR) THEN
202 PRINT *, "phyetat0: Le champ <TS"//str2//"> est absent"
203 stop 1
204 ENDIF
205 ierr = NF90_GET_VAR(ncid, varid, tsol(1, nsrf))
206 IF (ierr /= NF90_NOERR) THEN
207 PRINT *, "phyetat0: Lecture echouee pour <TS"//str2//">"
208 stop 1
209 ENDIF
210 xmin = 1.0E+20
211 xmax = -1.0E+20
212 DO i = 1, klon
213 xmin = MIN(tsol(i, nsrf), xmin)
214 xmax = MAX(tsol(i, nsrf), xmax)
215 ENDDO
216 PRINT *, 'Temperature du sol TS**:', nsrf, xmin, xmax
217 ENDDO
218 ELSE
219 PRINT *, 'phyetat0: Le champ <TS> est present'
220 PRINT *, ' J ignore donc les autres temperatures TS**'
221 call nf95_get_var(ncid, varid, tsol(:, 1))
222 xmin = 1.0E+20
223 xmax = -1.0E+20
224 DO i = 1, klon
225 xmin = MIN(tsol(i, 1), xmin)
226 xmax = MAX(tsol(i, 1), xmax)
227 ENDDO
228 PRINT *, 'Temperature du sol <TS>', xmin, xmax
229 DO nsrf = 2, nbsrf
230 DO i = 1, klon
231 tsol(i, nsrf) = tsol(i, 1)
232 ENDDO
233 ENDDO
234 ENDIF
235
236 ! Lecture des temperatures du sol profond:
237
238 DO nsrf = 1, nbsrf
239 DO isoil=1, nsoilmx
240 IF (isoil.GT.99 .AND. nsrf.GT.99) THEN
241 PRINT *, "Trop de couches ou sous-mailles"
242 stop 1
243 ENDIF
244 WRITE(str7, '(i2.2, "srf", i2.2)') isoil, nsrf
245 ierr = NF90_INQ_VARID(ncid, 'Tsoil'//str7, varid)
246 IF (ierr /= NF90_NOERR) THEN
247 PRINT *, "phyetat0: Le champ <Tsoil"//str7//"> est absent"
248 PRINT *, " Il prend donc la valeur de surface"
249 DO i=1, klon
250 tsoil(i, isoil, nsrf)=tsol(i, nsrf)
251 ENDDO
252 ELSE
253 ierr = NF90_GET_VAR(ncid, varid, tsoil(1, isoil, nsrf))
254 IF (ierr /= NF90_NOERR) THEN
255 PRINT *, "Lecture echouee pour <Tsoil"//str7//">"
256 stop 1
257 ENDIF
258 ENDIF
259 ENDDO
260 ENDDO
261
262 !IM "slab" ocean
263
264 ! Lecture de tslab (pour slab ocean seulement):
265
266 IF (ocean .eq. 'slab ') then
267 ierr = NF90_INQ_VARID(ncid, "TSLAB", varid)
268 IF (ierr /= NF90_NOERR) THEN
269 PRINT *, "phyetat0: Le champ <TSLAB> est absent"
270 stop 1
271 ENDIF
272 call nf95_get_var(ncid, varid, tslab)
273 xmin = 1.0E+20
274 xmax = -1.0E+20
275 DO i = 1, klon
276 xmin = MIN(tslab(i), xmin)
277 xmax = MAX(tslab(i), xmax)
278 ENDDO
279 PRINT *, 'Ecart de la SST tslab:', xmin, xmax
280
281 ! Lecture de seaice (pour slab ocean seulement):
282
283 ierr = NF90_INQ_VARID(ncid, "SEAICE", varid)
284 IF (ierr /= NF90_NOERR) THEN
285 PRINT *, "phyetat0: Le champ <SEAICE> est absent"
286 stop 1
287 ENDIF
288 call nf95_get_var(ncid, varid, seaice)
289 xmin = 1.0E+20
290 xmax = -1.0E+20
291 DO i = 1, klon
292 xmin = MIN(seaice(i), xmin)
293 xmax = MAX(seaice(i), xmax)
294 ENDDO
295 PRINT *, 'Masse de la glace de mer seaice:', xmin, xmax
296 ELSE
297 tslab = 0.
298 seaice = 0.
299 ENDIF
300
301 ! Lecture de l'humidite de l'air juste au dessus du sol:
302
303 ierr = NF90_INQ_VARID(ncid, "QS", varid)
304 IF (ierr /= NF90_NOERR) THEN
305 PRINT *, 'phyetat0: Le champ <QS> est absent'
306 PRINT *, ' Mais je vais essayer de lire QS**'
307 DO nsrf = 1, nbsrf
308 IF (nsrf.GT.99) THEN
309 PRINT *, "Trop de sous-mailles"
310 stop 1
311 ENDIF
312 WRITE(str2, '(i2.2)') nsrf
313 ierr = NF90_INQ_VARID(ncid, "QS"//str2, varid)
314 IF (ierr /= NF90_NOERR) THEN
315 PRINT *, "phyetat0: Le champ <QS"//str2//"> est absent"
316 stop 1
317 ENDIF
318 ierr = NF90_GET_VAR(ncid, varid, qsurf(1, nsrf))
319 IF (ierr /= NF90_NOERR) THEN
320 PRINT *, "phyetat0: Lecture echouee pour <QS"//str2//">"
321 stop 1
322 ENDIF
323 xmin = 1.0E+20
324 xmax = -1.0E+20
325 DO i = 1, klon
326 xmin = MIN(qsurf(i, nsrf), xmin)
327 xmax = MAX(qsurf(i, nsrf), xmax)
328 ENDDO
329 PRINT *, 'Humidite pres du sol QS**:', nsrf, xmin, xmax
330 ENDDO
331 ELSE
332 PRINT *, 'phyetat0: Le champ <QS> est present'
333 PRINT *, ' J ignore donc les autres humidites QS**'
334 call nf95_get_var(ncid, varid, qsurf(:, 1))
335 xmin = 1.0E+20
336 xmax = -1.0E+20
337 DO i = 1, klon
338 xmin = MIN(qsurf(i, 1), xmin)
339 xmax = MAX(qsurf(i, 1), xmax)
340 ENDDO
341 PRINT *, 'Humidite pres du sol <QS>', xmin, xmax
342 DO nsrf = 2, nbsrf
343 DO i = 1, klon
344 qsurf(i, nsrf) = qsurf(i, 1)
345 ENDDO
346 ENDDO
347 ENDIF
348
349 ! Eau dans le sol (pour le modele de sol "bucket")
350
351 ierr = NF90_INQ_VARID(ncid, "QSOL", varid)
352 IF (ierr == NF90_NOERR) THEN
353 call nf95_get_var(ncid, varid, qsol)
354 else
355 PRINT *, 'phyetat0: Le champ <QSOL> est absent'
356 PRINT *, ' Valeur par defaut nulle'
357 qsol = 0.
358 ENDIF
359 xmin = 1.0E+20
360 xmax = -1.0E+20
361 DO i = 1, klon
362 xmin = MIN(qsol(i), xmin)
363 xmax = MAX(qsol(i), xmax)
364 ENDDO
365 PRINT *, 'Eau dans le sol (mm) <QSOL>', xmin, xmax
366
367 ! Lecture de neige au sol:
368
369 ierr = NF90_INQ_VARID(ncid, "SNOW", varid)
370 IF (ierr /= NF90_NOERR) THEN
371 PRINT *, 'phyetat0: Le champ <SNOW> est absent'
372 PRINT *, ' Mais je vais essayer de lire SNOW**'
373 DO nsrf = 1, nbsrf
374 IF (nsrf.GT.99) THEN
375 PRINT *, "Trop de sous-mailles"
376 stop 1
377 ENDIF
378 WRITE(str2, '(i2.2)') nsrf
379 ierr = NF90_INQ_VARID(ncid, "SNOW"//str2, varid)
380 IF (ierr /= NF90_NOERR) THEN
381 PRINT *, "phyetat0: Le champ <SNOW"//str2//"> est absent"
382 stop 1
383 ENDIF
384 ierr = NF90_GET_VAR(ncid, varid, snow(1, nsrf))
385 IF (ierr /= NF90_NOERR) THEN
386 PRINT *, "phyetat0: Lecture echouee pour <SNOW"//str2//">"
387 stop 1
388 ENDIF
389 xmin = 1.0E+20
390 xmax = -1.0E+20
391 DO i = 1, klon
392 xmin = MIN(snow(i, nsrf), xmin)
393 xmax = MAX(snow(i, nsrf), xmax)
394 ENDDO
395 PRINT *, 'Neige du sol SNOW**:', nsrf, xmin, xmax
396 ENDDO
397 ELSE
398 PRINT *, 'phyetat0: Le champ <SNOW> est present'
399 PRINT *, ' J ignore donc les autres neiges SNOW**'
400 call nf95_get_var(ncid, varid, snow(:, 1))
401 xmin = 1.0E+20
402 xmax = -1.0E+20
403 DO i = 1, klon
404 xmin = MIN(snow(i, 1), xmin)
405 xmax = MAX(snow(i, 1), xmax)
406 ENDDO
407 PRINT *, 'Neige du sol <SNOW>', xmin, xmax
408 DO nsrf = 2, nbsrf
409 DO i = 1, klon
410 snow(i, nsrf) = snow(i, 1)
411 ENDDO
412 ENDDO
413 ENDIF
414
415 ! Lecture de albedo au sol:
416
417 ierr = NF90_INQ_VARID(ncid, "ALBE", varid)
418 IF (ierr /= NF90_NOERR) THEN
419 PRINT *, 'phyetat0: Le champ <ALBE> est absent'
420 PRINT *, ' Mais je vais essayer de lire ALBE**'
421 DO nsrf = 1, nbsrf
422 IF (nsrf.GT.99) THEN
423 PRINT *, "Trop de sous-mailles"
424 stop 1
425 ENDIF
426 WRITE(str2, '(i2.2)') nsrf
427 ierr = NF90_INQ_VARID(ncid, "ALBE"//str2, varid)
428 IF (ierr /= NF90_NOERR) THEN
429 PRINT *, "phyetat0: Le champ <ALBE"//str2//"> est absent"
430 stop 1
431 ENDIF
432 ierr = NF90_GET_VAR(ncid, varid, albe(1, nsrf))
433 IF (ierr /= NF90_NOERR) THEN
434 PRINT *, "phyetat0: Lecture echouee pour <ALBE"//str2//">"
435 stop 1
436 ENDIF
437 xmin = 1.0E+20
438 xmax = -1.0E+20
439 DO i = 1, klon
440 xmin = MIN(albe(i, nsrf), xmin)
441 xmax = MAX(albe(i, nsrf), xmax)
442 ENDDO
443 PRINT *, 'Albedo du sol ALBE**:', nsrf, xmin, xmax
444 ENDDO
445 ELSE
446 PRINT *, 'phyetat0: Le champ <ALBE> est present'
447 PRINT *, ' J ignore donc les autres ALBE**'
448 call nf95_get_var(ncid, varid, albe(:, 1))
449 xmin = 1.0E+20
450 xmax = -1.0E+20
451 DO i = 1, klon
452 xmin = MIN(albe(i, 1), xmin)
453 xmax = MAX(albe(i, 1), xmax)
454 ENDDO
455 PRINT *, 'Neige du sol <ALBE>', xmin, xmax
456 DO nsrf = 2, nbsrf
457 DO i = 1, klon
458 albe(i, nsrf) = albe(i, 1)
459 ENDDO
460 ENDDO
461 ENDIF
462
463 ! Lecture de albedo au sol LW:
464
465 ierr = NF90_INQ_VARID(ncid, "ALBLW", varid)
466 IF (ierr /= NF90_NOERR) THEN
467 PRINT *, 'phyetat0: Le champ <ALBLW> est absent'
468 ! PRINT *, ' Mais je vais essayer de lire ALBLW**'
469 PRINT *, ' Mais je vais prendre ALBE**'
470 DO nsrf = 1, nbsrf
471 DO i = 1, klon
472 alblw(i, nsrf) = albe(i, nsrf)
473 ENDDO
474 ENDDO
475 ELSE
476 PRINT *, 'phyetat0: Le champ <ALBLW> est present'
477 PRINT *, ' J ignore donc les autres ALBLW**'
478 call nf95_get_var(ncid, varid, alblw(:, 1))
479 xmin = 1.0E+20
480 xmax = -1.0E+20
481 DO i = 1, klon
482 xmin = MIN(alblw(i, 1), xmin)
483 xmax = MAX(alblw(i, 1), xmax)
484 ENDDO
485 PRINT *, 'Neige du sol <ALBLW>', xmin, xmax
486 DO nsrf = 2, nbsrf
487 DO i = 1, klon
488 alblw(i, nsrf) = alblw(i, 1)
489 ENDDO
490 ENDDO
491 ENDIF
492
493 ! Lecture de evaporation:
494
495 ierr = NF90_INQ_VARID(ncid, "EVAP", varid)
496 IF (ierr /= NF90_NOERR) THEN
497 PRINT *, 'phyetat0: Le champ <EVAP> est absent'
498 PRINT *, ' Mais je vais essayer de lire EVAP**'
499 DO nsrf = 1, nbsrf
500 IF (nsrf.GT.99) THEN
501 PRINT *, "Trop de sous-mailles"
502 stop 1
503 ENDIF
504 WRITE(str2, '(i2.2)') nsrf
505 ierr = NF90_INQ_VARID(ncid, "EVAP"//str2, varid)
506 IF (ierr /= NF90_NOERR) THEN
507 PRINT *, "phyetat0: Le champ <EVAP"//str2//"> est absent"
508 stop 1
509 ENDIF
510 ierr = NF90_GET_VAR(ncid, varid, evap(1, nsrf))
511 IF (ierr /= NF90_NOERR) THEN
512 PRINT *, "phyetat0: Lecture echouee pour <EVAP"//str2//">"
513 stop 1
514 ENDIF
515 xmin = 1.0E+20
516 xmax = -1.0E+20
517 DO i = 1, klon
518 xmin = MIN(evap(i, nsrf), xmin)
519 xmax = MAX(evap(i, nsrf), xmax)
520 ENDDO
521 PRINT *, 'evap du sol EVAP**:', nsrf, xmin, xmax
522 ENDDO
523 ELSE
524 PRINT *, 'phyetat0: Le champ <EVAP> est present'
525 PRINT *, ' J ignore donc les autres EVAP**'
526 call nf95_get_var(ncid, varid, evap(:, 1))
527 xmin = 1.0E+20
528 xmax = -1.0E+20
529 DO i = 1, klon
530 xmin = MIN(evap(i, 1), xmin)
531 xmax = MAX(evap(i, 1), xmax)
532 ENDDO
533 PRINT *, 'Evap du sol <EVAP>', xmin, xmax
534 DO nsrf = 2, nbsrf
535 DO i = 1, klon
536 evap(i, nsrf) = evap(i, 1)
537 ENDDO
538 ENDDO
539 ENDIF
540
541 ! Lecture precipitation liquide:
542
543 ierr = NF90_INQ_VARID(ncid, "rain_f", varid)
544 IF (ierr /= NF90_NOERR) THEN
545 PRINT *, 'phyetat0: Le champ <rain_f> est absent'
546 stop 1
547 ENDIF
548 ierr = NF90_GET_VAR(ncid, varid, rain_fall)
549 IF (ierr /= NF90_NOERR) THEN
550 PRINT *, 'phyetat0: Lecture echouee pour <rain_f>'
551 stop 1
552 ENDIF
553 xmin = 1.0E+20
554 xmax = -1.0E+20
555 DO i = 1, klon
556 xmin = MIN(rain_fall(i), xmin)
557 xmax = MAX(rain_fall(i), xmax)
558 ENDDO
559 PRINT *, 'Precipitation liquide rain_f:', xmin, xmax
560
561 ! Lecture precipitation solide:
562
563 ierr = NF90_INQ_VARID(ncid, "snow_f", varid)
564 IF (ierr /= NF90_NOERR) THEN
565 PRINT *, 'phyetat0: Le champ <snow_f> est absent'
566 stop 1
567 ENDIF
568 ierr = NF90_GET_VAR(ncid, varid, snow_fall)
569 IF (ierr /= NF90_NOERR) THEN
570 PRINT *, 'phyetat0: Lecture echouee pour <snow_f>'
571 stop 1
572 ENDIF
573 xmin = 1.0E+20
574 xmax = -1.0E+20
575 DO i = 1, klon
576 xmin = MIN(snow_fall(i), xmin)
577 xmax = MAX(snow_fall(i), xmax)
578 ENDDO
579 PRINT *, 'Precipitation solide snow_f:', xmin, xmax
580
581 ! Lecture rayonnement solaire au sol:
582
583 ierr = NF90_INQ_VARID(ncid, "solsw", varid)
584 IF (ierr /= NF90_NOERR) THEN
585 PRINT *, 'phyetat0: Le champ <solsw> est absent'
586 PRINT *, 'mis a zero'
587 solsw = 0.
588 ELSE
589 call nf95_get_var(ncid, varid, solsw)
590 ENDIF
591 xmin = 1.0E+20
592 xmax = -1.0E+20
593 DO i = 1, klon
594 xmin = MIN(solsw(i), xmin)
595 xmax = MAX(solsw(i), xmax)
596 ENDDO
597 PRINT *, 'Rayonnement solaire au sol solsw:', xmin, xmax
598
599 ! Lecture rayonnement IF au sol:
600
601 ierr = NF90_INQ_VARID(ncid, "sollw", varid)
602 IF (ierr /= NF90_NOERR) THEN
603 PRINT *, 'phyetat0: Le champ <sollw> est absent'
604 PRINT *, 'mis a zero'
605 sollw = 0.
606 ELSE
607 call nf95_get_var(ncid, varid, sollw)
608 ENDIF
609 xmin = 1.0E+20
610 xmax = -1.0E+20
611 DO i = 1, klon
612 xmin = MIN(sollw(i), xmin)
613 xmax = MAX(sollw(i), xmax)
614 ENDDO
615 PRINT *, 'Rayonnement IF au sol sollw:', xmin, xmax
616
617 ! Lecture derive des flux:
618
619 ierr = NF90_INQ_VARID(ncid, "fder", varid)
620 IF (ierr /= NF90_NOERR) THEN
621 PRINT *, 'phyetat0: Le champ <fder> est absent'
622 PRINT *, 'mis a zero'
623 fder = 0.
624 ELSE
625 call nf95_get_var(ncid, varid, fder)
626 ENDIF
627 xmin = 1.0E+20
628 xmax = -1.0E+20
629 DO i = 1, klon
630 xmin = MIN(fder(i), xmin)
631 xmax = MAX(fder(i), xmax)
632 ENDDO
633 PRINT *, 'Derive des flux fder:', xmin, xmax
634
635 ! Lecture du rayonnement net au sol:
636
637 ierr = NF90_INQ_VARID(ncid, "RADS", varid)
638 IF (ierr /= NF90_NOERR) THEN
639 PRINT *, 'phyetat0: Le champ <RADS> est absent'
640 stop 1
641 ENDIF
642 ierr = NF90_GET_VAR(ncid, varid, radsol)
643 IF (ierr /= NF90_NOERR) THEN
644 PRINT *, 'phyetat0: Lecture echouee pour <RADS>'
645 stop 1
646 ENDIF
647 xmin = 1.0E+20
648 xmax = -1.0E+20
649 DO i = 1, klon
650 xmin = MIN(radsol(i), xmin)
651 xmax = MAX(radsol(i), xmax)
652 ENDDO
653 PRINT *, 'Rayonnement net au sol radsol:', xmin, xmax
654
655 ! Lecture de la longueur de rugosite
656
657 ierr = NF90_INQ_VARID(ncid, "RUG", varid)
658 IF (ierr /= NF90_NOERR) THEN
659 PRINT *, 'phyetat0: Le champ <RUG> est absent'
660 PRINT *, ' Mais je vais essayer de lire RUG**'
661 DO nsrf = 1, nbsrf
662 IF (nsrf.GT.99) THEN
663 PRINT *, "Trop de sous-mailles"
664 stop 1
665 ENDIF
666 WRITE(str2, '(i2.2)') nsrf
667 ierr = NF90_INQ_VARID(ncid, "RUG"//str2, varid)
668 IF (ierr /= NF90_NOERR) THEN
669 PRINT *, "phyetat0: Le champ <RUG"//str2//"> est absent"
670 stop 1
671 ENDIF
672 ierr = NF90_GET_VAR(ncid, varid, frugs(1, nsrf))
673 IF (ierr /= NF90_NOERR) THEN
674 PRINT *, "phyetat0: Lecture echouee pour <RUG"//str2//">"
675 stop 1
676 ENDIF
677 xmin = 1.0E+20
678 xmax = -1.0E+20
679 DO i = 1, klon
680 xmin = MIN(frugs(i, nsrf), xmin)
681 xmax = MAX(frugs(i, nsrf), xmax)
682 ENDDO
683 PRINT *, 'rugosite du sol RUG**:', nsrf, xmin, xmax
684 ENDDO
685 ELSE
686 PRINT *, 'phyetat0: Le champ <RUG> est present'
687 PRINT *, ' J ignore donc les autres RUG**'
688 call nf95_get_var(ncid, varid, frugs(:, 1))
689 xmin = 1.0E+20
690 xmax = -1.0E+20
691 DO i = 1, klon
692 xmin = MIN(frugs(i, 1), xmin)
693 xmax = MAX(frugs(i, 1), xmax)
694 ENDDO
695 PRINT *, 'rugosite <RUG>', xmin, xmax
696 DO nsrf = 2, nbsrf
697 DO i = 1, klon
698 frugs(i, nsrf) = frugs(i, 1)
699 ENDDO
700 ENDDO
701 ENDIF
702
703 ! Lecture de l'age de la neige:
704
705 ierr = NF90_INQ_VARID(ncid, "AGESNO", varid)
706 IF (ierr /= NF90_NOERR) THEN
707 PRINT *, 'phyetat0: Le champ <AGESNO> est absent'
708 PRINT *, ' Mais je vais essayer de lire AGESNO**'
709 DO nsrf = 1, nbsrf
710 IF (nsrf.GT.99) THEN
711 PRINT *, "Trop de sous-mailles"
712 stop 1
713 ENDIF
714 WRITE(str2, '(i2.2)') nsrf
715 ierr = NF90_INQ_VARID(ncid, "AGESNO"//str2, varid)
716 IF (ierr /= NF90_NOERR) THEN
717 PRINT *, "phyetat0: Le champ <AGESNO"//str2//"> est absent"
718 agesno = 50.0
719 ENDIF
720 ierr = NF90_GET_VAR(ncid, varid, agesno(1, nsrf))
721 IF (ierr /= NF90_NOERR) THEN
722 PRINT *, "phyetat0: Lecture echouee pour <AGESNO"//str2//">"
723 stop 1
724 ENDIF
725 xmin = 1.0E+20
726 xmax = -1.0E+20
727 DO i = 1, klon
728 xmin = MIN(agesno(i, nsrf), xmin)
729 xmax = MAX(agesno(i, nsrf), xmax)
730 ENDDO
731 PRINT *, 'Age de la neige AGESNO**:', nsrf, xmin, xmax
732 ENDDO
733 ELSE
734 PRINT *, 'phyetat0: Le champ <AGESNO> est present'
735 PRINT *, ' J ignore donc les autres AGESNO**'
736 call nf95_get_var(ncid, varid, agesno(:, 1))
737 xmin = 1.0E+20
738 xmax = -1.0E+20
739 DO i = 1, klon
740 xmin = MIN(agesno(i, 1), xmin)
741 xmax = MAX(agesno(i, 1), xmax)
742 ENDDO
743 PRINT *, 'Age de la neige <AGESNO>', xmin, xmax
744 DO nsrf = 2, nbsrf
745 DO i = 1, klon
746 agesno(i, nsrf) = agesno(i, 1)
747 ENDDO
748 ENDDO
749 ENDIF
750
751 ierr = NF90_INQ_VARID(ncid, "ZMEA", varid)
752 IF (ierr /= NF90_NOERR) THEN
753 PRINT *, 'phyetat0: Le champ <ZMEA> est absent'
754 stop 1
755 ENDIF
756 ierr = NF90_GET_VAR(ncid, varid, zmea)
757 IF (ierr /= NF90_NOERR) THEN
758 PRINT *, 'phyetat0: Lecture echouee pour <ZMEA>'
759 stop 1
760 ENDIF
761 xmin = 1.0E+20
762 xmax = -1.0E+20
763 DO i = 1, klon
764 xmin = MIN(zmea(i), xmin)
765 xmax = MAX(zmea(i), xmax)
766 ENDDO
767 PRINT *, 'OROGRAPHIE SOUS-MAILLE zmea:', xmin, xmax
768
769 ierr = NF90_INQ_VARID(ncid, "ZSTD", varid)
770 IF (ierr /= NF90_NOERR) THEN
771 PRINT *, 'phyetat0: Le champ <ZSTD> est absent'
772 stop 1
773 ENDIF
774 ierr = NF90_GET_VAR(ncid, varid, zstd)
775 IF (ierr /= NF90_NOERR) THEN
776 PRINT *, 'phyetat0: Lecture echouee pour <ZSTD>'
777 stop 1
778 ENDIF
779 xmin = 1.0E+20
780 xmax = -1.0E+20
781 DO i = 1, klon
782 xmin = MIN(zstd(i), xmin)
783 xmax = MAX(zstd(i), xmax)
784 ENDDO
785 PRINT *, 'OROGRAPHIE SOUS-MAILLE zstd:', xmin, xmax
786
787 ierr = NF90_INQ_VARID(ncid, "ZSIG", varid)
788 IF (ierr /= NF90_NOERR) THEN
789 PRINT *, 'phyetat0: Le champ <ZSIG> est absent'
790 stop 1
791 ENDIF
792 ierr = NF90_GET_VAR(ncid, varid, zsig)
793 IF (ierr /= NF90_NOERR) THEN
794 PRINT *, 'phyetat0: Lecture echouee pour <ZSIG>'
795 stop 1
796 ENDIF
797 xmin = 1.0E+20
798 xmax = -1.0E+20
799 DO i = 1, klon
800 xmin = MIN(zsig(i), xmin)
801 xmax = MAX(zsig(i), xmax)
802 ENDDO
803 PRINT *, 'OROGRAPHIE SOUS-MAILLE zsig:', xmin, xmax
804
805 ierr = NF90_INQ_VARID(ncid, "ZGAM", varid)
806 IF (ierr /= NF90_NOERR) THEN
807 PRINT *, 'phyetat0: Le champ <ZGAM> est absent'
808 stop 1
809 ENDIF
810 ierr = NF90_GET_VAR(ncid, varid, zgam)
811 IF (ierr /= NF90_NOERR) THEN
812 PRINT *, 'phyetat0: Lecture echouee pour <ZGAM>'
813 stop 1
814 ENDIF
815 xmin = 1.0E+20
816 xmax = -1.0E+20
817 DO i = 1, klon
818 xmin = MIN(zgam(i), xmin)
819 xmax = MAX(zgam(i), xmax)
820 ENDDO
821 PRINT *, 'OROGRAPHIE SOUS-MAILLE zgam:', xmin, xmax
822
823 ierr = NF90_INQ_VARID(ncid, "ZTHE", varid)
824 IF (ierr /= NF90_NOERR) THEN
825 PRINT *, 'phyetat0: Le champ <ZTHE> est absent'
826 stop 1
827 ENDIF
828 ierr = NF90_GET_VAR(ncid, varid, zthe)
829 IF (ierr /= NF90_NOERR) THEN
830 PRINT *, 'phyetat0: Lecture echouee pour <ZTHE>'
831 stop 1
832 ENDIF
833 xmin = 1.0E+20
834 xmax = -1.0E+20
835 DO i = 1, klon
836 xmin = MIN(zthe(i), xmin)
837 xmax = MAX(zthe(i), xmax)
838 ENDDO
839 PRINT *, 'OROGRAPHIE SOUS-MAILLE zthe:', xmin, xmax
840
841 ierr = NF90_INQ_VARID(ncid, "ZPIC", varid)
842 IF (ierr /= NF90_NOERR) THEN
843 PRINT *, 'phyetat0: Le champ <ZPIC> est absent'
844 stop 1
845 ENDIF
846 ierr = NF90_GET_VAR(ncid, varid, zpic)
847 IF (ierr /= NF90_NOERR) THEN
848 PRINT *, 'phyetat0: Lecture echouee pour <ZPIC>'
849 stop 1
850 ENDIF
851 xmin = 1.0E+20
852 xmax = -1.0E+20
853 DO i = 1, klon
854 xmin = MIN(zpic(i), xmin)
855 xmax = MAX(zpic(i), xmax)
856 ENDDO
857 PRINT *, 'OROGRAPHIE SOUS-MAILLE zpic:', xmin, xmax
858
859 ierr = NF90_INQ_VARID(ncid, "ZVAL", varid)
860 IF (ierr /= NF90_NOERR) THEN
861 PRINT *, 'phyetat0: Le champ <ZVAL> est absent'
862 stop 1
863 ENDIF
864 ierr = NF90_GET_VAR(ncid, varid, zval)
865 IF (ierr /= NF90_NOERR) THEN
866 PRINT *, 'phyetat0: Lecture echouee pour <ZVAL>'
867 stop 1
868 ENDIF
869 xmin = 1.0E+20
870 xmax = -1.0E+20
871 DO i = 1, klon
872 xmin = MIN(zval(i), xmin)
873 xmax = MAX(zval(i), xmax)
874 ENDDO
875 PRINT *, 'OROGRAPHIE SOUS-MAILLE zval:', xmin, xmax
876
877 ancien_ok = .TRUE.
878
879 ierr = NF90_INQ_VARID(ncid, "TANCIEN", varid)
880 IF (ierr /= NF90_NOERR) THEN
881 PRINT *, "phyetat0: Le champ <TANCIEN> est absent"
882 PRINT *, "Depart legerement fausse. Mais je continue"
883 ancien_ok = .FALSE.
884 ELSE
885 call nf95_get_var(ncid, varid, t_ancien)
886 ENDIF
887
888 ierr = NF90_INQ_VARID(ncid, "QANCIEN", varid)
889 IF (ierr /= NF90_NOERR) THEN
890 PRINT *, "phyetat0: Le champ <QANCIEN> est absent"
891 PRINT *, "Depart legerement fausse. Mais je continue"
892 ancien_ok = .FALSE.
893 ELSE
894 call nf95_get_var(ncid, varid, q_ancien)
895 ENDIF
896
897 ierr = NF90_INQ_VARID(ncid, "CLWCON", varid)
898 IF (ierr /= NF90_NOERR) THEN
899 PRINT *, "phyetat0: Le champ CLWCON est absent"
900 PRINT *, "Depart legerement fausse. Mais je continue"
901 clwcon = 0.
902 ELSE
903 call nf95_get_var(ncid, varid, clwcon)
904 ENDIF
905 xmin = 1.0E+20
906 xmax = -1.0E+20
907 xmin = MINval(clwcon)
908 xmax = MAXval(clwcon)
909 PRINT *, 'Eau liquide convective (ecart-type) clwcon:', xmin, xmax
910
911 ierr = NF90_INQ_VARID(ncid, "RNEBCON", varid)
912 IF (ierr /= NF90_NOERR) THEN
913 PRINT *, "phyetat0: Le champ RNEBCON est absent"
914 PRINT *, "Depart legerement fausse. Mais je continue"
915 rnebcon = 0.
916 ELSE
917 call nf95_get_var(ncid, varid, rnebcon)
918 ENDIF
919 xmin = 1.0E+20
920 xmax = -1.0E+20
921 xmin = MINval(rnebcon)
922 xmax = MAXval(rnebcon)
923 PRINT *, 'Nebulosite convective (ecart-type) rnebcon:', xmin, xmax
924
925 ! Lecture ratqs
926
927 ierr = NF90_INQ_VARID(ncid, "RATQS", varid)
928 IF (ierr /= NF90_NOERR) THEN
929 PRINT *, "phyetat0: Le champ <RATQS> est absent"
930 PRINT *, "Depart legerement fausse. Mais je continue"
931 ratqs = 0.
932 ELSE
933 call nf95_get_var(ncid, varid, ratqs)
934 ENDIF
935 xmin = 1.0E+20
936 xmax = -1.0E+20
937 xmin = MINval(ratqs)
938 xmax = MAXval(ratqs)
939 PRINT *, '(ecart-type) ratqs:', xmin, xmax
940
941 ! Lecture run_off_lic_0
942
943 ierr = NF90_INQ_VARID(ncid, "RUNOFFLIC0", varid)
944 IF (ierr /= NF90_NOERR) THEN
945 PRINT *, "phyetat0: Le champ <RUNOFFLIC0> est absent"
946 PRINT *, "Depart legerement fausse. Mais je continue"
947 run_off_lic_0 = 0.
948 ELSE
949 call nf95_get_var(ncid, varid, run_off_lic_0)
950 ENDIF
951 xmin = 1.0E+20
952 xmax = -1.0E+20
953 xmin = MINval(run_off_lic_0)
954 xmax = MAXval(run_off_lic_0)
955 PRINT *, '(ecart-type) run_off_lic_0:', xmin, xmax
956
957 call NF95_CLOSE(ncid)
958
959 END SUBROUTINE phyetat0
960
961 end module phyetat0_m
  ViewVC Help
Powered by ViewVC 1.1.21