Context Navigation

calltherm_interface.F90 @ 288

Last change on this file since 288 was 284, checked in by acolaitis, 13 years ago

--- AC 07/09/2011 ---

Added new flag for the Richardson-based surface layer :

callrichsl, can be changed in callphys.def

One should always use the thermals model when using this surface layer model.
Somes cases (weakly unstable with low winds), when not using thermals, won't be well represented by the
Richardson surface layer. This stands for Mesoscale and Gcm but not for LES model.

Correct configs :

callrichsl = .true.
calltherm = .true.

callrichsl = .false.
calltherm = .false.

callrichsl = .false.
calltherm = .true.

Previously unstable config :

callrichsl = .true.
calltherm = .false.

To be able to run without thermals and with the new surface layer, a modification has been made to

physiq.F to account for gustiness in GCM and MESOSCALE for negative Richardson, so that :

callrichsl = .true.
calltherm = .false.

can now be used without problems, but is not recommended.

Consequently, callrichsl = .false. is now the default configuration for thermals.

We recall the available options in callphys.def for thermals :

outptherm = BOOLEAN (.false. by default) : outputs thermals related quantities (lots of diagfi)
nsplit_thermals = INTEGER (50 by default in gcm, 2 in mesoscale) : subtimestep for thermals model.

It is recommended to use at least 40 in the gcm, and at least 2 in the mesoscale.
The user can lower these values but should check it's log for anomalies or errors regarding
tracer transport in the thermals, or "granulosity" in the outputs for wmax, lmax and hfmax.

File size: 7.7 KB

Rev	Line
[161]	1	!
	2	! AC 2011-01-05
	3	!
[185]	4	SUBROUTINE calltherm_interface (firstcall, &
[161]	5	& long,lati,zzlev,zzlay, &
	6	& ptimestep,pu,pv,pt,pq,pdu,pdv,pdt,pdq,q2, &
[185]	7	& pplay,pplev,pphi,zpopsk, &
	8	& pdu_th,pdv_th,pdt_th,pdq_th,lmax_th,zmax_th,pbl_dtke,hfmax,wmax)
[161]	9
	10	USE ioipsl_getincom
	11
	12	implicit none
	13	#include "callkeys.h"
[185]	14	#include "dimensions.h"
	15	#include "dimphys.h"
	16
[161]	17	!--------------------------------------------------------
	18	! Variables d'entree
	19	!--------------------------------------------------------
	20
	21	REAL, INTENT(IN) :: ptimestep
[185]	22	REAL, INTENT(IN) :: pplev(ngridmx,nlayermx+1),pplay(ngridmx,nlayermx)
	23	REAL, INTENT(IN) :: pphi(ngridmx,nlayermx)
	24	REAL, INTENT(IN) :: pu(ngridmx,nlayermx),pv(ngridmx,nlayermx)
	25	REAL, INTENT(IN) :: pt(ngridmx,nlayermx),pq(ngridmx,nlayermx,nqmx)
	26	REAL, INTENT(IN) :: zzlay(ngridmx,nlayermx)
	27	REAL, INTENT(IN) :: zzlev(ngridmx,nlayermx+1)
[161]	28	LOGICAL, INTENT(IN) :: firstcall
[185]	29	REAL, INTENT(IN) :: pdu(ngridmx,nlayermx),pdv(ngridmx,nlayermx)
	30	REAL, INTENT(IN) :: pdq(ngridmx,nlayermx,nqmx),pdt(ngridmx,nlayermx)
	31	REAL, INTENT(IN) :: q2(ngridmx,nlayermx+1)
	32	REAL, INTENT(IN) :: long(ngridmx),lati(ngridmx)
	33	REAL, INTENT(IN) :: zpopsk(ngridmx,nlayermx)
[161]	34
	35	!--------------------------------------------------------
	36	! Variables de sortie (ou entree/sortie)
	37	!--------------------------------------------------------
	38
[185]	39	REAL pdu_th(ngridmx,nlayermx),pdv_th(ngridmx,nlayermx)
	40	REAL pdt_th(ngridmx,nlayermx),pdq_th(ngridmx,nlayermx,nqmx)
	41	INTEGER lmax_th(ngridmx)
	42	REAL zmax_th(ngridmx)
	43	REAL pbl_dtke(ngridmx,nlayermx+1)
[161]	44
	45	!--------------------------------------------------------
	46	! Variables du thermique
	47	!--------------------------------------------------------
