Context Navigation

mod_1D_cases_read.F90 @ 2221

Last change on this file since 2221 was 2191, checked in by fhourdin, 10 years ago
Ajout du cas 1D CINDY-DYNAMO, utilisant le nouveau format standard, amené à être étendu aux autres cas. Addition of the CINDY-DYNAMO 1D case, using the new standard format for 1D cases, that will be extended to all 1D cases. Catherine Rio
File size: 33.7 KB

Rev	Line
[2118]	1	MODULE mod_1D_cases_read
	2
	3	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
[2191]	4	!Declarations specifiques au cas standard
	5	character*80 :: fich_cas
	6	! Discrétisation
	7	integer nlev_cas, nt_cas
[2118]	8
	9
[2191]	10	integer year_ini_cas, day_ini_cas, mth_ini_cas
	11	real heure_ini_cas
	12	real day_ju_ini_cas ! Julian day of case first day
	13	parameter (year_ini_cas=2011)
	14	parameter (mth_ini_cas=10)
	15	parameter (day_ini_cas=1) ! 10 = 10Juil2006
	16	parameter (heure_ini_cas=0.) !0h en secondes
	17	real pdt_cas
	18	parameter (pdt_cas=3.*3600)
[2118]	19
[2191]	20	!CR ATTENTION TEST AMMA
	21	! parameter (year_ini_cas=2006)
	22	! parameter (mth_ini_cas=7)
	23	! parameter (day_ini_cas=10) ! 10 = 10Juil2006
	24	! parameter (heure_ini_cas=0.) !0h en secondes
	25	! parameter (pdt_cas=1800.)
[2118]	26
[2191]	27	!profils environnementaux
	28	real, allocatable:: plev_cas(:,:)
[2118]	29
[2191]	30	real, allocatable:: z_cas(:,:)
	31	real, allocatable:: t_cas(:,:),q_cas(:,:),rh_cas(:,:)
	32	real, allocatable:: th_cas(:,:),rv_cas(:,:)
	33	real, allocatable:: u_cas(:,:)
	34	real, allocatable:: v_cas(:,:)
[2118]	35
[2191]	36	!forcing
	37	real, allocatable:: ht_cas(:,:),vt_cas(:,:),dt_cas(:,:),dtrad_cas(:,:)
	38	real, allocatable:: hth_cas(:,:),vth_cas(:,:),dth_cas(:,:)
	39	real, allocatable:: hq_cas(:,:),vq_cas(:,:),dq_cas(:,:)
	40	real, allocatable:: hr_cas(:,:),vr_cas(:,:),dr_cas(:,:)
	41	real, allocatable:: hu_cas(:,:),vu_cas(:,:),du_cas(:,:)
	42	real, allocatable:: hv_cas(:,:),vv_cas(:,:),dv_cas(:,:)
	43	real, allocatable:: vitw_cas(:,:)
	44	real, allocatable:: ug_cas(:,:),vg_cas(:,:)
	45	real, allocatable:: lat_cas(:),sens_cas(:),ts_cas(:)
[2118]	46
	47	!champs interpoles
[2191]	48	real, allocatable:: plev_prof_cas(:)
	49	real, allocatable:: t_prof_cas(:)
	50	real, allocatable:: q_prof_cas(:)
	51	real, allocatable:: u_prof_cas(:)
	52	real, allocatable:: v_prof_cas(:)
[2118]	53
[2191]	54	real, allocatable:: vitw_prof_cas(:)
	55	real, allocatable:: ug_prof_cas(:)
	56	real, allocatable:: vg_prof_cas(:)
	57	real, allocatable:: ht_prof_cas(:)
	58	real, allocatable:: hq_prof_cas(:)
	59	real, allocatable:: vt_prof_cas(:)
	60	real, allocatable:: vq_prof_cas(:)
	61	real, allocatable:: dt_prof_cas(:)
	62	real, allocatable:: dtrad_prof_cas(:)
	63	real, allocatable:: dq_prof_cas(:)
	64	real, allocatable:: hu_prof_cas(:)
	65	real, allocatable:: hv_prof_cas(:)
	66	real, allocatable:: vu_prof_cas(:)
	67	real, allocatable:: vv_prof_cas(:)
	68	real, allocatable:: du_prof_cas(:)
	69	real, allocatable:: dv_prof_cas(:)
[2118]	70
[2191]	71
	72	real lat_prof_cas,sens_prof_cas,ts_prof_cas
	73
	74
	75
[2118]	76	CONTAINS
	77
[2191]	78	SUBROUTINE read_1D_cas
[2118]	79	implicit none
	80
	81	#include "netcdf.inc"
	82
	83	INTEGER nid,rid,ierr
[2191]	84	INTEGER ii,jj
[2118]	85
[2191]	86	fich_cas='setup/cas.nc'
	87	print*,'fich_cas ',fich_cas
	88	ierr = NF_OPEN(fich_cas,NF_NOWRITE,nid)
	89	print*,'fich_cas,NF_NOWRITE,nid ',fich_cas,NF_NOWRITE,nid
[2118]	90	if (ierr.NE.NF_NOERR) then
	91	write(,) 'ERROR: GROS Pb opening forcings nc file '
	92	write(,) NF_STRERROR(ierr)
	93	stop ""
	94	endif
	95	!.......................................................................
[2191]	96	ierr=NF_INQ_DIMID(nid,'lat',rid)
	97	IF (ierr.NE.NF_NOERR) THEN
	98	print*, 'Oh probleme lecture dimension lat'
	99	ENDIF
	100	ierr=NF_INQ_DIMLEN(nid,rid,ii)
	101	print*,'OK nid,rid,lat',nid,rid,ii
	102	!.......................................................................
	103	ierr=NF_INQ_DIMID(nid,'lon',rid)
	104	IF (ierr.NE.NF_NOERR) THEN
	105	print*, 'Oh probleme lecture dimension lon'
	106	ENDIF
	107	ierr=NF_INQ_DIMLEN(nid,rid,jj)
	108	print*,'OK nid,rid,lat',nid,rid,jj
	109	!.......................................................................
[2118]	110	ierr=NF_INQ_DIMID(nid,'lev',rid)
	111	IF (ierr.NE.NF_NOERR) THEN
	112	print*, 'Oh probleme lecture dimension zz'
	113	ENDIF
[2191]	114	ierr=NF_INQ_DIMLEN(nid,rid,nlev_cas)
	115	print*,'OK nid,rid,nlev_cas',nid,rid,nlev_cas
[2118]	116	!.......................................................................
	117	ierr=NF_INQ_DIMID(nid,'time',rid)
	118	print*,'nid,rid',nid,rid
