Context Navigation

source: LMDZ6/trunk/libf/phylmd/dyn1d/1D_interp_cases.h @ 4142

Last change on this file since 4142 was 4104, checked in by evignon, 3 years ago
modification des routines pour lecture de la v1 des cas 1D au format standard + corrections de bug pour le geopotentiel à la surface
Property svn:keywords set to `Id`
File size: 8.5 KB

Rev	Line
[2017]	1
[3780]	2	print*,'FORCING CASE forcing_case2'
[3594]	3	! print*, &
	4	! & '#### ITAP,day,day1,(day-day1)86400,(day-day1)86400/pdt_cas=', &
	5	! & daytime,day1,(daytime-day1)*86400., &
	6	! & (daytime-day1)*86400/pdt_cas
[2191]	7
	8	! time interpolation:
[3686]	9	CALL interp_case_time_std(daytime,day1,annee_ref &
	10	! & ,year_ini_cas,day_ju_ini_cas,nt_cas,pdt_cas,nlev_cas &
	11	& ,nt_cas,nlev_cas &
	12	& ,ts_cas,ps_cas,plev_cas,t_cas,th_cas,thv_cas,thl_cas,qv_cas,ql_cas,qi_cas &
	13	& ,u_cas,v_cas,ug_cas,vg_cas &
	14	& ,temp_nudg_cas,qv_nudg_cas,u_nudg_cas,v_nudg_cas &
[4104]	15	& ,invtau_temp_nudg_cas,invtau_qv_nudg_cas,invtau_u_nudg_cas,invtau_v_nudg_cas &
[3781]	16	& ,vitw_cas,omega_cas,tke_cas,du_cas,hu_cas,vu_cas &
[3686]	17	& ,dv_cas,hv_cas,vv_cas,dt_cas,ht_cas,vt_cas,dtrad_cas &
	18	& ,dq_cas,hq_cas,vq_cas,dth_cas,hth_cas,vth_cas,lat_cas,sens_cas,ustar_cas &
[3781]	19	& ,uw_cas,vw_cas,q1_cas,q2_cas,tkes_cas &
[2716]	20	!
[3686]	21	& ,ts_prof_cas,ps_prof_cas,plev_prof_cas,t_prof_cas,theta_prof_cas,thv_prof_cas &
	22	& ,thl_prof_cas,qv_prof_cas,ql_prof_cas,qi_prof_cas &
	23	& ,u_prof_cas,v_prof_cas,ug_prof_cas,vg_prof_cas &
	24	& ,temp_nudg_prof_cas,qv_nudg_prof_cas,u_nudg_prof_cas,v_nudg_prof_cas &
[4104]	25	& ,invtau_temp_nudg_prof_cas,invtau_qv_nudg_prof_cas,invtau_u_nudg_prof_cas,invtau_v_nudg_prof_cas &
[3781]	26	& ,vitw_prof_cas,omega_prof_cas,tke_prof_cas &
[3686]	27	& ,du_prof_cas,hu_prof_cas,vu_prof_cas &
	28	& ,dv_prof_cas,hv_prof_cas,vv_prof_cas,dt_prof_cas,ht_prof_cas,vt_prof_cas &
	29	& ,dtrad_prof_cas,dq_prof_cas,hq_prof_cas,vq_prof_cas &
	30	& ,dth_prof_cas,hth_prof_cas,vth_prof_cas,lat_prof_cas &
[3781]	31	& ,sens_prof_cas,ustar_prof_cas,uw_prof_cas,vw_prof_cas,q1_prof_cas,q2_prof_cas,tkes_prof_cas)
[3780]	32	! EV tg instead of ts_cur
	33	tg = ts_prof_cas
[2716]	34	! psurf=plev_prof_cas(1)
	35	psurf=ps_prof_cas
	36
	37	! vertical interpolation:
