33 USE yomhook
,ONLY : lhook, dr_hook
34 USE parkind1
,ONLY : jprb
41 INTEGER,
INTENT(IN) :: kgrid_par
42 INTEGER,
INTENT(IN) :: kl
43 INTEGER,
INTENT(IN) :: kc
44 REAL,
DIMENSION(KGRID_PAR),
INTENT(IN) :: pgrid_par
45 REAL,
DIMENSION(KL,KC),
INTENT(OUT) :: pcorner_lat
46 REAL,
DIMENSION(KL,KC),
INTENT(OUT) :: pcorner_lon
60 REAL,
DIMENSION(KL) :: zlon
61 REAL,
DIMENSION(KL) :: zlat
63 REAL(KIND=JPRB) :: zhook_handle
67 IF (lhook) CALL dr_hook(
'GET_MESH_CORNER_LONLAT_REG',0,zhook_handle)
75 CALL
abor1_sfx(
'GET_MESH_CORNER_LONLAT_REG: WRONG NUMBER OF POINT')
80 zdlon = (zlonmax-zlonmin) / float(ilon)
81 zdlat = (zlatmax-zlatmin) / float(ilat)
95 pcorner_lon(:,1) = zlon(:)-zdlon/2.
96 pcorner_lat(:,1) = zlat(:)-zdlat/2.
98 pcorner_lon(:,3) = zlon(:)+zdlon/2.
99 pcorner_lat(:,3) = zlat(:)+zdlat/2.
101 pcorner_lon(:,2) = pcorner_lon(:,3)
102 pcorner_lat(:,2) = pcorner_lat(:,1)
104 pcorner_lon(:,4) = pcorner_lon(:,1)
105 pcorner_lat(:,4) = pcorner_lat(:,3)
107 IF (lhook) CALL dr_hook(
'GET_MESH_CORNER_LONLAT_REG',1,zhook_handle)
subroutine abor1_sfx(YTEXT)
subroutine get_mesh_corner_lonlat_reg(KGRID_PAR, KL, KC, PGRID_PAR, PCORNER_LAT, PCORNER_LON)
subroutine get_gridtype_lonlat_reg(PGRID_PAR, PLONMIN, PLONMAX, PLATMIN, PLATMAX, KLON, KLAT, KL, PLON, PLAT)