source: trunk/libIGCM/AA_metrics_LMDZ @ 1159

#-Q- curie ######################
#-Q- curie ## CURIE   TGCC/CEA ##
#-Q- curie ######################
#-Q- curie #MSUB -r METRICS_LMDZ       # Nom du job
#-Q- curie #MSUB -eo
#-Q- curie #MSUB -n 1              # Reservation du processus
#-Q- curie #MSUB -T 14400          # Limite de temps elapsed du job
#-Q- curie #MSUB -q ::default_node::
#-Q- curie #MSUB -Q normal
#-Q- curie #MSUB -A ::default_project::
#-Q- curie set +x
#-Q- ada #!/bin/ksh
#-Q- ada #######################
#-Q- ada ## ADA         IDRIS ##
#-Q- ada #######################
#-Q- ada # @ job_type = serial
#-Q- ada # @ requirements = (Feature == "prepost")
#-Q- ada # Temps Elapsed max. d'une requete hh:mm:ss
#-Q- ada # @ wall_clock_limit = 4:00:00
#-Q- ada # Nom du travail LoadLeveler
#-Q- ada # @ job_name   = METRICS_LMDZ
#-Q- ada # Fichier de sortie standard du travail
#-Q- ada # @ output     = $(job_name).$(jobid)
#-Q- ada # Fichier de sortie d'erreur du travail
#-Q- ada # @ error      =  $(job_name).$(jobid)
#-Q- ada # pour recevoir un mail en cas de depassement du temps Elapsed (ou autre pb.)
#-Q- ada # @ notification = error
#-Q- ada # @ environment  = $POST_DIR ; $Script_Post_Output ; $SUBMIT_DIR ; $libIGCM ; $DEBUG_debug ; $BigBrother ; $R_INIT ; $R_BC ; $StandAlone ; $RESOL_ATM ; $RESOL_OCE ; $RESOL_ICE ; $RESOL_MBG ; $RESOL_SRF ; $RESOL_SBG ; $YEARS ; $DIR_OUT ; $MASTER ; $libIGCM_SX
#-Q- ada # @ queue
#-Q- lxiv8 ######################
#-Q- lxiv8 ## OBELIX      LSCE ##
#-Q- lxiv8 ######################
#-Q- lxiv8 #PBS -N METRICS_LMDZ
#-Q- lxiv8 #PBS -m a
#-Q- lxiv8 #PBS -j oe
#-Q- lxiv8 #PBS -q test
#-Q- lxiv8 #PBS -o METRICS_LMDZ.$$
#-Q- lxiv8 #PBS -S /bin/ksh
#-Q- default #!/bin/ksh
#-Q- default ##################
#-Q- default ## DEFAULT HOST ##
#-Q- default ##################

#**************************************************************
# Author: Sebastien Denvil
# Contact: Sebastien.Denvil__at__ipsl.jussieu.fr
# $Revision::                                          $ Revision of last commit
# $Author::                                            $ Author of last commit
# $Date::                                              $ Date of last commit
# IPSL (2006)
#  This software is governed by the CeCILL licence see libIGCM/libIGCM_CeCILL.LIC
#
#**************************************************************

#set -eu
#set -vx

date

#-Q- ada export OMP_NUM_THREADS=1

#D- Task type (computing or post-processing)
TaskType=post-processing

######################################################################

. ${libIGCM}/libIGCM_debug/libIGCM_debug.ksh
. ${libIGCM}/libIGCM_card/libIGCM_card.ksh
. ${libIGCM}/libIGCM_date/libIGCM_date.ksh
#-------
. ${libIGCM}/libIGCM_sys/libIGCM_sys.ksh
. ${libIGCM}/libIGCM_config/libIGCM_config.ksh
#-------
( ${DEBUG_debug} ) && IGCM_debug_Check
( ${DEBUG_debug} ) && IGCM_card_Check
( ${DEBUG_debug} ) && IGCM_date_Check

######################################################################

#set -vx

function METRICS_sed
{
# Usage : METRICS_sed metrics_template var_name myvalue
#         In file metrics_template modify the line containing var_name=xxx into var_name=myvalue
    IGCM_debug_PushStack "METRICS_sed"
    
    # Test if the fichier exist
    if [ ! -f ${1} ] ; then
        echo "WARNING : ${1} file does not exist. Following will not be done : METRICS_sed : ${1} ${2} ${3}"
        IGCM_debug_PopStack "METRICS_sed"
        return
    fi

    sed -e "s/^${2}\ *=.*/${2} = ${3}/" ${1} > ${1}.tmp
    RET=$?
    echo "METRICS_sed : ${1} ${2} ${3}"
    \mv ${1}.tmp ${1}

    IGCM_debug_PopStack "METRICS_sed"
    return $RET
}

