51 USE yomhook
,ONLY : lhook, dr_hook
52 USE parkind1
,ONLY : jprb
59 INTEGER,
INTENT(IN) :: kgrid_par
60 REAL,
DIMENSION(KGRID_PAR),
INTENT(IN) :: pgrid_par
61 LOGICAL,
DIMENSION(720,360),
INTENT(OUT) :: olatlonmask
67 REAL,
DIMENSION(720,360) :: zlon_mask
68 REAL,
DIMENSION(720,360) :: zlat_mask
71 REAL(KIND=JPRB) :: zhook_handle
74 IF (lhook) CALL dr_hook(
'LATLONMASK_CARTESIAN',0,zhook_handle)
79 olatlonmask(:,:) = .false.
87 zlon_mask(:,:)= spread( (/ ( jlon /2. - 0.25 , jlon=1,720 ) /) , dim=2, ncopies=360 )
88 zlat_mask(:,:)= spread( (/ ( (jlat-180)/2. - 0.25 , jlat=1,360 ) /) , dim=1, ncopies=720 )
93 zlon_mask(:,:)=zlon_mask(:,:)+nint((zlon0-zlon_mask(:,:))/360.)*360.
96 WHERE ( zlon_mask(:,:) -0.25 <= zlon0 .AND. zlon0 <= zlon_mask(:,:) +0.25 &
97 .AND. zlat_mask(:,:) -0.25 <= zlon0 .AND. zlat0 <= zlat_mask(:,:) +0.25 )
98 olatlonmask(:,:) = .true.
100 IF (lhook) CALL dr_hook(
'LATLONMASK_CARTESIAN',1,zhook_handle)
subroutine latlonmask_cartesian(KGRID_PAR, PGRID_PAR, OLATLONMASK)
subroutine get_gridtype_cartesian(PGRID_PAR, PLAT0, PLON0, KIMAX, KJMAX, PX, PY, PDX, PDY, KL)