[185]	48	REAL u_seri(ngridmx,nlayermx), v_seri(ngridmx,nlayermx)
	49	REAL t_seri(ngridmx,nlayermx)
	50	REAL d_t_ajs(ngridmx,nlayermx)
	51	REAL d_u_ajs(ngridmx,nlayermx), d_q_ajs(ngridmx,nlayermx,nqmx)
	52	REAL d_v_ajs(ngridmx,nlayermx)
	53	REAL fm_therm(ngridmx,nlayermx+1), entr_therm(ngridmx,nlayermx)
	54	REAL detr_therm(ngridmx,nlayermx)
	55	REAL zw2(ngridmx,nlayermx+1)
	56	REAL fraca(ngridmx,nlayermx+1)
	57	REAL ztla(ngridmx,nlayermx)
	58	REAL q_therm(ngridmx,nlayermx), pq_therm(ngridmx,nlayermx,nqmx)
	59	REAL dq_therm(ngridmx,nlayermx), dq_thermdown(ngridmx,nlayermx)
	60	REAL q2_therm(ngridmx,nlayermx), dq2_therm(ngridmx,nlayermx)
[161]	61
	62	LOGICAL qtransport_thermals,dtke_thermals
	63
	64	INTEGER l,ig,iq
	65
	66	! Variable de diagnostique : flux de chaleur vertical
	67
[185]	68	REAL heatFlux(ngridmx,nlayermx)
	69	REAL heatFlux_down(ngridmx,nlayermx)
	70	REAL buoyancyOut(ngridmx,nlayermx)
	71	REAL buoyancyEst(ngridmx,nlayermx)
	72	REAL hfmax(ngridmx),wmax(ngridmx)
[161]	73
	74	!---------------------------------------------------------
	75	!---------------------------------------------------------
	76	! **********************************************************************
	77	! Thermique
	78	! **********************************************************************
	79
	80	! Initialisation des sorties
	81
	82	lmax_th(:)=1
	83	pdu_th(:,:)=0.
	84	pdv_th(:,:)=0.
	85	pdt_th(:,:)=0.
	86	entr_therm(:,:)=0.
	87	detr_therm(:,:)=0.
	88	q2_therm(:,:)=0.
	89	dq2_therm(:,:)=0.
	90	dq_therm(:,:)=0.
	91	dq_thermdown(:,:)=0.
	92	ztla(:,:)=0.
	93	pbl_dtke(:,:)=0.
	94	fm_therm(:,:)=0.
	95	zw2(:,:)=0.
	96	fraca(:,:)=0.
	97	if (tracer) then
	98	pdq_th(:,:,:)=0.
	99	end if
	100
	101	! Dans le model terrestres, les seri sont des q+dq tendances déja cumulées. Il n'y a donc pas de
	102	! cumulage à l'intérieur de la routine comme dans le model martien. On le fait ici :
	103
	104	u_seri(:,:)=pu(:,:)+pdu(:,:)*ptimestep
	105	v_seri(:,:)=pv(:,:)+pdv(:,:)*ptimestep
	106	t_seri(:,:)=pt(:,:)+pdt(:,:)*ptimestep
	107
	108	pq_therm(:,:,:)=0.
[165]	109	qtransport_thermals=.true.
[161]	110	call getin("qtransport_thermals",qtransport_thermals)
	111	if(qtransport_thermals) then
	112	if(tracer) then
	113	pq_therm(:,:,:)=pq(:,:,:)+pdq(:,:,:)*ptimestep
	114	endif
	115	endif
	116
	117	d_t_ajs(:,:)=0.
	118	d_u_ajs(:,:)=0.
	119	d_v_ajs(:,:)=0.
	120	d_q_ajs(:,:,:)=0.
	121	heatFlux(:,:)=0.
	122	heatFlux_down(:,:)=0.
	123	buoyancyOut(:,:)=0.
	124	buoyancyEst(:,:)=0.
	125
[165]	126	dtke_thermals=.false.
[161]	127	call getin("dtke_thermals",dtke_thermals)
	128	if(dtke_thermals) then
	129
[185]	130	DO l=1,nlayermx
[161]	131	q2_therm(:,l)=0.5*(q2(:,l)+q2(:,l+1))
	132	ENDDO
	133	endif
	134
