7 hprogram,kgrid_par,klu,oread,ksize,pgrid_par,kresp,hdir)
50 USE yomhook
,ONLY : lhook, dr_hook
51 USE parkind1
,ONLY : jprb
62 CHARACTER(LEN=6),
INTENT(IN) :: hprogram
63 INTEGER,
INTENT(INOUT) :: kgrid_par
64 INTEGER,
INTENT(IN) :: klu
65 LOGICAL,
INTENT(IN) :: oread
66 INTEGER,
INTENT(IN) :: ksize
67 REAL,
DIMENSION(KSIZE),
INTENT(OUT) :: pgrid_par
68 INTEGER,
INTENT(OUT) :: kresp
69 CHARACTER(LEN=1),
INTENT(IN) :: hdir
81 INTEGER,
DIMENSION(:),
ALLOCATABLE :: inlopa
83 REAL,
DIMENSION(KLU) :: zlat
84 REAL,
DIMENSION(KLU) :: zlon
85 REAL,
DIMENSION(KLU) :: zlat_xy
86 REAL,
DIMENSION(KLU) :: zlon_xy
87 REAL,
DIMENSION(KLU) :: zmesh_size
89 REAL,
DIMENSION(KLU) :: zlatsup
90 REAL,
DIMENSION(KLU) :: zlonsup
91 REAL,
DIMENSION(KLU) :: zlatinf
92 REAL,
DIMENSION(KLU) :: zloninf
96 REAL,
DIMENSION(:),
POINTER :: zgrid_par=>null()
98 REAL(KIND=JPRB) :: zhook_handle
104 IF (lhook) CALL dr_hook(
'READ_GRIDTYPE_GAUSS',0,zhook_handle)
106 hprogram,
'LAPO',zlapo, kresp,hdir=hdir)
108 hprogram,
'LOPO',zlopo,kresp,hdir=hdir)
110 hprogram,
'CODIL',zcodil,kresp,hdir=hdir)
118 hprogram,
'NLATI',inlati,kresp,hdir=hdir)
119 ALLOCATE(inlopa(inlati))
122 hprogram,
'NLOPA',inlopa(:),kresp,hdir=hdir)
125 hprogram,
'NLOPA',inlopa(:),kresp,hdir=
'-')
128 hprogram,
'LATGAUSS',zlat(:),kresp,hdir=hdir)
130 hprogram,
'LONGAUSS',zlon(:),kresp,hdir=hdir)
132 hprogram,
'LAT_G_XY',zlat_xy(:),kresp,hdir=hdir)
134 hprogram,
'LON_G_XY',zlon_xy(:),kresp,hdir=hdir)
136 hprogram,
'MESHGAUSS',zmesh_size(:),kresp,hdir=hdir)
138 hprogram,
'LONINF',zloninf(:),kresp,hdir=hdir)
140 hprogram,
'LATINF',zlatinf(:),kresp,hdir=hdir)
142 hprogram,
'LONSUP',zlonsup(:),kresp,hdir=hdir)
144 hprogram,
'LATSUP',zlatsup(:),kresp,hdir=hdir)
152 klu,zlat,zlon,zlat_xy,zlon_xy,zmesh_size, &
153 zloninf,zlatinf,zlonsup,zlatsup )
158 IF (
SIZE(pgrid_par) /=
SIZE(zgrid_par))
THEN
160 WRITE(iluout,*)
'size of PGRID_PAR =',
SIZE(pgrid_par)
161 WRITE(iluout,*)
'size of ZGRID_PAR =',
SIZE(zgrid_par)
162 CALL
abor1_sfx(
'READ_GRIDTYPE_GAUSS: SIZE OF PGRID_PAR IS NOT CORRECT')
165 pgrid_par = zgrid_par
167 kgrid_par =
SIZE(zgrid_par)
170 DEALLOCATE(zgrid_par)
171 IF (lhook) CALL dr_hook(
'READ_GRIDTYPE_GAUSS',1,zhook_handle)
subroutine put_gridtype_gauss(PGRID_PAR, KNLATI, PLAPO, PLOPO, PCODIL, KNLOPA, KL, PLAT, PLON, PLAT_XY, PLON_XY, PMESH_SIZE, PLONINF, PLATINF, PLONSUP, PLATSUP)
subroutine read_gridtype_gauss(HPROGRAM, KGRID_PAR, KLU, OREAD, KSIZE, PGRID_PAR, KRESP, HDIR)
subroutine abor1_sfx(YTEXT)
subroutine get_luout(HPROGRAM, KLUOUT)