forked from DOI-USGS/volcano-ash3d-metreader
-
Notifications
You must be signed in to change notification settings - Fork 4
/
gen_GRIB_index.f90
33 lines (26 loc) · 886 Bytes
/
gen_GRIB_index.f90
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
program gen_grib_index
use MetReader
use eccodes
implicit none
integer :: nargs
integer :: iostatus
character (len=100) :: arg
character(len=130) :: grib_file
INTERFACE
subroutine MR_Set_Gen_Index_GRIB(grib_file)
character(len=130),intent(in) :: grib_file
end subroutine MR_Set_Gen_Index_GRIB
END INTERFACE
nargs = command_argument_count()
if (nargs.ne.1) then
!do io=1,MR_nio;if(VB(io).le.verbosity_error)then
! write(errlog(io),*)"MR ERROR: no grib file given"
write(6,*)"MR ERROR: no grib file given"
!endif;enddo
stop 1
else
call get_command_argument(number=1, value=arg, status=iostatus)
read(arg,*)grib_file
call MR_Set_Gen_Index_GRIB(grib_file)
endif
end program