14 USE yomhook
,ONLY : lhook, dr_hook
15 USE parkind1
,ONLY : jprb
21 TYPE(slt_t),
INTENT(INOUT) :: slt
23 CHARACTER(LEN=6),
INTENT(IN) :: hprogram
26 CHARACTER(LEN=4) :: crgunits
29 REAL(KIND=JPRB) :: zhook_handle
32 IF (lhook) CALL dr_hook(
'INIT_SLT',0,zhook_handle)
35 ALLOCATE(slt%XEMISRADIUS_SLT(nsltmde))
36 ALLOCATE(slt%XEMISSIG_SLT (nsltmde))
41 IF(cemisparam_slt.eq.
"Vig01")
THEN
43 xemisradius_ini_slt(:) = (/ 0.2, 2.0, 12. /)
44 xemissig_ini_slt(:) = (/ 1.9, 2.0, 3.00 /)
47 xemisradius_ini_slt(:) = 0.5*(/0.28, 2.25, 15.32/)
48 xemissig_ini_slt(:) = (/1.59, 2.00, 2.00 /)
53 jmode_idx = jorder_slt(jmode)
55 slt%XEMISSIG_SLT (jmode) = xemissig_ini_slt(jmode_idx)
56 slt%XEMISRADIUS_SLT(jmode) = xemisradius_ini_slt(jmode_idx)
58 IF (crgunits==
"MASS") &
59 slt%XEMISRADIUS_SLT(jmode) = slt%XEMISRADIUS_SLT(jmode) * exp(-3.d0 * (log(slt%XEMISSIG_SLT(jmode)))**2)
63 IF (lhook) CALL dr_hook(
'INIT_SLT',1,zhook_handle)
subroutine init_slt(SLT, HPROGRAM)