Context Navigation

phyetat0.F @ 996

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

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

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

File size: 51.6 KB

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

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

Download in other formats:

Original Format