Context Navigation

phyetat0.F @ 5451

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

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

Download in other formats:

Original Format