Description
This example shows: how to read levels from TEMP BUFR messages.
Example BUFR file: PraticaTemp.bufr
Source code
! !Copyright 2005-2015 ECMWF. ! ! This software is licensed under the terms of the Apache Licence Version 2.0 !which can be obtained at http://www.apache.org/licenses/LICENSE-2.0. ! ! In applying this licence, ECMWF does not waive the privileges and immunities granted to it by ! virtue of its status as an intergovernmental organisation nor does it submit to any jurisdiction. ! ! ! FORTRAN 90 Implementation: bufr_read_temp ! ! Description: how to read levels from TEMP BUFR messages. ! ! Please note that TEMP reports can be encoded in various ways in BUFR. Therefore the code ! below might not work directly for other types of TEMP messages than the one used in the ! example. It is advised to use bufr_dump first to understand the structure of these messages. program bufr_read_temp use eccodes implicit none integer :: ifile integer :: iret,ierr integer :: ibufr integer :: i, count=0 integer(kind=4),dimension(:), allocatable :: timePeriod,extendedVerticalSoundingSignificance integer(kind=4) :: blockNumber,stationNumber,numberOfLevels character(len=30) :: str real(kind=8),dimension(:), allocatable :: pressure,airTemperature,dewpointTemperature real(kind=8),dimension(:), allocatable :: geopotentialHeight,latitudeDisplacement,longitudeDisplacement real(kind=8),dimension(:), allocatable :: windDirection,windSpeed character(len=128) :: keyName call codes_open_file(ifile,'PraticaTemp.bufr','r') !call codes_open_file(ifile,'temp.bufr','r') ! the first bufr message is loaded from file ! ibufr is the bufr id to be used in subsequent calls call codes_bufr_new_from_file(ifile,ibufr,iret) do while (iret/=CODES_END_OF_FILE) write(*,*) 'message: ',count call codes_set(ibufr,'unpack',1) call codes_get(ibufr,'timePeriod',timePeriod) call codes_get(ibufr,'pressure',pressure) call codes_get(ibufr,'extendedVerticalSoundingSignificance',extendedVerticalSoundingSignificance) call codes_get(ibufr,'geopotentialHeight',geopotentialHeight) call codes_get(ibufr,'latitudeDisplacement',latitudeDisplacement) call codes_get(ibufr,'longitudeDisplacement',longitudeDisplacement) call codes_get(ibufr,'airTemperature',airTemperature) call codes_get(ibufr,'dewpointTemperature',dewpointTemperature) call codes_get(ibufr,'windDirection',windDirection) call codes_get(ibufr,'windSpeed',windSpeed) call codes_get(ibufr,'blockNumber',blockNumber) call codes_get(ibufr,'stationNumber',stationNumber) print *,'station',blockNumber,stationNumber print *,'timePeriod pressure geopotentialHeight latitudeDisplacement & &longitudeDisplacement airTemperature windDirection windSpeed significance' do i=1,size(windSpeed) write(*,'(I5,6X,F9.1,2X,F9.2,10X,F8.2,14X,F8.2,16X,F8.2,6X,F8.2,4X,F8.2,4X,I0)') timePeriod(i),pressure(i),& &geopotentialHeight(i),latitudeDisplacement(i),& &longitudeDisplacement(i),airTemperature(i),windDirection(i),windSpeed(i),extendedVerticalSoundingSignificance(i) enddo ! free arrays deallocate(timePeriod) deallocate(pressure) deallocate(geopotentialHeight) deallocate(latitudeDisplacement) deallocate(longitudeDisplacement) deallocate(airTemperature) deallocate(dewpointTemperature) deallocate(windDirection) deallocate(windSpeed) deallocate(extendedVerticalSoundingSignificance) ! release the bufr message call codes_release(ibufr) ! load the next bufr message call codes_bufr_new_from_file(ifile,ibufr,iret) count=count+1 end do ! close file call codes_close_file(ifile) end program bufr_read_temp
# # Copyright 2005-2015 ECMWF. # # This software is licensed under the terms of the Apache Licence Version 2.0 # which can be obtained at http://www.apache.org/licenses/LICENSE-2.0. # # In applying this licence, ECMWF does not waive the privileges and immunities # granted to it by virtue of its status as an intergovernmental organisation # nor does it submit to any jurisdiction. # # Python implementation: bufr_read_temp # # Description: how to read levels from TEMP BUFR # messages. # # # Please note that TEMP reports can be encoded in various ways in BUFR. # Therefore the code below might not work directly for other types of TEMP # messages than the one used in the example. It is advised to use bufr_dump to # understand the structure of the messages. # import traceback import sys from eccodes import * INPUT = 'PraticaTemp.bufr' VERBOSE = 1 # verbose error reporting def example(): # open bufr file f = open(INPUT) cnt = 0 # loop for the messages in the file while 1: # get handle for message bufr = codes_bufr_new_from_file(f) if bufr is None: break print "message: %s" % cnt # we need to instruct ecCodes to expand all the descriptors # i.e. unpack the data section codes_set(bufr, 'unpack', 1) #get all the timePeriods timePeriod = codes_get_array(bufr, "timePeriod") pressure = codes_get_array(bufr, "pressure") extendedVerticalSoundingSignificance = codes_get_array(bufr, "extendedVerticalSoundingSignificance") geopotentialHeight = codes_get_array(bufr, "geopotentialHeight") latitudeDisplacement = codes_get_array(bufr, "latitudeDisplacement") longitudeDisplacement = codes_get_array(bufr, "longitudeDisplacement") airTemperature = codes_get_array(bufr, "airTemperature") dewpointTemperature = codes_get_array(bufr, "dewpointTemperature") windDirection = codes_get_array(bufr, "windDirection") windSpeed = codes_get_array(bufr, "windSpeed") blockNumber = codes_get(bufr, "blockNumber") stationNumber = codes_get(bufr, "stationNumber") print 'station %d%d' % (blockNumber,stationNumber) print 'timePeriod pressure geopotentialHeight latitudeDisplacement longitudeDisplacement airTemperature windDirection windSpeed significance' for i in range(0,len(windSpeed)-1): print timePeriod[i],pressure[i],geopotentialHeight[i],latitudeDisplacement[i],longitudeDisplacement[i],airTemperature[i],windDirection[i],windSpeed[i],extendedVerticalSoundingSignificance[i] cnt += 1 # delete handle codes_release(bufr) # close the file f.close() def main(): try: example() except CodesInternalError, err: if VERBOSE: traceback.print_exc(file=sys.stderr) else: print >>sys.stderr, err.msg return 1 if __name__ == "__main__": sys.exit(main())