######################
## CURIE   TGCC/CEA ##
######################
#MSUB -r create_MM             # Nom du job               
#MSUB -eo
#MSUB -n 1              # Reservation du processus
#MSUB -T 86400          # Limite de temps elapsed du job
#MSUB -q xlarge
#MSUB -Q normal
#MSUB -A gen2211
set +x

#### Script permettant d'extraire depuis les fichiers journaliers de inca pack en frequence 1Y les moyennes mensuelles

EXPER=ESM1989.4
EXPERIMENT=historical
RUNDIR=$SCRATCHDIR/POST_climato
mkdir -p $RUNDIR
cd $RUNDIR

STORESIMU=$STOREDIR/IGCM_OUT/IPSLCM5CHT/PROD/${EXPERIMENT}/${EXPER}

mkdir -p ${STORESIMU}/CHM/Analyse/MM
#
for YEAR in   1999 2000
do

    if [[ ! -s ${EXPER}_${YEAR}0101_${YEAR}1231_1M_inca_avgr.nc ]]; then 
	ccc_hsm get ${STORESIMU}/CHM/Output/MO/${EXPER}_${YEAR}0101_${YEAR}1231_1M_inca_avgr.nc
	cp ${STORESIMU}/CHM/Output/MO/${EXPER}_${YEAR}0101_${YEAR}1231_1M_inca_avgr.nc .
    fi

    ncwa -O -a time_counter -d time_counter,0,30    ${EXPER}_${YEAR}0101_${YEAR}1231_1M_inca_avgr.nc ${EXPER}_${YEAR}_01_MM_inca_avgr.nc 
    ncwa -O -a time_counter -d time_counter,31,58   ${EXPER}_${YEAR}0101_${YEAR}1231_1M_inca_avgr.nc ${EXPER}_${YEAR}_02_MM_inca_avgr.nc 
    ncwa -O -a time_counter -d time_counter,59,89   ${EXPER}_${YEAR}0101_${YEAR}1231_1M_inca_avgr.nc ${EXPER}_${YEAR}_03_MM_inca_avgr.nc 
    ncwa -O -a time_counter -d time_counter,90,119  ${EXPER}_${YEAR}0101_${YEAR}1231_1M_inca_avgr.nc ${EXPER}_${YEAR}_04_MM_inca_avgr.nc 
    ncwa -O -a time_counter -d time_counter,120,150 ${EXPER}_${YEAR}0101_${YEAR}1231_1M_inca_avgr.nc ${EXPER}_${YEAR}_05_MM_inca_avgr.nc 
    ncwa -O -a time_counter -d time_counter,151,180 ${EXPER}_${YEAR}0101_${YEAR}1231_1M_inca_avgr.nc ${EXPER}_${YEAR}_06_MM_inca_avgr.nc 
    ncwa -O -a time_counter -d time_counter,181,211 ${EXPER}_${YEAR}0101_${YEAR}1231_1M_inca_avgr.nc ${EXPER}_${YEAR}_07_MM_inca_avgr.nc 
    ncwa -O -a time_counter -d time_counter,212,242 ${EXPER}_${YEAR}0101_${YEAR}1231_1M_inca_avgr.nc ${EXPER}_${YEAR}_08_MM_inca_avgr.nc 
    ncwa -O -a time_counter -d time_counter,243,272 ${EXPER}_${YEAR}0101_${YEAR}1231_1M_inca_avgr.nc ${EXPER}_${YEAR}_09_MM_inca_avgr.nc 
    ncwa -O -a time_counter -d time_counter,273,303 ${EXPER}_${YEAR}0101_${YEAR}1231_1M_inca_avgr.nc ${EXPER}_${YEAR}_10_MM_inca_avgr.nc 
    ncwa -O -a time_counter -d time_counter,304,333 ${EXPER}_${YEAR}0101_${YEAR}1231_1M_inca_avgr.nc ${EXPER}_${YEAR}_11_MM_inca_avgr.nc 
    ncwa -O -a time_counter -d time_counter,334,364 ${EXPER}_${YEAR}0101_${YEAR}1231_1M_inca_avgr.nc ${EXPER}_${YEAR}_12_MM_inca_avgr.nc 
    
    ncecat -F -O  ${EXPER}_${YEAR}_01_MM_inca_avgr.nc   ${EXPER}_${YEAR}_02_MM_inca_avgr.nc  ${EXPER}_${YEAR}_03_MM_inca_avgr.nc  ${EXPER}_${YEAR}_04_MM_inca_avgr.nc  ${EXPER}_${YEAR}_05_MM_inca_avgr.nc  ${EXPER}_${YEAR}_06_MM_inca_avgr.nc  ${EXPER}_${YEAR}_07_MM_inca_avgr.nc  ${EXPER}_${YEAR}_08_MM_inca_avgr.nc  ${EXPER}_${YEAR}_09_MM_inca_avgr.nc  ${EXPER}_${YEAR}_10_MM_inca_avgr.nc ${EXPER}_${YEAR}_11_MM_inca_avgr.nc  ${EXPER}_${YEAR}_12_MM_inca_avgr.nc  ${EXPER}_${YEAR}_allMM_inca_avgr.nc
#ncwa -O -a time_counter ${EXPER}_${YEAR}_allMM_inca_avgr.nc ${EXPER}_${YEAR}_allMM_inca_avgr.nc   #enleve dim degenerate

    ncra -F -O ${EXPER}_${YEAR}_allMM_inca_avgr.nc ${EXPER}_${YEAR}_AM_inca_avgr.nc

    ncrename -d record,time_counter ${EXPER}_${YEAR}_allMM_inca_avgr.nc
    mv ${EXPER}_${YEAR}_allMM_inca_avgr.nc ${EXPER}_${YEAR}_AM_inca_avgr.nc ${EXPER}_${YEAR}_*_MM_inca_avgr.nc  ${STORESIMU}/CHM/Analyse/MM/.
            
done
exit