47 character(len=255) :: filename
48 character(len=255) :: wfname1(maxwf)
49 integer :: metdata_format
50 integer :: gfileid, status, gribid, centrenum
52 if ( maxwf.le.0 )
then
53 print*,
'No wind file available'
54 metdata_format = unknown_metdata
58 filename = path(3)(1:length(3)) // trim(wfname(1))
62 call grib_open_file(gfileid, trim(filename),
'r', status)
63 if (status.ne. grib_success)
then
64 print *,
'Unable to open: ', filename
65 metdata_format = unknown_metdata
72 call grib_new_from_file(gfileid, gribid, status)
75 call grib_get(gribid,
'centre', centrenum)
77 if (centrenum .EQ. 7)
then
78 metdata_format = gfs_metdata
79 elseif (centrenum == 98)
then
80 metdata_format = ecmwf_metdata
82 metdata_format = unknown_metdata
85 call grib_close_file(gfileid)