#!/bin/sh
#to convert ostia from netcdf to grib, ref
#https://www.cpc.ncep.noaa.gov/products/wesley/wgrib2/netcdf_2_grib_ostia.html
#

#step1: ncdump -h xxnc
#step2: making a templaet file
#wgrib2 ~wd51we/grib2/examples/gep19.aec -d 1  -set_pdt +0 -set_var HGT  -new_grid_winds earth \
#   -new_grid latlon 180.025000:7200:0.050000 -89.975000:3600:0.05 ostia.template.tmp

# setting grid == 0 to save space
#wgrib2 ostia.template.tmp -rpn 0 -set_lev surface -set_ftime anl -grib_out ostia.template

#step3: wgrib2 convert script
#wgrib2=/export/cpc-lw7-webisuzaki/wd51we/bin/wgrib2.v2.0.9rc3.n4
wgrib2=/ncrc/home2/Wesley.Ebisuzaki//bin/wgrib2
wgrib2=wgrib2

indir=/lustre/f2/scratch/$USER/ostiadump
indir=./
# for date in `days 20240312..20230312`
for date in 20240101
do
  yr=`echo $date |cut -c1-4`
  mon=`echo $date |cut -c5-6`
  # file=$indir/$yr/$mon/${date}120000-UKMO-L4_GHRSST-SSTfnd-OSTIA-GLOB-v02.0-fv02.0.nc
  # new_file=$yr/${date}120000-UKMO-L4_GHRSST-SSTfnd-OSTIA-GLOB-v02.0-fv02.0.grb
  file=${date}120000-UKMO-L4_GHRSST-SSTfnd-OSTIA-GLOB-v02.0-fv02.0.nc
  new_file=${date}120000-UKMO-L4_GHRSST-SSTfnd-OSTIA-GLOB-v02.0-fv02.0.grb
  $wgrib2 ostia.template \
       -import_netcdf $file  \
       "sea_ice_fraction" "0:1:0:3600:0:7200" \
       -set_var ICEC -set_scaling -2 0 -set center 254 -set_grib_type c1 -set_date ${date}12 \
       -grib_out $new_file \
       -import_netcdf $file  \
       analysed_sst "0:1:0:3600:0:7200" \
       -set_var TMP -set center 74 -set center 74 -set_date ${date}12 \
       -set_scaling -2 0 -set_grib_type c3 -grib_out $new_file \
       -import_netcdf $file  \
       mask "0:1:0:3600:0:7200" \
       -set_var MASK -set center 74 -set_date ${date}12 \
       -set_scaling -2 0 -set_grib_type c1 -grib_out $new_file \
       -import_netcdf $file \
       analysis_error "0:1:0:3600:0:7200" \
       -set_var TMP -set table_4.3 7 -set center 74 -set_date ${date}12 \
       -set_scaling -2 0 -set_grib_type c3 -grib_out $new_file
done