[2191]	119	nt_cas=0
[2118]	120	IF (ierr.NE.NF_NOERR) THEN
	121	stop 'probleme lecture dimension sens'
	122	ENDIF
[2191]	123	ierr=NF_INQ_DIMLEN(nid,rid,nt_cas)
	124	print*,'nid,rid,nlev_cas',nid,rid,nt_cas
[2118]	125
	126	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
[2191]	127	!profils moyens:
	128	allocate(plev_cas(nlev_cas,nt_cas))
	129	allocate(z_cas(nlev_cas,nt_cas))
	130	allocate(t_cas(nlev_cas,nt_cas),q_cas(nlev_cas,nt_cas),rh_cas(nlev_cas,nt_cas))
	131	allocate(th_cas(nlev_cas,nt_cas),rv_cas(nlev_cas,nt_cas))
	132	allocate(u_cas(nlev_cas,nt_cas))
	133	allocate(v_cas(nlev_cas,nt_cas))
[2118]	134
[2191]	135	!forcing
	136	allocate(ht_cas(nlev_cas,nt_cas),vt_cas(nlev_cas,nt_cas),dt_cas(nlev_cas,nt_cas),dtrad_cas(nlev_cas,nt_cas))
	137	allocate(hq_cas(nlev_cas,nt_cas),vq_cas(nlev_cas,nt_cas),dq_cas(nlev_cas,nt_cas))
	138	allocate(hth_cas(nlev_cas,nt_cas),vth_cas(nlev_cas,nt_cas),dth_cas(nlev_cas,nt_cas))
	139	allocate(hr_cas(nlev_cas,nt_cas),vr_cas(nlev_cas,nt_cas),dr_cas(nlev_cas,nt_cas))
	140	allocate(hu_cas(nlev_cas,nt_cas),vu_cas(nlev_cas,nt_cas),du_cas(nlev_cas,nt_cas))
	141	allocate(hv_cas(nlev_cas,nt_cas),vv_cas(nlev_cas,nt_cas),dv_cas(nlev_cas,nt_cas))
	142	allocate(vitw_cas(nlev_cas,nt_cas))
	143	allocate(ug_cas(nlev_cas,nt_cas))
	144	allocate(vg_cas(nlev_cas,nt_cas))
	145	allocate(lat_cas(nt_cas),sens_cas(nt_cas),ts_cas(nt_cas))
[2118]	146
	147
	148	!champs interpoles
[2191]	149	allocate(plev_prof_cas(nlev_cas))
	150	allocate(t_prof_cas(nlev_cas))
	151	allocate(q_prof_cas(nlev_cas))
	152	allocate(u_prof_cas(nlev_cas))
	153	allocate(v_prof_cas(nlev_cas))
[2118]	154
[2191]	155	allocate(vitw_prof_cas(nlev_cas))
	156	allocate(ug_prof_cas(nlev_cas))
	157	allocate(vg_prof_cas(nlev_cas))
	158	allocate(ht_prof_cas(nlev_cas))
	159	allocate(hq_prof_cas(nlev_cas))
	160	allocate(hu_prof_cas(nlev_cas))
	161	allocate(hv_prof_cas(nlev_cas))
	162	allocate(vt_prof_cas(nlev_cas))
	163	allocate(vq_prof_cas(nlev_cas))
	164	allocate(vu_prof_cas(nlev_cas))
	165	allocate(vv_prof_cas(nlev_cas))
	166	allocate(dt_prof_cas(nlev_cas))
	167	allocate(dtrad_prof_cas(nlev_cas))
	168	allocate(dq_prof_cas(nlev_cas))
	169	allocate(du_prof_cas(nlev_cas))
	170	allocate(dv_prof_cas(nlev_cas))
	171
[2118]	172	print*,'Allocations OK'
[2191]	173	call read_cas(nid,nlev_cas,nt_cas &
	174	& ,z_cas,plev_cas,t_cas,q_cas,rh_cas,th_cas,rv_cas,u_cas,v_cas &
	175	& ,ug_cas,vg_cas,vitw_cas,du_cas,hu_cas,vu_cas,dv_cas,hv_cas,vv_cas &
	176	& ,dt_cas,dtrad_cas,ht_cas,vt_cas,dq_cas,hq_cas,vq_cas &
	177	& ,dth_cas,hth_cas,vth_cas,dr_cas,hr_cas,vr_cas,sens_cas,lat_cas,ts_cas)
[2118]	178
	179
[2191]	180	END SUBROUTINE read_1D_cas
[2118]	181
	182
[2191]	183
[2118]	184	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
	185	SUBROUTINE deallocate_1D_cases
[2191]	186	!profils environnementaux:
	187	deallocate(plev_cas)
[2118]	188
[2191]	189	deallocate(z_cas)
	190	deallocate(t_cas,q_cas,rh_cas)
	191	deallocate(th_cas,rv_cas)
	192	deallocate(u_cas)
	193	deallocate(v_cas)
[2118]	194
[2191]	195	!forcing
	196	deallocate(ht_cas,vt_cas,dt_cas,dtrad_cas)
	197	deallocate(hq_cas,vq_cas,dq_cas)
	198	deallocate(hth_cas,vth_cas,dth_cas)
	199	deallocate(hr_cas,vr_cas,dr_cas)
	200	deallocate(hu_cas,vu_cas,du_cas)
	201	deallocate(hv_cas,vv_cas,dv_cas)
	202	deallocate(vitw_cas)
	203	deallocate(ug_cas)
	204	deallocate(vg_cas)
	205	deallocate(lat_cas,sens_cas,ts_cas)
[2118]	206
	207	!champs interpoles
[2191]	208	deallocate(plev_prof_cas)
	209	deallocate(t_prof_cas)
	210	deallocate(q_prof_cas)
	211	deallocate(u_prof_cas)
	212	deallocate(v_prof_cas)
	213
	214	deallocate(vitw_prof_cas)
	215	deallocate(ug_prof_cas)
	216	deallocate(vg_prof_cas)
	217	deallocate(ht_prof_cas)
	218	deallocate(hq_prof_cas)
	219	deallocate(hu_prof_cas)
	220	deallocate(hv_prof_cas)
	221	deallocate(vt_prof_cas)
	222	deallocate(vq_prof_cas)
	223	deallocate(vu_prof_cas)
	224	deallocate(vv_prof_cas)
	225	deallocate(dt_prof_cas)
	226	deallocate(dtrad_prof_cas)
	227	deallocate(dq_prof_cas)
	228	deallocate(du_prof_cas)
	229	deallocate(dv_prof_cas)
	230	deallocate(t_prof_cas)
	231	deallocate(q_prof_cas)
	232	deallocate(u_prof_cas)
	233	deallocate(v_prof_cas)
	234
