7 hfiletype,hinterp_type,kni)
49 USE modd_grid_gauss, ONLY : xila1, xilo1, xila2, xilo2, ninla, ninlo, nilen, lrotpole, xcoef, xlap, xlop
51 USE yomhook
,ONLY : lhook, dr_hook
52 USE parkind1
,ONLY : jprb
61 CHARACTER(LEN=6),
INTENT(IN) :: hfiletype
62 CHARACTER(LEN=6),
INTENT(OUT) :: hinterp_type
63 INTEGER,
INTENT(OUT) :: kni
68 CHARACTER(LEN=12) :: yrecfm
73 REAL,
DIMENSION(:),
ALLOCATABLE :: zw
80 INTEGER,
DIMENSION(:),
ALLOCATABLE :: inlopa
82 REAL(KIND=JPRB) :: zhook_handle
85 IF (lhook) CALL dr_hook(
'PREP_GRID_GAUSS',0,zhook_handle)
93 hfiletype,yrecfm,zlapo,iresp)
96 hfiletype,yrecfm,zlopo,iresp)
99 hfiletype,yrecfm,zcodil,iresp)
108 hfiletype,yrecfm,inlati,iresp)
110 IF (
ALLOCATED(inlopa))
DEALLOCATE(inlopa)
111 ALLOCATE(inlopa(inlati))
112 IF (
ALLOCATED(ninlo))
DEALLOCATE(ninlo)
113 ALLOCATE(ninlo(inlati))
116 hfiletype,yrecfm,inlopa,iresp,hdir=
'A')
125 inlati2 = nint(
REAL(inlati)/2.0)
135 IF (zlapo>89.99 .AND. abs(zlopo)<0.00001)
THEN
141 xila1=90.0*(1.0-0.5/
REAL(inlati))
143 xila2=-90.0*(1.0-0.5/
REAL(inlati))
144 xilo2=360.0*(
REAL(inlopa(inlati2))-1.0)/
REAL(inlopa(inlati2))
146 hinterp_type =
'HORIBL'
148 IF (lhook) CALL dr_hook(
'PREP_GRID_GAUSS',1,zhook_handle)
subroutine prep_grid_gauss(HFILETYPE, HINTERP_TYPE, KNI)