Context Navigation

phyetat0.F @ 3733

Last change on this file since 3733 was 704, checked in by Laurent Fairhead, 18 years ago
Inclusion des modifs de Y. Meurdesoif pour la version V3 LF
Property svn:eol-style set to `native` Property svn:keywords set to `Author Date Id Revision`
File size: 50.3 KB

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

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

Download in other formats:

Original Format