[2118]	235	END SUBROUTINE deallocate_1D_cases
	236
	237
	238	END MODULE mod_1D_cases_read
	239	!=====================================================================
[2191]	240	subroutine read_cas(nid,nlevel,ntime &
	241	& ,zz,pp,temp,qv,rh,theta,rv,u,v,ug,vg,w, &
	242	& du,hu,vu,dv,hv,vv,dt,dtrad,ht,vt,dq,hq,vq, &
	243	& dth,hth,vth,dr,hr,vr,sens,flat,ts)
[2118]	244
[2191]	245	!program reading forcing of the case study
[2118]	246	implicit none
	247	#include "netcdf.inc"
	248
	249	integer ntime,nlevel
	250
[2191]	251	real zz(nlevel,ntime)
	252	real pp(nlevel,ntime)
	253	real temp(nlevel,ntime),qv(nlevel,ntime),rh(nlevel,ntime)
	254	real theta(nlevel,ntime),rv(nlevel,ntime)
	255	real u(nlevel,ntime)
	256	real v(nlevel,ntime)
	257	real ug(nlevel,ntime)
	258	real vg(nlevel,ntime)
	259	real w(nlevel,ntime)
	260	real du(nlevel,ntime),hu(nlevel,ntime),vu(nlevel,ntime)
	261	real dv(nlevel,ntime),hv(nlevel,ntime),vv(nlevel,ntime)
	262	real dt(nlevel,ntime),ht(nlevel,ntime),vt(nlevel,ntime)
	263	real dtrad(nlevel,ntime)
	264	real dq(nlevel,ntime),hq(nlevel,ntime),vq(nlevel,ntime)
	265	real dth(nlevel,ntime),hth(nlevel,ntime),vth(nlevel,ntime)
	266	real dr(nlevel,ntime),hr(nlevel,ntime),vr(nlevel,ntime)
	267	real flat(ntime),sens(ntime),ts(ntime)
[2118]	268
	269
	270	integer nid, ierr,rid
	271	integer nbvar3d
[2191]	272	parameter(nbvar3d=34)
[2118]	273	integer var3didin(nbvar3d)
	274
	275	ierr=NF_INQ_VARID(nid,"zz",var3didin(1))
	276	if(ierr/=NF_NOERR) then
	277	write(,) NF_STRERROR(ierr)
	278	stop 'lev'
	279	endif
[2191]	280
	281	ierr=NF_INQ_VARID(nid,"pp",var3didin(2))
	282	if(ierr/=NF_NOERR) then
	283	write(,) NF_STRERROR(ierr)
	284	stop 'plev'
	285	endif
[2118]	286
	287
[2191]	288	ierr=NF_INQ_VARID(nid,"temp",var3didin(3))
[2118]	289	if(ierr/=NF_NOERR) then
	290	write(,) NF_STRERROR(ierr)
	291	stop 'temp'
	292	endif
	293
[2191]	294	ierr=NF_INQ_VARID(nid,"qv",var3didin(4))
[2118]	295	if(ierr/=NF_NOERR) then
	296	write(,) NF_STRERROR(ierr)
	297	stop 'qv'
	298	endif
	299
[2191]	300	ierr=NF_INQ_VARID(nid,"rh",var3didin(5))
[2118]	301	if(ierr/=NF_NOERR) then
	302	write(,) NF_STRERROR(ierr)
[2191]	303	stop 'rh'
	304	endif
	305
	306	ierr=NF_INQ_VARID(nid,"theta",var3didin(6))
	307	if(ierr/=NF_NOERR) then
	308	write(,) NF_STRERROR(ierr)
	309	stop 'theta'
	310	endif
	311
	312	ierr=NF_INQ_VARID(nid,"rv",var3didin(7))
	313	if(ierr/=NF_NOERR) then
	314	write(,) NF_STRERROR(ierr)
	315	stop 'rv'
	316	endif
	317
	318
	319	ierr=NF_INQ_VARID(nid,"u",var3didin(8))
	320	if(ierr/=NF_NOERR) then
	321	write(,) NF_STRERROR(ierr)
[2118]	322	stop 'u'
	323	endif
	324
[2191]	325	ierr=NF_INQ_VARID(nid,"v",var3didin(9))
[2118]	326	if(ierr/=NF_NOERR) then
	327	write(,) NF_STRERROR(ierr)
	328	stop 'v'
	329	endif
	330
[2191]	331	ierr=NF_INQ_VARID(nid,"ug",var3didin(10))
[2118]	332	if(ierr/=NF_NOERR) then
	333	write(,) NF_STRERROR(ierr)
[2191]	334	stop 'ug'
[2118]	335	endif
	336
[2191]	337	ierr=NF_INQ_VARID(nid,"vg",var3didin(11))
[2118]	338	if(ierr/=NF_NOERR) then
	339	write(,) NF_STRERROR(ierr)
[2191]	340	stop 'vg'
[2118]	341	endif
	342
[2191]	343	ierr=NF_INQ_VARID(nid,"w",var3didin(12))
[2118]	344	if(ierr/=NF_NOERR) then
	345	write(,) NF_STRERROR(ierr)
[2191]	346	stop 'w'
[2118]	347	endif
[2191]	348
	349	ierr=NF_INQ_VARID(nid,"advu",var3didin(13))
	350	if(ierr/=NF_NOERR) then
	351	write(,) NF_STRERROR(ierr)
	352	stop 'advu'
	353	endif
	354
	355	ierr=NF_INQ_VARID(nid,"hu",var3didin(14))
	356	if(ierr/=NF_NOERR) then
	357	write(,) NF_STRERROR(ierr)
	358	stop 'hu'
	359	endif
	360
	361	ierr=NF_INQ_VARID(nid,"vu",var3didin(15))
	362	if(ierr/=NF_NOERR) then
	363	write(,) NF_STRERROR(ierr)
	364	stop 'vu'
	365	endif
	366
	367	ierr=NF_INQ_VARID(nid,"advv",var3didin(16))
	368	if(ierr/=NF_NOERR) then
	369	write(,) NF_STRERROR(ierr)
	370	stop 'advv'
	371	endif
	372
	373	ierr=NF_INQ_VARID(nid,"hv",var3didin(17))
	374	if(ierr/=NF_NOERR) then
	375	write(,) NF_STRERROR(ierr)
	376	stop 'hv'
	377	endif
	378
	379	ierr=NF_INQ_VARID(nid,"vv",var3didin(18))
	380	if(ierr/=NF_NOERR) then
	381	write(,) NF_STRERROR(ierr)
	382	stop 'vv'
	383	endif
	384
	385	ierr=NF_INQ_VARID(nid,"advT",var3didin(19))
	386	if(ierr/=NF_NOERR) then
	387	write(,) NF_STRERROR(ierr)
	388	stop 'advT'
	389	endif
	390
	391	ierr=NF_INQ_VARID(nid,"hT",var3didin(20))
	392	if(ierr/=NF_NOERR) then
	393	write(,) NF_STRERROR(ierr)
	394	stop 'hT'
	395	endif
	396
	397	ierr=NF_INQ_VARID(nid,"vT",var3didin(21))
	398	if(ierr/=NF_NOERR) then
	399	write(,) NF_STRERROR(ierr)
	400	stop 'vT'
	401	endif
	402
	403	ierr=NF_INQ_VARID(nid,"advq",var3didin(22))
	404	if(ierr/=NF_NOERR) then
	405	write(,) NF_STRERROR(ierr)
	406	stop 'advq'
	407	endif
