46 USE yomhook
,ONLY : lhook, dr_hook
47 USE parkind1
,ONLY : jprb
54 INTEGER,
INTENT(IN) :: kgrid_par
55 INTEGER,
INTENT(IN) :: kl
56 REAL,
DIMENSION(KGRID_PAR),
INTENT(IN) :: pgrid_par
57 REAL,
DIMENSION(KL),
INTENT(OUT) :: plat
58 REAL,
DIMENSION(KL),
INTENT(OUT) :: plon
59 REAL,
DIMENSION(KL),
INTENT(OUT) :: pmesh_size
60 REAL,
DIMENSION(KL),
INTENT(OUT) :: pdir
65 REAL,
DIMENSION(:),
ALLOCATABLE :: zx
66 REAL,
DIMENSION(:),
ALLOCATABLE :: zy
67 REAL,
DIMENSION(:),
ALLOCATABLE :: zdx
68 REAL,
DIMENSION(:),
ALLOCATABLE :: zdy
69 REAL,
DIMENSION(:),
ALLOCATABLE :: zdlat
70 REAL,
DIMENSION(:),
ALLOCATABLE :: zdlon
71 REAL(KIND=JPRB) :: zhook_handle
77 IF (lhook) CALL dr_hook(
'LATLON_GRIDTYPE_LONLATVAL',0,zhook_handle)
78 ALLOCATE(zx(
SIZE(plat)))
79 ALLOCATE(zy(
SIZE(plat)))
80 ALLOCATE(zdx(
SIZE(plat)))
81 ALLOCATE(zdy(
SIZE(plat)))
82 ALLOCATE(zdlon(
SIZE(plat)))
83 ALLOCATE(zdlat(
SIZE(plat)))
104 pmesh_size(:) = xradius**2 * xpi/180.*(zdlon(:)) &
105 * (sin((plat(:)+zdlat(:)/2.)*xpi/180.)-sin((plat(:)-zdlat(:)/2.)*xpi/180.))
125 IF (lhook) CALL dr_hook(
'LATLON_GRIDTYPE_LONLATVAL',1,zhook_handle)
subroutine get_gridtype_lonlatval(PGRID_PAR, KL, PX, PY, PDX, PDY)
subroutine latlon_gridtype_lonlatval(KGRID_PAR, KL, PGRID_PAR, PLAT, PLON, PMESH_SIZE, PDIR)
subroutine latlon_lonlatval(PX, PY, PLAT, PLON)