Changeset 370 for trunk/libIGCM/AA_create_se

Timestamp:

10/27/10 14:44:39 (14 years ago)

Author:

mafoipsl

Message:

Add AA_create_multi_se script to build SE_nnnY based on previous SE.
Use optionnal MultiSeasonalFrequency?[=50 by default] parameter in config.card Post Section.
Use climatological time axis for SE and SE_nnnY files.
Add DIR_OUT to allow atlas to work with SE_nnnY.
Add POST_DIR default value for interactive session.

File:

• trunk/libIGCM/AA_create_se (modified) (6 diffs)

trunk/libIGCM/AA_create_se

@@ -133 +133 @@
 #D- Flag to determine atlas job's output directory
 #D- Default : value from libIGCM_post.ksh if any
-POST_DIR=${POST_DIR:=${PBS_O_WORKDIR}}
+POST_DIR=${POST_DIR:=${PBS_O_WORKDIR:=$(pwd)}}
 
 #D- Increased verbosity (1, 2, 3)
@@ -397 +397 @@
           IGCM_sys_Rm tmpfile_ncra
 
+          # detect time counter : time or time_counter ie unlimited variable
+          var_unlim=$(ncdump -h ${liste_file_se_tmp[1]}|grep UNLIMITED| cut -d ' ' -f 1)
+
           for mois in  01 02 03 04 05 06 07 08 09 10 11 12 ; do
           #
@@ -402 +405 @@
           #
               OUT_SE[${nb}]=${config_UserChoices_JobName}_SE${mois}_${an_deb}_${an_fin}_${file}.nc
-              IGCM_sys_ncra -O -v ${list_var_final_ncra} -d time_counter,${nb},,12 ${liste_file_se_tmp[*]} ${OUT_SE[${nb}]}
+              IGCM_sys_ncra -O -v ${list_var_final_ncra} -d ${var_unlim},${nb},,12 ${liste_file_se_tmp[*]} ${OUT_SE[${nb}]}
           done
           IGCM_sys_Rm ${liste_file_se_tmp[*]}
@@ -459 +462 @@
               list_var_final_ncra=`cat tmpfile_ncra | sort | uniq -c | awk -v nbfile=$nbfile '{if ($1 == nbfile) {print $2}}' | paste -s -d ','`
               IGCM_sys_Rm tmpfile_ncra
+
               IGCM_sys_ncra -O -v ${list_var_final_ncra} ${liste_file_se_tmp[*]} ${OUT_SE[${nb}]}
 
@@ -550 +554 @@
         IGCM_sys_Rm ${OUT_SE[*]}
         #
-        if [ ${comp} = "CPL" ] ; then
-              ${ncatted} -a modulo,time,c,c," " ${RESULT_SE}
-        else
-              ${ncatted} -a modulo,time_counter,c,c," " ${RESULT_SE}
-        fi
+        # detect time counter : time or time_counter ie unlimited variable
+        var_unlim=$(ncdump -h ${RESULT_SE}|grep UNLIMITED| cut -d ' ' -f 1)
+
+        # add modulo attribute on time axis
+        ${ncatted} -a modulo,${var_unlim},c,c," " ${RESULT_SE}
         #
         eval IGCM_sys_Put_Out ${RESULT_SE} \${R_OUT_${comp}}/Analyse/SE/${RESULT_SE}
@@ -571 +575 @@
 
 YEARS=${an_deb}_${an_fin}
+DIR_OUT=SE
 
 export POST_DIR; export SUBMIT_DIR; export libIGCM; export R_INIT; export R_BC; export StandAlone
 export RESOL_ATM; export RESOL_OCE ; export RESOL_ICE ; export RESOL_MBG ; export RESOL_SRF ; export RESOL_SBG
-export R_SAVE; export config_UserChoices_JobName; export config_UserChoices_TagName; export YEARS; 
-listVarEnv="POST_DIR,SUBMIT_DIR,libIGCM,R_INIT,R_BC,StandAlone,RESOL_ATM,RESOL_OCE,RESOL_ICE,RESOL_MBG,RESOL_SRF,RESOL_SBG,R_SAVE,config_UserChoices_JobName,config_UserChoices_TagName,YEARS"
+export R_SAVE; export config_UserChoices_JobName; export config_UserChoices_TagName; export YEARS; export DIR_OUT 
+listVarEnv="POST_DIR,SUBMIT_DIR,libIGCM,R_INIT,R_BC,StandAlone,RESOL_ATM,RESOL_OCE,RESOL_ICE,RESOL_MBG,RESOL_SRF,RESOL_SBG,R_SAVE,config_UserChoices_JobName,config_UserChoices_TagName,YEARS,DIR_OUT"
 export listVarEnv;