[2118]	408
[2191]	409	ierr=NF_INQ_VARID(nid,"hq",var3didin(23))
[2118]	410	if(ierr/=NF_NOERR) then
	411	write(,) NF_STRERROR(ierr)
[2191]	412	stop 'hq'
	413	endif
	414
	415	ierr=NF_INQ_VARID(nid,"vq",var3didin(24))
	416	if(ierr/=NF_NOERR) then
	417	write(,) NF_STRERROR(ierr)
	418	stop 'vq'
	419	endif
	420
	421	ierr=NF_INQ_VARID(nid,"advth",var3didin(25))
	422	if(ierr/=NF_NOERR) then
	423	write(,) NF_STRERROR(ierr)
	424	stop 'advth'
	425	endif
	426
	427	ierr=NF_INQ_VARID(nid,"hth",var3didin(26))
	428	if(ierr/=NF_NOERR) then
	429	write(,) NF_STRERROR(ierr)
	430	stop 'hth'
	431	endif
	432
	433	ierr=NF_INQ_VARID(nid,"vth",var3didin(27))
	434	if(ierr/=NF_NOERR) then
	435	write(,) NF_STRERROR(ierr)
	436	stop 'vth'
	437	endif
	438
	439	ierr=NF_INQ_VARID(nid,"advr",var3didin(28))
	440	if(ierr/=NF_NOERR) then
	441	write(,) NF_STRERROR(ierr)
	442	stop 'advr'
	443	endif
	444
	445	ierr=NF_INQ_VARID(nid,"hr",var3didin(29))
	446	if(ierr/=NF_NOERR) then
	447	write(,) NF_STRERROR(ierr)
	448	stop 'hr'
	449	endif
	450
	451	ierr=NF_INQ_VARID(nid,"vr",var3didin(30))
	452	if(ierr/=NF_NOERR) then
	453	write(,) NF_STRERROR(ierr)
	454	stop 'vr'
	455	endif
	456
	457	ierr=NF_INQ_VARID(nid,"radT",var3didin(31))
	458	if(ierr/=NF_NOERR) then
	459	write(,) NF_STRERROR(ierr)
	460	stop 'radT'
	461	endif
	462
	463	ierr=NF_INQ_VARID(nid,"sens",var3didin(32))
	464	if(ierr/=NF_NOERR) then
	465	write(,) NF_STRERROR(ierr)
[2118]	466	stop 'sens'
	467	endif
	468
[2191]	469	ierr=NF_INQ_VARID(nid,"flat",var3didin(33))
[2118]	470	if(ierr/=NF_NOERR) then
	471	write(,) NF_STRERROR(ierr)
	472	stop 'flat'
	473	endif
	474
[2191]	475	ierr=NF_INQ_VARID(nid,"ts",var3didin(34))
[2118]	476	if(ierr/=NF_NOERR) then
	477	write(,) NF_STRERROR(ierr)
[2191]	478	stop 'ts'
	479	endif
[2118]	480
	481	#ifdef NC_DOUBLE
	482	ierr = NF_GET_VAR_DOUBLE(nid,var3didin(1),zz)
	483	#else
	484	ierr = NF_GET_VAR_REAL(nid,var3didin(1),zz)
	485	#endif
	486	if(ierr/=NF_NOERR) then
	487	write(,) NF_STRERROR(ierr)
	488	stop "getvarup"
	489	endif
	490	! write(,)'lecture z ok',zz
	491
	492	#ifdef NC_DOUBLE
[2191]	493	ierr = NF_GET_VAR_DOUBLE(nid,var3didin(2),pp)
[2118]	494	#else
[2191]	495	ierr = NF_GET_VAR_REAL(nid,var3didin(2),pp)
[2118]	496	#endif
	497	if(ierr/=NF_NOERR) then
	498	write(,) NF_STRERROR(ierr)
	499	stop "getvarup"
	500	endif
[2191]	501	! write(,)'lecture pp ok',pp
[2118]	502
[2191]	503
[2118]	504	#ifdef NC_DOUBLE
[2191]	505	ierr = NF_GET_VAR_DOUBLE(nid,var3didin(3),temp)
[2118]	506	#else
[2191]	507	ierr = NF_GET_VAR_REAL(nid,var3didin(3),temp)
[2118]	508	#endif
	509	if(ierr/=NF_NOERR) then
	510	write(,) NF_STRERROR(ierr)
	511	stop "getvarup"
	512	endif
[2191]	513	! write(,)'lecture T ok',temp
	514
	515	#ifdef NC_DOUBLE
	516	ierr = NF_GET_VAR_DOUBLE(nid,var3didin(4),qv)
	517	#else
	518	ierr = NF_GET_VAR_REAL(nid,var3didin(4),qv)
	519	#endif
	520	if(ierr/=NF_NOERR) then
	521	write(,) NF_STRERROR(ierr)
	522	stop "getvarup"
	523	endif
[2118]	524	! write(,)'lecture qv ok',qv
	525
	526	#ifdef NC_DOUBLE
[2191]	527	ierr = NF_GET_VAR_DOUBLE(nid,var3didin(5),rh)
[2118]	528	#else
[2191]	529	ierr = NF_GET_VAR_REAL(nid,var3didin(5),rh)
[2118]	530	#endif
	531	if(ierr/=NF_NOERR) then
	532	write(,) NF_STRERROR(ierr)
	533	stop "getvarup"
	534	endif
[2191]	535	! write(,)'lecture rh ok',rh
	536
	537	#ifdef NC_DOUBLE
	538	ierr = NF_GET_VAR_DOUBLE(nid,var3didin(6),theta)
	539	#else
	540	ierr = NF_GET_VAR_REAL(nid,var3didin(6),theta)
	541	#endif
	542	if(ierr/=NF_NOERR) then
	543	write(,) NF_STRERROR(ierr)
	544	stop "getvarup"
	545	endif
	546	! write(,)'lecture theta ok',theta
	547
	548	#ifdef NC_DOUBLE
	549	ierr = NF_GET_VAR_DOUBLE(nid,var3didin(7),rv)
	550	#else
	551	ierr = NF_GET_VAR_REAL(nid,var3didin(7),rv)
	552	#endif
	553	if(ierr/=NF_NOERR) then
	554	write(,) NF_STRERROR(ierr)
	555	stop "getvarup"
	556	endif
	557	! write(,)'lecture rv ok',rv
	558
	559	#ifdef NC_DOUBLE
	560	ierr = NF_GET_VAR_DOUBLE(nid,var3didin(8),u)
	561	#else
	562	ierr = NF_GET_VAR_REAL(nid,var3didin(8),u)
	563	#endif
	564	if(ierr/=NF_NOERR) then
	565	write(,) NF_STRERROR(ierr)
	566	stop "getvarup"
	567	endif
[2118]	568	! write(,)'lecture u ok',u
	569
	570	#ifdef NC_DOUBLE
[2191]	571	ierr = NF_GET_VAR_DOUBLE(nid,var3didin(9),v)
[2118]	572	#else
[2191]	573	ierr = NF_GET_VAR_REAL(nid,var3didin(9),v)
[2118]	574	#endif
	575	if(ierr/=NF_NOERR) then
	576	write(,) NF_STRERROR(ierr)
	577	stop "getvarup"
	578	endif
	579	! write(,)'lecture v ok',v
	580
	581	#ifdef NC_DOUBLE
[2191]	582	ierr = NF_GET_VAR_DOUBLE(nid,var3didin(10),ug)
[2118]	583	#else
[2191]	584	ierr = NF_GET_VAR_REAL(nid,var3didin(10),ug)
[2118]	585	#endif
	586	if(ierr/=NF_NOERR) then
	587	write(,) NF_STRERROR(ierr)
	588	stop "getvarup"
	589	endif
[2191]	590	! write(,)'lecture ug ok',ug
[2118]	591
	592	#ifdef NC_DOUBLE
[2191]	593	ierr = NF_GET_VAR_DOUBLE(nid,var3didin(11),vg)
[2118]	594	#else
[2191]	595	ierr = NF_GET_VAR_REAL(nid,var3didin(11),vg)
[2118]	596	#endif
	597	if(ierr/=NF_NOERR) then
	598	write(,) NF_STRERROR(ierr)
	599	stop "getvarup"
	600	endif
[2191]	601	! write(,)'lecture vg ok',vg
	602
	603	#ifdef NC_DOUBLE
	604	ierr = NF_GET_VAR_DOUBLE(nid,var3didin(12),w)
	605	#else
	606	ierr = NF_GET_VAR_REAL(nid,var3didin(12),w)
	607	#endif
	608	if(ierr/=NF_NOERR) then
	609	write(,) NF_STRERROR(ierr)
	610	stop "getvarup"
	611	endif
	612	! write(,)'lecture w ok',w
	613
	614	#ifdef NC_DOUBLE
	615	ierr = NF_GET_VAR_DOUBLE(nid,var3didin(13),du)
	616	#else
	617	ierr = NF_GET_VAR_REAL(nid,var3didin(13),du)
	618	#endif
	619	if(ierr/=NF_NOERR) then
	620	write(,) NF_STRERROR(ierr)
	621	stop "getvarup"
	622	endif
	623	! write(,)'lecture du ok',du
	624
	625	#ifdef NC_DOUBLE
	626	ierr = NF_GET_VAR_DOUBLE(nid,var3didin(14),hu)
	627	#else
	628	ierr = NF_GET_VAR_REAL(nid,var3didin(14),hu)
	629	#endif
	630	if(ierr/=NF_NOERR) then
	631	write(,) NF_STRERROR(ierr)
	632	stop "getvarup"
	633	endif
	634	! write(,)'lecture hu ok',hu
	635
	636	#ifdef NC_DOUBLE
	637	ierr = NF_GET_VAR_DOUBLE(nid,var3didin(15),vu)
	638	#else
	639	ierr = NF_GET_VAR_REAL(nid,var3didin(15),vu)
	640	#endif
	641	if(ierr/=NF_NOERR) then
	642	write(,) NF_STRERROR(ierr)
	643	stop "getvarup"
	644	endif
	645	! write(,)'lecture vu ok',vu
	646
	647	#ifdef NC_DOUBLE
	648	ierr = NF_GET_VAR_DOUBLE(nid,var3didin(16),dv)
	649	#else
	650	ierr = NF_GET_VAR_REAL(nid,var3didin(16),dv)
	651	#endif
	652	if(ierr/=NF_NOERR) then
	653	write(,) NF_STRERROR(ierr)
	654	stop "getvarup"
	655	endif
	656	! write(,)'lecture dv ok',dv
	657
	658	#ifdef NC_DOUBLE
	659	ierr = NF_GET_VAR_DOUBLE(nid,var3didin(17),hv)
	660	#else
	661	ierr = NF_GET_VAR_REAL(nid,var3didin(17),hv)
	662	#endif
	663	if(ierr/=NF_NOERR) then
	664	write(,) NF_STRERROR(ierr)
	665	stop "getvarup"
	666	endif
	667	! write(,)'lecture hv ok',hv
	668
	669	#ifdef NC_DOUBLE
	670	ierr = NF_GET_VAR_DOUBLE(nid,var3didin(18),vv)
	671	#else
	672	ierr = NF_GET_VAR_REAL(nid,var3didin(18),vv)
	673	#endif
	674	if(ierr/=NF_NOERR) then
	675	write(,) NF_STRERROR(ierr)
	676	stop "getvarup"
	677	endif
	678	! write(,)'lecture vv ok',vv
	679
	680	#ifdef NC_DOUBLE
	681	ierr = NF_GET_VAR_DOUBLE(nid,var3didin(19),dt)
	682	#else
	683	ierr = NF_GET_VAR_REAL(nid,var3didin(19),dt)
	684	#endif
	685	if(ierr/=NF_NOERR) then
	686	write(,) NF_STRERROR(ierr)
	687	stop "getvarup"
	688	endif
[2118]	689	! write(,)'lecture dt ok',dt
	690
	691	#ifdef NC_DOUBLE
