#################################################################
#  ##### #     #  ###  #####
#  #     ##   ## #   # #
#  ##### # # # # #   # #####
#      # #  #  # #   #     #     BUFR to GRIB converted
#      # #     # #   #     #
#  ##### #     #  ###  #####
##################################################################

# SMOS NRT BUFR product convertion in GRIB
# Patricia de Rosnay  (Patricia.Rosnay@ecmwf.int)
# October 2011
##################################################################

# This software converts SMOS NRT BUFR files to GRIB files, either on a regular or a Gaussian Reduced grid to be choosed by the user.
# The software relies on:
# 1- A Metview programme that converts BUFR to GEOPOINT and then GEOPOINT to GRIB on a regular Latitude-Longitude Grid
# 2- A scin command to interpolate the GRIB file on a regular Gaussian grid if requested. 

# Information on Gaussian reduced grid: http://www.ecmwf.int/products/data/technical/gaussian/index.html

# The control script is smos_bufr2grib.sh
# From the script the user controls:
#	 	the input file, incidence angle binning angle, output grid and resolution
# 		alias free FOV, or entire FOV
# The script launches the metview macro  smos_bufr2gribLL.met
#  In the default version the wetland filters out all pixels with a water fraction larger than 5%, this can be changed
#  The user can easily modify all the filters in smos_bufr2gribLL.met. Comments in the metview programme make it easy to understand
#  and to modifyi, so please do not hesitate.

# To convert the data over Ocean only or on both ocean and land, 
# change the the water fraction custom values, line 82 of smos_bufr2gribLL.met.


# To run the converter:
sh smos_bufr2grib.sh

###################################################################
# Requested tools:
# 1- Linux or Unix computer
# 2- Metview: http://www.ecmwf.int/publications/manuals/metview/index.html
# 3- SCIN: ECMWF interpolation software. Not yet available outside ECMWF,
#          so skip this part if not at ECMWF by puting 
grid_resol="regular"
# in smos_bufr2grib.sh


# Note on the SMOS NRT BUFR file content
# radiometric accuracy 012167
# snapshot accuracy: 012166 (Corbella eq)
# polarization: 002099 (0h, 1V)
# snapshot overallquality: 033028  (must be =1 but is always=3)
# SMOS info flag: 025174
# Faraday angle 025083
# TB 12080
# For more information see the SMOS BUFR specification document
# (SMOS_NRT_BUFR_ECMWF_v1.9.pdf)