Changeset 35

Timestamp:

11/26/07 16:20:37 (16 years ago)

Author:

pinsard

Message:

work in progress

Location:

trunk

Files:

• detectbarotropicmode.sh (modified) (7 diffs)
• getmodelout.sh (modified) (6 diffs)

trunk/detectbarotropicmode.sh

@@ -1 @@
-#!/bin/sh -x
+#!/bin/sh
 #
 # module :
-# build file for detection of barotrope mode
-#
-# usage :
-# $ detectbarotropemode.sh -l list
+# build file for detection of barotropic mode
 #
 # example :
-# $ detectbarotropemode.sh -l ${GEOMAG_OD}/list_G42
+# $ detectbarotropicmode.sh -l ${GEOMAG_OD}/list_G42 -r ORCA025 -exp G42 -g gridT \
+# -latmin -50 -latmax -30 -longmin 295 -longmax 335
 #
 # update :
@@ -20 +18 @@
 # 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
+command=detectbarotropicmode.sh
 #
-usage=" Usage : ${command} -f list"
+usage=" Usage : ${command} -l list -r orcares -exp drakkar_exp -g grid -latmin latmin -latmax latmax -longmin longmin -longmax longmax"
 #
 while [ ! -z "${1}" ] # ++ pb bash
@@ -43 +30 @@
     shift
  ;;
-
+ -r)
+    orcares=${2}
+    shift
+ ;;
+ -exp)
+    drakkar_exp=${2}
+    shift
+ ;;
+ -g)
+    grid=${2}
+    shift
+ ;;
+ -latmin)
+    latmin=${2}
+    shift
+ ;;
+ -latmax)
+    latmax=${2}
+    shift
+ ;;
+ -longmin)
+    longmin=${2}
+    shift
+ ;;
+ -longmax)
+    longmax=${2}
+    shift
+ ;;
  *) # other choice
     echo "${usage}"
@@ -61 +75 @@
 # ++ blindage en tout genre
 #
-# 
+case ${grid} in
+gridT)
+  echo "iii : grid=${orcares}"
+ ;;
+*)
+  echo "eee : pb \${grid} = ${grid}"
+  exit 1
+;;
+esac
+#
+case ${orcares} in
+ORCA025)
+  echo "iii : orcares=${orcares}"
+ ;;
+*)
+  echo "eee : pb \${orcares} = ${orcares}"
+  exit 1
+;;
+esac
+#
+case ${drakkar_exp} in
+G42)
+  echo "iii : drakkar_exp=${drakkar_exp}"
+ ;;
+*)
+  echo "eee : pb \${drakkar_exp} = ${drakkar_exp}"
+  exit 1
+;;
+esac
+#
 # ++ conversion of latitude and longitude in x and y indexes
-case in ${orcares} in
+case ${orcares} in
 ORCA025)
-   xmin=10
-   xmax=30
-   ymin=10
-   ymax=40
-;
+   xmin=529
+   xmax=689
+   ymin=266
+   ymax=372
+   echo "iii : forcage xmin xmax ymin ymax : ${xmin} ${xmax} ${ymin} ${ymax}"
+;;
 *)
   echo "eee : pb conversion lat,long to x,y"
@@ -80 +124 @@
 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}
+ echo "iii : extraction from ${file}" # ++ voir seulement si debug demandé
+ filessh=${GEOMAG_OD}/`basename ${file} .nc`_ssh.nc
+ rm -f ${filessh} 2> /dev/null
+ ncks --dimension x,${xmin},${xmax} --dimension y,${ymin},${ymax}  -v sossheig,nav_lon,nav_lat \
+ ${file} ${filessh}
  status_ncks=${?}
  if [ ${status_ncks} -ne 0 ]
@@ -92 +135 @@
     exit 1
  else
-    list_ssh="${list_ssh} ${file_ssh}"
+    list_ssh="${list_ssh} ${filessh}"
  fi
 done
+line=`head -n 1  ${list}`
+bfile=`basename ${line}`
+datemin=`echo ${bfile}  | awk -F "_" '{print $2}'`
+line=`tail -n 1  ${list}`
+bfile=`basename ${line}`
+datemax=`echo ${bfile}  | awk -F "_" '{print $2}'`
 #
