Context Navigation

phyetat0.F @ 977

Last change on this file since 977 was 977, checked in by Laurent Fairhead, 16 years ago
Erreur consecutive au renommage des variables LF
Property svn:eol-style set to `native` Property svn:keywords set to `Author Date Id Revision`
File size: 53.9 KB

Line
1	!
2	! $Header$
3	!
4	c
5	c
6	SUBROUTINE phyetat0 (fichnom,
7	. ocean_in, ok_veget_in,
8	. clesphy0,
9	. tabcntr0)
10
11	USE dimphy
12	USE mod_grid_phy_lmdz
13	USE mod_phys_lmdz_para
14	USE iophy
15	USE ocean_slab_mod, ONLY : ocean_slab_init
16	USE ocean_cpl_mod, ONLY : ocean_cpl_init
17	USE ocean_forced_mod, ONLY : ocean_forced_init
18	USE fonte_neige_mod, ONLY : fonte_neige_init
19	USE pbl_surface_mod, ONLY : pbl_surface_init
20	USE surface_data, ONLY : ocean, ok_veget
21	USE phys_state_var_mod
22
23	IMPLICIT none
24	c======================================================================
25	c Auteur(s) Z.X. Li (LMD/CNRS) date: 19930818
26	c Objet: Lecture de l'etat initial pour la physique
27	c======================================================================
28	#include "dimensions.h"
29	#include "netcdf.inc"
30	#include "indicesol.h"
31	#include "dimsoil.h"
32	#include "clesphys.h"
33	#include "temps.h"
34	#include "thermcell.h"
35	#include "compbl.h"
36	c======================================================================
37	CHARACTER() fichnom
38
39	c les variables globales lues dans le fichier restart
40	REAL rlat_glo(klon_glo), rlon_glo(klon_glo)
41	REAL pctsrf_glo(klon_glo, nbsrf)
42	REAL tsol_glo(klon_glo,nbsrf)
43	REAL alb1_glo(klon_glo,nbsrf)
44	REAL alb2_glo(klon_glo,nbsrf)
45	REAL rain_fall_glo(klon_glo)
46	REAL snow_fall_glo(klon_glo)
47	real solsw_glo(klon_glo)
48	real sollw_glo(klon_glo)
49	REAL radsol_glo(klon_glo)
50	REAL zmea_glo(klon_glo)
51	REAL zstd_glo(klon_glo)
52	REAL zsig_glo(klon_glo)
53	REAL zgam_glo(klon_glo)
54	REAL zthe_glo(klon_glo)
55	REAL zpic_glo(klon_glo)
56	REAL zval_glo(klon_glo)
57	REAL rugsrel_glo(klon_glo)
58	REAL t_ancien_glo(klon_glo,klev), q_ancien_glo(klon_glo,klev)
59	REAL clwcon_glo(klon_glo,klev)
60	REAL rnebcon_glo(klon_glo,klev)
61	REAL ratqs_glo(klon_glo,klev)
62	REAL pbl_tke_glo(klon_glo,klev+1,nbsrf)
63	REAL zmax0_glo(klon_glo), f0_glo(klon)
64	REAL ema_work1_glo(klon_glo, klev), ema_work2_glo(klon_glo, klev)
65	REAL wake_deltat_glo(klon,klev), wake_deltaq_glo(klon,klev)
66	REAL wake_s_glo(klon), wake_cstar_glo(klon), wake_fip_glo(klon)
67	REAL tsoil_p(klon,nsoilmx,nbsrf)
68	REAL tslab_p(klon), seaice_p(klon)
69	REAL qsurf_p(klon,nbsrf)
70	REAL qsol_p(klon)
71	REAL snow_p(klon,nbsrf)
72	REAL evap_p(klon,nbsrf)
73	real fder_p(klon)
74	REAL frugs_p(klon,nbsrf)
75	REAL agesno_p(klon,nbsrf)
76	REAL run_off_lic_0_p(klon)
77
78	LOGICAL,SAVE :: ancien_ok_glo
79	!$OMP THREADPRIVATE(ancien_ok_glo)
80
81	REAL zmasq_glo(klon_glo)
82	REAL tsoil(klon_glo,nsoilmx,nbsrf)
83	cIM "slab" ocean
84	REAL tslab(klon_glo), seaice(klon_glo)
85	REAL qsurf(klon_glo,nbsrf)
86	REAL qsol(klon_glo)
87	REAL snow(klon_glo,nbsrf)
88	REAL evap(klon_glo,nbsrf)
89	real fder(klon_glo)
90	REAL frugs(klon_glo,nbsrf)
91	REAL agesno(klon_glo,nbsrf)
92	REAL fractint(klon_glo)
93	REAL run_off_lic_0(klon_glo)
94
95	CHARACTER*6 ocean_in
96	LOGICAL ok_veget_in
97
98	INTEGER longcles
99	PARAMETER ( longcles = 20 )
100	REAL clesphy0( longcles )
101	c
102	REAL xmin, xmax
103	c
104	INTEGER nid, nvarid
105	INTEGER ierr, i, nsrf, isoil ,k
106	INTEGER length
107	PARAMETER (length=100)
108	REAL tab_cntrl(length), tabcntr0(length)
109	REAL,SAVE :: tab_cntrl_omp(length)
110	CHARACTER*7 str7
111	CHARACTER*2 str2
112
113	c FH1D
114	c real iolat(jjm+1)
115	real iolat(jjm+1-1/iim)
116	c
117	c Ouvrir le fichier contenant l'etat initial:
118	c
119
120	c$OMP MASTER
121	print *,'MASTER -x , omp_rank=',omp_rank
122	c$OMP END MASTER
123
124	c$OMP MASTER
125	IF (is_mpi_root) THEN
126	print*,'fichnom ',fichnom
127	ierr = NF_OPEN (fichnom, NF_NOWRITE,nid)
128	IF (ierr.NE.NF_NOERR) THEN
129	write(6,*)' Pb d''ouverture du fichier '//fichnom
130	write(6,*)' ierr = ', ierr
131	CALL ABORT
132	ENDIF
133	ENDIF
134	c$OMP END MASTER
135	c
136	c Lecture des parametres de controle:
137	c
138	c$OMP MASTER
139	IF (is_mpi_root) THEN
140
141	ierr = NF_INQ_VARID (nid, "controle", nvarid)
142	IF (ierr.NE.NF_NOERR) THEN
143	PRINT*, 'phyetat0: Le champ <controle> est absent'
144	CALL abort
145	ENDIF
146	#ifdef NC_DOUBLE
147	ierr = NF_GET_VAR_DOUBLE(nid, nvarid, tab_cntrl)
148	#else
149	ierr = NF_GET_VAR_REAL(nid, nvarid, tab_cntrl)
150	#endif
151	IF (ierr.NE.NF_NOERR) THEN
152	PRINT*, 'phyetat0: Lecture echouee pour <controle>'
153	CALL abort
154	ENDIF
155	ENDIF
156
157	c$OMP END MASTER
158
159	CALL bcast(tab_cntrl)
160
161	c
162	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
163	! FH 2008/05/09 On elimine toutes les clefs physiques dans la dynamique
164	! Les constantes de la physiques sont lues dans la physique seulement.
165	! Les egalites du type
166	! tab_cntrl( 5 )=clesphy0(1)
167	! sont remplacees par
168	! clesphy0(1)=tab_cntrl( 5 )
169	! On inverse aussi la logique.
170	! On remplit les tab_cntrl avec les parametres lus dans les .def
171	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
172
173	DO i = 1, length
174	tabcntr0( i ) = tab_cntrl( i )
175	ENDDO
176	c
177	tab_cntrl(1)=dtime
178	tab_cntrl(2)=radpas
179	co2_ppm_etat0 = tab_cntrl(3)
180	solaire_etat0 = tab_cntrl(4)
181	tab_cntrl(5)=iflag_con
182	tab_cntrl(6)=nbapp_rad
183
184	if (cycle_diurne) tab_cntrl( 7) =1.
185	if (soil_model) tab_cntrl( 8) =1.
186	if (new_oliq) tab_cntrl( 9) =1.
187	if (ok_orodr) tab_cntrl(10) =1.
188	if (ok_orolf) tab_cntrl(11) =1.
189	if (ok_limitvrai) tab_cntrl(12) =1.
190
191
192	itau_phy = tab_cntrl(15)
193
194
195
196	IF( clesphy0(1).NE.tab_cntrl( 5 ) ) THEN
197	clesphy0(1)=tab_cntrl( 5 )
198	ENDIF
199
200	IF( clesphy0(2).NE.tab_cntrl( 6 ) ) THEN
201	clesphy0(2)=tab_cntrl( 6 )
202	ENDIF
203
204	IF( clesphy0(3).NE.tab_cntrl( 7 ) ) THEN
205	clesphy0(3)=tab_cntrl( 7 )
206	ENDIF
207
208	IF( clesphy0(4).NE.tab_cntrl( 8 ) ) THEN
209	clesphy0(4)=tab_cntrl( 8 )
210	ENDIF
211
212	IF( clesphy0(5).NE.tab_cntrl( 9 ) ) THEN
213	clesphy0(5)=tab_cntrl( 9 )
214	ENDIF
215
216	IF( clesphy0(6).NE.tab_cntrl( 10 ) ) THEN
217	clesphy0(6)=tab_cntrl( 10 )
218	ENDIF
219
220	IF( clesphy0(7).NE.tab_cntrl( 11 ) ) THEN
221	clesphy0(7)=tab_cntrl( 11 )
222	ENDIF
223
224	IF( clesphy0(8).NE.tab_cntrl( 12 ) ) THEN
225	clesphy0(8)=tab_cntrl( 12 )
226	ENDIF
227
228
229	c
230	c Lecture des latitudes (coordonnees):
231	c
232	c$OMP MASTER
233	IF (is_mpi_root) THEN
234
235	ierr = NF_INQ_VARID (nid, "latitude", nvarid)
236	IF (ierr.NE.NF_NOERR) THEN
237	PRINT*, 'phyetat0: Le champ <latitude> est absent'
238	CALL abort
239	ENDIF
240	#ifdef NC_DOUBLE
241	ierr = NF_GET_VAR_DOUBLE(nid, nvarid, rlat_glo)
242	#else
243	ierr = NF_GET_VAR_REAL(nid, nvarid, rlat_glo)
244	#endif
245	IF (ierr.NE.NF_NOERR) THEN
246	PRINT*, 'phyetat0: Lecture echouee pour <latitude>'
247	CALL abort
248	ENDIF
249
250	c
251	c Lecture des longitudes (coordonnees):
252	c
253	ierr = NF_INQ_VARID (nid, "longitude", nvarid)
254	IF (ierr.NE.NF_NOERR) THEN
255	PRINT*, 'phyetat0: Le champ <longitude> est absent'
256	CALL abort
257	ENDIF
258	#ifdef NC_DOUBLE
259	ierr = NF_GET_VAR_DOUBLE(nid, nvarid, rlon_glo)
260	#else
261	ierr = NF_GET_VAR_REAL(nid, nvarid, rlon_glo)
262	#endif
263	IF (ierr.NE.NF_NOERR) THEN
264	PRINT*, 'phyetat0: Lecture echouee pour <latitude>'
265	CALL abort
266	ENDIF
267	C
268	C
269	C Lecture du masque terre mer
270	C
271
272	ierr = NF_INQ_VARID (nid, "masque", nvarid)
273	IF (ierr .EQ. NF_NOERR) THEN
274	#ifdef NC_DOUBLE
275	ierr = NF_GET_VAR_DOUBLE(nid, nvarid, zmasq_glo)
276	#else
277	ierr = NF_GET_VAR_REAL(nid, nvarid, zmasq_glo)
278	#endif
279	IF (ierr.NE.NF_NOERR) THEN
280	PRINT*, 'phyetat0: Lecture echouee pour <masque>'
281	CALL abort
282	ENDIF
283	else
284	PRINT*, 'phyetat0: Le champ <masque> est absent'
285	PRINT*, 'fichier startphy non compatible avec phyetat0'
286	C CALL abort
287	ENDIF
288
289
290	C Lecture des fractions pour chaque sous-surface
291	C
292	C initialisation des sous-surfaces
293	C
294	pctsrf = 0.
295	C
296	C fraction de terre
297	C
298
299	ierr = NF_INQ_VARID (nid, "FTER", nvarid)
300	IF (ierr .EQ. NF_NOERR) THEN
301	#ifdef NC_DOUBLE
302	ierr = NF_GET_VAR_DOUBLE(nid, nvarid,
303	. pctsrf_glo(1 : klon_glo,is_ter))
304	#else
305	ierr = NF_GET_VAR_REAL(nid, nvarid,
306	. pctsrf_glo(1 : klon_glo,is_ter))
307	#endif
308	IF (ierr.NE.NF_NOERR) THEN
309	PRINT*, 'phyetat0: Lecture echouee pour <FTER>'
310	CALL abort
311	ENDIF
312	else
313	PRINT*, 'phyetat0: Le champ <FTER> est absent'
314	c@$$ CALL abort
315	ENDIF
316
317	C
318	C fraction de glace de terre
319	C
320	ierr = NF_INQ_VARID (nid, "FLIC", nvarid)
321	IF (ierr .EQ. NF_NOERR) THEN
322	#ifdef NC_DOUBLE
323	ierr = NF_GET_VAR_DOUBLE(nid, nvarid,
324	. pctsrf_glo(1 : klon_glo,is_lic))
325	#else
326	ierr = NF_GET_VAR_REAL(nid, nvarid,
327	. pctsrf_glo(1 : klon_glo,is_lic))
328	#endif
329	IF (ierr.NE.NF_NOERR) THEN
330	PRINT*, 'phyetat0: Lecture echouee pour <FLIC>'
331	CALL abort
332	ENDIF
333	else
334	PRINT*, 'phyetat0: Le champ <FLIC> est absent'
335	c@$$ CALL abort
336	ENDIF
337	C
338	C fraction d'ocean
339	C
340	ierr = NF_INQ_VARID (nid, "FOCE", nvarid)
341	IF (ierr .EQ. NF_NOERR) THEN
342	#ifdef NC_DOUBLE
343	ierr = NF_GET_VAR_DOUBLE(nid, nvarid,
344	. pctsrf_glo(1 : klon_glo,is_oce))
345	#else
346	ierr = NF_GET_VAR_REAL(nid, nvarid,
347	. pctsrf_glo(1 : klon_glo,is_oce))
348	#endif
349	IF (ierr.NE.NF_NOERR) THEN
350	PRINT*, 'phyetat0: Lecture echouee pour <FOCE>'
351	CALL abort
352	ENDIF
353	else
354	PRINT*, 'phyetat0: Le champ <FOCE> est absent'
355	c@$$ CALL abort
356	ENDIF
357
358	C
359	C fraction glace de mer
360	C
361	ierr = NF_INQ_VARID (nid, "FSIC", nvarid)
362	IF (ierr .EQ. NF_NOERR) THEN
363	#ifdef NC_DOUBLE
364	ierr = NF_GET_VAR_DOUBLE(nid, nvarid,
365	. pctsrf_glo(1 : klon_glo,is_sic))
366	#else
367	ierr = NF_GET_VAR_REAL(nid, nvarid,
368	. pctsrf_glo(1 : klon_glo, is_sic))
369	#endif
370	IF (ierr.NE.NF_NOERR) THEN
371	PRINT*, 'phyetat0: Lecture echouee pour <FSIC>'
372	CALL abort
373	ENDIF
374	else
375	PRINT*, 'phyetat0: Le champ <FSIC> est absent'
376	c@$$ CALL abort
377	ENDIF
378
379	C
380	C Verification de l'adequation entre le masque et les sous-surfaces
381	C
382	fractint( 1 : klon_glo) = pctsrf_glo(1 : klon_glo, is_ter)
383	$ + pctsrf_glo(1 : klon_glo, is_lic)
384	DO i = 1 , klon_glo
385	IF ( abs(fractint(i) - zmasq_glo(i) ) .GT. EPSFRA ) THEN
386	WRITE(,) 'phyetat0: attention fraction terre pas ',
387	$ 'coherente ', i, zmasq_glo(i), pctsrf_glo(i, is_ter)
388	$ ,pctsrf_glo(i, is_lic)
389	ENDIF
390	END DO
391	fractint (1 : klon_glo) = pctsrf_glo(1 : klon_glo, is_oce)
392	$ + pctsrf_glo(1 : klon_glo, is_sic)
393	DO i = 1 , klon_glo
394	IF ( abs( fractint(i) - (1. - zmasq_glo(i))) .GT. EPSFRA ) THEN
395	WRITE(,) 'phyetat0 attention fraction ocean pas ',
396	$ 'coherente ', i, zmasq_glo(i) , pctsrf_glo(i, is_oce)
397	$ ,pctsrf_glo(i, is_sic)
398	ENDIF
399	END DO
400
401	C
402	c Lecture des temperatures du sol:
403	c
404
405	ierr = NF_INQ_VARID (nid, "TS", nvarid)
406	IF (ierr.NE.NF_NOERR) THEN
407	PRINT*, 'phyetat0: Le champ <TS> est absent'
408	PRINT, ' Mais je vais essayer de lire TS*'
409	DO nsrf = 1, nbsrf
410	IF (nsrf.GT.99) THEN
411	PRINT*, "Trop de sous-mailles"
412	CALL abort
413	ENDIF
414	WRITE(str2,'(i2.2)') nsrf
415	ierr = NF_INQ_VARID (nid, "TS"//str2, nvarid)
416	IF (ierr.NE.NF_NOERR) THEN
417	PRINT*, "phyetat0: Le champ <TS"//str2//"> est absent"
418	CALL abort
419	ENDIF
420	#ifdef NC_DOUBLE
421	ierr = NF_GET_VAR_DOUBLE(nid, nvarid, tsol_glo(1,nsrf))
422	#else
423	ierr = NF_GET_VAR_REAL(nid, nvarid, tsol_glo(1,nsrf))
424	#endif
425	IF (ierr.NE.NF_NOERR) THEN
426	PRINT*, "phyetat0: Lecture echouee pour <TS"//str2//">"
427	CALL abort
428	ENDIF
429
430	xmin = 1.0E+20
431	xmax = -1.0E+20
432	DO i = 1, klon_glo
433	xmin = MIN(tsol_glo(i,nsrf),xmin)
434	xmax = MAX(tsol_glo(i,nsrf),xmax)
435	ENDDO
436	PRINT,'Temperature du sol TS*:', nsrf, xmin, xmax
437	ENDDO
438	ELSE
439	PRINT*, 'phyetat0: Le champ <TS> est present'
440	PRINT, ' J ignore donc les autres temperatures TS*'
441	#ifdef NC_DOUBLE
442	ierr = NF_GET_VAR_DOUBLE(nid, nvarid, tsol_glo(1,1))
443	#else
444	ierr = NF_GET_VAR_REAL(nid, nvarid, tsol_glo(1,1))
445	#endif
446	IF (ierr.NE.NF_NOERR) THEN
447	PRINT*, "phyetat0: Lecture echouee pour <TS>"
448	CALL abort
449	ENDIF
450	xmin = 1.0E+20
451	xmax = -1.0E+20
452	DO i = 1, klon_glo
453	xmin = MIN(tsol_glo(i,1),xmin)
454	xmax = MAX(tsol_glo(i,1),xmax)
455	ENDDO
456	PRINT*,'Temperature du sol <TS>', xmin, xmax
457	DO nsrf = 2, nbsrf
458	DO i = 1, klon_glo
459	tsol_glo(i,nsrf) = tsol_glo(i,1)
460	ENDDO
461	ENDDO
462	ENDIF
463
464	c
465	c Lecture des temperatures du sol profond:
466	c
467	DO nsrf = 1, nbsrf
468	DO isoil=1, nsoilmx
469	IF (isoil.GT.99 .AND. nsrf.GT.99) THEN
470	PRINT*, "Trop de couches ou sous-mailles"
471	CALL abort
472	ENDIF
473	WRITE(str7,'(i2.2,"srf",i2.2)') isoil, nsrf
474	ierr = NF_INQ_VARID (nid, 'Tsoil'//str7, nvarid)
475	IF (ierr.NE.NF_NOERR) THEN
476	PRINT*, "phyetat0: Le champ <Tsoil"//str7//"> est absent"
477	PRINT*, " Il prend donc la valeur de surface"
478	DO i=1, klon_glo
479	tsoil(i,isoil,nsrf)=tsol_glo(i,nsrf)
480	ENDDO
481	ELSE
482	#ifdef NC_DOUBLE
483	ierr = NF_GET_VAR_DOUBLE(nid, nvarid, tsoil(1,isoil,nsrf))
484	#else
485	ierr = NF_GET_VAR_REAL(nid, nvarid, tsoil(1,isoil,nsrf))
486	#endif
487	IF (ierr.NE.NF_NOERR) THEN
488	PRINT*, "Lecture echouee pour <Tsoil"//str7//">"
489	CALL abort
490	ENDIF
491	ENDIF
492	ENDDO
493	ENDDO
494	c
495	cIM "slab" ocean
496	c
497	c Lecture de tslab (pour slab ocean seulement):
498	c
499	IF (ocean_in .eq. 'slab ') then
500	ierr = NF_INQ_VARID (nid, "TSLAB", nvarid)
501	IF (ierr.NE.NF_NOERR) THEN
502	PRINT*, "phyetat0: Le champ <TSLAB> est absent"
503	CALL abort
504	ENDIF
505	#ifdef NC_DOUBLE
506	ierr = NF_GET_VAR_DOUBLE(nid, nvarid, tslab)
507	#else
508	ierr = NF_GET_VAR_REAL(nid, nvarid, tslab)
509	#endif
510	IF (ierr.NE.NF_NOERR) THEN
511	PRINT*, "phyetat0: Lecture echouee pour <TSLAB>"
512	CALL abort
513	ENDIF
514	xmin = 1.0E+20
515	xmax = -1.0E+20
516	DO i = 1, klon_glo
517	xmin = MIN(tslab(i),xmin)
518	xmax = MAX(tslab(i),xmax)
519	ENDDO
520	PRINT*,'Min, Max tslab (utilise si OCEAN=slab )', xmin, xmax
521	c
522	c Lecture de seaice (pour slab ocean seulement):
523	c
524	ierr = NF_INQ_VARID (nid, "SEAICE", nvarid)
525	IF (ierr.NE.NF_NOERR) THEN
526	PRINT*, "phyetat0: Le champ <SEAICE> est absent"
527	CALL abort
528	ENDIF
529	#ifdef NC_DOUBLE
530	ierr = NF_GET_VAR_DOUBLE(nid, nvarid, seaice)
531	#else
532	ierr = NF_GET_VAR_REAL(nid, nvarid, seaice)
533	#endif
534	IF (ierr.NE.NF_NOERR) THEN
535	PRINT*, "phyetat0: Lecture echouee pour <SEAICE>"
536	CALL abort
537	ENDIF
538	xmin = 1.0E+20
539	xmax = -1.0E+20
540	DO i = 1, klon_glo
541	xmin = MIN(seaice(i),xmin)
542	xmax = MAX(seaice(i),xmax)
543	ENDDO
544	PRINT*,'Masse de la glace de mer (utilise si OCEAN=slab)',
545	$ xmin, xmax
546	ELSE
547	tslab = 0.
548	seaice = 0.
549	ENDIF
550	c
551	c Lecture de l'humidite de l'air juste au dessus du sol:
552	c
553	ierr = NF_INQ_VARID (nid, "QS", nvarid)
554	IF (ierr.NE.NF_NOERR) THEN
555	PRINT*, 'phyetat0: Le champ <QS> est absent'
556	PRINT, ' Mais je vais essayer de lire QS*'
557	DO nsrf = 1, nbsrf
558	IF (nsrf.GT.99) THEN
559	PRINT*, "Trop de sous-mailles"
560	CALL abort
561	ENDIF
562	WRITE(str2,'(i2.2)') nsrf
563	ierr = NF_INQ_VARID (nid, "QS"//str2, nvarid)
564	IF (ierr.NE.NF_NOERR) THEN
565	PRINT*, "phyetat0: Le champ <QS"//str2//"> est absent"
566	CALL abort
567	ENDIF
568	#ifdef NC_DOUBLE
569	ierr = NF_GET_VAR_DOUBLE(nid, nvarid, qsurf(1,nsrf))
570	#else
571	ierr = NF_GET_VAR_REAL(nid, nvarid, qsurf(1,nsrf))
572	#endif
573	IF (ierr.NE.NF_NOERR) THEN
574	PRINT*, "phyetat0: Lecture echouee pour <QS"//str2//">"
575	CALL abort
576	ENDIF
577	xmin = 1.0E+20
578	xmax = -1.0E+20
579	DO i = 1, klon_glo
580	xmin = MIN(qsurf(i,nsrf),xmin)
581	xmax = MAX(qsurf(i,nsrf),xmax)
582	ENDDO
583	PRINT,'Humidite pres du sol QS*:', nsrf, xmin, xmax
584	ENDDO
585	ELSE
586	PRINT*, 'phyetat0: Le champ <QS> est present'
587	PRINT, ' J ignore donc les autres humidites QS*'
588	#ifdef NC_DOUBLE
589	ierr = NF_GET_VAR_DOUBLE(nid, nvarid, qsurf(1,1))
590	#else
591	ierr = NF_GET_VAR_REAL(nid, nvarid, qsurf(1,1))
592	#endif
593	IF (ierr.NE.NF_NOERR) THEN
594	PRINT*, "phyetat0: Lecture echouee pour <QS>"
595	CALL abort
596	ENDIF
597	xmin = 1.0E+20
598	xmax = -1.0E+20
599	DO i = 1, klon_glo
600	xmin = MIN(qsurf(i,1),xmin)
601	xmax = MAX(qsurf(i,1),xmax)
602	ENDDO
603	PRINT*,'Humidite pres du sol <QS>', xmin, xmax
604	DO nsrf = 2, nbsrf
605	DO i = 1, klon_glo
606	qsurf(i,nsrf) = qsurf(i,1)
607	ENDDO
608	ENDDO
609	ENDIF
610	C
611	C Eau dans le sol (pour le modele de sol "bucket")
612	C
613	ierr = NF_INQ_VARID (nid, "QSOL", nvarid)
614	IF (ierr .EQ. NF_NOERR) THEN
615	#ifdef NC_DOUBLE
616	ierr = NF_GET_VAR_DOUBLE(nid, nvarid, qsol)
617	#else
618	ierr = NF_GET_VAR_REAL(nid, nvarid, qsol)
619	#endif
620	IF (ierr.NE.NF_NOERR) THEN
621	PRINT*, 'phyetat0: Lecture echouee pour <QSOL>'
622	CALL abort
623	ENDIF
624	else
625	PRINT*, 'phyetat0: Le champ <QSOL> est absent'
626	PRINT*, ' Valeur par defaut nulle'
627	qsol(:)=0.
628	c@$$ CALL abort
629	ENDIF
630	xmin = 1.0E+20
631	xmax = -1.0E+20
632	DO i = 1, klon_glo
633	xmin = MIN(qsol(i),xmin)
634	xmax = MAX(qsol(i),xmax)
635	ENDDO
636	PRINT*,'Eau dans le sol (mm) <QSOL>', xmin, xmax
637	c
638	c Lecture de neige au sol:
639	c
640	ierr = NF_INQ_VARID (nid, "SNOW", nvarid)
641	IF (ierr.NE.NF_NOERR) THEN
642	PRINT*, 'phyetat0: Le champ <SNOW> est absent'
643	PRINT, ' Mais je vais essayer de lire SNOW*'
644	DO nsrf = 1, nbsrf
645	IF (nsrf.GT.99) THEN
646	PRINT*, "Trop de sous-mailles"
647	CALL abort
648	ENDIF
649	WRITE(str2,'(i2.2)') nsrf
650	ierr = NF_INQ_VARID (nid, "SNOW"//str2, nvarid)
651	IF (ierr.NE.NF_NOERR) THEN
652	PRINT*, "phyetat0: Le champ <SNOW"//str2//"> est absent"
653	CALL abort
654	ENDIF
655	#ifdef NC_DOUBLE
656	ierr = NF_GET_VAR_DOUBLE(nid, nvarid, snow(1,nsrf))
657	#else
658	ierr = NF_GET_VAR_REAL(nid, nvarid, snow(1,nsrf))
659	#endif
660	IF (ierr.NE.NF_NOERR) THEN
661	PRINT*, "phyetat0: Lecture echouee pour <SNOW"//str2//">"
662	CALL abort
663	ENDIF
664	xmin = 1.0E+20
665	xmax = -1.0E+20
666	DO i = 1, klon_glo
667	xmin = MIN(snow(i,nsrf),xmin)
668	xmax = MAX(snow(i,nsrf),xmax)
669	ENDDO
670	PRINT,'Neige du sol SNOW*:', nsrf, xmin, xmax
671	ENDDO
672	ELSE
673	PRINT*, 'phyetat0: Le champ <SNOW> est present'
674	PRINT, ' J ignore donc les autres neiges SNOW*'
675	#ifdef NC_DOUBLE
676	ierr = NF_GET_VAR_DOUBLE(nid, nvarid, snow(1,1))
677	#else
678	ierr = NF_GET_VAR_REAL(nid, nvarid, snow(1,1))
679	#endif
680	IF (ierr.NE.NF_NOERR) THEN
681	PRINT*, "phyetat0: Lecture echouee pour <SNOW>"
682	CALL abort
683	ENDIF
684	xmin = 1.0E+20
685	xmax = -1.0E+20
686	DO i = 1, klon_glo
687	xmin = MIN(snow(i,1),xmin)
688	xmax = MAX(snow(i,1),xmax)
689	ENDDO
690	PRINT*,'Neige du sol <SNOW>', xmin, xmax
691	DO nsrf = 2, nbsrf
692	DO i = 1, klon_glo
693	snow(i,nsrf) = snow(i,1)
694	ENDDO
695	ENDDO
696	ENDIF
697	c
698	c Lecture de albedo de l'interval visible au sol:
699	c
700	ierr = NF_INQ_VARID (nid, "ALBE", nvarid)
701	IF (ierr.NE.NF_NOERR) THEN
702	PRINT*, 'phyetat0: Le champ <ALBE> est absent'
703	PRINT, ' Mais je vais essayer de lire ALBE*'
704	DO nsrf = 1, nbsrf
705	IF (nsrf.GT.99) THEN
706	PRINT*, "Trop de sous-mailles"
707	CALL abort
708	ENDIF
709	WRITE(str2,'(i2.2)') nsrf
710	ierr = NF_INQ_VARID (nid, "ALBE"//str2, nvarid)
711	IF (ierr.NE.NF_NOERR) THEN
712	PRINT*, "phyetat0: Le champ <ALBE"//str2//"> est absent"
713	CALL abort
714	ENDIF
715	#ifdef NC_DOUBLE
716	ierr = NF_GET_VAR_DOUBLE(nid, nvarid, alb1_glo(1,nsrf))
717	#else
718	ierr = NF_GET_VAR_REAL(nid, nvarid, alb1_glo(1,nsrf))
719	#endif
720	IF (ierr.NE.NF_NOERR) THEN
721	PRINT*, "phyetat0: Lecture echouee pour <ALBE"//str2//">"
722	CALL abort
723	ENDIF
724	xmin = 1.0E+20
725	xmax = -1.0E+20
726	DO i = 1, klon_glo
727	xmin = MIN(alb1_glo(i,nsrf),xmin)
728	xmax = MAX(alb1_glo(i,nsrf),xmax)
729	ENDDO
730	PRINT,'Albedo du sol ALBE*:', nsrf, xmin, xmax
731	ENDDO
732	ELSE
733	PRINT*, 'phyetat0: Le champ <ALBE> est present'
734	PRINT, ' J ignore donc les autres ALBE*'
735	#ifdef NC_DOUBLE
736	ierr = NF_GET_VAR_DOUBLE(nid, nvarid, alb1_glo(1,1))
737	#else
738	ierr = NF_GET_VAR_REAL(nid, nvarid, alb1_glo(1,1))
739	#endif
740	IF (ierr.NE.NF_NOERR) THEN
741	PRINT*, "phyetat0: Lecture echouee pour <ALBE>"
742	CALL abort
743	ENDIF
744	xmin = 1.0E+20
745	xmax = -1.0E+20
746	DO i = 1, klon_glo
747	xmin = MIN(alb1_glo(i,1),xmin)
748	xmax = MAX(alb1_glo(i,1),xmax)
749	ENDDO
750	PRINT*,'Neige du sol <ALBE>', xmin, xmax
751	DO nsrf = 2, nbsrf
752	DO i = 1, klon_glo
753	alb1_glo(i,nsrf) = alb1_glo(i,1)
754	ENDDO
755	ENDDO
756	ENDIF
757
758	c
759	c Lecture de albedo au sol dans l'interval proche infra-rouge:
760	c
761	ierr = NF_INQ_VARID (nid, "ALBLW", nvarid)
762	IF (ierr.NE.NF_NOERR) THEN
763	PRINT*, 'phyetat0: Le champ <ALBLW> est absent'
764	c PRINT, ' Mais je vais essayer de lire ALBLW*'
765	PRINT, ' Mais je vais prendre ALBE*'
766	DO nsrf = 1, nbsrf
767	DO i = 1, klon_glo
768	alb2_glo(i,nsrf) = alb1_glo(i,nsrf)
769	ENDDO
770	ENDDO
771	ELSE
772	PRINT*, 'phyetat0: Le champ <ALBLW> est present'
773	PRINT, ' J ignore donc les autres ALBLW*'
774	#ifdef NC_DOUBLE
775	ierr = NF_GET_VAR_DOUBLE(nid, nvarid, alb2_glo(1,1))
776	#else
777	ierr = NF_GET_VAR_REAL(nid, nvarid, alb2_glo(1,1))
778	#endif
779	IF (ierr.NE.NF_NOERR) THEN
780	PRINT*, "phyetat0: Lecture echouee pour <ALBLW>"
781	CALL abort
782	ENDIF
783	xmin = 1.0E+20
784	xmax = -1.0E+20
785	DO i = 1, klon_glo
786	xmin = MIN(alb2_glo(i,1),xmin)
787	xmax = MAX(alb2_glo(i,1),xmax)
788	ENDDO
789	PRINT*,'Neige du sol <ALBLW>', xmin, xmax
790	DO nsrf = 2, nbsrf
791	DO i = 1, klon_glo
792	alb2_glo(i,nsrf) = alb2_glo(i,1)
793	ENDDO
794	ENDDO
795	ENDIF
796	c
797	c Lecture de evaporation:
798	c
799	ierr = NF_INQ_VARID (nid, "EVAP", nvarid)
800	IF (ierr.NE.NF_NOERR) THEN
801	PRINT*, 'phyetat0: Le champ <EVAP> est absent'
802	PRINT, ' Mais je vais essayer de lire EVAP*'
803	DO nsrf = 1, nbsrf
804	IF (nsrf.GT.99) THEN
805	PRINT*, "Trop de sous-mailles"
806	CALL abort
807	ENDIF
808	WRITE(str2,'(i2.2)') nsrf
809	ierr = NF_INQ_VARID (nid, "EVAP"//str2, nvarid)
810	IF (ierr.NE.NF_NOERR) THEN
811	PRINT*, "phyetat0: Le champ <EVAP"//str2//"> est absent"
812	CALL abort
813	ENDIF
814	#ifdef NC_DOUBLE
815	ierr = NF_GET_VAR_DOUBLE(nid, nvarid, evap(1,nsrf))
816	#else
817	ierr = NF_GET_VAR_REAL(nid, nvarid, evap(1,nsrf))
818	#endif
819	IF (ierr.NE.NF_NOERR) THEN
820	PRINT*, "phyetat0: Lecture echouee pour <EVAP"//str2//">"
821	CALL abort
822	ENDIF
823	xmin = 1.0E+20
824	xmax = -1.0E+20
825	DO i = 1, klon_glo
826	xmin = MIN(evap(i,nsrf),xmin)
827	xmax = MAX(evap(i,nsrf),xmax)
828	ENDDO
829	PRINT,'evap du sol EVAP*:', nsrf, xmin, xmax
830	ENDDO
831	ELSE
832	PRINT*, 'phyetat0: Le champ <EVAP> est present'
833	PRINT, ' J ignore donc les autres EVAP*'
834	#ifdef NC_DOUBLE
835	ierr = NF_GET_VAR_DOUBLE(nid, nvarid, evap(1,1))
836	#else
837	ierr = NF_GET_VAR_REAL(nid, nvarid, evap(1,1))
838	#endif
839	IF (ierr.NE.NF_NOERR) THEN
840	PRINT*, "phyetat0: Lecture echouee pour <EVAP>"
841	CALL abort
842	ENDIF
843	xmin = 1.0E+20
844	xmax = -1.0E+20
845	DO i = 1, klon_glo
846	xmin = MIN(evap(i,1),xmin)
847	xmax = MAX(evap(i,1),xmax)
848	ENDDO
849	PRINT*,'Evap du sol <EVAP>', xmin, xmax
850	DO nsrf = 2, nbsrf
851	DO i = 1, klon_glo
852	evap(i,nsrf) = evap(i,1)
853	ENDDO
854	ENDDO
855	ENDIF
856	c
857	c Lecture precipitation liquide:
858	c
859	ierr = NF_INQ_VARID (nid, "rain_f", nvarid)
860	IF (ierr.NE.NF_NOERR) THEN
861	PRINT*, 'phyetat0: Le champ <rain_f> est absent'
862	CALL abort
863	ENDIF
864	#ifdef NC_DOUBLE
865	ierr = NF_GET_VAR_DOUBLE(nid, nvarid, rain_fall_glo)
866	#else
867	ierr = NF_GET_VAR_REAL(nid, nvarid, rain_fall_glo)
868	#endif
869	IF (ierr.NE.NF_NOERR) THEN
870	PRINT*, 'phyetat0: Lecture echouee pour <rain_f>'
871	CALL abort
872	ENDIF
873	xmin = 1.0E+20
874	xmax = -1.0E+20
875	DO i = 1, klon_glo
876	xmin = MIN(rain_fall_glo(i),xmin)
877	xmax = MAX(rain_fall_glo(i),xmax)
878	ENDDO
879	PRINT*,'Precipitation liquide rain_f:', xmin, xmax
880	c
881	c Lecture precipitation solide:
882	c
883	ierr = NF_INQ_VARID (nid, "snow_f", nvarid)
884	IF (ierr.NE.NF_NOERR) THEN
885	PRINT*, 'phyetat0: Le champ <snow_f> est absent'
886	CALL abort
887	ENDIF
888	#ifdef NC_DOUBLE
889	ierr = NF_GET_VAR_DOUBLE(nid, nvarid, snow_fall_glo)
890	#else
891	ierr = NF_GET_VAR_REAL(nid, nvarid, snow_fall_glo)
892	#endif
893	IF (ierr.NE.NF_NOERR) THEN
894	PRINT*, 'phyetat0: Lecture echouee pour <snow_f>'
895	CALL abort
896	ENDIF
897	xmin = 1.0E+20
898	xmax = -1.0E+20
899	DO i = 1, klon_glo
900	xmin = MIN(snow_fall_glo(i),xmin)
901	xmax = MAX(snow_fall_glo(i),xmax)
902	ENDDO
903	PRINT*,'Precipitation solide snow_f:', xmin, xmax
904	c
905	c Lecture rayonnement solaire au sol:
906	c
907	ierr = NF_INQ_VARID (nid, "solsw", nvarid)
908	IF (ierr.NE.NF_NOERR) THEN
909	PRINT*, 'phyetat0: Le champ <solsw> est absent'
910	PRINT*, 'mis a zero'
911	solsw_glo = 0.
912	ELSE
913	#ifdef NC_DOUBLE
914	ierr = NF_GET_VAR_DOUBLE(nid, nvarid, solsw_glo)
915	#else
916	ierr = NF_GET_VAR_REAL(nid, nvarid, solsw_glo)
917	#endif
918	IF (ierr.NE.NF_NOERR) THEN
919	PRINT*, 'phyetat0: Lecture echouee pour <solsw>'
920	CALL abort
921	ENDIF
922	ENDIF
923	xmin = 1.0E+20
924	xmax = -1.0E+20
925	DO i = 1, klon_glo
926	xmin = MIN(solsw_glo(i),xmin)
927	xmax = MAX(solsw_glo(i),xmax)
928	ENDDO
929	PRINT*,'Rayonnement solaire au sol solsw:', xmin, xmax
930	c
931	c Lecture rayonnement IF au sol:
932	c
933	ierr = NF_INQ_VARID (nid, "sollw", nvarid)
934	IF (ierr.NE.NF_NOERR) THEN
935	PRINT*, 'phyetat0: Le champ <sollw> est absent'
936	PRINT*, 'mis a zero'
937	sollw_glo = 0.
938	ELSE
939	#ifdef NC_DOUBLE
940	ierr = NF_GET_VAR_DOUBLE(nid, nvarid, sollw_glo)
941	#else
942	ierr = NF_GET_VAR_REAL(nid, nvarid, sollw_glo)
943	#endif
944	IF (ierr.NE.NF_NOERR) THEN
945	PRINT*, 'phyetat0: Lecture echouee pour <sollw>'
946	CALL abort
947	ENDIF
948	ENDIF
949	xmin = 1.0E+20
950	xmax = -1.0E+20
951	DO i = 1, klon_glo
952	xmin = MIN(sollw_glo(i),xmin)
953	xmax = MAX(sollw_glo(i),xmax)
954	ENDDO
955	PRINT*,'Rayonnement IF au sol sollw:', xmin, xmax
956
957	ENDIF ! is_mpi_root
958	c$OMP END MASTER
959
960
961	c$OMP MASTER
962	IF (is_mpi_root) THEN
963	c
964	c Lecture derive des flux:
965	c
966	ierr = NF_INQ_VARID (nid, "fder", nvarid)
967	IF (ierr.NE.NF_NOERR) THEN
968	PRINT*, 'phyetat0: Le champ <fder> est absent'
969	PRINT*, 'mis a zero'
970	fder = 0.
971	ELSE
972	#ifdef NC_DOUBLE
973	ierr = NF_GET_VAR_DOUBLE(nid, nvarid, fder)
974	#else
975	ierr = NF_GET_VAR_REAL(nid, nvarid, fder)
976	#endif
977	IF (ierr.NE.NF_NOERR) THEN
978	PRINT*, 'phyetat0: Lecture echouee pour <fder>'
979	CALL abort
980	ENDIF
981	ENDIF
982	xmin = 1.0E+20
983	xmax = -1.0E+20
984	DO i = 1, klon_glo
985	xmin = MIN(fder(i),xmin)
986	xmax = MAX(fder(i),xmax)
987	ENDDO
988	PRINT*,'Derive des flux fder:', xmin, xmax
989
990	c
991	c Lecture du rayonnement net au sol:
992	c
993	ierr = NF_INQ_VARID (nid, "RADS", nvarid)
994	IF (ierr.NE.NF_NOERR) THEN
995	PRINT*, 'phyetat0: Le champ <RADS> est absent'
996	CALL abort
997	ENDIF
998	#ifdef NC_DOUBLE
999	ierr = NF_GET_VAR_DOUBLE(nid, nvarid, radsol_glo)
1000	#else
1001	ierr = NF_GET_VAR_REAL(nid, nvarid, radsol_glo)
1002	#endif
1003	IF (ierr.NE.NF_NOERR) THEN
1004	PRINT*, 'phyetat0: Lecture echouee pour <RADS>'
1005	CALL abort
1006	ENDIF
1007	xmin = 1.0E+20
1008	xmax = -1.0E+20
1009	DO i = 1, klon_glo
1010	xmin = MIN(radsol_glo(i),xmin)
1011	xmax = MAX(radsol_glo(i),xmax)
1012	ENDDO
1013	PRINT*,'Rayonnement net au sol radsol:', xmin, xmax
1014	c
1015	c Lecture de la longueur de rugosite
1016	c
1017	c
1018	ierr = NF_INQ_VARID (nid, "RUG", nvarid)
1019	IF (ierr.NE.NF_NOERR) THEN
1020	PRINT*, 'phyetat0: Le champ <RUG> est absent'
1021	PRINT, ' Mais je vais essayer de lire RUG*'
1022	DO nsrf = 1, nbsrf
1023	IF (nsrf.GT.99) THEN
1024	PRINT*, "Trop de sous-mailles"
1025	CALL abort
1026	ENDIF
1027	WRITE(str2,'(i2.2)') nsrf
1028	ierr = NF_INQ_VARID (nid, "RUG"//str2, nvarid)
1029	IF (ierr.NE.NF_NOERR) THEN
1030	PRINT*, "phyetat0: Le champ <RUG"//str2//"> est absent"
1031	CALL abort
1032	ENDIF
1033	#ifdef NC_DOUBLE
1034	ierr = NF_GET_VAR_DOUBLE(nid, nvarid, frugs(1,nsrf))
1035	#else
1036	ierr = NF_GET_VAR_REAL(nid, nvarid, frugs(1,nsrf))
1037	#endif
1038	IF (ierr.NE.NF_NOERR) THEN
1039	PRINT*, "phyetat0: Lecture echouee pour <RUG"//str2//">"
1040	CALL abort
1041	ENDIF
1042	xmin = 1.0E+20
1043	xmax = -1.0E+20
1044	DO i = 1, klon_glo
1045	xmin = MIN(frugs(i,nsrf),xmin)
1046	xmax = MAX(frugs(i,nsrf),xmax)
1047	ENDDO
1048	PRINT,'rugosite du sol RUG*:', nsrf, xmin, xmax
1049	ENDDO
1050	ELSE
1051	PRINT*, 'phyetat0: Le champ <RUG> est present'
1052	PRINT, ' J ignore donc les autres RUG*'
1053	#ifdef NC_DOUBLE
1054	ierr = NF_GET_VAR_DOUBLE(nid, nvarid, frugs(1,1))
1055	#else
1056	ierr = NF_GET_VAR_REAL(nid, nvarid, frugs(1,1))
1057	#endif
1058	IF (ierr.NE.NF_NOERR) THEN
1059	PRINT*, "phyetat0: Lecture echouee pour <RUG>"
1060	CALL abort
1061	ENDIF
1062	xmin = 1.0E+20
1063	xmax = -1.0E+20
1064	DO i = 1, klon_glo
1065	xmin = MIN(frugs(i,1),xmin)
1066	xmax = MAX(frugs(i,1),xmax)
1067	ENDDO
1068	PRINT*,'rugosite <RUG>', xmin, xmax
1069	DO nsrf = 2, nbsrf
1070	DO i = 1, klon_glo
1071	frugs(i,nsrf) = frugs(i,1)
1072	ENDDO
1073	ENDDO
1074	ENDIF
1075
1076	c
1077	c Lecture de l'age de la neige:
1078	c
1079	ierr = NF_INQ_VARID (nid, "AGESNO", nvarid)
1080	IF (ierr.NE.NF_NOERR) THEN
1081	PRINT*, 'phyetat0: Le champ <AGESNO> est absent'
1082	PRINT, ' Mais je vais essayer de lire AGESNO*'
1083	DO nsrf = 1, nbsrf
1084	IF (nsrf.GT.99) THEN
1085	PRINT*, "Trop de sous-mailles"
1086	CALL abort
1087	ENDIF
1088	WRITE(str2,'(i2.2)') nsrf
1089	ierr = NF_INQ_VARID (nid, "AGESNO"//str2, nvarid)
1090	IF (ierr.NE.NF_NOERR) THEN
1091	PRINT*, "phyetat0: Le champ <AGESNO"//str2//"> est absent"
1092	agesno = 50.0
1093	ENDIF
1094	#ifdef NC_DOUBLE
1095	ierr = NF_GET_VAR_DOUBLE(nid, nvarid, agesno(1,nsrf))
1096	#else
1097	ierr = NF_GET_VAR_REAL(nid, nvarid, agesno(1,nsrf))
1098	#endif
1099	IF (ierr.NE.NF_NOERR) THEN
1100	PRINT*, "phyetat0: Lecture echouee pour <AGESNO"//str2//">"
1101	CALL abort
1102	ENDIF
1103	xmin = 1.0E+20
1104	xmax = -1.0E+20
1105	DO i = 1, klon_glo
1106	xmin = MIN(agesno(i,nsrf),xmin)
1107	xmax = MAX(agesno(i,nsrf),xmax)
1108	ENDDO
1109	PRINT,'Age de la neige AGESNO*:', nsrf, xmin, xmax
1110	ENDDO
1111	ELSE
1112	PRINT*, 'phyetat0: Le champ <AGESNO> est present'
1113	PRINT, ' J ignore donc les autres AGESNO*'
1114	#ifdef NC_DOUBLE
1115	ierr = NF_GET_VAR_DOUBLE(nid, nvarid, agesno(1,1))
1116	#else
1117	ierr = NF_GET_VAR_REAL(nid, nvarid, agesno(1,1))
1118	#endif
1119	IF (ierr.NE.NF_NOERR) THEN
1120	PRINT*, "phyetat0: Lecture echouee pour <AGESNO>"
1121	CALL abort
1122	ENDIF
1123	xmin = 1.0E+20
1124	xmax = -1.0E+20
1125	DO i = 1, klon_glo
1126	xmin = MIN(agesno(i,1),xmin)
1127	xmax = MAX(agesno(i,1),xmax)
1128	ENDDO
1129	PRINT*,'Age de la neige <AGESNO>', xmin, xmax
1130	DO nsrf = 2, nbsrf
1131	DO i = 1, klon_glo
1132	agesno(i,nsrf) = agesno(i,1)
1133	ENDDO
1134	ENDDO
1135	ENDIF
1136
1137	c
1138	ierr = NF_INQ_VARID (nid, "ZMEA", nvarid)
1139	IF (ierr.NE.NF_NOERR) THEN
1140	PRINT*, 'phyetat0: Le champ <ZMEA> est absent'
1141	CALL abort
1142	ENDIF
1143	#ifdef NC_DOUBLE
1144	ierr = NF_GET_VAR_DOUBLE(nid, nvarid, zmea_glo)
1145	#else
1146	ierr = NF_GET_VAR_REAL(nid, nvarid, zmea_glo)
1147	#endif
1148	IF (ierr.NE.NF_NOERR) THEN
1149	PRINT*, 'phyetat0: Lecture echouee pour <ZMEA>'
1150	CALL abort
1151	ENDIF
1152	xmin = 1.0E+20
1153	xmax = -1.0E+20
1154	DO i = 1, klon_glo
1155	xmin = MIN(zmea_glo(i),xmin)
1156	xmax = MAX(zmea_glo(i),xmax)
1157	ENDDO
1158	PRINT*,'OROGRAPHIE SOUS-MAILLE zmea:', xmin, xmax
1159	c
1160	c
1161	ierr = NF_INQ_VARID (nid, "ZSTD", nvarid)
1162	IF (ierr.NE.NF_NOERR) THEN
1163	PRINT*, 'phyetat0: Le champ <ZSTD> est absent'
1164	CALL abort
1165	ENDIF
1166	#ifdef NC_DOUBLE
1167	ierr = NF_GET_VAR_DOUBLE(nid, nvarid, zstd_glo)
1168	#else
1169	ierr = NF_GET_VAR_REAL(nid, nvarid, zstd_glo)
1170	#endif
1171	IF (ierr.NE.NF_NOERR) THEN
1172	PRINT*, 'phyetat0: Lecture echouee pour <ZSTD>'
1173	CALL abort
1174	ENDIF
1175	xmin = 1.0E+20
1176	xmax = -1.0E+20
1177	DO i = 1, klon_glo
1178	xmin = MIN(zstd_glo(i),xmin)
1179	xmax = MAX(zstd_glo(i),xmax)
1180	ENDDO
1181	PRINT*,'OROGRAPHIE SOUS-MAILLE zstd:', xmin, xmax
1182	c
1183	c
1184	ierr = NF_INQ_VARID (nid, "ZSIG", nvarid)
1185	IF (ierr.NE.NF_NOERR) THEN
1186	PRINT*, 'phyetat0: Le champ <ZSIG> est absent'
1187	CALL abort
1188	ENDIF
1189	#ifdef NC_DOUBLE
1190	ierr = NF_GET_VAR_DOUBLE(nid, nvarid, zsig_glo)
1191	#else
1192	ierr = NF_GET_VAR_REAL(nid, nvarid, zsig_glo)
1193	#endif
1194	IF (ierr.NE.NF_NOERR) THEN
1195	PRINT*, 'phyetat0: Lecture echouee pour <ZSIG>'
1196	CALL abort
1197	ENDIF
1198	xmin = 1.0E+20
1199	xmax = -1.0E+20
1200	DO i = 1, klon_glo
1201	xmin = MIN(zsig_glo(i),xmin)
1202	xmax = MAX(zsig_glo(i),xmax)
1203	ENDDO
1204	PRINT*,'OROGRAPHIE SOUS-MAILLE zsig:', xmin, xmax
1205	c
1206	c
1207	ierr = NF_INQ_VARID (nid, "ZGAM", nvarid)
1208	IF (ierr.NE.NF_NOERR) THEN
1209	PRINT*, 'phyetat0: Le champ <ZGAM> est absent'
1210	CALL abort
1211	ENDIF
1212	#ifdef NC_DOUBLE
1213	ierr = NF_GET_VAR_DOUBLE(nid, nvarid, zgam_glo)
1214	#else
1215	ierr = NF_GET_VAR_REAL(nid, nvarid, zgam_glo)
1216	#endif
1217	IF (ierr.NE.NF_NOERR) THEN
1218	PRINT*, 'phyetat0: Lecture echouee pour <ZGAM>'
1219	CALL abort
1220	ENDIF
1221	xmin = 1.0E+20
1222	xmax = -1.0E+20
1223	DO i = 1, klon_glo
1224	xmin = MIN(zgam_glo(i),xmin)
1225	xmax = MAX(zgam_glo(i),xmax)
1226	ENDDO
1227	PRINT*,'OROGRAPHIE SOUS-MAILLE zgam:', xmin, xmax
1228	c
1229	c
1230	ierr = NF_INQ_VARID (nid, "ZTHE", nvarid)
1231	IF (ierr.NE.NF_NOERR) THEN
1232	PRINT*, 'phyetat0: Le champ <ZTHE> est absent'
1233	CALL abort
1234	ENDIF
1235	#ifdef NC_DOUBLE
1236	ierr = NF_GET_VAR_DOUBLE(nid, nvarid, zthe_glo)
1237	#else
1238	ierr = NF_GET_VAR_REAL(nid, nvarid, zthe_glo)
1239	#endif
1240	IF (ierr.NE.NF_NOERR) THEN
1241	PRINT*, 'phyetat0: Lecture echouee pour <ZTHE>'
1242	CALL abort
1243	ENDIF
1244	xmin = 1.0E+20
1245	xmax = -1.0E+20
1246	DO i = 1, klon_glo
1247	xmin = MIN(zthe_glo(i),xmin)
1248	xmax = MAX(zthe_glo(i),xmax)
1249	ENDDO
1250	PRINT*,'OROGRAPHIE SOUS-MAILLE zthe:', xmin, xmax
1251	c
1252	c
1253	ierr = NF_INQ_VARID (nid, "ZPIC", nvarid)
1254	IF (ierr.NE.NF_NOERR) THEN
1255	PRINT*, 'phyetat0: Le champ <ZPIC> est absent'
1256	CALL abort
1257	ENDIF
1258	#ifdef NC_DOUBLE
1259	ierr = NF_GET_VAR_DOUBLE(nid, nvarid, zpic_glo)
1260	#else
1261	ierr = NF_GET_VAR_REAL(nid, nvarid, zpic_glo)
1262	#endif
1263	IF (ierr.NE.NF_NOERR) THEN
1264	PRINT*, 'phyetat0: Lecture echouee pour <ZPIC>'
1265	CALL abort
1266	ENDIF
1267	xmin = 1.0E+20
1268	xmax = -1.0E+20
1269	DO i = 1, klon_glo
1270	xmin = MIN(zpic_glo(i),xmin)
1271	xmax = MAX(zpic_glo(i),xmax)
1272	ENDDO
1273	PRINT*,'OROGRAPHIE SOUS-MAILLE zpic:', xmin, xmax
1274	c
1275	ierr = NF_INQ_VARID (nid, "ZVAL", nvarid)
1276	IF (ierr.NE.NF_NOERR) THEN
1277	PRINT*, 'phyetat0: Le champ <ZVAL> est absent'
1278	CALL abort
1279	ENDIF
1280	#ifdef NC_DOUBLE
1281	ierr = NF_GET_VAR_DOUBLE(nid, nvarid, zval_glo)
1282	#else
1283	ierr = NF_GET_VAR_REAL(nid, nvarid, zval_glo)
1284	#endif
1285	IF (ierr.NE.NF_NOERR) THEN
1286	PRINT*, 'phyetat0: Lecture echouee pour <ZVAL>'
1287	CALL abort
1288	ENDIF
1289	xmin = 1.0E+20
1290	xmax = -1.0E+20
1291	DO i = 1, klon_glo
1292	xmin = MIN(zval_glo(i),xmin)
1293	xmax = MAX(zval_glo(i),xmax)
1294	ENDDO
1295	PRINT*,'OROGRAPHIE SOUS-MAILLE zval:', xmin, xmax
1296	c
1297	c
1298	ierr = NF_INQ_VARID (nid, "RUGSREL", nvarid)
1299	IF (ierr.NE.NF_NOERR) THEN
1300	PRINT*, 'phyetat0: Le champ <RUGSREL> est absent'
1301	CALL abort
1302	ENDIF
1303	#ifdef NC_DOUBLE
1304	ierr = NF_GET_VAR_DOUBLE(nid, nvarid, rugsrel_glo)
1305	#else
1306	ierr = NF_GET_VAR_REAL(nid, nvarid, rugsrel_glo)
1307	#endif
1308	IF (ierr.NE.NF_NOERR) THEN
1309	PRINT*, 'phyetat0: Lecture echouee pour <RUGSREL>'
1310	CALL abort
1311	ENDIF
1312	xmin = 1.0E+20
1313	xmax = -1.0E+20
1314	DO i = 1, klon_glo
1315	xmin = MIN(rugsrel_glo(i),xmin)
1316	xmax = MAX(rugsrel_glo(i),xmax)
1317	ENDDO
1318	PRINT*,'Rugosite relief (ecart-type) rugsrel:', xmin, xmax
1319	c
1320	c
1321	ancien_ok_glo = .TRUE.
1322	c
1323	ierr = NF_INQ_VARID (nid, "TANCIEN", nvarid)
1324	IF (ierr.NE.NF_NOERR) THEN
1325	PRINT*, "phyetat0: Le champ <TANCIEN> est absent"
1326	PRINT*, "Depart legerement fausse. Mais je continue"
1327	ancien_ok_glo = .FALSE.
1328	ELSE
1329	#ifdef NC_DOUBLE
1330	ierr = NF_GET_VAR_DOUBLE(nid, nvarid, t_ancien_glo)
1331	#else
1332	ierr = NF_GET_VAR_REAL(nid, nvarid, t_ancien_glo)
1333	#endif
1334	IF (ierr.NE.NF_NOERR) THEN
1335	PRINT*, "phyetat0: Lecture echouee pour <TANCIEN>"
1336	CALL abort
1337	ENDIF
1338	ENDIF
1339	c
1340	ierr = NF_INQ_VARID (nid, "QANCIEN", nvarid)
1341	IF (ierr.NE.NF_NOERR) THEN
1342	PRINT*, "phyetat0: Le champ <QANCIEN> est absent"
1343	PRINT*, "Depart legerement fausse. Mais je continue"
1344	ancien_ok_glo = .FALSE.
1345	ELSE
1346	#ifdef NC_DOUBLE
1347	ierr = NF_GET_VAR_DOUBLE(nid, nvarid, q_ancien_glo)
1348	#else
1349	ierr = NF_GET_VAR_REAL(nid, nvarid, q_ancien_glo)
1350	#endif
1351	IF (ierr.NE.NF_NOERR) THEN
1352	PRINT*, "phyetat0: Lecture echouee pour <QANCIEN>"
1353	CALL abort
1354	ENDIF
1355	ENDIF
1356	c
1357	clwcon_glo=0.
1358	ierr = NF_INQ_VARID (nid, "CLWCON", nvarid)
1359	IF (ierr.NE.NF_NOERR) THEN
1360	PRINT*, "phyetat0: Le champ CLWCON est absent"
1361	PRINT*, "Depart legerement fausse. Mais je continue"
1362	c clwcon_glo = 0.
1363	ELSE
1364	#ifdef NC_DOUBLE
1365	ierr = NF_GET_VAR_DOUBLE(nid, nvarid, clwcon_glo)
1366	#else
1367	ierr = NF_GET_VAR_REAL(nid, nvarid, clwcon_glo)
1368	#endif
1369	IF (ierr.NE.NF_NOERR) THEN
1370	PRINT*, "phyetat0: Lecture echouee pour <CLWCON>"
1371	CALL abort
1372	ENDIF
1373	ENDIF
1374	xmin = 1.0E+20
1375	xmax = -1.0E+20
1376	xmin = MINval(clwcon_glo)
1377	xmax = MAXval(clwcon_glo)
1378	PRINT*,'Eau liquide convective (ecart-type) clwcon:', xmin, xmax
1379	c
1380	rnebcon_glo = 0.
1381	ierr = NF_INQ_VARID (nid, "RNEBCON", nvarid)
1382	IF (ierr.NE.NF_NOERR) THEN
1383	PRINT*, "phyetat0: Le champ RNEBCON est absent"
1384	PRINT*, "Depart legerement fausse. Mais je continue"
1385	c rnebcon_glo = 0.
1386	ELSE
1387	#ifdef NC_DOUBLE
1388	ierr = NF_GET_VAR_DOUBLE(nid, nvarid, rnebcon_glo)
1389	#else
1390	ierr = NF_GET_VAR_REAL(nid, nvarid, rnebcon_glo)
1391	#endif
1392	IF (ierr.NE.NF_NOERR) THEN
1393	PRINT*, "phyetat0: Lecture echouee pour <RNEBCON>"
1394	CALL abort
1395	ENDIF
1396	ENDIF
1397	xmin = 1.0E+20
1398	xmax = -1.0E+20
1399	xmin = MINval(rnebcon_glo)
1400	xmax = MAXval(rnebcon_glo)
1401	PRINT*,'Nebulosite convective (ecart-type) rnebcon:', xmin, xmax
1402
1403	c
1404	c Lecture ratqs
1405	c
1406	ratqs_glo=0.
1407	ierr = NF_INQ_VARID (nid, "RATQS", nvarid)
1408	IF (ierr.NE.NF_NOERR) THEN
1409	PRINT*, "phyetat0: Le champ <RATQS> est absent"
1410	PRINT*, "Depart legerement fausse. Mais je continue"
1411	ratqs_glo = 0.
1412	ELSE
1413	#ifdef NC_DOUBLE
1414	ierr = NF_GET_VAR_DOUBLE(nid, nvarid, ratqs_glo)
1415	#else
1416	ierr = NF_GET_VAR_REAL(nid, nvarid, ratqs_glo)
1417	#endif
1418	IF (ierr.NE.NF_NOERR) THEN
1419	PRINT*, "phyetat0: Lecture echouee pour <RATQS>"
1420	CALL abort
1421	ENDIF
1422	ENDIF
1423	xmin = 1.0E+20
1424	xmax = -1.0E+20
1425	xmin = MINval(ratqs_glo)
1426	xmax = MAXval(ratqs_glo)
1427	PRINT*,'(ecart-type) ratqs:', xmin, xmax
1428	c
1429	c Lecture run_off_lic_0
1430	c
1431	ierr = NF_INQ_VARID (nid, "RUNOFFLIC0", nvarid)
1432	IF (ierr.NE.NF_NOERR) THEN
1433	PRINT*, "phyetat0: Le champ <RUNOFFLIC0> est absent"
1434	PRINT*, "Depart legerement fausse. Mais je continue"
1435	run_off_lic_0 = 0.
1436	ELSE
1437	#ifdef NC_DOUBLE
1438	ierr = NF_GET_VAR_DOUBLE(nid, nvarid, run_off_lic_0)
1439	#else
1440	ierr = NF_GET_VAR_REAL(nid, nvarid, run_off_lic_0)
1441	#endif
1442	IF (ierr.NE.NF_NOERR) THEN
1443	PRINT*, "phyetat0: Lecture echouee pour <RUNOFFLIC0>"
1444	CALL abort
1445	ENDIF
1446	ENDIF
1447	xmin = 1.0E+20
1448	xmax = -1.0E+20
1449	xmin = MINval(run_off_lic_0)
1450	xmax = MAXval(run_off_lic_0)
1451	PRINT*,'(ecart-type) run_off_lic_0:', xmin, xmax
1452
1453
1454	c Lecture de l'energie cinetique turbulente
1455	c
1456
1457	IF (iflag_pbl>1) then
1458	PRINT*, 'phyetat0: Le champ <TKE> est absent'
1459	PRINT, ' Mais je vais essayer de lire TKE*'
1460	DO nsrf = 1, nbsrf
1461	IF (nsrf.GT.99) THEN
1462	PRINT*, "Trop de sous-mailles"
1463	CALL abort
1464	ENDIF
1465	WRITE(str2,'(i2.2)') nsrf
1466	ierr = NF_INQ_VARID (nid, "TKE"//str2, nvarid)
1467	IF (ierr.NE.NF_NOERR) THEN
1468	PRINT*, "WARNING phyetat0: <TKE"//str2//"> est absent"
1469	pbl_tke_glo(:,:,nsrf)=1.e-8
1470	ELSE
1471	#ifdef NC_DOUBLE
1472	ierr = NF_GET_VAR_DOUBLE(nid, nvarid,
1473	& pbl_tke_glo(1,1,nsrf))
1474	#else
1475	ierr = NF_GET_VAR_REAL(nid, nvarid, pbl_tke_glo(1,1,nsrf))
1476	#endif
1477	IF (ierr.NE.NF_NOERR) THEN
1478	PRINT*, "WARNING phyetat0: echec lect <TKE"//str2//">"
1479	CALL abort
1480	ENDIF
1481	ENDIF
1482
1483	xmin = 1.0E+20
1484	xmax = -1.0E+20
1485	DO k = 1, klev
1486	DO i = 1, klon_glo
1487	xmin = MIN(pbl_tke_glo(i,k,nsrf),xmin)
1488	xmax = MAX(pbl_tke_glo(i,k,nsrf),xmax)
1489	ENDDO
1490	ENDDO
1491	PRINT,'Temperature du sol TKE*:', nsrf, xmin, xmax
1492	ENDDO
1493	ENDIF
1494	c
1495	c zmax0
1496	ierr = NF_INQ_VARID (nid, "ZMAX0", nvarid)
1497	IF (ierr.NE.NF_NOERR) THEN
1498	PRINT*, "phyetat0: Le champ <ZMAX0> est absent"
1499	PRINT*, "Depart legerement fausse. Mais je continue"
1500	zmax0_glo=40.
1501	ELSE
1502	#ifdef NC_DOUBLE
1503	ierr = NF_GET_VAR_DOUBLE(nid, nvarid, zmax0_glo)
1504	#else
1505	ierr = NF_GET_VAR_REAL(nid, nvarid, zmax0_glo)
1506	#endif
1507	IF (ierr.NE.NF_NOERR) THEN
1508	PRINT*, "phyetat0: Lecture echouee pour <ZMAX0>"
1509	CALL abort
1510	ENDIF
1511	ENDIF
1512	xmin = 1.0E+20
1513	xmax = -1.0E+20
1514	xmin = MINval(zmax0_glo)
1515	xmax = MAXval(zmax0_glo)
1516	PRINT*,'(ecart-type) zmax0:', xmin, xmax
1517	c
1518	c f0(ig)=1.e-5
1519	c f0
1520	ierr = NF_INQ_VARID (nid, "f0", nvarid)
1521	IF (ierr.NE.NF_NOERR) THEN
1522	PRINT*, "phyetat0: Le champ <f0> est absent"
1523	PRINT*, "Depart legerement fausse. Mais je continue"
1524	f0_glo=1.e-5
1525	ELSE
1526	#ifdef NC_DOUBLE
1527	ierr = NF_GET_VAR_DOUBLE(nid, nvarid, f0_glo)
1528	#else
1529	ierr = NF_GET_VAR_REAL(nid, nvarid, f0_glo)
1530	#endif
1531	IF (ierr.NE.NF_NOERR) THEN
1532	PRINT*, "phyetat0: Lecture echouee pour <f0>"
1533	CALL abort
1534	ENDIF
1535	ENDIF
1536	xmin = 1.0E+20
1537	xmax = -1.0E+20
1538	xmin = MINval(f0_glo)
1539	xmax = MAXval(f0_glo)
1540	PRINT*,'(ecart-type) f0:', xmin, xmax
1541	c
1542	c ema_work1
1543	c
1544	ierr = NF_INQ_VARID (nid, "EMA_WORK1", nvarid)
1545	IF (ierr.NE.NF_NOERR) THEN
1546	PRINT*, "phyetat0: Le champ <EMA_WORK1> est absent"
1547	PRINT*, "Depart legerement fausse. Mais je continue"
1548	ema_work1_glo=0.
1549	ELSE
1550	#ifdef NC_DOUBLE
1551	ierr = NF_GET_VAR_DOUBLE(nid, nvarid, ema_work1_glo)
1552	#else
1553	ierr = NF_GET_VAR_REAL(nid, nvarid, ema_work1_glo)
1554	#endif
1555	IF (ierr.NE.NF_NOERR) THEN
1556	PRINT*, "phyetat0: Lecture echouee pour <EMA_WORK1>"
1557	CALL abort
1558	ENDIF
1559	xmin = 1.0E+20
1560	xmax = -1.0E+20
1561	DO k = 1, klev
1562	DO i = 1, klon
1563	xmin = MIN(ema_work1_glo(i,k),xmin)
1564	xmax = MAX(ema_work1_glo(i,k),xmax)
1565	ENDDO
1566	ENDDO
1567	PRINT*,'ema_work1:', xmin, xmax
1568	ENDIF
1569	c
1570	c ema_work2
1571	c
1572	ierr = NF_INQ_VARID (nid, "EMA_WORK2", nvarid)
1573	IF (ierr.NE.NF_NOERR) THEN
1574	PRINT*, "phyetat0: Le champ <EMA_WORK2> est absent"
1575	PRINT*, "Depart legerement fausse. Mais je continue"
1576	ema_work2_glo=0.
1577	ELSE
1578	#ifdef NC_DOUBLE
1579	ierr = NF_GET_VAR_DOUBLE(nid, nvarid, ema_work2_glo)
1580	#else
1581	ierr = NF_GET_VAR_REAL(nid, nvarid, ema_work2_glo)
1582	#endif
1583	IF (ierr.NE.NF_NOERR) THEN
1584	PRINT*, "phyetat0: Lecture echouee pour <EMA_WORK2>"
1585	CALL abort
1586	ENDIF
1587	xmin = 1.0E+20
1588	xmax = -1.0E+20
1589	DO k = 1, klev
1590	DO i = 1, klon
1591	xmin = MIN(ema_work2_glo(i,k),xmin)
1592	xmax = MAX(ema_work2_glo(i,k),xmax)
1593	ENDDO
1594	ENDDO
1595	PRINT*,'ema_work2:', xmin, xmax
1596	ENDIF
1597	c
1598	c wake_deltat
1599	c
1600	ierr = NF_INQ_VARID (nid, "WAKE_DELTAT", nvarid)
1601	IF (ierr.NE.NF_NOERR) THEN
1602	PRINT*, "phyetat0: Le champ <WAKE_DELTAT> est absent"
1603	PRINT*, "Depart legerement fausse. Mais je continue"
1604	wake_deltat_glo=0.
1605	ELSE
1606	#ifdef NC_DOUBLE
1607	ierr = NF_GET_VAR_DOUBLE(nid, nvarid, wake_deltat_glo)
1608	#else
1609	ierr = NF_GET_VAR_REAL(nid, nvarid, wake_deltat_glo)
1610	#endif
1611	IF (ierr.NE.NF_NOERR) THEN
1612	PRINT*, "phyetat0: Lecture echouee pour <WAKE_DELTAT>"
1613	CALL abort
1614	ENDIF
1615	xmin = 1.0E+20
1616	xmax = -1.0E+20
1617	DO k = 1, klev
1618	DO i = 1, klon
1619	xmin = MIN(wake_deltat_glo(i,k),xmin)
1620	xmax = MAX(wake_deltat_glo(i,k),xmax)
1621	ENDDO
1622	ENDDO
1623	PRINT*,'wake_deltat:', xmin, xmax
1624	ENDIF
1625	c
1626	c wake_deltaq
1627	c
1628	ierr = NF_INQ_VARID (nid, "WAKE_DELTAQ", nvarid)
1629	IF (ierr.NE.NF_NOERR) THEN
1630	PRINT*, "phyetat0: Le champ <WAKE_DELTAQ> est absent"
1631	PRINT*, "Depart legerement fausse. Mais je continue"
1632	wake_deltaq_glo=0.
1633	ELSE
1634	#ifdef NC_DOUBLE
1635	ierr = NF_GET_VAR_DOUBLE(nid, nvarid, wake_deltaq_glo)
1636	#else
1637	ierr = NF_GET_VAR_REAL(nid, nvarid, wake_deltaq_glo)
1638	#endif
1639	IF (ierr.NE.NF_NOERR) THEN
1640	PRINT*, "phyetat0: Lecture echouee pour <WAKE_DELTAQ>"
1641	CALL abort
1642	ENDIF
1643	xmin = 1.0E+20
1644	xmax = -1.0E+20
1645	DO k = 1, klev
1646	DO i = 1, klon
1647	xmin = MIN(wake_deltaq_glo(i,k),xmin)
1648	xmax = MAX(wake_deltaq_glo(i,k),xmax)
1649	ENDDO
1650	ENDDO
1651	PRINT*,'wake_deltaq:', xmin, xmax
1652	ENDIF
1653	c
1654	c wake_s
1655	c
1656	ierr = NF_INQ_VARID (nid, "WAKE_S", nvarid)
1657	IF (ierr.NE.NF_NOERR) THEN
1658	PRINT*, "phyetat0: Le champ <WAKE_S> est absent"
1659	PRINT*, "Depart legerement fausse. Mais je continue"
1660	wake_s_glo=0.
1661	ELSE
1662	#ifdef NC_DOUBLE
1663	ierr = NF_GET_VAR_DOUBLE(nid, nvarid, wake_s_glo)
1664	#else
1665	ierr = NF_GET_VAR_REAL(nid, nvarid, wake_s_glo)
1666	#endif
1667	IF (ierr.NE.NF_NOERR) THEN
1668	PRINT*, "phyetat0: Lecture echouee pour <WAKE_S>"
1669	CALL abort
1670	ENDIF
1671	ENDIF
1672	xmin = 1.0E+20
1673	xmax = -1.0E+20
1674	xmin = MINval(wake_s_glo)
1675	xmax = MAXval(wake_s_glo)
1676	PRINT*,'(ecart-type) wake_s:', xmin, xmax
1677	c
1678	c wake_cstar
1679	c
1680	ierr = NF_INQ_VARID (nid, "WAKE_CSTAR", nvarid)
1681	IF (ierr.NE.NF_NOERR) THEN
1682	PRINT*, "phyetat0: Le champ <WAKE_CSTAR> est absent"
1683	PRINT*, "Depart legerement fausse. Mais je continue"
1684	wake_cstar_glo=0.
1685	ELSE
1686	#ifdef NC_DOUBLE
1687	ierr = NF_GET_VAR_DOUBLE(nid, nvarid, wake_cstar_glo)
1688	#else
1689	ierr = NF_GET_VAR_REAL(nid, nvarid, wake_cstar_glo)
1690	#endif
1691	IF (ierr.NE.NF_NOERR) THEN
1692	PRINT*, "phyetat0: Lecture echouee pour <WAKE_CSTAR>"
1693	CALL abort
1694	ENDIF
1695	ENDIF
1696	xmin = 1.0E+20
1697	xmax = -1.0E+20
1698	xmin = MINval(wake_cstar_glo)
1699	xmax = MAXval(wake_cstar_glo)
1700	PRINT*,'(ecart-type) wake_cstar:', xmin, xmax
1701	c
1702	c wake_fip
1703	c
1704	ierr = NF_INQ_VARID (nid, "WAKE_FIP", nvarid)
1705	IF (ierr.NE.NF_NOERR) THEN
1706	PRINT*, "phyetat0: Le champ <WAKE_FIP> est absent"
1707	PRINT*, "Depart legerement fausse. Mais je continue"
1708	wake_fip_glo=0.
1709	ELSE
1710	#ifdef NC_DOUBLE
1711	ierr = NF_GET_VAR_DOUBLE(nid, nvarid, wake_fip_glo)
1712	#else
1713	ierr = NF_GET_VAR_REAL(nid, nvarid, wake_fip_glo)
1714	#endif
1715	IF (ierr.NE.NF_NOERR) THEN
1716	PRINT*, "phyetat0: Lecture echouee pour <WAKE_FIP>"
1717	CALL abort
1718	ENDIF
1719	ENDIF
1720	xmin = 1.0E+20
1721	xmax = -1.0E+20
1722	xmin = MINval(wake_fip_glo)
1723	xmax = MAXval(wake_fip_glo)
1724	PRINT*,'(ecart-type) wake_fip:', xmin, xmax
1725	c
1726	c Fermer le fichier:
1727	c
1728	ierr = NF_CLOSE(nid)
1729	ENDIF ! is_mpi_root
1730	c
1731	c$OMP END MASTER
1732
1733	c$OMP MASTER
1734	cym en attendant mieux
1735	iolat(1)=rlat_glo(1)
1736
1737	!FH1D
1738	!iolat(jjm+1)=rlat(klon_glo)
1739	iolat(jjm+1-1/iim)=rlat_glo(klon_glo)
1740	if (iim.gt.1) then
1741	do i=2,jjm
1742	iolat(i)=rlat_glo(2+(i-2)*iim)
1743	enddo
1744	endif
1745
1746	CALL bcast_mpi(iolat)
1747	CALL bcast_mpi(rlon_glo)
1748
1749	!FH1D
1750	! call init_iophy(iolat,rlon(2:iim+1))
1751	call init_iophy(iolat,rlon_glo(2-1/iim:iim+1-1/iim))
1752
1753	c$OMP END MASTER
1754
1755	call Scatter( rlat_glo,rlat)
1756	call Scatter( rlon_glo,rlon)
1757	call Scatter( tsol_glo,ftsol)
1758	IF (iflag_pbl>1) then
1759	call Scatter( pbl_tke_glo,pbl_tke)
1760	endif
1761	call Scatter( zmax0_glo,zmax0)
1762	call Scatter( f0_glo,f0)
1763	call Scatter( ema_work1_glo, ema_work1)
1764	call Scatter( ema_work2_glo, ema_work2)
1765	call Scatter( wake_deltat_glo, wake_deltat)
1766	call Scatter( wake_deltaq_glo, wake_deltaq)
1767	call Scatter( wake_s_glo, wake_s)
1768	call Scatter( wake_cstar_glo, wake_cstar)
1769	call Scatter( wake_fip_glo, wake_fip)
1770	call Scatter( tsoil,tsoil_p)
1771	call Scatter( tslab,tslab_p)
1772	call Scatter( seaice,seaice_p)
1773	call Scatter( qsurf,qsurf_p)
1774	call Scatter( qsol,qsol_p)
1775	call Scatter( snow,snow_p)
1776	call Scatter( alb1_glo,falb1)
1777	call Scatter( alb2_glo,falb2)
1778	call Scatter( evap,evap_p)
1779	call Scatter( radsol_glo,radsol)
1780	call Scatter( rain_fall_glo,rain_fall)
1781	call Scatter( snow_fall_glo,snow_fall)
1782	call Scatter( sollw_glo,sollw)
1783	call Scatter( solsw_glo,solsw)
1784	call Scatter( fder,fder_p)
1785	call Scatter( frugs,frugs_p)
1786	call Scatter( agesno,agesno_p)
1787	call Scatter( zmea_glo,zmea)
1788	call Scatter( zstd_glo,zstd)
1789	call Scatter( zsig_glo,zsig)
1790	call Scatter( zgam_glo,zgam)
1791	call Scatter( zthe_glo,zthe)
1792	call Scatter( zpic_glo,zpic)
1793	call Scatter( zval_glo,zval)
1794	call Scatter( rugsrel_glo,rugoro)
1795	call Scatter( pctsrf_glo,pctsrf)
1796	call Scatter( run_off_lic_0,run_off_lic_0_p)
1797	call Scatter( t_ancien_glo,t_ancien)
1798	call Scatter( q_ancien_glo,q_ancien)
1799	call Scatter( rnebcon_glo,rnebcon)
1800	call Scatter( clwcon_glo,clwcon)
1801	call Scatter( ratqs_glo,ratqs)
1802	call Scatter( zmasq_glo,zmasq)
1803
1804	c
1805	c Initilalize variables in module surface_data
1806	c
1807	ok_veget = ok_veget_in
1808	ocean = ocean_in
1809	c
1810	c Initialize module pbl_surface_mod
1811	c
1812	CALL pbl_surface_init(qsol_p, fder_p, snow_p, qsurf_p,
1813	$ evap_p, frugs_p, agesno_p, tsoil_p)
1814
1815	c Initialize ocean module according to ocean type
1816	IF ( ocean == 'slab' ) THEN
1817	c initilalize module ocean_slab_init
1818	CALL ocean_slab_init(dtime, tslab_p, seaice_p, pctsrf)
1819	ELSEIF ( ocean == 'couple' ) THEN
1820	c initilalize module ocean_cpl_init
1821	CALL ocean_cpl_init(dtime, rlon, rlat)
1822	ELSE
1823	c initilalize module ocean_forced_init
1824	CALL ocean_forced_init
1825	ENDIF
1826	c
1827	c Initilialize module fonte_neige_mod
1828	c
1829	CALL fonte_neige_init(run_off_lic_0_p)
1830
1831
1832	RETURN
1833	END

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

Download in other formats:

Original Format