4A1-GRIB API库235.pdf

GRIB API Fortran 90 - C - Python interfaces part 2 Dominique Lucas - Christian Weihrauch User Support Slide 1 com-grib 2013 - Fortran/C/Python APIs - II ©ecmwf Slide 1 Content GRIB API indexed access Encoding a loaded GRIB message C API Python API Slide 2 com-grib 2013 - Fortran/C/Python APIs - II ©ecmwf Slide 2 Input arguments Example – grib_get Output arguments ! Load all the GRIB messages contained in file.grib1 call grib_open_file(ifile, file.grib1,r) Loop on all the messages in a file. n=1 A new grib message is loaded call grib_new_from_file(ifile,igrib(n), iret) from file. igrib(n) is the grib id to LOOP: do while (iret /= GRIB_END_OF_FILE) be used in subsequent calls n=n+1; call grib_new_from_file(ifile,igrib(n), iret) end do LOOP ! Decode/encode data from the loaded message call grib_get( igrib(indx) , “dataDate”, date) Values is declared as call grib_get(igrib(indx), “typeOfLevel”, typeOfLevel) real, dimension(:), allocatable:: values call grib_get(igrib(indx), “level”, level) call grib_get_size(igrib(indx), “values”, nb_values); allocate(values(nb_values)) call grib_get(igrib(indx), “values”, values) print*, date, levelType, level, values(1), values(nb_values) ! Release do i=1,n call grib_release(igrib(n)) Slide 3 end do deallocate(values) call grib_close_file(ifile) com-grib 2013 - Fort