1 | #!/bin/ksh |
---|
2 | |
---|
3 | #************************************************************** |
---|
4 | # Author: Martial Mancip |
---|
5 | # Contact: Martial.Mancip_ipsl.jussieu.fr |
---|
6 | # $Date$ |
---|
7 | # $Author$ |
---|
8 | # $Revision$ |
---|
9 | # IPSL (2006) |
---|
10 | # This software is governed by the CeCILL licence see libIGCM/libIGCM_CeCILL.LIC |
---|
11 | # |
---|
12 | #************************************************************** |
---|
13 | |
---|
14 | #!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! |
---|
15 | # We keep only first time axis as time_counter !!!! |
---|
16 | #!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! |
---|
17 | |
---|
18 | function IGCM_Patch_20070220_histcom_time_axis { |
---|
19 | IGCM_debug_PushStack "IGCM_Patch_histcom_time_axis" $@ |
---|
20 | |
---|
21 | typeset file AXISt AXISt_i VARName var ATTN ATTNV VAR_AXISt VAR_AXISt_i |
---|
22 | |
---|
23 | ncdump -h ${1} > ${1}_ncdump-h |
---|
24 | set +A AXISt $(cat ${1}_ncdump-h | grep 't_ave_.*(time_counter)' | sed -e 's/[[:space:]]float t_ave_\([0-9]*\).*/\1/') |
---|
25 | # liste => tous les axes de temps vont s'appeler time_counter !! |
---|
26 | if [ $? -eq 0 ] ; then |
---|
27 | |
---|
28 | file=t_$1 |
---|
29 | IGCM_sys_Mv ${1} ${file} |
---|
30 | IGCM_sys_Chmod 644 ${file} |
---|
31 | |
---|
32 | if [ ${#AXISt[*]} -gt 1 ] ; then |
---|
33 | echo "WARNING : for file ${1} we have more than one time axis !" |
---|
34 | echo "We choose first one ! = t_ave_${AXISt[0]}." |
---|
35 | fi |
---|
36 | #!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! |
---|
37 | # We keep only first time axis as time_counter !!!! |
---|
38 | #!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! |
---|
39 | |
---|
40 | AXISt_i=${AXISt[0]} |
---|
41 | # Rename time_counter variable |
---|
42 | IGCM_sys_ncrename -O --hst -v t_ave_${AXISt_i},time_counter ${file} |
---|
43 | ncdump -h ${file} > ${file}_ncdump-h |
---|
44 | |
---|
45 | # We have also to rename variables indexed by this time value |
---|
46 | # liste => pb avec deux variables indéxées par deux valeurs de temps différentes => même nom au final! |
---|
47 | set +A VAR_AXISt $( cat ${file}_ncdump-h | grep "[[:space:]]*[a-z]* [a-zA-Z_]*${AXISt_i}(.*) ;" | sed -e "s/[[:space:]]*[a-z]* \([a-zA-Z_]*\)_${AXISt_i}(.*) ;/\1/" ) |
---|
48 | for VAR_AXISt_i in ${VAR_AXISt[*]} ; do |
---|
49 | IGCM_sys_ncrename -O --hst -v ${VAR_AXISt_i}_${AXISt_i},${VAR_AXISt_i} ${file} |
---|
50 | done |
---|
51 | |
---|
52 | # A-t-on besoin de renommer les attributs des autres variables ?? |
---|
53 | # Rename relative associate string in other variables |
---|
54 | |
---|
55 | set +A VARName $( cat ${file}_ncdump-h | grep "[0-9a-zA-Z_]*.*:associate.*" | grep "t_ave_${AXISt_i}" | sed -e "s/[[:space:]]*\(.*\):associate.*/\1/" ) |
---|
56 | for var in ${VARName[*]} ; do |
---|
57 | # Old attribute value |
---|
58 | ATTN=$( cat ${file}_ncdump-h | grep "${var}:associate.*" | grep "t_ave_${AXISt_i}" | sed -e "s/.*:associate = \"\(.*\)\" ;/\1/" ) |
---|
59 | # New attribute value |
---|
60 | ATTNV=$( echo ${ATTN} | sed -e "s/t_ave_${AXISt_i}/time_counter/" ) #| sed -e 's/[[:space:]]/\\ /g') |
---|
61 | # change it ! |
---|
62 | ${ncatted} -O --hst -a associate,${var},m,c,"${ATTNV}" ${file} |
---|
63 | # This is not an IGCM function because of a |
---|
64 | # Problem with ksh and string passed in this IGCM_sys_ncatted function to ncatted ! |
---|
65 | done |
---|
66 | IGCM_debug_Verif_Exit_Post |
---|
67 | IGCM_sys_Mv ${file} ${1} |
---|
68 | IGCM_sys_Rm ${file}_ncdump-h |
---|
69 | fi |
---|
70 | IGCM_sys_Rm ${1}_ncdump-h |
---|
71 | |
---|
72 | IGCM_debug_PopStack "IGCM_Patch_histcom_time_axis" |
---|
73 | } |
---|