[3781]	38	CALL interp2_case_vertical_std(play,plev,nlev_cas,plev_prof_cas &
[3686]	39	& ,t_prof_cas,theta_prof_cas,thv_prof_cas,thl_prof_cas &
	40	& ,qv_prof_cas,ql_prof_cas,qi_prof_cas,u_prof_cas,v_prof_cas &
	41	& ,ug_prof_cas,vg_prof_cas &
	42	& ,temp_nudg_prof_cas,qv_nudg_prof_cas,u_nudg_prof_cas,v_nudg_prof_cas &
[4104]	43	& ,invtau_temp_nudg_prof_cas,invtau_qv_nudg_prof_cas,invtau_u_nudg_prof_cas,invtau_v_nudg_prof_cas &
	44
[3781]	45	& ,vitw_prof_cas,omega_prof_cas,tke_prof_cas &
[3686]	46	& ,du_prof_cas,hu_prof_cas,vu_prof_cas,dv_prof_cas,hv_prof_cas,vv_prof_cas &
	47	& ,dt_prof_cas,ht_prof_cas,vt_prof_cas,dtrad_prof_cas,dq_prof_cas,hq_prof_cas,vq_prof_cas &
	48	& ,dth_prof_cas,hth_prof_cas,vth_prof_cas &
[2716]	49	!
[3686]	50	& ,t_mod_cas,theta_mod_cas,thv_mod_cas,thl_mod_cas,qv_mod_cas,ql_mod_cas,qi_mod_cas &
	51	& ,u_mod_cas,v_mod_cas,ug_mod_cas,vg_mod_cas &
	52	& ,temp_nudg_mod_cas,qv_nudg_mod_cas,u_nudg_mod_cas,v_nudg_mod_cas &
[4104]	53	& ,invtau_temp_nudg_mod_cas,invtau_qv_nudg_mod_cas,invtau_u_nudg_mod_cas,invtau_v_nudg_mod_cas &
[3781]	54	& ,w_mod_cas,omega_mod_cas,tke_mod_cas &
[3686]	55	& ,du_mod_cas,hu_mod_cas,vu_mod_cas,dv_mod_cas,hv_mod_cas,vv_mod_cas &
	56	& ,dt_mod_cas,ht_mod_cas,vt_mod_cas,dtrad_mod_cas,dq_mod_cas,hq_mod_cas,vq_mod_cas &
	57	& ,dth_mod_cas,hth_mod_cas,vth_mod_cas,mxcalc)
[2716]	58
	59
	60	DO l=1,llm
	61	teta(l)=temp(l)(100000./play(l))*(rd/rcpd)
	62	ENDDO
	63	!calcul de l'advection verticale a partir du omega
	64	!Calcul des gradients verticaux
	65	!initialisation
	66	d_t_z(:)=0.
	67	d_th_z(:)=0.
	68	d_q_z(:)=0.
[2920]	69	d_u_z(:)=0.
	70	d_v_z(:)=0.
[2716]	71	d_t_dyn_z(:)=0.
	72	d_th_dyn_z(:)=0.
	73	d_q_dyn_z(:)=0.
[2920]	74	d_u_dyn_z(:)=0.
	75	d_v_dyn_z(:)=0.
[3541]	76	if (1==0) then
	77	DO l=2,llm-1
	78	d_t_z(l)=(temp(l+1)-temp(l-1))/(play(l+1)-play(l-1))
	79	d_th_z(l)=(teta(l+1)-teta(l-1))/(play(l+1)-play(l-1))
	80	d_q_z(l)=(q(l+1,1)-q(l-1,1))/(play(l+1)-play(l-1))
	81	d_u_z(l)=(u(l+1)-u(l-1))/(play(l+1)-play(l-1))
	82	d_v_z(l)=(v(l+1)-v(l-1))/(play(l+1)-play(l-1))
	83	ENDDO
	84	else
	85	DO l=2,llm-1
[3686]	86	IF (omega(l)>0.) THEN
[3541]	87	d_t_z(l)=(temp(l+1)-temp(l))/(play(l+1)-play(l))
	88	d_th_z(l)=(teta(l+1)-teta(l))/(play(l+1)-play(l))
	89	d_q_z(l)=(q(l+1,1)-q(l,1))/(play(l+1)-play(l))
	90	d_u_z(l)=(u(l+1)-u(l))/(play(l+1)-play(l))
	91	d_v_z(l)=(v(l+1)-v(l))/(play(l+1)-play(l))
[3686]	92	ELSE
[3541]	93	d_t_z(l)=(temp(l-1)-temp(l))/(play(l-1)-play(l))
	94	d_th_z(l)=(teta(l-1)-teta(l))/(play(l-1)-play(l))
	95	d_q_z(l)=(q(l-1,1)-q(l,1))/(play(l-1)-play(l))
	96	d_u_z(l)=(u(l-1)-u(l))/(play(l-1)-play(l))
	97	d_v_z(l)=(v(l-1)-v(l))/(play(l-1)-play(l))
[3686]	98	ENDIF
	99	ENDDO
[3541]	100	endif
[2716]	101	d_t_z(1)=d_t_z(2)
[3541]	102	d_t_z(1)=d_t_z(2)
[2716]	103	d_th_z(1)=d_th_z(2)
	104	d_q_z(1)=d_q_z(2)
