#!/bin/bash
#
#
# Check with cdo info that the mean fields of LWnet are negative.
#
#
EXP="ETHZ_Avg IPSL_Avg IPSL_Alt WFDE5_CRU_GPCC"

v="LWnet"
for e in ${EXP} ; do
    if [ ! -e ${e}/${v}_${e}_8914_orig.nc ] ; then
	mv ${e}/${v}_${e}_8914.nc ${e}/${v}_${e}_8914_orig.nc
    else
	/bin/rm -rf ${e}/${v}_${e}_8914.nc
    fi
    ##    cdo mulc,-1.0 ${e}/${v}_${e}_8914_orig.nc ${e}/${v}_${e}_8914.nc
    /bin/rm -f ${e}/tmp.nc
    ncap2 -s "LWnet=LWnet*-1.0" ${e}/${v}_${e}_8914_orig.nc ${e}/tmp.nc
    ncatted -h -a _FillValue,${v},o,f,-1.e+33 ${e}/tmp.nc ${e}/${v}_${e}_8914.nc
    /bin/rm -f ${e}/tmp.nc
done

EXP="ETHZ_Avg IPSL_Avg IPSL_Alt"
for e in ${EXP} ; do
    echo "=====>> $e $v <<====="
    ncdump -h ${e}/${v}_${e}_8914.nc | grep UNLIM
    /bin/rm -f ../MeteoCat/${e}/TS_DA/MeteoCat_${e}_19890101_20140101_1D_${v}.nc
    cp ${e}/${v}_${e}_8914.nc ../MeteoCat/${e}/TS_DA/MeteoCat_${e}_19890101_20140101_1D_${v}.nc
done

EXP="WFDE5_CRU_GPCC"
for e in ${EXP} ; do
    echo "=====>> $e $v <<====="
    ncdump -h ${e}/${v}_${e}_8914.nc | grep UNLIM
    /bin/rm -f ../MeteoCat/${e}/TS_DA/MeteoCat_${e}_19890101_20140101_1D_${v}.nc
    ncrename -v lon,LON -v lat,LAT ${e}/${v}_${e}_8914.nc \
	     ../MeteoCat/${e}/TS_DA/MeteoCat_${e}_19890101_20140101_1D_${v}.nc
done