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

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

Parent Directory Parent Directory | Revision Log Revision Log

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