38 USE yomhook
,ONLY : lhook, dr_hook
39 USE parkind1
,ONLY : jprb
46 INTEGER,
INTENT(IN) :: kgrid_par
47 INTEGER,
INTENT(IN) :: kl
48 INTEGER,
INTENT(IN) :: knear_nbr
49 REAL,
DIMENSION(KGRID_PAR),
INTENT(IN) :: pgrid_par
50 INTEGER,
DIMENSION(:,:),
POINTER :: knear
62 REAL(KIND=JPRB) :: zhook_handle
65 IF (lhook) CALL dr_hook(
'GET_NEAR_MESHES_LONLAT_ROT',0,zhook_handle)
70 idist = int(sqrt(float(knear_nbr)))
72 IF (ilon*ilat==kl)
THEN
76 jl = jlon + ilon * (jlat-1)
78 DO jx=-(idist-1)/2,idist/2
79 DO jy=-(idist-1)/2,idist/2
80 IF (jlon+jx>0 .AND. jlon+jx<ilon+1 .AND. jlat+jy>0 .AND. jlat+jy<ilat+1)
THEN
82 knear(jl,icount) = (jlon+jx) + ilon * (jlat+jy-1)
89 IF (lhook) CALL dr_hook(
'GET_NEAR_MESHES_LONLAT_ROT',1,zhook_handle)
subroutine get_gridtype_lonlat_rot(PGRID_PAR, PWEST, PSOUTH, PDLON, PDLAT, PPOLON, PPOLAT, KLON, KLAT, KL, PLON, PLAT)
subroutine get_near_meshes_lonlat_rot(KGRID_PAR, KL, PGRID_PAR, KNEAR_NBR, KNEAR)