#!/bin/bash
set -e
#author:heyc
# get the bias between obs(initial state) and simulation
# Revised:
# 2008.12.31
# 2009.01.02
# 2011.05.31
# 2012.08.20
# 2012.10.01

#~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
# -name of experiment
ExpID=H04S5

# varname
export VAR=sss
#export VAR=sst

# value range, numeric or don't set
case $VAR in
	"sss")
		export minn=-1.0
		export maxn=1.0
		export cnsp=0.1
		;;
	"sst")
		export minn=-5.0
		export maxn=5.0
		export cnsp=0.5
		;;
	*)
		echo 'get "nicecon"...'
esac


ROOT="~"
PATHIN=${ROOT}/output/$ExpID/gen2reg
PATHOUT=~/output/$ExpID/sim-clm
PATHOUT=`echo ${PATHOUT} | tr "." "_"`
mkdir -p ${PATHOUT}/
#~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
echo "ExpID:	$ExpID"
echo "VAR:	$VAR"
echo "PREFIX:	$PREFIX"
echo "ista:	$ista"
echo "iend:	$iend"
echo "imon:	$imon"
echo "minn:	$minn"
echo "maxn:	$maxn"
echo "cnsp:	$cnsp"
#~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
ln -sf ~/tools/data/grid_mnp2.nc grid.nc
ln -sf ~/output/${ExpID}/ave2d/${ExpID}_ave2d_1948-2007.nc in.nc
ln -sf ~/tools/grids/MNP2/Levitus/Relax/PHC3.0/20120307/${VAR}_clim_phc3.0.nc clim.nc

# -plot diff of simulated vs inital value
ncl micom_sim-clm_phc30_annual.ncl

# postproc figures and data
#epstopdf figure.eps && rm -f figure.eps  &&  putn "figure.*"
epstopdf figure.eps && rm -f figure.eps
mv figure.pdf H04S5_ave2d_${VAR}_1948-2007.01_sim-phc30_annual.pdf
FIGURE=${ExpID}_2d_ave_sim-clm

# files history logs
SNAME=`echo $0 | cut -d '/' -f2`
echo -e "`date +%Y%m%d%H%M`:" >> ${PATHOUT}/history.log
echo -e "\t`pwd`/${SNAME}\t--->\t${FIGURE}.*" >> ${PATHOUT}/history.log