Context Navigation

phyetat0.F @ 978

Last change on this file since 978 was 978, checked in by Laurent Fairhead, 16 years ago
Encore quelques bugs consecutifs au renommage des variables YM LF
Property svn:eol-style set to `native` Property svn:keywords set to `Author Date Id Revision`
File size: 54.0 KB

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

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

Download in other formats:

Original Format