#*********************************************************
#to add lon/lat dimension to the RegCM4 output file
#by Thanh NGO-DUC (ngoducthanh@vnu.edu.vn)
#at VNU Hanoi University of Science
#on Aug 2012
#system requirements: NCO should be already installed
#*********************************************************
module load cdo
module load nco
FILEIN=$1
FILEOUT=$2
cp $FILEIN $FILEOUT
# change the lat, lon dimesion names
ncrename -O -d y,lat -d x,lon $FILEOUT
# extract the 2D lon,lat variables
ncks -O -v xlon,xlat $FILEOUT tmp_lonlat.nc
# change the name of the 2D lat, lon variables
#ncrename -v lon,lon_tmp $FILEOUT
#ncrename -v lat,lat_tmp $FILEOUT
# change the 2D lon, lat variables to 1D lon, lat variables
ncwa -O -a lon tmp_lonlat.nc tmp_lat.nc
ncwa -O -a lat tmp_lonlat.nc tmp_lon.nc
ncrename -v xlat,lat tmp_lat.nc
ncrename -v xlon,lon tmp_lon.nc
# add the new 1D lon, lat variables to $FILEOUT
ncks -A -v lat tmp_lat.nc $FILEOUT
ncks -A -v lon tmp_lon.nc $FILEOUT
# remove the previous 2D lon,lat variables
ncks -O -x -v xlat,xlon,iy,jx $FILEOUT $FILEOUT
# remove the intermediate files
rm tmp_lat.nc tmp_lon.nc tmp_lonlat.nc
#END