Changeset 33 for trunk/detectbarotopemode.sh

Timestamp:

11/26/07 12:39:53 (17 years ago)

Author:

pinsard

Message:

introduce a shell to get outputs of model in GEOMAG environement

File:

• trunk/detectbarotopemode.sh (modified) (1 diff)

trunk/detectbarotopemode.sh

@@ -2 +2 @@
 #
 # module :
-# build file for detection of baratrope mode
+# build file for detection of barotrope mode
 #
 # usage :
-# $ ./detectbaratropemode.sh ++
+# $ detectbarotropemode.sh -l list
+#
+# example :
+# $ detectbarotropemode.sh -l ${GEOMAG_OD}/list_G42
 #
 # update :
 # $Id$
 # pas fini ++
+# reee522 2007-11-26T11:33:42Z rhodes (IRIX64)
+# externalization of  copy of data file in getmodelout.sh
+# reee522 2007-11-26T10:55:44Z rhodes (IRIX64)
+# add ncks, ncrcat and some informations inspired by
+# http://dods.ipsl.jussieu.fr/projets_ES/misc/myls
 # fplod 2007-11-23T17:02:56Z aedon.locean-ipsl.upmc.fr (Darwin)
 # creation
 # 
+command=detectbarotopemode.sh
+#
 # ++ simulation parametres
 yearmin=8
 yearmax=10
+latmin=-50
+latmax=-30
+longmin=295
+longmax=335
 grid=gridT
+orcares=ORCA025
 DRAKKAR_EXP=G42
 #
+usage=" Usage : ${command} -f list"
+#
+while [ ! -z "${1}" ] # ++ pb bash
+do
+ case ${1} in
+ -l)
+    list=${2}
+    shift
+ ;;
+
+ *) # other choice
+    echo "${usage}"
+ ;;
+ esac
+ shift # next flag
+done
+#
+set -u
+#
+# check GEOMAG environement
+if [ ! -d ${GEOMAG_LOG} ]
+then
+ echo " eee : \${GEOMAG_LOG} not found"
+ exit 1
+fi
 # ++ blindage en tout genre
 #
 # 
-#++list=/tmp/${LOGNAME}/listfiles${$}
-list=./listfiles${$}
+# ++ conversion of latitude and longitude in x and y indexes
+case in ${orcares} in
+ORCA025)
+   xmin=10
+   xmax=30
+   ymin=10
+   ymax=40
+;
+*)
+  echo "eee : pb conversion lat,long to x,y"
+  exit 1
+;;
+esac
 #
-# pour info liste des fichiers à prendre en compte
-iyear=${yearmin}
-while [ ${iyear} -lt ${yearmax} ]
-do 
-  cyear=`printf "%4.4d" ${iyear}`
-  ls /u/rech/cli/rcli002/ORCA025/ORCA025-${DRAKKAR_EXP}-S/${cyear}/ORCA025-${DRAKKAR_EXP}_y${cyear}m??d??_${grid}.nc >> ${list}
-
-  # ++ calcul taille cf myls dans faq es dods
-  # iyear=$(( ${iyear} + 1 )) ++ pas ok sur rhodes
-  # iyear=$(expr  ${iyear} + 1 )  ++ pas ok sur rhodes
-  iyear=`expr  ${iyear} + 1 ` 
+# extract ssh and create a new set of files
+list_ssh=""
+for file in  `cat ${list}`
+do
+ fileb=`basename ${file}`
+ filessh=`basename ${file} .nc`_ssh.nc
+ rm -f ${GEOMAG_OD}/${filessh} 2> /dev/null
+ ncks --dimension x,${xmin},${xmax} --dimension y,${ymin},${ymax}  \
+       -v sossheig \
+ ${GEOMAG_ID}/${fileb} ${GEOMAG_OD}/${filessh}
+ status_ncks=${?}
+ if [ ${status_ncks} -ne 0 ]
+ then
+    echo "eee : pb with ncks ${file}"
+    exit 1
+ else
+    list_ssh="${list_ssh} ${file_ssh}"
+ fi
 done
 #
-# put files in ${GEOMAG_ID}
-# ++ ln cp ou scp 
-# ++ -f ou test si deja existant 
-for file in  `cat ${list}`
-do
-  ln -sf ${file} ${GEOMAG_ID}
-done 
-#
-#
-# extract ssh and create a new set of files
-for file in  `cat ${list}`
-do
- filessh=`basename ${file} .nc`_ssh.nc
- # ++ test si exist deja
- ncks -v sossheig ${GEOMAG_ID}/${file} ${GEOMAG_OD}/${filessh} #++ lat long et sortie x, y, t
- #++ test status
- ls -l  ${GEOMAG_ID}/${file} ${GEOMAG_OD}/${filessh}
- read a
-done
-#
+datemin=19980101#++
+datemax=19991231#++
 # concatenation ++
-# ++ ncrcat
-
+varlist=`cat ${list_ssh}`
+filetot=ssh_zone_${datemin}_${datemax}_${orcares}_${DRAKKAR_EXP}.nc
+ncrcat ${varlist} ${filetot}
+status_ncrcat=${?}
+if [ ${status_ncrcat} -ne 0 ]
+then
+   echo "eee : pb with ncrcat"
+   exit 1
+fi
 # 
 # menage
 rm ${list}
+rm ${list_ssh}
 #
 # end