5 SUBROUTINE add_noise(PADDTIMECORR,PASSIM_WINH,PWHITE_NOISE,PRED_NOISE)
7 USE yomhook
, ONLY : lhook,dr_hook
8 USE parkind1
, ONLY : jprb
12 REAL,
INTENT(IN) :: paddtimecorr
13 REAL,
INTENT(IN) :: pwhite_noise
14 REAL,
INTENT(IN) :: passim_winh
15 REAL,
INTENT(INOUT) :: pred_noise
18 REAL(KIND=JPRB) :: zhook_handle
20 IF (lhook) CALL dr_hook(
'ADD_NOISE',0,zhook_handle)
27 zalpha = 1./(1. + ((passim_winh/24.0)/paddtimecorr))
29 zmu = sqrt(1. - zalpha**2)*pwhite_noise
31 pred_noise = zalpha*pred_noise + zmu
35 IF (lhook) CALL dr_hook(
'ADD_NOISE',1,zhook_handle)
subroutine add_noise(PADDTIMECORR, PASSIM_WINH, PWHITE_NOISE, PRED_NOISE)