Context Navigation

phyetat0.F @ 996

Last change on this file since 996 was 996, checked in by lsce, 16 years ago

Modifications liées au calcul des nouveau sous-fractions
Nettoyage de ocean slab : il reste uniquement la version avec glace de mer forcé
Nouveaux variables pour distiguer la version et type d'ocean : type_ocean=force/slab/couple, version_ocean=opa8/nemo pour couplé ou version_ocean=sicOBS pour slab

Property svn:eol-style set to native
Property svn:keywords set to Author Date Id Revision

File size: 51.6 KB

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

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

Download in other formats:

Original Format