Context Navigation

phyetat0.F @ 1068

Last change on this file since 1068 was 956, checked in by lmdzadmin, 17 years ago

Nettoyage du controle des parametres physiques. FH

Les parametres cycle_diurne, soil_model, new_oliq, ok_orodr, ok_orolf, ok_limitvrai, nbapp_rad et iflag_con
sont maintenant geres par la physique uniquement.
ecritphy est elimine.
dimphy.F90 et clesphys.h ne sont plus utilises par le code dynamique.
Le test academique obtenu en compilant avec
makegcm -p nophys gcm
fonctionne. FH
IM

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

File size: 45.4 KB

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

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

Download in other formats:

Original Format