#!/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