/[lmdze]/trunk/Sources/phylmd/phyetat0.f
ViewVC logotype

Contents of /trunk/Sources/phylmd/phyetat0.f

Parent Directory Parent Directory | Revision Log Revision Log

Revision 43 - (show annotations)
Fri Apr 8 12:43:31 2011 UTC (13 years, 1 month ago) by guez
Original Path: trunk/libf/phylmd/phyetat0.f90
File size: 29762 byte(s) 
"start_init_phys" is now called directly by "etat0" instead of through
"start_init_dyn". "qsol_2d" is no longer a variable of module
"start_init_phys_m", it is an argument of
"start_init_phys". "start_init_dyn" now receives "tsol_2d" from
"etat0".

Split file "vlspltqs.f" into "vlspltqs.f90", "vlxqs.f90" and
""vlyqs.f90".

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