41 USE modd_pgdwork, ONLY : nsize, xsumval, xsumval2, xsumval3, &
42 xmean_work, xstd_work, xskew_work, &
46 USE yomhook
,ONLY : lhook, dr_hook
47 USE parkind1
,ONLY : jprb
57 REAL,
DIMENSION(SIZE(NSIZE)) :: zsize
60 REAL(KIND=JPRB) :: zhook_handle
63 IF (lhook) CALL dr_hook(
'AVERAGE2_CTI',0,zhook_handle)
64 zsize(:)=
REAL(nsize(:))
73 xmean_work(:)=xsumval(:)/zsize(:)
80 WHERE (xmax_work(:)-xmin_work(:)>=1.0)
81 xstd_work(:)=sqrt( max(0.,xsumval2(:)/nsize(:) - xmean_work(:)*xmean_work(:)) )
91 WHERE(xstd_work(:)>0.0)
93 xskew_work(:)=xsumval3(:)-zsize(:)*xmean_work(:)*xmean_work(:)*xmean_work(:) &
94 -3.0*zsize(:)*xmean_work(:)*xstd_work(:)*xstd_work(:)
96 xskew_work(:)=xskew_work(:)/(zsize(:)*xstd_work(:)*xstd_work(:)*xstd_work(:))
103 IF (lhook) CALL dr_hook(
'AVERAGE2_CTI',1,zhook_handle)