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
73 INTEGER :: jloop, jover
74 INTEGER :: icoverclass
75 REAL(KIND=JPRB) :: zhook_handle
82 IF (lhook) CALL dr_hook(
'AVERAGE1_COVER',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 icoverclass = nint(pvalue(jloop))
112 IF (icoverclass<1 .OR. icoverclass >
SIZE(xsumcover,2) ) cycle
117 nsize(iindex(jover,jloop))=nsize(iindex(jover,jloop))+1
122 xsumcover(iindex(jover,jloop),icoverclass)=xsumcover(iindex(jover,jloop),icoverclass)+1.
128 IF (lhook) CALL dr_hook(
'AVERAGE1_COVER',1,zhook_handle)
subroutine get_mesh_index(KLUOUT, KNBLINES, PLAT, PLON, KINDEX, PVALUE, PNODATA, KSSO, KISSOX, KISSOY)
subroutine average1_cover(KLUOUT, KNBLINES, PLAT, PLON, PVALUE, PNODATA)