#!/bin/ksh

#**************************************************************
# Author: Marie-Alice Foujols
# Contact: Marie-Alice.Foujols__at__ipsl.jussieu.fr
# $Revision::                                          $ Revision of last commit
# $Author::                                            $ Author of last commit
# $Date::                                              $ Date of last commit
# IPSL (2009)
#  This software is governed by the CeCILL licence see libIGCM/libIGCM_CeCILL.LIC
#
#**************************************************************

#------------------------------------------------------------
# ncks CPL files for 1D frequency, ncra for 1M averaged files
#------------------------------------------------------------

function IGCM_Patch_20140801_ncks_cdo {
    IGCM_debug_PushStack "IGCM_Patch_ncks_ncra" $@

    PREFIX=TEMPO
    if [ -f FIRST_TIME_CPL ] ; then
       echo FIRST_TIME_CPL exists
    else
       # create FIRST_TIME_CPL file
       date > FIRST_TIME_CPL

       # Get date from directory name
       DateDeb=$( basename $(pwd) | awk -F _ '{ print $2}' )
       DateDeb=$( IGCM_date_ConvertFormatToHuman ${DateDeb} )
       
       rm -f ${PREFIX}_*_cpl_*.nc 

       list_file_atm=$( find . \( -name "S*lmdz.x*.nc" -o -name "C*lmdz.x*.nc" -o -name "S*LMDZ*.nc" -o -name "C*LMDZ*.nc" \) )

       for file in ${list_file_atm} 
       do
           IGCM_sys_ncks -A ${file} ${PREFIX}_1D_cpl_atm.nc  > /dev/null 2>&1
       done
       rm -f S*lmdz.x*.nc C*lmdz.x*.nc S*LMDZ*.nc C*LMDZ*.nc

       ## Correctly defines the time axis for CDO
       # DateDeb is first date of the file, not start date of run
       # (according to time var in OASIS files)
       #
       ncatted \
	   -a axis,time,c,c,"T"                                  \
	   -a long_name,time,c,c,"Time axis"                     \
	   -a title,time,c,c,"Time"                              \
	   -a calendar,time,c,c,"noleap"                         \
	   -a units,time,c,c,"seconds since ${DateDeb} 00:00:00" \
	   -a time_origin,time,c,c,"${DateDeb} 00:00:00"         \
	   ${PREFIX}_1D_cpl_atm.nc 

       for file in O*oceanx*.nc
       do
           IGCM_sys_ncks -A ${file} ${PREFIX}_1D_cpl_oce.nc   > /dev/null 2>&1
       done
       rm -f O*oceanx*.nc

       ncatted \
	   -a axis,time,c,c,"T"                                  \
	   -a long_name,time,c,c,"Time axis"                     \
	   -a title,time,c,c,"Time"                              \
	   -a calendar,time,c,c,"noleap"                         \
	   -a units,time,c,c,"seconds since ${DateDeb} 00:00:00" \
	   -a time_origin,time,c,c,"${DateDeb} 00:00:00"         \
	   ${PREFIX}_1D_cpl_oce.nc 
    fi

    # Monthly average
    IGCM_sys_cdo monavg ${PREFIX}_1D_cpl_atm.nc  ${PREFIX}_1M_cpl_atm.nc
    IGCM_sys_cdo monavg ${PREFIX}_1D_cpl_oce.nc  ${PREFIX}_1M_cpl_oce.nc
    
    IGCM_debug_PopStack "IGCM_Patch_ncks_ncra"
}