1 | |
---|
2 | print*,'FORCING CASE forcing_case2' |
---|
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 |
---|
7 | |
---|
8 | ! time interpolation: |
---|
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 & |
---|
15 | & ,invtau_temp_nudg_cas,invtau_qv_nudg_cas,invtau_u_nudg_cas,invtau_v_nudg_cas & |
---|
16 | & ,vitw_cas,omega_cas,tke_cas,du_cas,hu_cas,vu_cas & |
---|
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 & |
---|
19 | & ,uw_cas,vw_cas,q1_cas,q2_cas,tkes_cas & |
---|
20 | ! |
---|
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 & |
---|
25 | & ,invtau_temp_nudg_prof_cas,invtau_qv_nudg_prof_cas,invtau_u_nudg_prof_cas,invtau_v_nudg_prof_cas & |
---|
26 | & ,vitw_prof_cas,omega_prof_cas,tke_prof_cas & |
---|
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 & |
---|
31 | & ,sens_prof_cas,ustar_prof_cas,uw_prof_cas,vw_prof_cas,q1_prof_cas,q2_prof_cas,tkes_prof_cas) |
---|
32 | ! EV tg instead of ts_cur |
---|
33 | tg = ts_prof_cas |
---|
34 | ! psurf=plev_prof_cas(1) |
---|
35 | psurf=ps_prof_cas |
---|
36 | |
---|
37 | ! vertical interpolation: |
---|
38 | CALL interp2_case_vertical_std(play,plev,nlev_cas,plev_prof_cas & |
---|
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 & |
---|
43 | & ,invtau_temp_nudg_prof_cas,invtau_qv_nudg_prof_cas,invtau_u_nudg_prof_cas,invtau_v_nudg_prof_cas & |
---|
44 | |
---|
45 | & ,vitw_prof_cas,omega_prof_cas,tke_prof_cas & |
---|
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 & |
---|
49 | ! |
---|
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 & |
---|
53 | & ,invtau_temp_nudg_mod_cas,invtau_qv_nudg_mod_cas,invtau_u_nudg_mod_cas,invtau_v_nudg_mod_cas & |
---|
54 | & ,w_mod_cas,omega_mod_cas,tke_mod_cas & |
---|
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) |
---|
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. |
---|
69 | d_u_z(:)=0. |
---|
70 | d_v_z(:)=0. |
---|
71 | d_t_dyn_z(:)=0. |
---|
72 | d_th_dyn_z(:)=0. |
---|
73 | d_q_dyn_z(:)=0. |
---|
74 | d_u_dyn_z(:)=0. |
---|
75 | d_v_dyn_z(:)=0. |
---|
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 |
---|
86 | IF (omega(l)>0.) THEN |
---|
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)) |
---|
92 | ELSE |
---|
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)) |
---|
98 | ENDIF |
---|
99 | ENDDO |
---|
100 | endif |
---|
101 | d_t_z(1)=d_t_z(2) |
---|
102 | d_t_z(1)=d_t_z(2) |
---|
103 | d_th_z(1)=d_th_z(2) |
---|
104 | d_q_z(1)=d_q_z(2) |
---|
105 | d_u_z(1)=d_u_z(2) |
---|
106 | d_v_z(1)=d_v_z(2) |
---|
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) |
---|
110 | d_u_z(llm)=d_u_z(llm-1) |
---|
111 | d_v_z(llm)=d_v_z(llm-1) |
---|
112 | |
---|
113 | ! TRAVAIL : PRENDRE DES NOTATIONS COHERENTES POUR W |
---|
114 | do l = 1, llm |
---|
115 | ! Modif w_mod_cas -> omega_mod_cas (MM+MPL 20170309) |
---|
116 | print*, l, llm |
---|
117 | print*, play(l), temp(l) |
---|
118 | omega(l) = -w_mod_cas(l)*play(l)*rg/(rd*temp(l)) |
---|
119 | enddo |
---|
120 | |
---|
121 | !Calcul de l advection verticale |
---|
122 | ! Modif w_mod_cas -> omega_mod_cas (MM+MPL 20170310) |
---|
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(:) |
---|
128 | |
---|
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 |
---|
136 | |
---|
137 | do l = 1, llm |
---|
138 | |
---|
139 | !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! |
---|
140 | ! Modif w_mod_cas -> omega_mod_cas (MM+MPL 20170309) |
---|
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 | |
---|
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 |
---|
148 | |
---|
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) |
---|
154 | |
---|
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 | !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! |
---|
165 | |
---|
166 | if (trad.eq.1) then |
---|
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 |
---|