[2191]	692	ierr = NF_GET_VAR_DOUBLE(nid,var3didin(20),ht)
[2118]	693	#else
[2191]	694	ierr = NF_GET_VAR_REAL(nid,var3didin(20),ht)
[2118]	695	#endif
	696	if(ierr/=NF_NOERR) then
	697	write(,) NF_STRERROR(ierr)
	698	stop "getvarup"
	699	endif
[2191]	700	! write(,)'lecture ht ok',ht
	701
	702	#ifdef NC_DOUBLE
	703	ierr = NF_GET_VAR_DOUBLE(nid,var3didin(21),vt)
	704	#else
	705	ierr = NF_GET_VAR_REAL(nid,var3didin(21),vt)
	706	#endif
	707	if(ierr/=NF_NOERR) then
	708	write(,) NF_STRERROR(ierr)
	709	stop "getvarup"
	710	endif
	711	! write(,)'lecture vt ok',vt
	712
	713	#ifdef NC_DOUBLE
	714	ierr = NF_GET_VAR_DOUBLE(nid,var3didin(22),dq)
	715	#else
	716	ierr = NF_GET_VAR_REAL(nid,var3didin(22),dq)
	717	#endif
	718	if(ierr/=NF_NOERR) then
	719	write(,) NF_STRERROR(ierr)
	720	stop "getvarup"
	721	endif
[2118]	722	! write(,)'lecture dq ok',dq
	723
	724	#ifdef NC_DOUBLE
[2191]	725	ierr = NF_GET_VAR_DOUBLE(nid,var3didin(23),hq)
[2118]	726	#else
[2191]	727	ierr = NF_GET_VAR_REAL(nid,var3didin(23),hq)
[2118]	728	#endif
	729	if(ierr/=NF_NOERR) then
	730	write(,) NF_STRERROR(ierr)
	731	stop "getvarup"
	732	endif
[2191]	733	! write(,)'lecture hq ok',hq
	734
	735	#ifdef NC_DOUBLE
	736	ierr = NF_GET_VAR_DOUBLE(nid,var3didin(24),vq)
	737	#else
	738	ierr = NF_GET_VAR_REAL(nid,var3didin(24),vq)
	739	#endif
	740	if(ierr/=NF_NOERR) then
	741	write(,) NF_STRERROR(ierr)
	742	stop "getvarup"
	743	endif
	744	! write(,)'lecture vq ok',vq
	745
	746	#ifdef NC_DOUBLE
	747	ierr = NF_GET_VAR_DOUBLE(nid,var3didin(25),dth)
	748	#else
	749	ierr = NF_GET_VAR_REAL(nid,var3didin(25),dth)
	750	#endif
	751	if(ierr/=NF_NOERR) then
	752	write(,) NF_STRERROR(ierr)
	753	stop "getvarup"
	754	endif
	755	! write(,)'lecture dth ok',dth
	756
	757	#ifdef NC_DOUBLE
	758	ierr = NF_GET_VAR_DOUBLE(nid,var3didin(26),hth)
	759	#else
	760	ierr = NF_GET_VAR_REAL(nid,var3didin(26),hth)
	761	#endif
	762	if(ierr/=NF_NOERR) then
	763	write(,) NF_STRERROR(ierr)
	764	stop "getvarup"
	765	endif
	766	! write(,)'lecture hth ok',hth
	767
	768	#ifdef NC_DOUBLE
	769	ierr = NF_GET_VAR_DOUBLE(nid,var3didin(27),vth)
	770	#else
	771	ierr = NF_GET_VAR_REAL(nid,var3didin(27),vth)
	772	#endif
	773	if(ierr/=NF_NOERR) then
	774	write(,) NF_STRERROR(ierr)
	775	stop "getvarup"
	776	endif
	777	! write(,)'lecture vth ok',vth
	778
	779	#ifdef NC_DOUBLE
	780	ierr = NF_GET_VAR_DOUBLE(nid,var3didin(28),dr)
	781	#else
	782	ierr = NF_GET_VAR_REAL(nid,var3didin(28),dr)
	783	#endif
	784	if(ierr/=NF_NOERR) then
	785	write(,) NF_STRERROR(ierr)
	786	stop "getvarup"
	787	endif
	788	! write(,)'lecture dr ok',dr
	789
	790	#ifdef NC_DOUBLE
	791	ierr = NF_GET_VAR_DOUBLE(nid,var3didin(29),hr)
	792	#else
	793	ierr = NF_GET_VAR_REAL(nid,var3didin(29),hr)
	794	#endif
	795	if(ierr/=NF_NOERR) then
	796	write(,) NF_STRERROR(ierr)
	797	stop "getvarup"
	798	endif
	799	! write(,)'lecture hr ok',hr
	800
	801	#ifdef NC_DOUBLE
	802	ierr = NF_GET_VAR_DOUBLE(nid,var3didin(30),vr)
	803	#else
	804	ierr = NF_GET_VAR_REAL(nid,var3didin(30),vr)
	805	#endif
	806	if(ierr/=NF_NOERR) then
	807	write(,) NF_STRERROR(ierr)
	808	stop "getvarup"
	809	endif
	810	! write(,)'lecture vr ok',vr
	811
	812	#ifdef NC_DOUBLE
	813	ierr = NF_GET_VAR_DOUBLE(nid,var3didin(31),dtrad)
	814	#else
	815	ierr = NF_GET_VAR_REAL(nid,var3didin(31),dtrad)
	816	#endif
	817	if(ierr/=NF_NOERR) then
	818	write(,) NF_STRERROR(ierr)
	819	stop "getvarup"
	820	endif
	821	! write(,)'lecture dtrad ok',dtrad
	822
	823	#ifdef NC_DOUBLE
	824	ierr = NF_GET_VAR_DOUBLE(nid,var3didin(32),sens)
	825	#else
	826	ierr = NF_GET_VAR_REAL(nid,var3didin(32),sens)
	827	#endif
	828	if(ierr/=NF_NOERR) then
	829	write(,) NF_STRERROR(ierr)
	830	stop "getvarup"
	831	endif
[2118]	832	! write(,)'lecture sens ok',sens
	833
	834	#ifdef NC_DOUBLE
[2191]	835	ierr = NF_GET_VAR_DOUBLE(nid,var3didin(33),flat)
[2118]	836	#else
[2191]	837	ierr = NF_GET_VAR_REAL(nid,var3didin(33),flat)
[2118]	838	#endif
	839	if(ierr/=NF_NOERR) then
	840	write(,) NF_STRERROR(ierr)
	841	stop "getvarup"
	842	endif
	843	! write(,)'lecture flat ok',flat
	844
	845	#ifdef NC_DOUBLE
