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

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

Parent Directory Parent Directory | Revision Log Revision Log

Revision 15 - (show annotations)
Fri Aug 1 15:24:12 2008 UTC (15 years, 8 months ago) by guez
File size: 32735 byte(s) 
-- Minor modification of input/output:

Added variable "Sigma_O3_Royer" to "histday.nc". "ecrit_day" is not
modified in "physiq". Removed variables "pyu1", "pyv1", "ftsol1",
"ftsol2", "ftsol3", "ftsol4", "psrf1", "psrf2", "psrf3", "psrf4"
"mfu", "mfd", "en_u", "en_d", "de_d", "de_u", "coefh" from
"histrac.nc".

Variable "raz_date" of module "conf_gcm_m" has logical type instead of
integer type.

-- Should not change any result at run time:

Modified calls to "IOIPSL_Lionel" procedures because the interfaces of
these procedures have been simplified.

Changed name of variable in module "start_init_orog_m": "masque" to
"mask".

Created a module containing procedure "phyredem".

Removed arguments "punjours", "pdayref" and "ptimestep" of procedure
"iniphysiq".

Renamed procedure "gr_phy_write" to "gr_phy_write_2d". Created
procedure "gr_phy_write_3d".

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