Context Navigation

phyetat0.F @ 719

Last change on this file since 719 was 689, checked in by lmdzadmin, 19 years ago
Commentaires variables slab IM
Property svn:eol-style set to `native` Property svn:keywords set to `Author Date Id Revision`
File size: 40.7 KB

Line
1	!
2	! $Header$
3	!
4	c
5	c
6	SUBROUTINE phyetat0 (fichnom,dtime,co2_ppm_etat0,solaire_etat0,
7	. rlat,rlon, pctsrf, tsol,tsoil,
8	cIM "slab" ocean
9	. ocean, tslab,seaice,
10	. qsurf,qsol,snow,
11	. albe, alblw, evap, rain_fall, snow_fall, solsw, sollw,
12	. fder,radsol,frugs,agesno,clesphy0,
13	. zmea,zstd,zsig,zgam,zthe,zpic,zval,rugsrel,tabcntr0,
14	. t_ancien,q_ancien,ancien_ok, rnebcon, ratqs,clwcon,
15	. run_off_lic_0)
16	IMPLICIT none
17	c======================================================================
18	c Auteur(s) Z.X. Li (LMD/CNRS) date: 19930818
19	c Objet: Lecture de l'etat initial pour la physique
20	c======================================================================
21	#include "dimensions.h"
22	#include "dimphy.h"
23	#include "netcdf.inc"
24	#include "indicesol.h"
25	#include "dimsoil.h"
26	#include "clesphys.h"
27	#include "temps.h"
28	c======================================================================
29	CHARACTER() fichnom
30	REAL dtime
31	INTEGER radpas
32	REAL rlat(klon), rlon(klon)
33	REAL co2_ppm_etat0
34	REAL solaire_etat0
35	REAL tsol(klon,nbsrf)
36	REAL tsoil(klon,nsoilmx,nbsrf)
37	cIM "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 zstd(klon)
55	REAL zsig(klon)
56	REAL zgam(klon)
57	REAL zthe(klon)
58	REAL zpic(klon)
59	REAL zval(klon)
60	REAL rugsrel(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*6 ocean
70
71	INTEGER longcles
72	PARAMETER ( longcles = 20 )
73	REAL clesphy0( longcles )
74	c
75	REAL xmin, xmax
76	c
77	INTEGER nid, nvarid
78	INTEGER ierr, i, nsrf, isoil
79	INTEGER length
80	PARAMETER (length=100)
81	REAL tab_cntrl(length), tabcntr0(length)
82	CHARACTER*7 str7
83	CHARACTER*2 str2
84	c
85	c Ouvrir le fichier contenant l'etat initial:
86	c
87	print*,'fichnom',fichnom
88	ierr = NF_OPEN (fichnom, NF_NOWRITE,nid)
89	IF (ierr.NE.NF_NOERR) THEN
90	write(6,*)' Pb d''ouverture du fichier '//fichnom
91	write(6,*)' ierr = ', ierr
92	CALL ABORT
93	ENDIF
94	c
95	c Lecture des parametres de controle:
96	c
97	ierr = NF_INQ_VARID (nid, "controle", nvarid)
98	IF (ierr.NE.NF_NOERR) THEN
99	PRINT*, 'phyetat0: Le champ <controle> est absent'
100	CALL abort
101	ENDIF
102	#ifdef NC_DOUBLE
103	ierr = NF_GET_VAR_DOUBLE(nid, nvarid, tab_cntrl)
104	#else
105	ierr = NF_GET_VAR_REAL(nid, nvarid, tab_cntrl)
106	#endif
107	IF (ierr.NE.NF_NOERR) THEN
108	PRINT*, 'phyetat0: Lecture echouee pour <controle>'
109	CALL abort
110	ELSE
111	c
112	DO i = 1, length
113	tabcntr0( i ) = tab_cntrl( i )
114	ENDDO
115	c
116	cycle_diurne = .FALSE.
117	soil_model = .FALSE.
118	new_oliq = .FALSE.
119	ok_orodr = .FALSE.
120	ok_orolf = .FALSE.
121	ok_limitvrai = .FALSE.
122
123
124	IF( clesphy0(1).NE.tab_cntrl( 5 ) ) THEN
125	tab_cntrl( 5 ) = clesphy0(1)
126	ENDIF
127
128	IF( clesphy0(2).NE.tab_cntrl( 6 ) ) THEN
129	tab_cntrl( 6 ) = clesphy0(2)
130	ENDIF
131
132	IF( clesphy0(3).NE.tab_cntrl( 7 ) ) THEN
133	tab_cntrl( 7 ) = clesphy0(3)
134	ENDIF
135
136	IF( clesphy0(4).NE.tab_cntrl( 8 ) ) THEN
137	tab_cntrl( 8 ) = clesphy0(4)
138	ENDIF
139
140	IF( clesphy0(5).NE.tab_cntrl( 9 ) ) THEN
141	tab_cntrl( 9 ) = clesphy0( 5 )
142	ENDIF
143
144	IF( clesphy0(6).NE.tab_cntrl( 10 ) ) THEN
145	tab_cntrl( 10 ) = clesphy0( 6 )
146	ENDIF
147
148	IF( clesphy0(7).NE.tab_cntrl( 11 ) ) THEN
149	tab_cntrl( 11 ) = clesphy0( 7 )
150	ENDIF
151
152	IF( clesphy0(8).NE.tab_cntrl( 12 ) ) THEN
153	tab_cntrl( 12 ) = clesphy0( 8 )
154	ENDIF
155
156
157	dtime = tab_cntrl(1)
158	radpas = tab_cntrl(2)
159	co2_ppm_etat0 = tab_cntrl(3)
160	solaire_etat0 = tab_cntrl(4)
161	iflag_con = tab_cntrl(5)
162	nbapp_rad = tab_cntrl(6)
163
164
165	cycle_diurne = .FALSE.
166	soil_model = .FALSE.
167	new_oliq = .FALSE.
168	ok_orodr = .FALSE.
169	ok_orolf = .FALSE.
170	ok_limitvrai = .FALSE.
171
172	IF( tab_cntrl( 7) .EQ. 1. ) cycle_diurne = .TRUE.
173	IF( tab_cntrl( 8) .EQ. 1. ) soil_model = .TRUE.
174	IF( tab_cntrl( 9) .EQ. 1. ) new_oliq = .TRUE.
175	IF( tab_cntrl(10) .EQ. 1. ) ok_orodr = .TRUE.
176	IF( tab_cntrl(11) .EQ. 1. ) ok_orolf = .TRUE.
177	IF( tab_cntrl(12) .EQ. 1. ) ok_limitvrai = .TRUE.
178
179	ENDIF
180
181	itau_phy = tab_cntrl(15)
182
183	c
184	c Lecture des latitudes (coordonnees):
185	c
186	ierr = NF_INQ_VARID (nid, "latitude", nvarid)
187	IF (ierr.NE.NF_NOERR) THEN
188	PRINT*, 'phyetat0: Le champ <latitude> est absent'
189	CALL abort
190	ENDIF
191	#ifdef NC_DOUBLE
192	ierr = NF_GET_VAR_DOUBLE(nid, nvarid, rlat)
193	#else
194	ierr = NF_GET_VAR_REAL(nid, nvarid, rlat)
195	#endif
196	IF (ierr.NE.NF_NOERR) THEN
197	PRINT*, 'phyetat0: Lecture echouee pour <latitude>'
198	CALL abort
199	ENDIF
200	c
201	c Lecture des longitudes (coordonnees):
202	c
203	ierr = NF_INQ_VARID (nid, "longitude", nvarid)
204	IF (ierr.NE.NF_NOERR) THEN
205	PRINT*, 'phyetat0: Le champ <longitude> est absent'
206	CALL abort
207	ENDIF
208	#ifdef NC_DOUBLE
209	ierr = NF_GET_VAR_DOUBLE(nid, nvarid, rlon)
210	#else
211	ierr = NF_GET_VAR_REAL(nid, nvarid, rlon)
212	#endif
213	IF (ierr.NE.NF_NOERR) THEN
214	PRINT*, 'phyetat0: Lecture echouee pour <latitude>'
215	CALL abort
216	ENDIF
217	C
218	C
219	C Lecture du masque terre mer
220	C
221	ierr = NF_INQ_VARID (nid, "masque", nvarid)
222	IF (ierr .EQ. NF_NOERR) THEN
223	#ifdef NC_DOUBLE
224	ierr = NF_GET_VAR_DOUBLE(nid, nvarid, zmasq)
225	#else
226	ierr = NF_GET_VAR_REAL(nid, nvarid, zmasq)
227	#endif
228	IF (ierr.NE.NF_NOERR) THEN
229	PRINT*, 'phyetat0: Lecture echouee pour <masque>'
230	CALL abort
231	ENDIF
232	else
233	PRINT*, 'phyetat0: Le champ <masque> est absent'
234	PRINT*, 'fichier startphy non compatible avec phyetat0'
235	C CALL abort
236	ENDIF
237	C Lecture des fractions pour chaque sous-surface
238	C
239	C initialisation des sous-surfaces
240	C
241	pctsrf = 0.
242	C
243	C fraction de terre
244	C
245	ierr = NF_INQ_VARID (nid, "FTER", nvarid)
246	IF (ierr .EQ. NF_NOERR) THEN
247	#ifdef NC_DOUBLE
248	ierr = NF_GET_VAR_DOUBLE(nid, nvarid, pctsrf(1 : klon,is_ter))
249	#else
250	ierr = NF_GET_VAR_REAL(nid, nvarid, pctsrf(1 : klon,is_ter))
251	#endif
252	IF (ierr.NE.NF_NOERR) THEN
253	PRINT*, 'phyetat0: Lecture echouee pour <FTER>'
254	CALL abort
255	ENDIF
256	else
257	PRINT*, 'phyetat0: Le champ <FTER> est absent'
258	c$$$ CALL abort
259	ENDIF
260	C
261	C fraction de glace de terre
262	C
263	ierr = NF_INQ_VARID (nid, "FLIC", nvarid)
264	IF (ierr .EQ. NF_NOERR) THEN
265	#ifdef NC_DOUBLE
266	ierr = NF_GET_VAR_DOUBLE(nid, nvarid, pctsrf(1 : klon,is_lic))
267	#else
268	ierr = NF_GET_VAR_REAL(nid, nvarid, pctsrf(1 : klon,is_lic))
269	#endif
270	IF (ierr.NE.NF_NOERR) THEN
271	PRINT*, 'phyetat0: Lecture echouee pour <FLIC>'
272	CALL abort
273	ENDIF
274	else
275	PRINT*, 'phyetat0: Le champ <FLIC> est absent'
276	c$$$ CALL abort
277	ENDIF
278	C
279	C fraction d'ocean
280	C
281	ierr = NF_INQ_VARID (nid, "FOCE", nvarid)
282	IF (ierr .EQ. NF_NOERR) THEN
283	#ifdef NC_DOUBLE
284	ierr = NF_GET_VAR_DOUBLE(nid, nvarid, pctsrf(1 : klon,is_oce))
285	#else
286	ierr = NF_GET_VAR_REAL(nid, nvarid, pctsrf(1 : klon,is_oce))
287	#endif
288	IF (ierr.NE.NF_NOERR) THEN
289	PRINT*, 'phyetat0: Lecture echouee pour <FOCE>'
290	CALL abort
291	ENDIF
292	else
293	PRINT*, 'phyetat0: Le champ <FOCE> est absent'
294	c$$$ CALL abort
295	ENDIF
296	C
297	C fraction glace de mer
298	C
299	ierr = NF_INQ_VARID (nid, "FSIC", nvarid)
300	IF (ierr .EQ. NF_NOERR) THEN
301	#ifdef NC_DOUBLE
302	ierr = NF_GET_VAR_DOUBLE(nid, nvarid, pctsrf(1 : klon,is_sic))
303	#else
304	ierr = NF_GET_VAR_REAL(nid, nvarid, pctsrf(1 : klon, is_sic))
305	#endif
306	IF (ierr.NE.NF_NOERR) THEN
307	PRINT*, 'phyetat0: Lecture echouee pour <FSIC>'
308	CALL abort
309	ENDIF
310	else
311	PRINT*, 'phyetat0: Le champ <FSIC> est absent'
312	c$$$ CALL abort
313	ENDIF
314	C
315	C Verification de l'adequation entre le masque et les sous-surfaces
316	C
317	fractint( 1 : klon) = pctsrf(1 : klon, is_ter)
318	$ + pctsrf(1 : klon, is_lic)
319	DO i = 1 , klon
320	IF ( abs(fractint(i) - zmasq(i) ) .GT. EPSFRA ) THEN
321	WRITE(,) 'phyetat0: attention fraction terre pas ',
322	$ 'coherente ', i, zmasq(i), pctsrf(i, is_ter)
323	$ ,pctsrf(i, is_lic)
324	ENDIF
325	END DO
326	fractint (1 : klon) = pctsrf(1 : klon, is_oce)
327	$ + pctsrf(1 : klon, is_sic)
328	DO i = 1 , klon
329	IF ( abs( fractint(i) - (1. - zmasq(i))) .GT. EPSFRA ) THEN
330	WRITE(,) 'phyetat0 attention fraction ocean pas ',
331	$ 'coherente ', i, zmasq(i) , pctsrf(i, is_oce)
332	$ ,pctsrf(i, is_sic)
333	ENDIF
334	END DO
335	C
336	c Lecture des temperatures du sol:
337	c
338	ierr = NF_INQ_VARID (nid, "TS", nvarid)
339	IF (ierr.NE.NF_NOERR) THEN
340	PRINT*, 'phyetat0: Le champ <TS> est absent'
341	PRINT, ' Mais je vais essayer de lire TS*'
342	DO nsrf = 1, nbsrf
343	IF (nsrf.GT.99) THEN
344	PRINT*, "Trop de sous-mailles"
345	CALL abort
346	ENDIF
347	WRITE(str2,'(i2.2)') nsrf
348	ierr = NF_INQ_VARID (nid, "TS"//str2, nvarid)
349	IF (ierr.NE.NF_NOERR) THEN
350	PRINT*, "phyetat0: Le champ <TS"//str2//"> est absent"
351	CALL abort
352	ENDIF
353	#ifdef NC_DOUBLE
354	ierr = NF_GET_VAR_DOUBLE(nid, nvarid, tsol(1,nsrf))
355	#else
356	ierr = NF_GET_VAR_REAL(nid, nvarid, tsol(1,nsrf))
357	#endif
358	IF (ierr.NE.NF_NOERR) THEN
359	PRINT*, "phyetat0: Lecture echouee pour <TS"//str2//">"
360	CALL abort
361	ENDIF
362	xmin = 1.0E+20
363	xmax = -1.0E+20
364	DO i = 1, klon
365	xmin = MIN(tsol(i,nsrf),xmin)
366	xmax = MAX(tsol(i,nsrf),xmax)
367	ENDDO
368	PRINT,'Temperature du sol TS*:', nsrf, xmin, xmax
369	ENDDO
370	ELSE
371	PRINT*, 'phyetat0: Le champ <TS> est present'
372	PRINT, ' J ignore donc les autres temperatures TS*'
373	#ifdef NC_DOUBLE
374	ierr = NF_GET_VAR_DOUBLE(nid, nvarid, tsol(1,1))
375	#else
376	ierr = NF_GET_VAR_REAL(nid, nvarid, tsol(1,1))
377	#endif
378	IF (ierr.NE.NF_NOERR) THEN
379	PRINT*, "phyetat0: Lecture echouee pour <TS>"
380	CALL abort
381	ENDIF
382	xmin = 1.0E+20
383	xmax = -1.0E+20
384	DO i = 1, klon
385	xmin = MIN(tsol(i,1),xmin)
386	xmax = MAX(tsol(i,1),xmax)
387	ENDDO
388	PRINT*,'Temperature du sol <TS>', xmin, xmax
389	DO nsrf = 2, nbsrf
390	DO i = 1, klon
391	tsol(i,nsrf) = tsol(i,1)
392	ENDDO
393	ENDDO
394	ENDIF
395	c
396	c Lecture des temperatures du sol profond:
397	c
398	DO nsrf = 1, nbsrf
399	DO isoil=1, nsoilmx
400	IF (isoil.GT.99 .AND. nsrf.GT.99) THEN
401	PRINT*, "Trop de couches ou sous-mailles"
402	CALL abort
403	ENDIF
404	WRITE(str7,'(i2.2,"srf",i2.2)') isoil, nsrf
405	ierr = NF_INQ_VARID (nid, 'Tsoil'//str7, nvarid)
406	IF (ierr.NE.NF_NOERR) THEN
407	PRINT*, "phyetat0: Le champ <Tsoil"//str7//"> est absent"
408	PRINT*, " Il prend donc la valeur de surface"
409	DO i=1, klon
410	tsoil(i,isoil,nsrf)=tsol(i,nsrf)
411	ENDDO
412	ELSE
413	#ifdef NC_DOUBLE
414	ierr = NF_GET_VAR_DOUBLE(nid, nvarid, tsoil(1,isoil,nsrf))
415	#else
416	ierr = NF_GET_VAR_REAL(nid, nvarid, tsoil(1,isoil,nsrf))
417	#endif
418	IF (ierr.NE.NF_NOERR) THEN
419	PRINT*, "Lecture echouee pour <Tsoil"//str7//">"
420	CALL abort
421	ENDIF
422	ENDIF
423	ENDDO
424	ENDDO
425	c
426	cIM "slab" ocean
427	c
428	c Lecture de tslab (pour slab ocean seulement):
429	c
430	IF (ocean .eq. 'slab ') then
431	ierr = NF_INQ_VARID (nid, "TSLAB", nvarid)
432	IF (ierr.NE.NF_NOERR) THEN
433	PRINT*, "phyetat0: Le champ <TSLAB> est absent"
434	CALL abort
435	ENDIF
436	#ifdef NC_DOUBLE
437	ierr = NF_GET_VAR_DOUBLE(nid, nvarid, tslab)
438	#else
439	ierr = NF_GET_VAR_REAL(nid, nvarid, tslab)
440	#endif
441	IF (ierr.NE.NF_NOERR) THEN
442	PRINT*, "phyetat0: Lecture echouee pour <TSLAB>"
443	CALL abort
444	ENDIF
445	xmin = 1.0E+20
446	xmax = -1.0E+20
447	DO i = 1, klon
448	xmin = MIN(tslab(i),xmin)
449	xmax = MAX(tslab(i),xmax)
450	ENDDO
451	PRINT*,'Min, Max tslab (utilise si OCEAN=slab )', xmin, xmax
452	c
453	c Lecture de seaice (pour slab ocean seulement):
454	c
455	ierr = NF_INQ_VARID (nid, "SEAICE", nvarid)
456	IF (ierr.NE.NF_NOERR) THEN
457	PRINT*, "phyetat0: Le champ <SEAICE> est absent"
458	CALL abort
459	ENDIF
460	#ifdef NC_DOUBLE
461	ierr = NF_GET_VAR_DOUBLE(nid, nvarid, seaice)
462	#else
463	ierr = NF_GET_VAR_REAL(nid, nvarid, seaice)
464	#endif
465	IF (ierr.NE.NF_NOERR) THEN
466	PRINT*, "phyetat0: Lecture echouee pour <SEAICE>"
467	CALL abort
468	ENDIF
469	xmin = 1.0E+20
470	xmax = -1.0E+20
471	DO i = 1, klon
472	xmin = MIN(seaice(i),xmin)
473	xmax = MAX(seaice(i),xmax)
474	ENDDO
475	PRINT*,'Masse de la glace de mer (utilise si OCEAN=slab)',
476	$ xmin, xmax
477	ELSE
478	tslab = 0.
479	seaice = 0.
480	ENDIF
481	c
482	c Lecture de l'humidite de l'air juste au dessus du sol:
483	c
484	ierr = NF_INQ_VARID (nid, "QS", nvarid)
485	IF (ierr.NE.NF_NOERR) THEN
486	PRINT*, 'phyetat0: Le champ <QS> est absent'
487	PRINT, ' Mais je vais essayer de lire QS*'
488	DO nsrf = 1, nbsrf
489	IF (nsrf.GT.99) THEN
490	PRINT*, "Trop de sous-mailles"
491	CALL abort
492	ENDIF
493	WRITE(str2,'(i2.2)') nsrf
494	ierr = NF_INQ_VARID (nid, "QS"//str2, nvarid)
495	IF (ierr.NE.NF_NOERR) THEN
496	PRINT*, "phyetat0: Le champ <QS"//str2//"> est absent"
497	CALL abort
498	ENDIF
499	#ifdef NC_DOUBLE
500	ierr = NF_GET_VAR_DOUBLE(nid, nvarid, qsurf(1,nsrf))
501	#else
502	ierr = NF_GET_VAR_REAL(nid, nvarid, qsurf(1,nsrf))
503	#endif
504	IF (ierr.NE.NF_NOERR) THEN
505	PRINT*, "phyetat0: Lecture echouee pour <QS"//str2//">"
506	CALL abort
507	ENDIF
508	xmin = 1.0E+20
509	xmax = -1.0E+20
510	DO i = 1, klon
511	xmin = MIN(qsurf(i,nsrf),xmin)
512	xmax = MAX(qsurf(i,nsrf),xmax)
513	ENDDO
514	PRINT,'Humidite pres du sol QS*:', nsrf, xmin, xmax
515	ENDDO
516	ELSE
517	PRINT*, 'phyetat0: Le champ <QS> est present'
518	PRINT, ' J ignore donc les autres humidites QS*'
519	#ifdef NC_DOUBLE
520	ierr = NF_GET_VAR_DOUBLE(nid, nvarid, qsurf(1,1))
521	#else
522	ierr = NF_GET_VAR_REAL(nid, nvarid, qsurf(1,1))
523	#endif
524	IF (ierr.NE.NF_NOERR) THEN
525	PRINT*, "phyetat0: Lecture echouee pour <QS>"
526	CALL abort
527	ENDIF
528	xmin = 1.0E+20
529	xmax = -1.0E+20
530	DO i = 1, klon
531	xmin = MIN(qsurf(i,1),xmin)
532	xmax = MAX(qsurf(i,1),xmax)
533	ENDDO
534	PRINT*,'Humidite pres du sol <QS>', xmin, xmax
535	DO nsrf = 2, nbsrf
536	DO i = 1, klon
537	qsurf(i,nsrf) = qsurf(i,1)
538	ENDDO
539	ENDDO
540	ENDIF
541	C
542	C Eau dans le sol (pour le modele de sol "bucket")
543	C
544	ierr = NF_INQ_VARID (nid, "QSOL", nvarid)
545	IF (ierr .EQ. NF_NOERR) THEN
546	#ifdef NC_DOUBLE
547	ierr = NF_GET_VAR_DOUBLE(nid, nvarid, qsol)
548	#else
549	ierr = NF_GET_VAR_REAL(nid, nvarid, qsol)
550	#endif
551	IF (ierr.NE.NF_NOERR) THEN
552	PRINT*, 'phyetat0: Lecture echouee pour <QSOL>'
553	CALL abort
554	ENDIF
555	else
556	PRINT*, 'phyetat0: Le champ <QSOL> est absent'
557	PRINT*, ' Valeur par defaut nulle'
558	qsol(:)=0.
559	c$$$ CALL abort
560	ENDIF
561	xmin = 1.0E+20
562	xmax = -1.0E+20
563	DO i = 1, klon
564	xmin = MIN(qsol(i),xmin)
565	xmax = MAX(qsol(i),xmax)
566	ENDDO
567	PRINT*,'Eau dans le sol (mm) <QSOL>', xmin, xmax
568	c
569	c Lecture de neige au sol:
570	c
571	ierr = NF_INQ_VARID (nid, "SNOW", nvarid)
572	IF (ierr.NE.NF_NOERR) THEN
573	PRINT*, 'phyetat0: Le champ <SNOW> est absent'
574	PRINT, ' Mais je vais essayer de lire SNOW*'
575	DO nsrf = 1, nbsrf
576	IF (nsrf.GT.99) THEN
577	PRINT*, "Trop de sous-mailles"
578	CALL abort
579	ENDIF
580	WRITE(str2,'(i2.2)') nsrf
581	ierr = NF_INQ_VARID (nid, "SNOW"//str2, nvarid)
582	IF (ierr.NE.NF_NOERR) THEN
583	PRINT*, "phyetat0: Le champ <SNOW"//str2//"> est absent"
584	CALL abort
585	ENDIF
586	#ifdef NC_DOUBLE
587	ierr = NF_GET_VAR_DOUBLE(nid, nvarid, snow(1,nsrf))
588	#else
589	ierr = NF_GET_VAR_REAL(nid, nvarid, snow(1,nsrf))
590	#endif
591	IF (ierr.NE.NF_NOERR) THEN
592	PRINT*, "phyetat0: Lecture echouee pour <SNOW"//str2//">"
593	CALL abort
594	ENDIF
595	xmin = 1.0E+20
596	xmax = -1.0E+20
597	DO i = 1, klon
598	xmin = MIN(snow(i,nsrf),xmin)
599	xmax = MAX(snow(i,nsrf),xmax)
600	ENDDO
601	PRINT,'Neige du sol SNOW*:', nsrf, xmin, xmax
602	ENDDO
603	ELSE
604	PRINT*, 'phyetat0: Le champ <SNOW> est present'
605	PRINT, ' J ignore donc les autres neiges SNOW*'
606	#ifdef NC_DOUBLE
607	ierr = NF_GET_VAR_DOUBLE(nid, nvarid, snow(1,1))
608	#else
609	ierr = NF_GET_VAR_REAL(nid, nvarid, snow(1,1))
610	#endif
611	IF (ierr.NE.NF_NOERR) THEN
612	PRINT*, "phyetat0: Lecture echouee pour <SNOW>"
613	CALL abort
614	ENDIF
615	xmin = 1.0E+20
616	xmax = -1.0E+20
617	DO i = 1, klon
618	xmin = MIN(snow(i,1),xmin)
619	xmax = MAX(snow(i,1),xmax)
620	ENDDO
621	PRINT*,'Neige du sol <SNOW>', xmin, xmax
622	DO nsrf = 2, nbsrf
623	DO i = 1, klon
624	snow(i,nsrf) = snow(i,1)
625	ENDDO
626	ENDDO
627	ENDIF
628	c
629	c Lecture de albedo au sol:
630	c
631	ierr = NF_INQ_VARID (nid, "ALBE", nvarid)
632	IF (ierr.NE.NF_NOERR) THEN
633	PRINT*, 'phyetat0: Le champ <ALBE> est absent'
634	PRINT, ' Mais je vais essayer de lire ALBE*'
635	DO nsrf = 1, nbsrf
636	IF (nsrf.GT.99) THEN
637	PRINT*, "Trop de sous-mailles"
638	CALL abort
639	ENDIF
640	WRITE(str2,'(i2.2)') nsrf
641	ierr = NF_INQ_VARID (nid, "ALBE"//str2, nvarid)
642	IF (ierr.NE.NF_NOERR) THEN
643	PRINT*, "phyetat0: Le champ <ALBE"//str2//"> est absent"
644	CALL abort
645	ENDIF
646	#ifdef NC_DOUBLE
647	ierr = NF_GET_VAR_DOUBLE(nid, nvarid, albe(1,nsrf))
648	#else
649	ierr = NF_GET_VAR_REAL(nid, nvarid, albe(1,nsrf))
650	#endif
651	IF (ierr.NE.NF_NOERR) THEN
652	PRINT*, "phyetat0: Lecture echouee pour <ALBE"//str2//">"
653	CALL abort
654	ENDIF
655	xmin = 1.0E+20
656	xmax = -1.0E+20
657	DO i = 1, klon
658	xmin = MIN(albe(i,nsrf),xmin)
659	xmax = MAX(albe(i,nsrf),xmax)
660	ENDDO
661	PRINT,'Albedo du sol ALBE*:', nsrf, xmin, xmax
662	ENDDO
663	ELSE
664	PRINT*, 'phyetat0: Le champ <ALBE> est present'
665	PRINT, ' J ignore donc les autres ALBE*'
666	#ifdef NC_DOUBLE
667	ierr = NF_GET_VAR_DOUBLE(nid, nvarid, albe(1,1))
668	#else
669	ierr = NF_GET_VAR_REAL(nid, nvarid, albe(1,1))
670	#endif
671	IF (ierr.NE.NF_NOERR) THEN
672	PRINT*, "phyetat0: Lecture echouee pour <ALBE>"
673	CALL abort
674	ENDIF
675	xmin = 1.0E+20
676	xmax = -1.0E+20
677	DO i = 1, klon
678	xmin = MIN(albe(i,1),xmin)
679	xmax = MAX(albe(i,1),xmax)
680	ENDDO
681	PRINT*,'Neige du sol <ALBE>', xmin, xmax
682	DO nsrf = 2, nbsrf
683	DO i = 1, klon
684	albe(i,nsrf) = albe(i,1)
685	ENDDO
686	ENDDO
687	ENDIF
688
689	c
690	c Lecture de albedo au sol LW:
691	c
692	ierr = NF_INQ_VARID (nid, "ALBLW", nvarid)
693	IF (ierr.NE.NF_NOERR) THEN
694	PRINT*, 'phyetat0: Le champ <ALBLW> est absent'
695	c PRINT, ' Mais je vais essayer de lire ALBLW*'
696	PRINT, ' Mais je vais prendre ALBE*'
697	DO nsrf = 1, nbsrf
698	DO i = 1, klon
699	alblw(i,nsrf) = albe(i,nsrf)
700	ENDDO
701	ENDDO
702	ELSE
703	PRINT*, 'phyetat0: Le champ <ALBLW> est present'
704	PRINT, ' J ignore donc les autres ALBLW*'
705	#ifdef NC_DOUBLE
706	ierr = NF_GET_VAR_DOUBLE(nid, nvarid, alblw(1,1))
707	#else
708	ierr = NF_GET_VAR_REAL(nid, nvarid, alblw(1,1))
709	#endif
710	IF (ierr.NE.NF_NOERR) THEN
711	PRINT*, "phyetat0: Lecture echouee pour <ALBLW>"
712	CALL abort
713	ENDIF
714	xmin = 1.0E+20
715	xmax = -1.0E+20
716	DO i = 1, klon
717	xmin = MIN(alblw(i,1),xmin)
718	xmax = MAX(alblw(i,1),xmax)
719	ENDDO
720	PRINT*,'Neige du sol <ALBLW>', xmin, xmax
721	DO nsrf = 2, nbsrf
722	DO i = 1, klon
723	alblw(i,nsrf) = alblw(i,1)
724	ENDDO
725	ENDDO
726	ENDIF
727	c
728	c Lecture de evaporation:
729	c
730	ierr = NF_INQ_VARID (nid, "EVAP", nvarid)
731	IF (ierr.NE.NF_NOERR) THEN
732	PRINT*, 'phyetat0: Le champ <EVAP> est absent'
733	PRINT, ' Mais je vais essayer de lire EVAP*'
734	DO nsrf = 1, nbsrf
735	IF (nsrf.GT.99) THEN
736	PRINT*, "Trop de sous-mailles"
737	CALL abort
738	ENDIF
739	WRITE(str2,'(i2.2)') nsrf
740	ierr = NF_INQ_VARID (nid, "EVAP"//str2, nvarid)
741	IF (ierr.NE.NF_NOERR) THEN
742	PRINT*, "phyetat0: Le champ <EVAP"//str2//"> est absent"
743	CALL abort
744	ENDIF
745	#ifdef NC_DOUBLE
746	ierr = NF_GET_VAR_DOUBLE(nid, nvarid, evap(1,nsrf))
747	#else
748	ierr = NF_GET_VAR_REAL(nid, nvarid, evap(1,nsrf))
749	#endif
750	IF (ierr.NE.NF_NOERR) THEN
751	PRINT*, "phyetat0: Lecture echouee pour <EVAP"//str2//">"
752	CALL abort
753	ENDIF
754	xmin = 1.0E+20
755	xmax = -1.0E+20
756	DO i = 1, klon
757	xmin = MIN(evap(i,nsrf),xmin)
758	xmax = MAX(evap(i,nsrf),xmax)
759	ENDDO
760	PRINT,'evap du sol EVAP*:', nsrf, xmin, xmax
761	ENDDO
762	ELSE
763	PRINT*, 'phyetat0: Le champ <EVAP> est present'
764	PRINT, ' J ignore donc les autres EVAP*'
765	#ifdef NC_DOUBLE
766	ierr = NF_GET_VAR_DOUBLE(nid, nvarid, evap(1,1))
767	#else
768	ierr = NF_GET_VAR_REAL(nid, nvarid, evap(1,1))
769	#endif
770	IF (ierr.NE.NF_NOERR) THEN
771	PRINT*, "phyetat0: Lecture echouee pour <EVAP>"
772	CALL abort
773	ENDIF
774	xmin = 1.0E+20
775	xmax = -1.0E+20
776	DO i = 1, klon
777	xmin = MIN(evap(i,1),xmin)
778	xmax = MAX(evap(i,1),xmax)
779	ENDDO
780	PRINT*,'Evap du sol <EVAP>', xmin, xmax
781	DO nsrf = 2, nbsrf
782	DO i = 1, klon
783	evap(i,nsrf) = evap(i,1)
784	ENDDO
785	ENDDO
786	ENDIF
787	c
788	c Lecture precipitation liquide:
789	c
790	ierr = NF_INQ_VARID (nid, "rain_f", nvarid)
791	IF (ierr.NE.NF_NOERR) THEN
792	PRINT*, 'phyetat0: Le champ <rain_f> est absent'
793	CALL abort
794	ENDIF
795	#ifdef NC_DOUBLE
796	ierr = NF_GET_VAR_DOUBLE(nid, nvarid, rain_fall)
797	#else
798	ierr = NF_GET_VAR_REAL(nid, nvarid, rain_fall)
799	#endif
800	IF (ierr.NE.NF_NOERR) THEN
801	PRINT*, 'phyetat0: Lecture echouee pour <rain_f>'
802	CALL abort
803	ENDIF
804	xmin = 1.0E+20
805	xmax = -1.0E+20
806	DO i = 1, klon
807	xmin = MIN(rain_fall(i),xmin)
808	xmax = MAX(rain_fall(i),xmax)
809	ENDDO
810	PRINT*,'Precipitation liquide rain_f:', xmin, xmax
811	c
812	c Lecture precipitation solide:
813	c
814	ierr = NF_INQ_VARID (nid, "snow_f", nvarid)
815	IF (ierr.NE.NF_NOERR) THEN
816	PRINT*, 'phyetat0: Le champ <snow_f> est absent'
817	CALL abort
818	ENDIF
819	#ifdef NC_DOUBLE
820	ierr = NF_GET_VAR_DOUBLE(nid, nvarid, snow_fall)
821	#else
822	ierr = NF_GET_VAR_REAL(nid, nvarid, snow_fall)
823	#endif
824	IF (ierr.NE.NF_NOERR) THEN
825	PRINT*, 'phyetat0: Lecture echouee pour <snow_f>'
826	CALL abort
827	ENDIF
828	xmin = 1.0E+20
829	xmax = -1.0E+20
830	DO i = 1, klon
831	xmin = MIN(snow_fall(i),xmin)
832	xmax = MAX(snow_fall(i),xmax)
833	ENDDO
834	PRINT*,'Precipitation solide snow_f:', xmin, xmax
835	c
836	c Lecture rayonnement solaire au sol:
837	c
838	ierr = NF_INQ_VARID (nid, "solsw", nvarid)
839	IF (ierr.NE.NF_NOERR) THEN
840	PRINT*, 'phyetat0: Le champ <solsw> est absent'
841	PRINT*, 'mis a zero'
842	solsw = 0.
843	ELSE
844	#ifdef NC_DOUBLE
845	ierr = NF_GET_VAR_DOUBLE(nid, nvarid, solsw)
846	#else
847	ierr = NF_GET_VAR_REAL(nid, nvarid, solsw)
848	#endif
849	IF (ierr.NE.NF_NOERR) THEN
850	PRINT*, 'phyetat0: Lecture echouee pour <solsw>'
851	CALL abort
852	ENDIF
853	ENDIF
854	xmin = 1.0E+20
855	xmax = -1.0E+20
856	DO i = 1, klon
857	xmin = MIN(solsw(i),xmin)
858	xmax = MAX(solsw(i),xmax)
859	ENDDO
860	PRINT*,'Rayonnement solaire au sol solsw:', xmin, xmax
861	c
862	c Lecture rayonnement IF au sol:
863	c
864	ierr = NF_INQ_VARID (nid, "sollw", nvarid)
865	IF (ierr.NE.NF_NOERR) THEN
866	PRINT*, 'phyetat0: Le champ <sollw> est absent'
867	PRINT*, 'mis a zero'
868	sollw = 0.
869	ELSE
870	#ifdef NC_DOUBLE
871	ierr = NF_GET_VAR_DOUBLE(nid, nvarid, sollw)
872	#else
873	ierr = NF_GET_VAR_REAL(nid, nvarid, sollw)
874	#endif
875	IF (ierr.NE.NF_NOERR) THEN
876	PRINT*, 'phyetat0: Lecture echouee pour <sollw>'
877	CALL abort
878	ENDIF
879	ENDIF
880	xmin = 1.0E+20
881	xmax = -1.0E+20
882	DO i = 1, klon
883	xmin = MIN(sollw(i),xmin)
884	xmax = MAX(sollw(i),xmax)
885	ENDDO
886	PRINT*,'Rayonnement IF au sol sollw:', xmin, xmax
887
888	c
889	c Lecture derive des flux:
890	c
891	ierr = NF_INQ_VARID (nid, "fder", nvarid)
892	IF (ierr.NE.NF_NOERR) THEN
893	PRINT*, 'phyetat0: Le champ <fder> est absent'
894	PRINT*, 'mis a zero'
895	fder = 0.
896	ELSE
897	#ifdef NC_DOUBLE
898	ierr = NF_GET_VAR_DOUBLE(nid, nvarid, fder)
899	#else
900	ierr = NF_GET_VAR_REAL(nid, nvarid, fder)
901	#endif
902	IF (ierr.NE.NF_NOERR) THEN
903	PRINT*, 'phyetat0: Lecture echouee pour <fder>'
904	CALL abort
905	ENDIF
906	ENDIF
907	xmin = 1.0E+20
908	xmax = -1.0E+20
909	DO i = 1, klon
910	xmin = MIN(fder(i),xmin)
911	xmax = MAX(fder(i),xmax)
912	ENDDO
913	PRINT*,'Derive des flux fder:', xmin, xmax
914
915	c
916	c Lecture du rayonnement net au sol:
917	c
918	ierr = NF_INQ_VARID (nid, "RADS", nvarid)
919	IF (ierr.NE.NF_NOERR) THEN
920	PRINT*, 'phyetat0: Le champ <RADS> est absent'
921	CALL abort
922	ENDIF
923	#ifdef NC_DOUBLE
924	ierr = NF_GET_VAR_DOUBLE(nid, nvarid, radsol)
925	#else
926	ierr = NF_GET_VAR_REAL(nid, nvarid, radsol)
927	#endif
928	IF (ierr.NE.NF_NOERR) THEN
929	PRINT*, 'phyetat0: Lecture echouee pour <RADS>'
930	CALL abort
931	ENDIF
932	xmin = 1.0E+20
933	xmax = -1.0E+20
934	DO i = 1, klon
935	xmin = MIN(radsol(i),xmin)
936	xmax = MAX(radsol(i),xmax)
937	ENDDO
938	PRINT*,'Rayonnement net au sol radsol:', xmin, xmax
939	c
940	c Lecture de la longueur de rugosite
941	c
942	c
943	ierr = NF_INQ_VARID (nid, "RUG", nvarid)
944	IF (ierr.NE.NF_NOERR) THEN
945	PRINT*, 'phyetat0: Le champ <RUG> est absent'
946	PRINT, ' Mais je vais essayer de lire RUG*'
947	DO nsrf = 1, nbsrf
948	IF (nsrf.GT.99) THEN
949	PRINT*, "Trop de sous-mailles"
950	CALL abort
951	ENDIF
952	WRITE(str2,'(i2.2)') nsrf
953	ierr = NF_INQ_VARID (nid, "RUG"//str2, nvarid)
954	IF (ierr.NE.NF_NOERR) THEN
955	PRINT*, "phyetat0: Le champ <RUG"//str2//"> est absent"
956	CALL abort
957	ENDIF
958	#ifdef NC_DOUBLE
959	ierr = NF_GET_VAR_DOUBLE(nid, nvarid, frugs(1,nsrf))
960	#else
961	ierr = NF_GET_VAR_REAL(nid, nvarid, frugs(1,nsrf))
962	#endif
963	IF (ierr.NE.NF_NOERR) THEN
964	PRINT*, "phyetat0: Lecture echouee pour <RUG"//str2//">"
965	CALL abort
966	ENDIF
967	xmin = 1.0E+20
968	xmax = -1.0E+20
969	DO i = 1, klon
970	xmin = MIN(frugs(i,nsrf),xmin)
971	xmax = MAX(frugs(i,nsrf),xmax)
972	ENDDO
973	PRINT,'rugosite du sol RUG*:', nsrf, xmin, xmax
974	ENDDO
975	ELSE
976	PRINT*, 'phyetat0: Le champ <RUG> est present'
977	PRINT, ' J ignore donc les autres RUG*'
978	#ifdef NC_DOUBLE
979	ierr = NF_GET_VAR_DOUBLE(nid, nvarid, frugs(1,1))
980	#else
981	ierr = NF_GET_VAR_REAL(nid, nvarid, frugs(1,1))
982	#endif
983	IF (ierr.NE.NF_NOERR) THEN
984	PRINT*, "phyetat0: Lecture echouee pour <RUG>"
985	CALL abort
986	ENDIF
987	xmin = 1.0E+20
988	xmax = -1.0E+20
989	DO i = 1, klon
990	xmin = MIN(frugs(i,1),xmin)
991	xmax = MAX(frugs(i,1),xmax)
992	ENDDO
993	PRINT*,'rugosite <RUG>', xmin, xmax
994	DO nsrf = 2, nbsrf
995	DO i = 1, klon
996	frugs(i,nsrf) = frugs(i,1)
997	ENDDO
998	ENDDO
999	ENDIF
1000
1001	c
1002	c Lecture de l'age de la neige:
1003	c
1004	ierr = NF_INQ_VARID (nid, "AGESNO", nvarid)
1005	IF (ierr.NE.NF_NOERR) THEN
1006	PRINT*, 'phyetat0: Le champ <AGESNO> est absent'
1007	PRINT, ' Mais je vais essayer de lire AGESNO*'
1008	DO nsrf = 1, nbsrf
1009	IF (nsrf.GT.99) THEN
1010	PRINT*, "Trop de sous-mailles"
1011	CALL abort
1012	ENDIF
1013	WRITE(str2,'(i2.2)') nsrf
1014	ierr = NF_INQ_VARID (nid, "AGESNO"//str2, nvarid)
1015	IF (ierr.NE.NF_NOERR) THEN
1016	PRINT*, "phyetat0: Le champ <AGESNO"//str2//"> est absent"
1017	agesno = 50.0
1018	ENDIF
1019	#ifdef NC_DOUBLE
1020	ierr = NF_GET_VAR_DOUBLE(nid, nvarid, agesno(1,nsrf))
1021	#else
1022	ierr = NF_GET_VAR_REAL(nid, nvarid, agesno(1,nsrf))
1023	#endif
1024	IF (ierr.NE.NF_NOERR) THEN
1025	PRINT*, "phyetat0: Lecture echouee pour <AGESNO"//str2//">"
1026	CALL abort
1027	ENDIF
1028	xmin = 1.0E+20
1029	xmax = -1.0E+20
1030	DO i = 1, klon
1031	xmin = MIN(agesno(i,nsrf),xmin)
1032	xmax = MAX(agesno(i,nsrf),xmax)
1033	ENDDO
1034	PRINT,'Age de la neige AGESNO*:', nsrf, xmin, xmax
1035	ENDDO
1036	ELSE
1037	PRINT*, 'phyetat0: Le champ <AGESNO> est present'
1038	PRINT, ' J ignore donc les autres AGESNO*'
1039	#ifdef NC_DOUBLE
1040	ierr = NF_GET_VAR_DOUBLE(nid, nvarid, agesno(1,1))
1041	#else
1042	ierr = NF_GET_VAR_REAL(nid, nvarid, agesno(1,1))
1043	#endif
1044	IF (ierr.NE.NF_NOERR) THEN
1045	PRINT*, "phyetat0: Lecture echouee pour <AGESNO>"
1046	CALL abort
1047	ENDIF
1048	xmin = 1.0E+20
1049	xmax = -1.0E+20
1050	DO i = 1, klon
1051	xmin = MIN(agesno(i,1),xmin)
1052	xmax = MAX(agesno(i,1),xmax)
1053	ENDDO
1054	PRINT*,'Age de la neige <AGESNO>', xmin, xmax
1055	DO nsrf = 2, nbsrf
1056	DO i = 1, klon
1057	agesno(i,nsrf) = agesno(i,1)
1058	ENDDO
1059	ENDDO
1060	ENDIF
1061
1062	c
1063	ierr = NF_INQ_VARID (nid, "ZMEA", nvarid)
1064	IF (ierr.NE.NF_NOERR) THEN
1065	PRINT*, 'phyetat0: Le champ <ZMEA> est absent'
1066	CALL abort
1067	ENDIF
1068	#ifdef NC_DOUBLE
1069	ierr = NF_GET_VAR_DOUBLE(nid, nvarid, zmea)
1070	#else
1071	ierr = NF_GET_VAR_REAL(nid, nvarid, zmea)
1072	#endif
1073	IF (ierr.NE.NF_NOERR) THEN
1074	PRINT*, 'phyetat0: Lecture echouee pour <ZMEA>'
1075	CALL abort
1076	ENDIF
1077	xmin = 1.0E+20
1078	xmax = -1.0E+20
1079	DO i = 1, klon
1080	xmin = MIN(zmea(i),xmin)
1081	xmax = MAX(zmea(i),xmax)
1082	ENDDO
1083	PRINT*,'OROGRAPHIE SOUS-MAILLE zmea:', xmin, xmax
1084	c
1085	c
1086	ierr = NF_INQ_VARID (nid, "ZSTD", nvarid)
1087	IF (ierr.NE.NF_NOERR) THEN
1088	PRINT*, 'phyetat0: Le champ <ZSTD> est absent'
1089	CALL abort
1090	ENDIF
1091	#ifdef NC_DOUBLE
1092	ierr = NF_GET_VAR_DOUBLE(nid, nvarid, zstd)
1093	#else
1094	ierr = NF_GET_VAR_REAL(nid, nvarid, zstd)
1095	#endif
1096	IF (ierr.NE.NF_NOERR) THEN
1097	PRINT*, 'phyetat0: Lecture echouee pour <ZSTD>'
1098	CALL abort
1099	ENDIF
1100	xmin = 1.0E+20
1101	xmax = -1.0E+20
1102	DO i = 1, klon
1103	xmin = MIN(zstd(i),xmin)
1104	xmax = MAX(zstd(i),xmax)
1105	ENDDO
1106	PRINT*,'OROGRAPHIE SOUS-MAILLE zstd:', xmin, xmax
1107	c
1108	c
1109	ierr = NF_INQ_VARID (nid, "ZSIG", nvarid)
1110	IF (ierr.NE.NF_NOERR) THEN
1111	PRINT*, 'phyetat0: Le champ <ZSIG> est absent'
1112	CALL abort
1113	ENDIF
1114	#ifdef NC_DOUBLE
1115	ierr = NF_GET_VAR_DOUBLE(nid, nvarid, zsig)
1116	#else
1117	ierr = NF_GET_VAR_REAL(nid, nvarid, zsig)
1118	#endif
1119	IF (ierr.NE.NF_NOERR) THEN
1120	PRINT*, 'phyetat0: Lecture echouee pour <ZSIG>'
1121	CALL abort
1122	ENDIF
1123	xmin = 1.0E+20
1124	xmax = -1.0E+20
1125	DO i = 1, klon
1126	xmin = MIN(zsig(i),xmin)
1127	xmax = MAX(zsig(i),xmax)
1128	ENDDO
1129	PRINT*,'OROGRAPHIE SOUS-MAILLE zsig:', xmin, xmax
1130	c
1131	c
1132	ierr = NF_INQ_VARID (nid, "ZGAM", nvarid)
1133	IF (ierr.NE.NF_NOERR) THEN
1134	PRINT*, 'phyetat0: Le champ <ZGAM> est absent'
1135	CALL abort
1136	ENDIF
1137	#ifdef NC_DOUBLE
1138	ierr = NF_GET_VAR_DOUBLE(nid, nvarid, zgam)
1139	#else
1140	ierr = NF_GET_VAR_REAL(nid, nvarid, zgam)
1141	#endif
1142	IF (ierr.NE.NF_NOERR) THEN
1143	PRINT*, 'phyetat0: Lecture echouee pour <ZGAM>'
1144	CALL abort
1145	ENDIF
1146	xmin = 1.0E+20
1147	xmax = -1.0E+20
1148	DO i = 1, klon
1149	xmin = MIN(zgam(i),xmin)
1150	xmax = MAX(zgam(i),xmax)
1151	ENDDO
1152	PRINT*,'OROGRAPHIE SOUS-MAILLE zgam:', xmin, xmax
1153	c
1154	c
1155	ierr = NF_INQ_VARID (nid, "ZTHE", nvarid)
1156	IF (ierr.NE.NF_NOERR) THEN
1157	PRINT*, 'phyetat0: Le champ <ZTHE> est absent'
1158	CALL abort
1159	ENDIF
1160	#ifdef NC_DOUBLE
1161	ierr = NF_GET_VAR_DOUBLE(nid, nvarid, zthe)
1162	#else
1163	ierr = NF_GET_VAR_REAL(nid, nvarid, zthe)
1164	#endif
1165	IF (ierr.NE.NF_NOERR) THEN
1166	PRINT*, 'phyetat0: Lecture echouee pour <ZTHE>'
1167	CALL abort
1168	ENDIF
1169	xmin = 1.0E+20
1170	xmax = -1.0E+20
1171	DO i = 1, klon
1172	xmin = MIN(zthe(i),xmin)
1173	xmax = MAX(zthe(i),xmax)
1174	ENDDO
1175	PRINT*,'OROGRAPHIE SOUS-MAILLE zthe:', xmin, xmax
1176	c
1177	c
1178	ierr = NF_INQ_VARID (nid, "ZPIC", nvarid)
1179	IF (ierr.NE.NF_NOERR) THEN
1180	PRINT*, 'phyetat0: Le champ <ZPIC> est absent'
1181	CALL abort
1182	ENDIF
1183	#ifdef NC_DOUBLE
1184	ierr = NF_GET_VAR_DOUBLE(nid, nvarid, zpic)
1185	#else
1186	ierr = NF_GET_VAR_REAL(nid, nvarid, zpic)
1187	#endif
1188	IF (ierr.NE.NF_NOERR) THEN
1189	PRINT*, 'phyetat0: Lecture echouee pour <ZPIC>'
1190	CALL abort
1191	ENDIF
1192	xmin = 1.0E+20
1193	xmax = -1.0E+20
1194	DO i = 1, klon
1195	xmin = MIN(zpic(i),xmin)
1196	xmax = MAX(zpic(i),xmax)
1197	ENDDO
1198	PRINT*,'OROGRAPHIE SOUS-MAILLE zpic:', xmin, xmax
1199	c
1200	ierr = NF_INQ_VARID (nid, "ZVAL", nvarid)
1201	IF (ierr.NE.NF_NOERR) THEN
1202	PRINT*, 'phyetat0: Le champ <ZVAL> est absent'
1203	CALL abort
1204	ENDIF
1205	#ifdef NC_DOUBLE
1206	ierr = NF_GET_VAR_DOUBLE(nid, nvarid, zval)
1207	#else
1208	ierr = NF_GET_VAR_REAL(nid, nvarid, zval)
1209	#endif
1210	IF (ierr.NE.NF_NOERR) THEN
1211	PRINT*, 'phyetat0: Lecture echouee pour <ZVAL>'
1212	CALL abort
1213	ENDIF
1214	xmin = 1.0E+20
1215	xmax = -1.0E+20
1216	DO i = 1, klon
1217	xmin = MIN(zval(i),xmin)
1218	xmax = MAX(zval(i),xmax)
1219	ENDDO
1220	PRINT*,'OROGRAPHIE SOUS-MAILLE zval:', xmin, xmax
1221	c
1222	c
1223	ierr = NF_INQ_VARID (nid, "RUGSREL", nvarid)
1224	IF (ierr.NE.NF_NOERR) THEN
1225	PRINT*, 'phyetat0: Le champ <RUGSREL> est absent'
1226	CALL abort
1227	ENDIF
1228	#ifdef NC_DOUBLE
1229	ierr = NF_GET_VAR_DOUBLE(nid, nvarid, rugsrel)
1230	#else
1231	ierr = NF_GET_VAR_REAL(nid, nvarid, rugsrel)
1232	#endif
1233	IF (ierr.NE.NF_NOERR) THEN
1234	PRINT*, 'phyetat0: Lecture echouee pour <RUGSREL>'
1235	CALL abort
1236	ENDIF
1237	xmin = 1.0E+20
1238	xmax = -1.0E+20
1239	DO i = 1, klon
1240	xmin = MIN(rugsrel(i),xmin)
1241	xmax = MAX(rugsrel(i),xmax)
1242	ENDDO
1243	PRINT*,'Rugosite relief (ecart-type) rugsrel:', xmin, xmax
1244	c
1245	c
1246	ancien_ok = .TRUE.
1247	c
1248	ierr = NF_INQ_VARID (nid, "TANCIEN", nvarid)
1249	IF (ierr.NE.NF_NOERR) THEN
1250	PRINT*, "phyetat0: Le champ <TANCIEN> est absent"
1251	PRINT*, "Depart legerement fausse. Mais je continue"
1252	ancien_ok = .FALSE.
1253	ELSE
1254	#ifdef NC_DOUBLE
1255	ierr = NF_GET_VAR_DOUBLE(nid, nvarid, t_ancien)
1256	#else
1257	ierr = NF_GET_VAR_REAL(nid, nvarid, t_ancien)
1258	#endif
1259	IF (ierr.NE.NF_NOERR) THEN
1260	PRINT*, "phyetat0: Lecture echouee pour <TANCIEN>"
1261	CALL abort
1262	ENDIF
1263	ENDIF
1264	c
1265	ierr = NF_INQ_VARID (nid, "QANCIEN", nvarid)
1266	IF (ierr.NE.NF_NOERR) THEN
1267	PRINT*, "phyetat0: Le champ <QANCIEN> est absent"
1268	PRINT*, "Depart legerement fausse. Mais je continue"
1269	ancien_ok = .FALSE.
1270	ELSE
1271	#ifdef NC_DOUBLE
1272	ierr = NF_GET_VAR_DOUBLE(nid, nvarid, q_ancien)
1273	#else
1274	ierr = NF_GET_VAR_REAL(nid, nvarid, q_ancien)
1275	#endif
1276	IF (ierr.NE.NF_NOERR) THEN
1277	PRINT*, "phyetat0: Lecture echouee pour <QANCIEN>"
1278	CALL abort
1279	ENDIF
1280	ENDIF
1281	c
1282	ierr = NF_INQ_VARID (nid, "CLWCON", nvarid)
1283	IF (ierr.NE.NF_NOERR) THEN
1284	PRINT*, "phyetat0: Le champ CLWCON est absent"
1285	PRINT*, "Depart legerement fausse. Mais je continue"
1286	clwcon = 0.
1287	ELSE
1288	#ifdef NC_DOUBLE
1289	ierr = NF_GET_VAR_DOUBLE(nid, nvarid, clwcon)
1290	#else
1291	ierr = NF_GET_VAR_REAL(nid, nvarid, clwcon)
1292	#endif
1293	IF (ierr.NE.NF_NOERR) THEN
1294	PRINT*, "phyetat0: Lecture echouee pour <CLWCON>"
1295	CALL abort
1296	ENDIF
1297	ENDIF
1298	xmin = 1.0E+20
1299	xmax = -1.0E+20
1300	xmin = MINval(clwcon)
1301	xmax = MAXval(clwcon)
1302	PRINT*,'Eau liquide convective (ecart-type) clwcon:', xmin, xmax
1303	c
1304	ierr = NF_INQ_VARID (nid, "RNEBCON", nvarid)
1305	IF (ierr.NE.NF_NOERR) THEN
1306	PRINT*, "phyetat0: Le champ RNEBCON est absent"
1307	PRINT*, "Depart legerement fausse. Mais je continue"
1308	rnebcon = 0.
1309	ELSE
1310	#ifdef NC_DOUBLE
1311	ierr = NF_GET_VAR_DOUBLE(nid, nvarid, rnebcon)
1312	#else
1313	ierr = NF_GET_VAR_REAL(nid, nvarid, rnebcon)
1314	#endif
1315	IF (ierr.NE.NF_NOERR) THEN
1316	PRINT*, "phyetat0: Lecture echouee pour <RNEBCON>"
1317	CALL abort
1318	ENDIF
1319	ENDIF
1320	xmin = 1.0E+20
1321	xmax = -1.0E+20
1322	xmin = MINval(rnebcon)
1323	xmax = MAXval(rnebcon)
1324	PRINT*,'Nebulosite convective (ecart-type) rnebcon:', xmin, xmax
1325
1326	c
1327	ierr = NF_INQ_VARID (nid, "QANCIEN", nvarid)
1328	IF (ierr.NE.NF_NOERR) THEN
1329	PRINT*, "phyetat0: Le champ <QANCIEN> est absent"
1330	PRINT*, "Depart legerement fausse. Mais je continue"
1331	ancien_ok = .FALSE.
1332	ELSE
1333	#ifdef NC_DOUBLE
1334	ierr = NF_GET_VAR_DOUBLE(nid, nvarid, q_ancien)
1335	#else
1336	ierr = NF_GET_VAR_REAL(nid, nvarid, q_ancien)
1337	#endif
1338	IF (ierr.NE.NF_NOERR) THEN
1339	PRINT*, "phyetat0: Lecture echouee pour <QANCIEN>"
1340	CALL abort
1341	ENDIF
1342	ENDIF
1343	c
1344	c Lecture ratqs
1345	c
1346	ierr = NF_INQ_VARID (nid, "RATQS", nvarid)
1347	IF (ierr.NE.NF_NOERR) THEN
1348	PRINT*, "phyetat0: Le champ <RATQS> est absent"
1349	PRINT*, "Depart legerement fausse. Mais je continue"
1350	ratqs = 0.
1351	ELSE
1352	#ifdef NC_DOUBLE
1353	ierr = NF_GET_VAR_DOUBLE(nid, nvarid, ratqs)
1354	#else
1355	ierr = NF_GET_VAR_REAL(nid, nvarid, ratqs)
1356	#endif
1357	IF (ierr.NE.NF_NOERR) THEN
1358	PRINT*, "phyetat0: Lecture echouee pour <RATQS>"
1359	CALL abort
1360	ENDIF
1361	ENDIF
1362	xmin = 1.0E+20
1363	xmax = -1.0E+20
1364	xmin = MINval(ratqs)
1365	xmax = MAXval(ratqs)
1366	PRINT*,'(ecart-type) ratqs:', xmin, xmax
1367	c
1368	c Lecture run_off_lic_0
1369	c
1370	ierr = NF_INQ_VARID (nid, "RUNOFFLIC0", nvarid)
1371	IF (ierr.NE.NF_NOERR) THEN
1372	PRINT*, "phyetat0: Le champ <RUNOFFLIC0> est absent"
1373	PRINT*, "Depart legerement fausse. Mais je continue"
1374	run_off_lic_0 = 0.
1375	ELSE
1376	#ifdef NC_DOUBLE
1377	ierr = NF_GET_VAR_DOUBLE(nid, nvarid, run_off_lic_0)
1378	#else
1379	ierr = NF_GET_VAR_REAL(nid, nvarid, run_off_lic_0)
1380	#endif
1381	IF (ierr.NE.NF_NOERR) THEN
1382	PRINT*, "phyetat0: Lecture echouee pour <RUNOFFLIC0>"
1383	CALL abort
1384	ENDIF
1385	ENDIF
1386	xmin = 1.0E+20
1387	xmax = -1.0E+20
1388	xmin = MINval(run_off_lic_0)
1389	xmax = MAXval(run_off_lic_0)
1390	PRINT*,'(ecart-type) run_off_lic_0:', xmin, xmax
1391	c
1392	c Fermer le fichier:
1393	c
1394	ierr = NF_CLOSE(nid)
1395	c
1396	RETURN
1397	END

Note: See TracBrowser for help on using the repository browser.

Download in other formats:

Original Format