[2920]	105	d_u_z(1)=d_u_z(2)
	106	d_v_z(1)=d_v_z(2)
[2716]	107	d_t_z(llm)=d_t_z(llm-1)
	108	d_th_z(llm)=d_th_z(llm-1)
	109	d_q_z(llm)=d_q_z(llm-1)
[2920]	110	d_u_z(llm)=d_u_z(llm-1)
	111	d_v_z(llm)=d_v_z(llm-1)
[2716]	112
[3541]	113	! TRAVAIL : PRENDRE DES NOTATIONS COHERENTES POUR W
	114	do l = 1, llm
	115	! Modif w_mod_cas -> omega_mod_cas (MM+MPL 20170309)
[3781]	116	print*, l, llm
	117	print*, play(l), temp(l)
[3541]	118	omega(l) = -w_mod_cas(l)play(l)rg/(rd*temp(l))
	119	enddo
	120
[2716]	121	!Calcul de l advection verticale
[2920]	122	! Modif w_mod_cas -> omega_mod_cas (MM+MPL 20170310)
[3541]	123	d_t_dyn_z(:)=omega(:)*d_t_z(:)
	124	d_th_dyn_z(:)=omega(:)*d_th_z(:)
	125	d_q_dyn_z(:)=omega(:)*d_q_z(:)
	126	d_u_dyn_z(:)=omega(:)*d_u_z(:)
	127	d_v_dyn_z(:)=omega(:)*d_v_z(:)
[2716]	128
[2920]	129	!geostrophic wind
	130	if (forc_geo.eq.1) then
	131	do l=1,llm
	132	ug(l) = ug_mod_cas(l)
	133	vg(l) = vg_mod_cas(l)
	134	enddo
	135	endif
[2716]	136
	137	do l = 1, llm
[3592]	138
	139	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
[2920]	140	! Modif w_mod_cas -> omega_mod_cas (MM+MPL 20170309)
[3592]	141	!!! omega2(l)= omega(l)/rg*airefi ! flxmass_w calcule comme ds physiq
	142	omega(l) = omega_mod_cas(l)
	143	omega2(l)= omega_mod_cas(l)/rg*airefi ! flxmass_w calcule comme ds physiq
	144
[3593]	145	! On effectue la somme du forcage total et de la decomposition
	146	! horizontal/vertical en supposant que soit l'un soit l'autre
	147	! sont remplis mais jamais les deux
[2716]	148
[3593]	149	d_t_adv(l) = dt_mod_cas(l)+ht_mod_cas(l)+vt_mod_cas(l)
	150	d_q_adv(l,1) = dq_mod_cas(l)+hq_mod_cas(l)+vq_mod_cas(l)
	151	d_q_adv(l,2) = 0.0
	152	d_u_adv(l) = du_mod_cas(l)+hu_mod_cas(l)+vu_mod_cas(l)
	153	d_v_adv(l) = dv_mod_cas(l)+hv_mod_cas(l)+vv_mod_cas(l)
[2716]	154
[3592]	155	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
	156	!! CONSERVE EN ATTENDANT QUE LE CAS EN QUESTION FONCTIONNE EN STD !!
	157	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
	158	!if (forc_w==1) then
	159	! d_q_adv(l,1)=d_q_adv(l,1)-d_q_dyn_z(l)
	160	! d_t_adv(l)=d_t_adv(l)-d_t_dyn_z(l)
	161	! d_v_adv(l)=d_v_adv(l)-d_v_dyn_z(l)
	162	! d_u_adv(l)=d_u_adv(l)-d_u_dyn_z(l)
	163	! endif
	164	!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!
[2716]	165
[3686]	166	if (trad.eq.1) then
[2716]	167	tend_rayo=1
	168	dt_cooling(l) = dtrad_mod_cas(l)
	169	! print *,'dt_cooling=',dt_cooling(l)
	170	else
	171	dt_cooling(l) = 0.0
	172	endif
	173	enddo
	174
	175	! Faut-il multiplier par -1 ? (MPL 20160713)
	176	IF(ok_flux_surf) THEN
	177	fsens=-1.*sens_prof_cas
	178	flat=-1.*lat_prof_cas
	179	print *,'1D_interp: sens,flat',fsens,flat
	180	ENDIF
	181	!
	182	IF (ok_prescr_ust) THEN
	183	ust=ustar_prof_cas
	184	print *,'ust=',ust
	185	ENDIF

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

Download in other formats: