grib2check.F90

Go to the documentation of this file.

subroutine grib2check(igrib,isec1,conversion_factor)
#if defined WITH_CTBTO_PATCHES
  use grib_api
  use par_mod
  use com_mod

  implicit none
  integer :: gribver,parcat,parnum,typsurf,valsurf,discipl,parid
  integer :: iret,igrib
  integer :: isec1(*)
  real    :: conversion_factor

  !HSO  grib api error messages
  character(len=24) :: griberrormsg = 'Error reading grib file'
  character(len=20) :: gribfunction = 'readwind'

  !print*,'GRiB Edition 2'
  !read the grib2 identifiers
  call grib_get_int(igrib,'discipline',discipl,iret)
  call grib_check(iret,gribfunction,griberrormsg)
  call grib_get_int(igrib,'parameterCategory',parcat,iret)
  call grib_check(iret,gribfunction,griberrormsg)
  call grib_get_int(igrib,'parameterNumber',parnum,iret)
  call grib_check(iret,gribfunction,griberrormsg)
  call grib_get_int(igrib,'typeOfFirstFixedSurface',typsurf,iret)
  call grib_check(iret,gribfunction,griberrormsg)
  call grib_get_int(igrib,'level',valsurf,iret)
  call grib_check(iret,gribfunction,griberrormsg)
  call grib_get_int(igrib,'paramId',parid,iret)
  call grib_check(iret,gribfunction,griberrormsg)

  !print*,discipl,parId,parCat,parNum,typSurf,valSurf

  !convert to grib1 identifiers
  isec1(6)=-1
  isec1(7)=-1
  isec1(8)=-1
  isec1(8)=valsurf     ! level


  conversion_factor=1.
  isec1(6)=parid

  if(parid .eq. 3063 ) then ! acpcp [kg/m2] convert into m
    conversion_factor=1000.
  endif
  if(parid .eq. 0 ) then !  snow depth [kg/m2] convert into m
    if ((parcat.eq.1).and.(parnum.eq.11).and.(typsurf.eq.1)) then ! SD
      isec1(6)=141         ! indicatorOfParameter    
      conversion_factor=1000.
    else
     print*,'***ERROR: undefined GRiB2 message found!',discipl, &
         parcat,parnum,typsurf
    
    endif
  endif

  

#endif
return
end subroutine grib2check