32 nvalnbr, nvalcount, xvallist
36 USE yomhook
,ONLY : lhook, dr_hook
37 USE parkind1
,ONLY : jprb
44 REAL,
DIMENSION(:),
INTENT(INOUT) :: ppgdarray
45 REAL(KIND=JPRB) :: zhook_handle
56 IF (lhook) CALL dr_hook(
'AVERAGE2_MESH',0,zhook_handle)
61 ppgdarray(:)=xsumval(:)/nsize(:)
66 ppgdarray(:)=nsize(:)/xsumval(:)
71 ppgdarray(:)=xcdref/exp(sqrt(nsize(:)/xsumval(:)))
75 DO jloop=1,
SIZE(nsize)
76 IF(nsize(jloop)==0) cycle
80 DO jval=1,nvalnbr(jloop)
81 IF (nvalcount(jloop,jval)>imax)
THEN
82 imax=nvalcount(jloop,jval)
87 ppgdarray(jloop)=xvallist(jloop,ival)
91 IF (lhook) CALL dr_hook(
'AVERAGE2_MESH',1,zhook_handle)
subroutine average2_mesh(PPGDARRAY)