[185]	135	CALL calltherm_mars(ptimestep,zzlev,zzlay &
[161]	136	& ,pplay,pplev,pphi &
	137	& ,u_seri,v_seri,t_seri,pq_therm, q2_therm &
	138	& ,d_u_ajs,d_v_ajs,d_t_ajs,d_q_ajs, dq2_therm &
	139	& ,fm_therm,entr_therm,detr_therm &
[185]	140	& ,lmax_th,zmax_th &
[161]	141	& ,zw2,fraca &
	142	& ,zpopsk,ztla,heatFlux,heatFlux_down &
[173]	143	& ,buoyancyOut,buoyancyEst,hfmax,wmax)
[161]	144
	145	! Accumulation des tendances. On n'accumule pas les quantités de traceurs car celle ci n'a pas du changer
	146	! étant donné qu'on ne prends en compte que q_seri de la vap d'eau = 0
	147
	148	! INCREMENTATION : les d_u_ sont des tendances alors que les pdu sont des dérivees, attention !
	149
	150	pdu_th(:,:)=d_u_ajs(:,:)/ptimestep
	151	pdv_th(:,:)=d_v_ajs(:,:)/ptimestep
	152	pdt_th(:,:)=d_t_ajs(:,:)/ptimestep
	153	if(qtransport_thermals) then
	154	if(tracer) then
	155	pdq_th(:,:,:)=d_q_ajs(:,:,:)/ptimestep
	156	endif
	157	endif
	158
	159
[185]	160	DO l=2,nlayermx
[161]	161	pbl_dtke(:,l)=0.5*(dq2_therm(:,l-1)+dq2_therm(:,l))/ptimestep
	162	ENDDO
	163
	164	pbl_dtke(:,1)=0.5*dq2_therm(:,1)/ptimestep
[185]	165	pbl_dtke(:,nlayermx+1)=0.
[161]	166	!! DIAGNOSTICS
	167
	168	if(outptherm) then
[185]	169	if (ngridmx .eq. 1) then
	170	call WRITEDIAGFI(ngridmx,'entr_therm','entrainement thermique',&
[161]	171	& 'kg/m-2',1,entr_therm)
[185]	172	call WRITEDIAGFI(ngridmx,'detr_therm','detrainement thermique',&
[161]	173	& 'kg/m-2',1,detr_therm)
[185]	174	call WRITEDIAGFI(ngridmx,'fm_therm','flux masse thermique',&
[161]	175	& 'kg/m-2',1,fm_therm)
[185]	176	call WRITEDIAGFI(ngridmx,'zw2','vitesse verticale thermique',&
[161]	177	& 'm/s',1,zw2)
[185]	178	call WRITEDIAGFI(ngridmx,'heatFlux_up','heatFlux_updraft',&
[161]	179	& 'SI',1,heatFlux)
[185]	180	call WRITEDIAGFI(ngridmx,'heatFlux_down','heatFlux_downdraft',&
[161]	181	& 'SI',1,heatFlux_down)
[185]	182	call WRITEDIAGFI(ngridmx,'fraca','fraction coverage',&
[161]	183	& 'percent',1,fraca)
[185]	184	call WRITEDIAGFI(ngridmx,'buoyancyOut','buoyancyOut',&
[161]	185	& 'm.s-2',1,buoyancyOut)
[185]	186	call WRITEDIAGFI(ngridmx,'buoyancyEst','buoyancyEst',&
[161]	187	& 'm.s-2',1,buoyancyEst)
[185]	188	call WRITEDIAGFI(ngridmx,'d_t_th', &
[161]	189	& 'tendance temp TH','K',1,d_t_ajs)
[185]	190	call WRITEDIAGFI(ngridmx,'zmax', &
	191	& 'pbl height','m',0,zmax_th)
[161]	192	else
	193
[185]	194	call WRITEDIAGFI(ngridmx,'entr_therm','entrainement thermique',&
[161]	195	& 'kg/m-2',3,entr_therm)
[185]	196	call WRITEDIAGFI(ngridmx,'detr_therm','detrainement thermique',&
[161]	197	& 'kg/m-2',3,detr_therm)
[185]	198	call WRITEDIAGFI(ngridmx,'fm_therm','flux masse thermique',&
[161]	199	& 'kg/m-2',3,fm_therm)
[185]	200	call WRITEDIAGFI(ngridmx,'zw2','vitesse verticale thermique',&
[161]	201	& 'm/s',3,zw2)
[185]	202	call WRITEDIAGFI(ngridmx,'heatFlux','heatFlux',&
[161]	203	& 'SI',3,heatFlux)
[185]	204	call WRITEDIAGFI(ngridmx,'buoyancyOut','buoyancyOut',&
[161]	205	& 'SI',3,buoyancyOut)
[185]	206	call WRITEDIAGFI(ngridmx,'d_t_th', &
[161]	207	& 'tendance temp TH','K',3,d_t_ajs)
	208
	209	endif
	210	endif
	211
	212	END

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

Download in other formats:

Original Format