-datemin=19980101#++
-datemax=19991231#++
-# concatenation ++
-varlist=`cat ${list_ssh}`
-filetot=ssh_zone_${datemin}_${datemax}_${orcares}_${DRAKKAR_EXP}.nc
-ncrcat ${varlist} ${filetot}
+# concatenation
+filetot=${GEOMAG_OD}/${orcares}-${drakkar_exp}_${datemin}_${datemax}_${grid}_ssh.nc
+rm -f  ${filetot} 2> /dev/null
+ncrcat ${list_ssh} ${filetot}
 status_ncrcat=${?}
 if [ ${status_ncrcat} -ne 0 ]
@@ -107 +154 @@
    echo "eee : pb with ncrcat"
    exit 1
+else
+   echo "iii : result in ${filetot}"
 fi
-# 
-# menage
-rm ${list}
+#
+# ++ global attributes
+# ++ convention
+# ++ nedtcf output or ascii
+# clean
 rm ${list_ssh}
 #
 # end
-exit 1
+exit 0

trunk/getmodelout.sh

@@ -1 @@
-#!/bin/sh -x
+#!/bin/sh
 #
 # module :
@@ -5 +5 @@
 #
 # usage :
-# $ getmodelout.sh -o list ++
+# $ getmodelout.sh -r orcares -exp drakkar_exp -g grid -yearmin yyyy -yearmax yyyy -o list
 #
 # example :
-# $ getmodelout.sh -o ${GEOMAG_OD}/list_G42
+# $ getmodelout.sh -r ORCA025 -exp G42 -g gridT -yearmin 9 -yearmax 10 -o ${GEOMAG_OD}/list_G42
 #
 # update :
@@ -18 +18 @@
 command=getmodelout.sh
 #
-# ++ simulation parametres
-yearmin=8
-yearmax=10
-grid=gridT
-orcares=ORCA025
-DRAKKAR_EXP=G42
-#
-usage=" Usage : ${command} -o list"
+usage=" Usage : ${command} -r orcares -exp drakkar_exp -g grid -o list"
 #
 while [ ! -z "${1}" ] # ++ pb bash
 do
  case ${1} in
+ -r)
+    orcares=${2}
+    shift
+ ;;
+ -exp)
+    drakkar_exp=${2}
+    shift
+ ;;
+ -g)
+    grid=${2}
+    shift
+ ;;
+ -yearmin) 
+   yearmin=${2}
+   shift
+ ;;
+ -yearmax) 
+   yearmax=${2}
+   shift
+ ;;
  -o) 
    list=${2}
@@ -49 +62 @@
  exit 1
 fi
+#
+if [ -f ${list} ]
+then
+  echo "eee : ${list} exist"
+  exit 1
+fi
 # ++ blindage en tout genre
 #
+case ${orcares} in
+ORCA025)
+   dirref=/u/rech/cli/rcli002/ # sur l'idris
+;;
+*)
+   echo "eee : pb \${orcares} = ${orcares}"
+;;
+esac
+#
+listfull=${GEOMAG_LOG}/lisfull${$}
+#
 iyear=${yearmin}
-while [ ${iyear} -lt ${yearmax} ]
+while [ ${iyear} -le ${yearmax} ]
 do 
   cyear=`printf "%4.4d" ${iyear}`
-  ls /u/rech/cli/rcli002/${orcares}/${orcares}-${DRAKKAR_EXP}-S/${cyear}/${orcares}-${DRAKKAR_EXP}_y${cyear}m??d??_${grid}.nc >> ${list}
+  ls ${dirref}/${orcares}/${orcares}-${drakkar_exp}-S/${cyear}/${orcares}-${drakkar_exp}_y${cyear}m??d??_${grid}.nc >> ${listfull}
   # ++ calcul taille cf myls dans faq es dods
-  
   # iyear=$(( ${iyear} + 1 )) ++ pas ok sur rhodes
   # iyear=$(expr  ${iyear} + 1 )  ++ pas ok sur rhodes
@@ -65 +94 @@
 # print total size of model files
 # ++ pb syntaxe unitsop="1024**2" # Mo
-#ls -l `cat ${list}` | \
+#ls -l `cat ${listfull}` | \
 #awk 'NF>4 {print $5, $9}' | \
 #awk 'sum=$1/'$unitsop'+sum {printf("%08.3f '$units' -- %s\n",$1/'$unitsop', $2)} 
@@ -72 +101 @@
 # ++ ln cp ou scp 
 # ++ -f ou test si deja existant 
-for file in  `cat ${list}`
+for file in  `cat ${listfull}`
 do
   ln -sf ${file} ${GEOMAG_ID}
+  echo ${GEOMAG_ID}/`basename ${file}` >> ${list}
 done 
 #