47 USE yomhook
,ONLY : lhook, dr_hook
48 USE parkind1
,ONLY : jprb
55 INTEGER,
INTENT(IN) :: kgrid_par
56 INTEGER,
INTENT(IN) :: kl
57 REAL,
DIMENSION(KGRID_PAR),
INTENT(IN) :: pgrid_par
58 REAL,
DIMENSION(KL),
INTENT(OUT) :: plat
59 REAL,
DIMENSION(KL),
INTENT(OUT) :: plon
60 REAL,
DIMENSION(KL),
INTENT(OUT) :: pmesh_size
61 REAL,
DIMENSION(KL),
INTENT(OUT) :: pdir
66 REAL,
DIMENSION(:),
ALLOCATABLE :: zx
67 REAL,
DIMENSION(:),
ALLOCATABLE :: zy
68 REAL,
DIMENSION(:),
ALLOCATABLE :: zmap
69 REAL,
DIMENSION(:),
ALLOCATABLE :: zdx
70 REAL,
DIMENSION(:),
ALLOCATABLE :: zdy
71 REAL,
DIMENSION(:),
ALLOCATABLE :: zydelty
72 REAL,
DIMENSION(:),
ALLOCATABLE :: zlatdy
73 REAL,
DIMENSION(:),
ALLOCATABLE :: zlondy
76 REAL(KIND=JPRB) :: zhook_handle
82 IF (lhook) CALL dr_hook(
'LATLON_GRIDTYPE_IGN',0,zhook_handle)
83 ALLOCATE(zx(
SIZE(plat)))
84 ALLOCATE(zy(
SIZE(plat)))
85 ALLOCATE(zdx(
SIZE(plat)))
86 ALLOCATE(zdy(
SIZE(plat)))
88 CALL
get_gridtype_ign(pgrid_par,klambert=ilambert,px=zx,py=zy,pdx=zdx,pdy=zdy )
105 ALLOCATE(zmap(
SIZE(plat)))
112 pmesh_size(:) = zdx(:) * zdy(:) / zmap(:)**2
120 ALLOCATE(zydelty(
SIZE(plat)))
121 ALLOCATE(zlatdy(
SIZE(plat)))
122 ALLOCATE(zlondy(
SIZE(plat)))
124 CALL
latlon_ign(ilambert,zx,zydelty,zlatdy,zlondy)
126 pdir(:)= atan( (xa(ilambert)*(zlondy(:)-plon(:))*xpi/180.) / xdelty) * xpi/180.
137 IF (lhook) CALL dr_hook(
'LATLON_GRIDTYPE_IGN',1,zhook_handle)
subroutine latlon_gridtype_ign(KGRID_PAR, KL, PGRID_PAR, PLAT, PLON, PMESH_SIZE, PDIR)
subroutine get_gridtype_ign(PGRID_PAR, KLAMBERT, KL, PX, PY, PDX, PDY, KDIMX, KDIMY, PXALL, PYALL)
subroutine latlon_ign(KLAMBERT, PX, PY, PLAT, PLON)
subroutine map_factor_ign(KLAMBERT, PX, PY, PMAP)