Context Navigation

phyetat0.F @ 1955

Last change on this file since 1955 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

Rev	Line
[524]	1	!
	2	! $Header$
	3	!
	4	c
	5	c
	6	SUBROUTINE phyetat0 (fichnom,dtime,co2_ppm_etat0,solaire_etat0,
[704]	7	. rlat_p,rlon_p, pctsrf_p, tsol_p,tsoil_p,
[644]	8	cIM "slab" ocean
[704]	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
[524]	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"
[704]	25	cym#include "dimphy.h"
[524]	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
[704]	35	REAL rlat_p(klon2), rlon_p(klon2)
[524]	36	REAL co2_ppm_etat0
	37	REAL solaire_etat0
[704]	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)
[524]	76	REAL tsol(klon,nbsrf)
	77	REAL tsoil(klon,nsoilmx,nbsrf)
[644]	78	cIM "slab" ocean
	79	REAL tslab(klon), seaice(klon)
[524]	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)
[704]	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)
[524]	110
[704]	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
[524]	147
[651]	148	CHARACTER*6 ocean
	149
[524]	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)
[704]	161	REAL,SAVE :: tab_cntrl_omp(length)
[524]	162	CHARACTER*7 str7
	163	CHARACTER*2 str2
[704]	164	real iolat(jjm+1)
[524]	165	c
	166	c Ouvrir le fichier contenant l'etat initial:
	167	c
[704]	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
[524]	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
[704]	181	c$OMP END MASTER
[524]	182	c
	183	c Lecture des parametres de controle:
	184	c
[704]	185	c$OMP MASTER
[524]	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
[704]	192	ierr = NF_GET_VAR_DOUBLE(nid, nvarid, tab_cntrl_omp)
[524]	193	#else
[704]	194	ierr = NF_GET_VAR_REAL(nid, nvarid, tab_cntrl_omp)
[524]	195	#endif
	196	IF (ierr.NE.NF_NOERR) THEN
	197	PRINT*, 'phyetat0: Lecture echouee pour <controle>'
	198	CALL abort
[704]	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
[524]	207	c
[704]	208
[524]	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
[704]	282	c$OMP MASTER
[524]	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
[704]	297
[524]	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
[704]	319
[524]	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
[704]	336
	337
[524]	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
[704]	346
[524]	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'
[704]	360	c@$$ CALL abort
[524]	361	ENDIF
[704]	362
[524]	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'
[704]	379	c@$$ CALL abort
[524]	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'
[704]	397	c@$$ CALL abort
[524]	398	ENDIF
[704]	399
[524]	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'
[704]	416	c@$$ CALL abort
[524]	417	ENDIF
[704]	418
[524]	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
[704]	440
[524]	441	C
	442	c Lecture des temperatures du sol:
	443	c
[704]	444
[524]	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
[704]	469
[524]	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
[704]	503
[524]	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
[644]	535	cIM "slab" ocean
[524]	536	c
[644]	537	c Lecture de tslab (pour slab ocean seulement):
	538	c
[651]	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
[524]	545	#ifdef NC_DOUBLE
[651]	546	ierr = NF_GET_VAR_DOUBLE(nid, nvarid, tslab)
[524]	547	#else
[651]	548	ierr = NF_GET_VAR_REAL(nid, nvarid, tslab)
[524]	549	#endif
[651]	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
[689]	560	PRINT*,'Min, Max tslab (utilise si OCEAN=slab )', xmin, xmax
[524]	561	c
[644]	562	c Lecture de seaice (pour slab ocean seulement):
	563	c
[651]	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
[644]	569	#ifdef NC_DOUBLE
[651]	570	ierr = NF_GET_VAR_DOUBLE(nid, nvarid, seaice)
[644]	571	#else
[651]	572	ierr = NF_GET_VAR_REAL(nid, nvarid, seaice)
[644]	573	#endif
[651]	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
[689]	584	PRINT*,'Masse de la glace de mer (utilise si OCEAN=slab)',
	585	$ xmin, xmax
[651]	586	ELSE
	587	tslab = 0.
	588	seaice = 0.
[644]	589	ENDIF
	590	c
[524]	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.
[704]	668	c@$$ CALL abort
[524]	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
[704]	997	c$OMP END MASTER
	998
	999
	1000	c$OMP MASTER
	1001
[524]	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
[704]	1396	clwcon=0.
[524]	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
[704]	1419	rnebcon=0.
[524]	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
[704]	1462	ratqs=0.
[524]	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
[704]	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
[524]	1679	RETURN
	1680	END

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

Download in other formats:

Original Format