#!/bin/bash

## Just run code in recipes, use it carefully
## example: use cdo to make cache file. Since cdo is quicker than ncl.
#;DIAG_NORCPM; RUNTHESECODES: echo '    no codes run here.'

#!/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 = curvilinear'   > $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