#!/bin/bash
gridfile='/nird/home/pgchiu/scratch/norcpm_verify/grid_tnx1v4_20170622.nc'
griddes=griddes.txt
## for grid file of BLOM/MICOM
xsize=$(ncdump -h $gridfile | grep 'x =' | sed 's/[ \t]*x = \([0-9]*\) .*/\1/')
ysize=$(ncdump -h $gridfile | grep 'y =' | sed 's/[ \t]*y = \([0-9]*\) .*/\1/')
gridsize=$(($xsize * $ysize))
ncdump -v plon $gridfile | sed -e'1,/data:/d' -e's/^.*=//' -e's/;//' -e's/}//'> plon.txt
ncdump -v plat $gridfile | sed -e'1,/data:/d' -e's/^.*=//' -e's/;//' -e's/}//'> plat.txt
ncdump -v pclon $gridfile | sed -e'1,/data:/d' -e's/^.*=//' -e's/;//' -e's/}//'> pclon.txt
ncdump -v pclat $gridfile | sed -e'1,/data:/d' -e's/^.*=//' -e's/;//' -e's/}//'> pclat.txt

echo 'gridtype = unstructured'   > $griddes
echo "gridsize = $gridsize"  >> $griddes
echo "xsize    = $xsize " >> $griddes
echo "ysize    = $ysize " >> $griddes
echo "xvals    = " >> $griddes
cat plon.txt >> $griddes
echo "yvals    = " >> $griddes
cat plat.txt >> $griddes

echo 'nvertex  = 4' >> $griddes
echo 'xbounds  = ' >> $griddes
cat pclon.txt >> $griddes
echo 'ybounds  = ' >> $griddes
cat pclat.txt >> $griddes