Context Navigation

phyetat0.F @ 947

Last change on this file since 947 was 937, checked in by lmdzadmin, 17 years ago
Ajout variables convection (ema_work1, ema_work2) dans startphy.nc IM
Property svn:eol-style set to `native` Property svn:keywords set to `Author Date Id Revision`
File size: 50.2 KB

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

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

Download in other formats:

Original Format