Changeset 693

Timestamp:

06/27/12 18:12:50 (12 years ago)

Author:

sdipsl

Message:

• add a function dedicated to rebuild CFMIP file station. Regular rebuild can't do it.

File:

• trunk/libIGCM/libIGCM_sys/libIGCM_sys_mercure.ksh (modified) (1 diff)

trunk/libIGCM/libIGCM_sys/libIGCM_sys_mercure.ksh

@@ -1484 +1484 @@
 }
 
+function IGCM_sys_rebuild_station {
+  IGCM_debug_PushStack "IGCM_sys_rebuild_station" -- $@
+  typeset i list_opt file_in file_out prefix_invert
+  if ( $DEBUG_sys ) ; then
+    echo "IGCM_sys_rebuild_station :" $@
+  fi
+  list_opt=$@
+
+  # Invert Axis : t,x -> x,t 
+  #               t,pres,x -> x,t,pres
+  # So that we can concatenate along x
+  i=0
+  for file_in in ${list_opt} ; do
+    (( i = i + 1))
+    [ ${i} = 1 ] && file_out=${file_in} && continue
+    prefix_invert=$( basename ${file_in} .nc )
+    IGCM_sys_ncpdq -a x,time_counter -a x,time_counter,presnivs ${file_in} ${prefix_invert}_xt.nc
+    list_invert[${#list_invert[*]}]=${prefix_invert}_xt.nc
+  done
+
+  # Concatenate
+  IGCM_sys_ncrcat ${list_invert} histstn_xt.nc
+
+  # Re-ivert file
+  IGCM_sys_ncpdq -a time_counter,x -a time_counter,presnivs,x histstn_xt.nc ${file_out}
+
+  IGCM_debug_PopStack "IGCM_sys_rebuild_station"
+}
+
 ############################################################
 # Check of space available on temporary filesytems