6 SUBROUTINE average1_cti(KLUOUT,KNBLINES,PLAT,PLON,PVALUE,PNODATA)
43 USE modd_pgdwork, ONLY : xsumval, xsumval2, xsumval3, nsize, &
46 USE modi_get_mesh_index
50 USE yomhook
,ONLY : lhook, dr_hook
51 USE parkind1
,ONLY : jprb
58 INTEGER,
INTENT(IN) :: kluout
59 INTEGER,
INTENT(IN) :: knblines
60 REAL,
DIMENSION(:),
INTENT(IN) :: plat
61 REAL,
DIMENSION(:),
INTENT(IN) :: plon
62 REAL,
DIMENSION(:),
INTENT(IN) :: pvalue
63 REAL,
OPTIONAL,
INTENT(IN) :: pnodata
68 INTEGER,
DIMENSION(NOVMX,SIZE(PLAT)) :: iindex
71 REAL,
DIMENSION(SIZE(PLAT)) :: zvalue
74 INTEGER :: jloop, jover
75 REAL(KIND=JPRB) :: zhook_handle
82 IF (lhook) CALL dr_hook(
'AVERAGE1_CTI',0,zhook_handle)
84 IF (present(pnodata))
THEN
87 CALL
get_mesh_index(kluout,knblines,plat,plon,iindex,zvalue,znodata)
98 DO jloop = 1 ,
SIZE(plat)
105 IF (iindex(jover,jloop)==0) cycle bloop
110 nsize(iindex(jover,jloop))=nsize(iindex(jover,jloop))+1
115 xsumval(iindex(jover,jloop))=xsumval(iindex(jover,jloop))+pvalue(jloop)
120 xsumval2(iindex(jover,jloop))=xsumval2(iindex(jover,jloop))+pvalue(jloop)**2
126 xsumval3(iindex(jover,jloop))=xsumval3(iindex(jover,jloop))+pvalue(jloop)**3
132 xmax_work(iindex(jover,jloop))=max(xmax_work(iindex(jover,jloop)),pvalue(jloop))
138 xmin_work(iindex(jover,jloop))=min(xmin_work(iindex(jover,jloop)),pvalue(jloop))
145 IF (lhook) CALL dr_hook(
'AVERAGE1_CTI',1,zhook_handle)
subroutine get_mesh_index(KLUOUT, KNBLINES, PLAT, PLON, KINDEX, PVALUE, PNODATA, KSSO, KISSOX, KISSOY)
subroutine average1_cti(KLUOUT, KNBLINES, PLAT, PLON, PVALUE, PNODATA)