一、安装wgrib2

前一篇文章已经详细描述。

二、fortran文件 test.f90

Program read_grib2

use wgrib2api

Implicit None

real, allocatable :: grid(:,:), lat(:,:), lon(:,:)

character (len=300) :: metadata, grid_info

character (len=200) :: file, inv, var, str_date

integer :: date

integer :: nx, ny, iret

Character (Len=99) :: invline

Character (Len=400) :: gridline

file = 'multi_reanal.glo_30m_ext.wind.198001.grb2'

inv = 'multi_reanal.glo_30m_ext.wind.198001.inv'

! make inv file, save in memory file #0

iret = grb2_mk_inv(file, inv)

If (iret/=0) Stop 1

! search using variable and regular date YYYYMMDDHH

var = ':UGRD:'

str_date = ':3 hour fcst:'

! iret = grb2_inq(file, inv, var, strdate, data2=grid)

iret = grb2_inq(file, inv, var, ':3 hour fcst:', data2=grid, lat=lat, lon=lon, desc=invline, grid_desc=gridline)

If (iret/=1) Stop 4

Write (*, *) 'SP(1,1) =', grid(1, 1), ' lat/lon=', lat(1, 1), lon(1, 1)

Write (*, *) 'inventory=', trim(invline)

Write (*, *) 'grid=', trim(gridline)

End Program read_grib2

三、编译文件compile.intel

#!/bin/sh

# compile for intel ifort

set -x

prog=$1

LIB="-L/$HOME/Downloads/grib2/lib -lwgrib2 -ldl "

LIB+="-L/opt/jasper-1.900.1_intel/lib -lpng -ljasper "

LIB+="-L/opt/zlib-1.2.11_intel/lib -lz "

LIB+="-L/opt/hdf5-1.10.2_intel/lib -lhdf5_hl -lhdf5 -ldl "

LIB+="-L/opt/netcdf-c-4.6.1_intel/lib -lnetcdf "

LIB+="-L/opt/netcdf-f-4.4.4_intel/lib -lnetcdff "

MOD="-I/$HOME/Downloads/grib2/lib "

MOD+="-I/opt/jasper-1.900.1_intel/include "

MOD+="-I/opt/hdf5-1.10.2_intel/include "

MOD+="-I/opt/netcdf-c-4.6.1_intel/include "

MOD+="-I/opt/netcdf-f-4.4.4_intel/include "

ifort -no-wrap-margin -o $prog $prog.f90 ${MOD} ${LIB}

exit

四、编译

./compile_intel test

五、执行

./test

结果为

SP(1,1) = 9.9990003E+20 lat/lon= -90.00000 0.0000000E+00

inventory=D=19800101000000:UGRD:surface:3 hour fcst:

grid=grid_template=0:winds(N/S):

lat-lon grid:(720 x 361) units 1e-06 input WE:NS output WE:SN res 48

lat 90.000000 to -90.000000 by 0.500000

lon 0.000000 to 359.500000 by 0.500000 #points=259920

fortran读取grib2文件