IGCM_sys_MkdirWork ${RUN_DIR_PATH}
IGCM_sys_Cd ${RUN_DIR_PATH}

# ------------------------------------------------------------------
# Test if all was right before proceeding further
# ------------------------------------------------------------------
IGCM_debug_Verif_Exit_Post

if [ ${StandAlone} = true ] ; then
    CARD_DIR=${SUBMIT_DIR}
else
    CARD_DIR=${RUN_DIR_PATH}
    IGCM_sys_Get_Master ${SUBMIT_DIR}/config.card ${RUN_DIR_PATH}
    IGCM_sys_Get_Master ${SUBMIT_DIR}/run.card    ${RUN_DIR_PATH}
    IGCM_sys_Get_Master ${SUBMIT_DIR}/COMP        ${RUN_DIR_PATH}
    IGCM_sys_Get_Master ${SUBMIT_DIR}/POST        ${RUN_DIR_PATH}
fi

#==================================
# First of all
#
# Read libIGCM compatibility version in config.card
# Read UserChoices section
# Read Ensemble section
# Read Post section
# Define all netcdf output directories
#==================================
IGCM_config_CommonConfiguration ${CARD_DIR}/config.card

#==================================
# Then set up several directories
#==================================
R_METRICS=METRICS/SE_${YEARS}
R_OUT_ATM=${R_SAVE}/ATM
DIR_OUT=${DIR_OUT:=SE}
RUN_DIR=${RUN_DIR_PATH}/METRICS_LMDZ_${YEARS}
IGCM_sys_MkdirArchive ${R_FIGR}
IGCM_sys_MkdirWork ${RUN_DIR}
IGCM_sys_Cd ${RUN_DIR}

#==================================
# Get the file
#==================================
IGCM_sys_Get ${R_OUT_ATM}/Analyse/${DIR_OUT}/${config_UserChoices_JobName}_SE_${YEARS}_1M_histmth.nc ${RUN_DIR}

#==================================
# SET UP THE PATH (will be done in libIGCM_sys later)
#==================================
export PATH=${PCMDI_MP}/work_install_v1/PCMDI_METRICS/bin:${PATH}

#==================================
# Appropriate configuration
#==================================
cp ${PCMDI_MP}/work_install_v1/PCMDI_METRICS/doc/input_parameters_sample_for_libIGCM_template.py metrics_template.py

#==================================
# Configure template
#==================================
IGCM_comp_modifyDefFile force metrics_template.py path_and_filename \'${config_UserChoices_JobName}_SE_${YEARS}_1M_histmth.nc\'
IGCM_comp_modifyDefFile force metrics_template.py attributes_provided \'IPSL_Extended\'
IGCM_comp_modifyDefFile force metrics_template.py case_id \'SE_${YEARS}\'
IGCM_comp_modifyDefFile force metrics_template.py obs_data_path \'${PCMDI_MP}/references-for-metrics/obs\'
IGCM_comp_modifyDefFile force metrics_template.py metrics_output_path \'${RUN_DIR}/METRICS\'
IGCM_comp_modifyDefFile force metrics_template.py model_clims_interpolated_output \'${RUN_DIR}\'

#==================================
# Compute metrics
#==================================
python ${PCMDI_MP}/work_install_v1/PCMDI_METRICS/bin/pcmdi_metrics_driver.py -p metrics_template.py

#==================================
# Save files
#==================================
IGCM_sys_Put_Dir METRICS ${R_FIGR}
IGCM_sys_Rm -Rf METRICS
IGCM_sys_Rm ${config_UserChoices_JobName}_SE_${YEARS}_1M_histmth.nc

#==================================
# DODS copy
#==================================
if [ X${config_Post_DodsCopy} = XTRUE ]  ; then
  IGCM_sys_TestDirArchive ${R_FIGR}/${R_METRICS}
  [ $? = 0 ] && IGCM_sys_Put_Dods ${R_METRICS}
fi
# Clean RUN_DIR_PATH (necessary for cesium and titane only)
IGCM_sys_RmRunDir -Rf ${RUN_DIR_PATH}