source: CONFIG/LMDZORINCA/branches/LMDZORINCA_BCOV/DRIVER/stomate.driver @ 1836

#!/bin/ksh

#D- Driver du script pour ORCHIDEE

#-----------------------------------------------------------------
function SBG_Initialize
{
    IGCM_debug_PushStack "SBG_Initialize"

    RESOL_SBG=ALL
    
    IGCM_debug_PopStack "SBG_Initialize"
}

#-----------------------------------------------------------------
function SBG_Update
{
    IGCM_debug_PushStack "SBG_Update"

    typeset STOMATE_WRITE_STEP

    case ${config_SBG_WriteFrequency} in
        *Y|*y) 
            WriteInYears=$( echo ${config_SBG_WriteFrequency} | awk -F '[yY]' '{print $1}' )
            PeriodLengthInYears=$( echo ${config_UserChoices_PeriodLength} | awk -F '[yY]' '{print $1}' )
            (( STOMATE_WRITE_STEP = PeriodLengthInDays * WriteInYears / PeriodLengthInYears )) ;;
        1M|1m) 
            case ${config_UserChoices_PeriodLength} in
            *Y|*y)
                STOMATE_WRITE_STEP=-1.
                ;;
            *M|*m)
                STOMATE_WRITE_STEP=-1.
                ;;
            *)
                (( STOMATE_WRITE_STEP = $( IGCM_date_DaysInMonth $year $month ) )) ;;
            esac
            ;;
        *M|*m) 
            WriteInMonths=$( echo ${config_SBG_WriteFrequency} | awk -F '[mM]' '{print $1}' )
            case ${config_UserChoices_PeriodLength} in
            *Y|*y)
                PeriodLengthInYears=$( echo ${config_UserChoices_PeriodLength} | awk -F '[yY]' '{print $1}' )
                (( STOMATE_WRITE_STEP = 30 ))
                ;;
            *M|*m)
                PeriodLengthInMonths=$( echo ${config_UserChoices_PeriodLength} | awk -F '[mM]' '{print $1}' )
                (( STOMATE_WRITE_STEP = PeriodLengthInDays * WriteInMonths / PeriodLengthInMonths ))
                ;;
            *)
                (( STOMATE_WRITE_STEP = $( IGCM_date_DaysInMonth $year $month ) )) ;;
            esac
            ;;
        5D|5d) 
            (( STOMATE_WRITE_STEP = 5 )) ;;
        1D|1d) 
            (( STOMATE_WRITE_STEP = 1 )) ;;
        *s)
            WriteInSeconds=$( echo ${config_SBG_WriteFrequency} | awk -F '[s]' '{print $1}' )
            (( STOMATE_WRITE_STEP = 1 )) ;;
        *) 
            IGCM_debug_Exit "SBG_Update " ${config_SBG_WriteFrequency} " invalid WriteFrequency : choose in 1Y, 1M, 5D, 1D." 
            IGCM_debug_Verif_Exit ;;
    esac

    ORCHIDEE_sed STOMATE_OK_STOMATE y

    ORCHIDEE_sed STOMATE_HIST_DT ${STOMATE_WRITE_STEP}
    ORCHIDEE_sed STOMATE_HISTLEVEL ${stomate_UserChoices_stomate_LEVEL}

    if ( [ ${CumulPeriod} -eq 1 ] && [ "${config_SBG_Restart}" = "n" ] ) ; then
        echo "STOMATE : without restart"
        ORCHIDEE_sed STOMATE_RESTART_FILEIN NONE
    else
        ORCHIDEE_sed STOMATE_RESTART_FILEIN stomate_rest_in.nc
    fi

    IGCM_debug_PopStack "SBG_Update"
}

#-----------------------------------------------------------------
function SBG_Finalize
{
#set -vx
    IGCM_debug_PushStack "SBG_Finalize"

    if [ $( IGCM_date_DaysBetweenGregorianDate ${PeriodDateEnd} ${year}1230 ) -ge 0 ] ; then
        IGCM_sys_Put_Out stomate_Cforcing.nc ${R_OUT_SBG_R}/${config_UserChoices_JobName}_${PeriodDateEnd}_stomate_Cforcing.nc
        IGCM_sys_Put_Out stomate_forcing.nc  ${R_OUT_SBG_R}/${config_UserChoices_JobName}_${PeriodDateEnd}_stomate_forcing.nc
#       rm -f stomate_Cforcing.nc
#       rm -f stomate_forcing.nc
    fi

    IGCM_debug_PopStack "SBG_Finalize"
}