Context Navigation

calltherm_interface.F90 @ 689

Last change on this file since 689 was 660, checked in by acolaitis, 13 years ago

Thermals: updated computation of wstar to take into account vertical turbulent heat flux from diffusion. Without that, wstar is underestimated. We use the sensible heat flux and tendancies computed in vdifc to compute the vertical turbulent heat flux contribution from the diffusion, and then compute the maximum in the PBL. The resulting total heat flux profile and wstar correspond very well to LES results.

File size: 19.6 KB

Line
1	!
2	! AC 2011-01-05
3	!
4	SUBROUTINE calltherm_interface (firstcall, &
5	& zzlev,zzlay, &
6	& ptimestep,pu,pv,pt,pq,pdu,pdv,pdt,pdq,q2, &
7	& pplay,pplev,pphi,zpopsk, &
8	& pdu_th,pdv_th,pdt_th,pdq_th,lmax,zmaxth,pbl_dtke, &
9	& pdhdif,hfmax,wstar,sensibFlux)
10
11	USE ioipsl_getincom
12
13	implicit none
14	#include "callkeys.h"
15	#include "dimensions.h"
16	#include "dimphys.h"
17	#include "comcstfi.h"
18	#include "tracer.h"
19
20	!--------------------------------------------------------
21	! Input Variables
22	!--------------------------------------------------------
23
24	! REAL, INTENT(IN) :: long(ngridmx),lati(ngridmx)
25	REAL, INTENT(IN) :: ptimestep
26	REAL, INTENT(IN) :: pplev(ngridmx,nlayermx+1),pplay(ngridmx,nlayermx)
27	REAL, INTENT(IN) :: pphi(ngridmx,nlayermx)
28	REAL, INTENT(IN) :: pu(ngridmx,nlayermx),pv(ngridmx,nlayermx)
29	REAL, INTENT(IN) :: pt(ngridmx,nlayermx),pq(ngridmx,nlayermx,nqmx)
30	REAL, INTENT(IN) :: zzlay(ngridmx,nlayermx)
31	REAL, INTENT(IN) :: zzlev(ngridmx,nlayermx+1)
32	LOGICAL, INTENT(IN) :: firstcall
33	REAL, INTENT(IN) :: pdu(ngridmx,nlayermx),pdv(ngridmx,nlayermx)
34	REAL, INTENT(IN) :: pdq(ngridmx,nlayermx,nqmx),pdt(ngridmx,nlayermx)
35	REAL, INTENT(IN) :: q2(ngridmx,nlayermx+1)
36	REAL, INTENT(IN) :: zpopsk(ngridmx,nlayermx)
37	REAL, INTENT(IN) :: pdhdif(ngridmx,nlayermx)
38	REAL, INTENT(IN) :: sensibFlux(ngridmx)
39
40	!--------------------------------------------------------
41	! Output Variables
42	!--------------------------------------------------------
43
44	REAL, INTENT(OUT) :: pdu_th(ngridmx,nlayermx)
45	REAL, INTENT(OUT) :: pdv_th(ngridmx,nlayermx)
46	REAL, INTENT(OUT) :: pdt_th(ngridmx,nlayermx)
47	REAL, INTENT(OUT) :: pdq_th(ngridmx,nlayermx,nqmx)
48	INTEGER, INTENT(OUT) :: lmax(ngridmx)
49	REAL, INTENT(OUT) :: zmaxth(ngridmx)
50	REAL, INTENT(OUT) :: pbl_dtke(ngridmx,nlayermx+1)
51	REAL, INTENT(OUT) :: wstar(ngridmx)
52
53	!--------------------------------------------------------
54	! Thermals local variables
55	!--------------------------------------------------------
56	REAL zu(ngridmx,nlayermx), zv(ngridmx,nlayermx)
57	REAL zt(ngridmx,nlayermx)
58	REAL d_t_ajs(ngridmx,nlayermx)
59	REAL d_u_ajs(ngridmx,nlayermx), d_q_ajs(ngridmx,nlayermx,nqmx)
60	REAL d_v_ajs(ngridmx,nlayermx)
61	REAL fm_therm(ngridmx,nlayermx+1), entr_therm(ngridmx,nlayermx)
62	REAL detr_therm(ngridmx,nlayermx),detrmod(ngridmx,nlayermx)
63	REAL zw2(ngridmx,nlayermx+1)
64	REAL fraca(ngridmx,nlayermx+1),zfraca(ngridmx,nlayermx+1)
65	REAL ztla(ngridmx,nlayermx)
66	REAL q_therm(ngridmx,nlayermx), pq_therm(ngridmx,nlayermx,nqmx)
67	REAL q2_therm(ngridmx,nlayermx), dq2_therm(ngridmx,nlayermx)
68	REAL lmax_real(ngridmx)
69	REAL masse(ngridmx,nlayermx)
70	LOGICAL qtransport_thermals,dtke_thermals
71	INTEGER l,ig,iq,ii(1),k
72	CHARACTER (LEN=20) modname
73
74	!--------------------------------------------------------
75	! Local variables for sub-timestep
76	!--------------------------------------------------------
77
78	REAL d_t_the(ngridmx,nlayermx), d_q_the(ngridmx,nlayermx,nqmx)
79	REAL d_u_the(ngridmx,nlayermx),d_v_the(ngridmx,nlayermx)
80	REAL dq2_the(ngridmx,nlayermx)
81	INTEGER isplit
82	INTEGER,SAVE :: nsplit_thermals
83	REAL, SAVE :: r_aspect_thermals
84	REAL fact
85	REAL zfm_therm(ngridmx,nlayermx+1),zdt
86	REAL zentr_therm(ngridmx,nlayermx),zdetr_therm(ngridmx,nlayermx)
87	REAL zheatFlux(ngridmx,nlayermx)
88	REAL zheatFlux_down(ngridmx,nlayermx)
89	REAL zbuoyancyOut(ngridmx,nlayermx)
90	REAL zbuoyancyEst(ngridmx,nlayermx)
91	REAL zzw2(ngridmx,nlayermx+1)
92	REAL zmax(ngridmx)
93	INTEGER ndt,zlmax
94
95	!--------------------------------------------------------
96	! Diagnostics
97	!--------------------------------------------------------
98
99	REAL heatFlux(ngridmx,nlayermx)
100	REAL heatFlux_down(ngridmx,nlayermx)
101	REAL buoyancyOut(ngridmx,nlayermx)
102	REAL buoyancyEst(ngridmx,nlayermx)
103	REAL hfmax(ngridmx),wmax(ngridmx)
104	REAL pbl_teta(ngridmx),dteta(ngridmx,nlayermx)
105	REAL rpdhd(ngridmx,nlayermx)
106	REAL wtdif(ngridmx,nlayermx),rho(ngridmx,nlayermx)
107	REAL wtth(ngridmx,nlayermx)
108
109	!--------------------------------------------------------
110	! Theta_m
111	!--------------------------------------------------------
112
113	INTEGER ico2
114	SAVE ico2
115
116	! **********************************************************************
117	! Initialization
118	! **********************************************************************
119
120	lmax(:)=0
121	pdu_th(:,:)=0.
122	pdv_th(:,:)=0.
123	pdt_th(:,:)=0.
124	entr_therm(:,:)=0.
125	detr_therm(:,:)=0.
126	q2_therm(:,:)=0.
127	dq2_therm(:,:)=0.
128	ztla(:,:)=0.
129	pbl_dtke(:,:)=0.
130	fm_therm(:,:)=0.
131	zw2(:,:)=0.
132	fraca(:,:)=0.
133	zfraca(:,:)=0.
134	if (tracer) then
135	pdq_th(:,:,:)=0.
136	end if
137	d_t_ajs(:,:)=0.
138	d_u_ajs(:,:)=0.
139	d_v_ajs(:,:)=0.
140	d_q_ajs(:,:,:)=0.
141	heatFlux(:,:)=0.
142	heatFlux_down(:,:)=0.
143	buoyancyOut(:,:)=0.
144	buoyancyEst(:,:)=0.
145	zmaxth(:)=0.
146	lmax_real(:)=0.
147
148
149	! **********************************************************************
150	! Preparing inputs for the thermals
151	! **********************************************************************
152
153	zu(:,:)=pu(:,:)+pdu(:,:)*ptimestep
154	zv(:,:)=pv(:,:)+pdv(:,:)*ptimestep
155	zt(:,:)=pt(:,:)+pdt(:,:)*ptimestep
156
157	pq_therm(:,:,:)=0.
158	qtransport_thermals=.true. !! default setting
159	!call getin("qtransport_thermals",qtransport_thermals)
160
161	if(qtransport_thermals) then
162	if(tracer) then
163	pq_therm(:,:,:)=pq(:,:,:)+pdq(:,:,:)*ptimestep
164	endif
165	endif
166
167	dtke_thermals=.false. !! default setting
168	call getin("dtke_thermals",dtke_thermals)
169	IF(dtke_thermals) THEN
170	DO l=1,nlayermx
171	q2_therm(:,l)=0.5*(q2(:,l)+q2(:,l+1))
172	ENDDO
173	ENDIF
174
175	! **********************************************************************
176	! Polar night mixing : theta_m
177	! **********************************************************************
178
179	if(firstcall) then
180	ico2=0
181	if (tracer) then
182	! Prepare Special treatment if one of the tracers is CO2 gas
183	do iq=1,nqmx
184	if (noms(iq).eq."co2") then
185	ico2=iq
186	end if
187	enddo
188	endif
189	endif !of if firstcall
190
191
192	! **********************************************************************
193	! **********************************************************************
194	! **********************************************************************
195	! CALLTHERM
196	! **********************************************************************
197	! **********************************************************************
198	! **********************************************************************
199
200	! r_aspect_thermals ! Mainly control the shape of the temperature profile
201	! in the surface layer. Decreasing it goes toward
202	! a convective-adjustment like profile.
203	! nsplit_thermals ! Sub-timestep for the thermals. Very dependant on the
204	! chosen timestep for the radiative transfer.
205	! It is recommended to run with 96 timestep per day and
206	! iradia = 1., configuration in which thermals can run
207	! very well with a sub-timestep of 10.
208	IF (firstcall) THEN
209	#ifdef MESOSCALE
210	!! valid for timesteps < 200s
211	nsplit_thermals=4
212	r_aspect_thermals=0.7
213	#else
214	IF ((ptimestep .le. 3699.*24./96.) .and. (iradia .eq. 1)) THEN
215	nsplit_thermals=10
216	ELSE
217	nsplit_thermals=35
218	ENDIF
219	r_aspect_thermals=1.
220	#endif
221	call getin("nsplit_thermals",nsplit_thermals)
222	call getin("r_aspect_thermals",r_aspect_thermals)
223	ENDIF
224
225	! **********************************************************************
226	! SUB-TIMESTEP LOOP
227	! **********************************************************************
228
229	zdt=ptimestep/REAL(nsplit_thermals)
230
231	DO isplit=1,nsplit_thermals
232
233	! Initialization of intermediary variables
234
235	! zfm_therm(:,:)=0. !init is done inside
236	! zentr_therm(:,:)=0.
237	! zdetr_therm(:,:)=0.
238	! zheatFlux(:,:)=0.
239	! zheatFlux_down(:,:)=0.
240	! zbuoyancyOut(:,:)=0.
241	! zbuoyancyEst(:,:)=0.
242	zzw2(:,:)=0.
243	zmax(:)=0.
244	lmax(:)=0
245	! d_t_the(:,:)=0. !init is done inside
246
247	! d_u_the(:,:)=0. !transported outside
248	! d_v_the(:,:)=0.
249	dq2_the(:,:)=0.
250
251	if (nqmx .ne. 0 .and. ico2 .ne. 0) then
252	d_q_the(:,:,ico2)=0.
253	endif
254
255	CALL thermcell_main_mars(zdt &
256	& ,pplay,pplev,pphi,zzlev,zzlay &
257	& ,zu,zv,zt,pq_therm,q2_therm &
258	& ,d_u_the,d_v_the,d_t_the,d_q_the,dq2_the &
259	& ,zfm_therm,zentr_therm,zdetr_therm,lmax,zmax &
260	& ,r_aspect_thermals &
261	& ,zzw2,fraca,zpopsk &
262	& ,ztla,zheatFlux,zheatFlux_down &
263	& ,zbuoyancyOut,zbuoyancyEst)
264
265	fact=1./REAL(nsplit_thermals)
266
267	d_t_the(:,:)=d_t_the(:,:)ptimestepfact
268	! d_u_the(:,:)=d_u_the(:,:)ptimestepfact
269	! d_v_the(:,:)=d_v_the(:,:)ptimestepfact
270	dq2_the(:,:)=dq2_the(:,:)*fact
271	if (ico2 .ne. 0) then
272	d_q_the(:,:,ico2)=d_q_the(:,:,ico2)ptimestepfact
273	endif
274
275	zmaxth(:)=zmaxth(:)+zmax(:)*fact
276	lmax_real(:)=lmax_real(:)+float(lmax(:))*fact
277	fm_therm(:,:)=fm_therm(:,:) &
278	& +zfm_therm(:,:)*fact
279	entr_therm(:,:)=entr_therm(:,:) &
280	& +zentr_therm(:,:)*fact
281	detr_therm(:,:)=detr_therm(:,:) &
282	& +zdetr_therm(:,:)*fact
283	zfraca(:,:)=zfraca(:,:) + fraca(:,:)*fact
284
285	heatFlux(:,:)=heatFlux(:,:) &
286	& +zheatFlux(:,:)*fact
287	heatFlux_down(:,:)=heatFlux_down(:,:) &
288	& +zheatFlux_down(:,:)*fact
289	buoyancyOut(:,:)=buoyancyOut(:,:) &
290	& +zbuoyancyOut(:,:)*fact
291	buoyancyEst(:,:)=buoyancyEst(:,:) &
292	& +zbuoyancyEst(:,:)*fact
293
294
295	zw2(:,:)=zw2(:,:) + zzw2(:,:)*fact
296
297	! accumulation de la tendance
298
299	d_t_ajs(:,:)=d_t_ajs(:,:)+d_t_the(:,:)
300	! d_u_ajs(:,:)=d_u_ajs(:,:)+d_u_the(:,:)
301	! d_v_ajs(:,:)=d_v_ajs(:,:)+d_v_the(:,:)
302	if (ico2 .ne. 0) then
303	d_q_ajs(:,:,ico2)=d_q_ajs(:,:,ico2)+d_q_the(:,:,ico2)
304	endif
305	! dq2_therm(:,:)=dq2_therm(:,:)+dq2_the(:,:)
306	! incrementation des variables meteo
307
308	zt(:,:) = zt(:,:) + d_t_the(:,:)
309	! zu(:,:) = zu(:,:) + d_u_the(:,:)
310	! zv(:,:) = zv(:,:) + d_v_the(:,:)
311	if (ico2 .ne. 0) then
312	pq_therm(:,:,ico2) = &
313	& pq_therm(:,:,ico2) + d_q_the(:,:,ico2)
314	endif
315	! q2_therm(:,:) = q2_therm(:,:) + dq2_therm(:,:)
316
317
318	ENDDO ! isplit
319	!****************************************************************
320
321	lmax(:)=nint(lmax_real(:))
322	zlmax=MAXVAL(lmax(:))+2
323	if (zlmax .ge. nlayermx) then
324	print*,'thermals have reached last layer of the model'
325	print*,'this is not good !'
326	endif
327
328
329	! Now that we have computed total entrainment and detrainment, we can
330	! advect u, v, and q in thermals. (theta already advected). We can do
331	! that separatly because u,v,and q are not used in thermcell_main for
332	! any thermals-related computation : they are purely passive.
333
334	! mass of cells
335	do l=1,nlayermx
336	masse(:,l)=(pplev(:,l)-pplev(:,l+1))/g
337	enddo
338
339	detrmod(:,:)=0.
340	do l=1,zlmax
341	do ig=1,ngridmx
342	detrmod(ig,l)=fm_therm(ig,l)-fm_therm(ig,l+1) &
343	& +entr_therm(ig,l)
344	if (detrmod(ig,l).lt.0.) then
345	entr_therm(ig,l)=entr_therm(ig,l)-detrmod(ig,l)
346	detrmod(ig,l)=0.
347	endif
348	enddo
349	enddo
350	ndt=10
351	call thermcell_dqup(ngridmx,nlayermx,ptimestep &
352	& ,fm_therm,entr_therm,detrmod, &
353	& masse,zu,d_u_ajs,ndt,zlmax)
354
355	call thermcell_dqup(ngridmx,nlayermx,ptimestep &
356	& ,fm_therm,entr_therm,detrmod, &
357	& masse,zv,d_v_ajs,ndt,zlmax)
358
359	if (nqmx .ne. 0.) then
360	DO iq=1,nqmx
361	if (iq .ne. ico2) then
362	call thermcell_dqup(ngridmx,nlayermx,ptimestep &
363	& ,fm_therm,entr_therm,detrmod, &
364	& masse,pq_therm(:,:,iq),d_q_ajs(:,:,iq),ndt,zlmax)
365	endif
366	ENDDO
367	endif
368
369	if (dtke_thermals) then
370	detrmod(:,:)=0.
371	ndt=10
372	do l=1,zlmax
373	do ig=1,ngridmx
374	detrmod(ig,l)=fm_therm(ig,l)-fm_therm(ig,l+1) &
375	& +entr_therm(ig,l)
376	if (detrmod(ig,l).lt.0.) then
377	entr_therm(ig,l)=entr_therm(ig,l)-detrmod(ig,l)
378	detrmod(ig,l)=0.
379	endif
380	enddo
381	enddo
382	call thermcell_dqup(ngridmx,nlayermx,ptimestep &
383	& ,fm_therm,entr_therm,detrmod, &
384	& masse,q2_therm,dq2_therm,ndt,zlmax)
385	endif
386
387	DO ig=1,ngridmx
388	wmax(ig)=MAXVAL(zw2(ig,:))
389	ENDDO
390
391	! **********************************************************************
392	! **********************************************************************
393	! **********************************************************************
394	! CALLTHERM END
395	! **********************************************************************
396	! **********************************************************************
397	! **********************************************************************
398
399
400	! **********************************************************************
401	! Preparing outputs
402	! **********************************************************************
403
404	do l=1,zlmax
405	pdu_th(:,l)=d_u_ajs(:,l)
406	pdv_th(:,l)=d_v_ajs(:,l)
407	enddo
408
409	if(qtransport_thermals) then
410	if(tracer) then
411	do iq=1,nqmx
412	if (iq .ne. ico2) then
413	do l=1,zlmax
414	pdq_th(:,l,iq)=d_q_ajs(:,l,iq)
415	enddo
416	else
417	do l=1,zlmax
418	pdq_th(:,l,iq)=d_q_ajs(:,l,iq)/ptimestep
419	enddo
420	endif
421	enddo
422	endif
423	endif
424
425	IF(dtke_thermals) THEN
426	DO l=2,nlayermx
427	pbl_dtke(:,l)=0.5*(dq2_therm(:,l-1)+dq2_therm(:,l))
428	ENDDO
429
430	pbl_dtke(:,1)=0.5*dq2_therm(:,1)
431	pbl_dtke(:,nlayermx+1)=0.
432	ENDIF
433
434	do l=1,zlmax
435	pdt_th(:,l)=d_t_ajs(:,l)/ptimestep
436	enddo
437
438
439	! **********************************************************************
440	! Compute the free convection velocity scale for vdifc
441	! **********************************************************************
442
443
444	! Potential temperature gradient
445
446	dteta(:,nlayermx)=0.
447	DO l=1,nlayermx-1
448	DO ig=1, ngridmx
449	dteta(ig,l) = ((zt(ig,l+1)-zt(ig,l))/zpopsk(ig,l)) &
450	& /(zzlay(ig,l+1)-zzlay(ig,l))
451	ENDDO
452	ENDDO
453
454	! Computation of the pbl mixed layer temperature
455
456	DO ig=1, ngridmx
457	ii=MINLOC(abs(dteta(ig,1:lmax(ig))))
458	pbl_teta(ig) = zt(ig,ii(1))/zpopsk(ig,ii(1))
459	ENDDO
460
461	! we must add the heat flux from the diffusion scheme to hfmax
462
463	! compute rho as it is after the diffusion
464
465	rho(:,:)=pplay(:,:) &
466	& /(r(pt(:,:)+pdhdif(:,:)zpopsk(:,:)*ptimestep))
467
468	! integrate -rho*pdhdif
469
470	rpdhd(:,:)=0.
471
472	DO ig=1,ngridmx
473	DO l=1,lmax(ig)
474	rpdhd(ig,l)=0.
475	DO k=1,l
476	rpdhd(ig,l)=rpdhd(ig,l)-rho(ig,k)pdhdif(ig,k) &
477	& (zzlev(ig,k+1)-zzlev(ig,k))
478	ENDDO
479	rpdhd(ig,l)=rpdhd(ig,l)-sensibFlux(ig)/cpp
480	ENDDO
481	ENDDO
482
483	! compute w'teta' from diffusion
484
485	wtdif(:,:)=rpdhd(:,:)/rho(:,:)
486
487	! compute rho as it is after the thermals
488
489	rho(:,:)=pplay(:,:)/(r*(zt(:,:)))
490	! integrate -rho*pdhdif
491
492	DO ig=1,ngridmx
493	DO l=1,lmax(ig)
494	rpdhd(ig,l)=0.
495	DO k=1,l
496	rpdhd(ig,l)=rpdhd(ig,l)-rho(ig,k)(pdt_th(ig,k)/zpopsk(ig,k)) &
497	& (zzlev(ig,k+1)-zzlev(ig,k))
498	ENDDO
499	rpdhd(ig,l)=rpdhd(ig,l)+ &
500	& rho(ig,1)*(heatFlux(ig,1)+heatFlux_down(ig,1))
501	ENDDO
502	ENDDO
503	rpdhd(:,nlayermx)=0.
504
505	! compute w'teta' from thermals
506
507	wtth(:,:)=rpdhd(:,:)/rho(:,:)
508
509	! We get the max heat flux from thermals and add the contribution from the diffusion
510
511	DO ig=1,ngridmx
512	hfmax(ig)=MAXVAL(wtth(ig,:)+wtdif(ig,:))
513	ENDDO
514	! We follow Spiga et. al 2010 (QJRMS)
515	! ------------
516
517	DO ig=1, ngridmx
518	IF (zmax(ig) .gt. 0.) THEN
519	wstar(ig)=(gzmaxth(ig)hfmax(ig)/pbl_teta(ig))**(1./3.)
520	ELSE
521	wstar(ig)=0.
522	ENDIF
523	ENDDO
524
525
526
527	! **********************************************************************
528	! Diagnostics
529	! **********************************************************************
530
531	if(outptherm) then
532	if (ngridmx .eq. 1) then
533	call WRITEDIAGFI(ngridmx,'entr_therm','entrainement thermique',&
534	& 'kg/m-2',1,entr_therm)
535	call WRITEDIAGFI(ngridmx,'detr_therm','detrainement thermique',&
536	& 'kg/m-2',1,detr_therm)
537	call WRITEDIAGFI(ngridmx,'fm_therm','flux masse thermique',&
538	& 'kg/m-2',1,fm_therm)
539	call WRITEDIAGFI(ngridmx,'zw2','vitesse verticale thermique',&
540	& 'm/s',1,zw2)
541	call WRITEDIAGFI(ngridmx,'heatFlux_up','heatFlux_updraft',&
542	& 'SI',1,heatFlux)
543	call WRITEDIAGFI(ngridmx,'heatFlux_down','heatFlux_downdraft',&
544	& 'SI',1,heatFlux_down)
545	call WRITEDIAGFI(ngridmx,'fraca','fraction coverage',&
546	& 'percent',1,fraca)
547	call WRITEDIAGFI(ngridmx,'buoyancyOut','buoyancyOut',&
548	& 'm.s-2',1,buoyancyOut)
549	call WRITEDIAGFI(ngridmx,'buoyancyEst','buoyancyEst',&
550	& 'm.s-2',1,buoyancyEst)
551	call WRITEDIAGFI(ngridmx,'d_t_th', &
552	& 'tendance temp TH','K',1,d_t_ajs)
553	call WRITEDIAGFI(ngridmx,'d_q_th', &
554	& 'tendance traceur TH','kg/kg',1,d_q_ajs)
555	call WRITEDIAGFI(ngridmx,'zmax', &
556	& 'pbl height','m',0,zmaxth)
557	call WRITEDIAGFI(ngridmx,'d_u_th', &
558	& 'tendance moment','m/s',1,pdu_th)
559	call WRITEDIAGFI(ngridmx,'wtdif', &
560	& 'heat flux from diffusion','K.m/s',1,wtdif)
561	call WRITEDIAGFI(ngridmx,'wtth', &
562	& 'heat flux from thermals','K.m/s',1,wtth)
563	call WRITEDIAGFI(ngridmx,'wttot', &
564	& 'heat flux PBL','K.m/s',1,wtdif(:,:)+wtth(:,:))
565
566	else
567
568	call WRITEDIAGFI(ngridmx,'entr_therm','entrainement thermique',&
569	& 'kg/m-2',3,entr_therm)
570	call WRITEDIAGFI(ngridmx,'detr_therm','detrainement thermique',&
571	& 'kg/m-2',3,detr_therm)
572	call WRITEDIAGFI(ngridmx,'fm_therm','flux masse thermique',&
573	& 'kg/m-2',3,fm_therm)
574	call WRITEDIAGFI(ngridmx,'zw2','vitesse verticale thermique',&
575	& 'm/s',3,zw2)
576	call WRITEDIAGFI(ngridmx,'heatFlux','heatFlux',&
577	& 'SI',3,heatFlux)
578	call WRITEDIAGFI(ngridmx,'buoyancyOut','buoyancyOut',&
579	& 'SI',3,buoyancyOut)
580	call WRITEDIAGFI(ngridmx,'d_t_th', &
581	& 'tendance temp TH','K',3,d_t_ajs)
582
583	endif
584	endif
585
586	END

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

Download in other formats:

Original Format