[2191]	846	ierr = NF_GET_VAR_DOUBLE(nid,var3didin(34),ts)
[2118]	847	#else
[2191]	848	ierr = NF_GET_VAR_REAL(nid,var3didin(34),ts)
[2118]	849	#endif
	850	if(ierr/=NF_NOERR) then
	851	write(,) NF_STRERROR(ierr)
	852	stop "getvarup"
	853	endif
[2191]	854	! write(,)'lecture ts ok',ts
[2118]	855
	856	return
[2191]	857	end subroutine read_cas
[2118]	858	!======================================================================
[2191]	859	SUBROUTINE interp_case_time(day,day1,annee_ref &
	860	& ,year_ini_cas,day_ini_cas,nt_cas,pdt_cas,nlev_cas &
	861	& ,ts_cas,plev_cas,t_cas,q_cas,u_cas,v_cas &
	862	& ,ug_cas,vg_cas,vitw_cas,du_cas,hu_cas,vu_cas &
	863	& ,dv_cas,hv_cas,vv_cas,dt_cas,ht_cas,vt_cas,dtrad_cas &
	864	& ,dq_cas,hq_cas,vq_cas,lat_cas,sens_cas &
	865	& ,ts_prof_cas,plev_prof_cas,t_prof_cas,q_prof_cas &
	866	& ,u_prof_cas,v_prof_cas,ug_prof_cas,vg_prof_cas &
	867	& ,vitw_prof_cas,du_prof_cas,hu_prof_cas,vu_prof_cas &
	868	& ,dv_prof_cas,hv_prof_cas,vv_prof_cas,dt_prof_cas &
	869	& ,ht_prof_cas,vt_prof_cas,dtrad_prof_cas,dq_prof_cas &
	870	& ,hq_prof_cas,vq_prof_cas,lat_prof_cas,sens_prof_cas)
	871
	872
[2118]	873	implicit none
	874
	875	!---------------------------------------------------------------------------------------
	876	! Time interpolation of a 2D field to the timestep corresponding to day
	877	!
	878	! day: current julian day (e.g. 717538.2)
	879	! day1: first day of the simulation
[2191]	880	! nt_cas: total nb of data in the forcing
	881	! pdt_cas: total time interval (in sec) between 2 forcing data
[2118]	882	!---------------------------------------------------------------------------------------
	883
	884	#include "compar1d.h"
	885
	886	! inputs:
	887	integer annee_ref
[2191]	888	integer nt_cas,nlev_cas
	889	integer year_ini_cas
	890	real day_ini_cas
	891	real day, day1,pdt_cas
	892	real ts_cas(nt_cas)
	893	real plev_cas(nlev_cas,nt_cas)
	894	real t_cas(nlev_cas,nt_cas),q_cas(nlev_cas,nt_cas)
	895	real u_cas(nlev_cas,nt_cas),v_cas(nlev_cas,nt_cas)
	896	real ug_cas(nlev_cas,nt_cas),vg_cas(nlev_cas,nt_cas)
	897	real vitw_cas(nlev_cas,nt_cas)
	898	real du_cas(nlev_cas,nt_cas),hu_cas(nlev_cas,nt_cas),vu_cas(nlev_cas,nt_cas)
	899	real dv_cas(nlev_cas,nt_cas),hv_cas(nlev_cas,nt_cas),vv_cas(nlev_cas,nt_cas)
	900	real dt_cas(nlev_cas,nt_cas),ht_cas(nlev_cas,nt_cas),vt_cas(nlev_cas,nt_cas)
	901	real dtrad_cas(nlev_cas,nt_cas)
	902	real dq_cas(nlev_cas,nt_cas),hq_cas(nlev_cas,nt_cas),vq_cas(nlev_cas,nt_cas)
	903	real lat_cas(nt_cas)
	904	real sens_cas(nt_cas)
	905
[2118]	906	! outputs:
[2191]	907	real plev_prof_cas(nlev_cas)
	908	real t_prof_cas(nlev_cas),q_prof_cas(nlev_cas)
	909	real u_prof_cas(nlev_cas),v_prof_cas(nlev_cas)
	910	real ug_prof_cas(nlev_cas),vg_prof_cas(nlev_cas)
	911	real vitw_prof_cas(nlev_cas)
	912	real du_prof_cas(nlev_cas),hu_prof_cas(nlev_cas),vu_prof_cas(nlev_cas)
	913	real dv_prof_cas(nlev_cas),hv_prof_cas(nlev_cas),vv_prof_cas(nlev_cas)
	914	real dt_prof_cas(nlev_cas),ht_prof_cas(nlev_cas),vt_prof_cas(nlev_cas)
	915	real dtrad_prof_cas(nlev_cas)
	916	real dq_prof_cas(nlev_cas),hq_prof_cas(nlev_cas),vq_prof_cas(nlev_cas)
	917	real lat_prof_cas,sens_prof_cas,ts_prof_cas
[2118]	918	! local:
[2191]	919	integer it_cas1, it_cas2,k
	920	real timeit,time_cas1,time_cas2,frac
[2118]	921
	922
[2191]	923	print*,'Check time',day1,day_ini_cas,day_ini_cas+90
	924
	925	if ((forcing_type.eq.10).and.(1.eq.1)) then
	926	! Check that initial day of the simulation consistent with the case:
	927	if (annee_ref.ne.2011) then
	928	print*,'Pour CINDY, annee_ref doit etre 2011'
	929	print*,'Changer annee_ref dans run.def'
	930	stop
	931	endif
	932	if (annee_ref.eq.2011 .and. day1.lt.day_ini_cas) then
	933	print*,'CINDY a débuté le 1 octobre 2011',day1,day_ini_cas
	934	print*,'Changer dayref dans run.def'
	935	stop
	936	endif
	937	if (annee_ref.eq.2011 .and. day1.gt.day_ini_cas+90) then
	938	print*,'CINDY a fini le 31 decembre'
	939	print*,'Changer dayref ou nday dans run.def',day1,day_ini_cas+90
	940	stop
	941	endif
	942	endif
	943	!CR:ATTENTION TEST AMMA
	944	if ((forcing_type.eq.10).and.(1.eq.0)) then
[2118]	945	! Check that initial day of the simulation consistent with AMMA case:
	946	if (annee_ref.ne.2006) then
	947	print*,'Pour AMMA, annee_ref doit etre 2006'
	948	print*,'Changer annee_ref dans run.def'
	949	stop
	950	endif
[2191]	951	if (annee_ref.eq.2006 .and. day1.lt.day_ini_cas) then
	952	print*,'AMMA a débuté le 10 juillet 2006',day1,day_ini_cas
[2118]	953	print*,'Changer dayref dans run.def'
	954	stop
	955	endif
[2191]	956	if (annee_ref.eq.2006 .and. day1.gt.day_ini_cas+1) then
[2118]	957	print*,'AMMA a fini le 11 juillet'
	958	print*,'Changer dayref ou nday dans run.def'
	959	stop
	960	endif
	961	endif
	962
[2191]	963	! Determine timestep relative to the 1st day:
[2118]	964	! timeit=(day-day1)*86400.
	965	! if (annee_ref.eq.1992) then
[2191]	966	! timeit=(day-day_ini_cas)*86400.
[2118]	967	! else
	968	! timeit=(day+61.-1.)*86400. ! 61 days between Nov01 and Dec31 1992
	969	! endif
[2191]	970	timeit=(day-day_ini_cas)*86400
[2118]	971
	972	! Determine the closest observation times:
[2191]	973	! it_cas1=INT(timeit/pdt_cas)+1
	974	! it_cas2=it_cas1 + 1
	975	! time_cas1=(it_cas1-1)*pdt_cas
	976	! time_cas2=(it_cas2-1)*pdt_cas
[2118]	977
[2191]	978	it_cas1=INT(timeit/pdt_cas)+1
	979	IF (it_cas1 .EQ. nt_cas) THEN
	980	it_cas2=it_cas1
[2118]	981	ELSE
[2191]	982	it_cas2=it_cas1 + 1
[2118]	983	ENDIF
[2191]	984	time_cas1=(it_cas1-1)*pdt_cas
	985	time_cas2=(it_cas2-1)*pdt_cas
[2118]	986
[2191]	987	if (it_cas1 .gt. nt_cas) then
	988	write(,) 'PB-stop: day, it_cas1, it_cas2, timeit: ' &
	989	& ,day,day_ini_cas,it_cas1,it_cas2,timeit/86400.
[2118]	990	stop
	991	endif
	992
	993	! time interpolation:
[2191]	994	frac=(time_cas2-timeit)/(time_cas2-time_cas1)
[2118]	995	frac=max(frac,0.0)
	996
[2191]	997	lat_prof_cas = lat_cas(it_cas2) &
	998	& -frac*(lat_cas(it_cas2)-lat_cas(it_cas1))
	999	sens_prof_cas = sens_cas(it_cas2) &
	1000	& -frac*(sens_cas(it_cas2)-sens_cas(it_cas1))
	1001	ts_prof_cas = ts_cas(it_cas2) &
	1002	& -frac*(ts_cas(it_cas2)-ts_cas(it_cas1))
[2118]	1003
[2191]	1004	do k=1,nlev_cas
	1005	plev_prof_cas(k) = plev_cas(k,it_cas2) &
	1006	& -frac*(plev_cas(k,it_cas2)-plev_cas(k,it_cas1))
	1007	t_prof_cas(k) = t_cas(k,it_cas2) &
	1008	& -frac*(t_cas(k,it_cas2)-t_cas(k,it_cas1))
	1009	q_prof_cas(k) = q_cas(k,it_cas2) &
	1010	& -frac*(q_cas(k,it_cas2)-q_cas(k,it_cas1))
	1011	u_prof_cas(k) = u_cas(k,it_cas2) &
	1012	& -frac*(u_cas(k,it_cas2)-u_cas(k,it_cas1))
	1013	v_prof_cas(k) = v_cas(k,it_cas2) &
	1014	& -frac*(v_cas(k,it_cas2)-v_cas(k,it_cas1))
	1015	ug_prof_cas(k) = ug_cas(k,it_cas2) &
	1016	& -frac*(ug_cas(k,it_cas2)-ug_cas(k,it_cas1))
	1017	vg_prof_cas(k) = vg_cas(k,it_cas2) &
	1018	& -frac*(vg_cas(k,it_cas2)-vg_cas(k,it_cas1))
	1019	vitw_prof_cas(k) = vitw_cas(k,it_cas2) &
	1020	& -frac*(vitw_cas(k,it_cas2)-vitw_cas(k,it_cas1))
	1021	du_prof_cas(k) = du_cas(k,it_cas2) &
	1022	& -frac*(du_cas(k,it_cas2)-du_cas(k,it_cas1))
	1023	hu_prof_cas(k) = hu_cas(k,it_cas2) &
	1024	& -frac*(hu_cas(k,it_cas2)-hu_cas(k,it_cas1))
	1025	vu_prof_cas(k) = vu_cas(k,it_cas2) &
	1026	& -frac*(vu_cas(k,it_cas2)-vu_cas(k,it_cas1))
	1027	dv_prof_cas(k) = dv_cas(k,it_cas2) &
	1028	& -frac*(dv_cas(k,it_cas2)-dv_cas(k,it_cas1))
	1029	hv_prof_cas(k) = hv_cas(k,it_cas2) &
	1030	& -frac*(hv_cas(k,it_cas2)-hv_cas(k,it_cas1))
	1031	vv_prof_cas(k) = vv_cas(k,it_cas2) &
	1032	& -frac*(vv_cas(k,it_cas2)-vv_cas(k,it_cas1))
	1033	dt_prof_cas(k) = dt_cas(k,it_cas2) &
	1034	& -frac*(dt_cas(k,it_cas2)-dt_cas(k,it_cas1))
	1035	ht_prof_cas(k) = ht_cas(k,it_cas2) &
	1036	& -frac*(ht_cas(k,it_cas2)-ht_cas(k,it_cas1))
	1037	vt_prof_cas(k) = vt_cas(k,it_cas2) &
	1038	& -frac*(vt_cas(k,it_cas2)-vt_cas(k,it_cas1))
	1039	dtrad_prof_cas(k) = dtrad_cas(k,it_cas2) &
	1040	& -frac*(dtrad_cas(k,it_cas2)-dtrad_cas(k,it_cas1))
	1041	dq_prof_cas(k) = dq_cas(k,it_cas2) &
	1042	& -frac*(dq_cas(k,it_cas2)-dq_cas(k,it_cas1))
	1043	hq_prof_cas(k) = hq_cas(k,it_cas2) &
	1044	& -frac*(hq_cas(k,it_cas2)-hq_cas(k,it_cas1))
	1045	vq_prof_cas(k) = vq_cas(k,it_cas2) &
	1046	& -frac*(vq_cas(k,it_cas2)-vq_cas(k,it_cas1))
[2118]	1047	enddo
	1048
	1049	return
	1050	END
	1051
[2191]	1052	!**********************************************************************************************

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

Download in other formats:

Original Format