source: trunk/UTIL/NCL/visu-ncl

Last change on this file was 808, checked in by slebonnois, 12 years ago

SL: Many changes for VENUS (related to newstart) and TITAN (related to clouds). Please read DOC/chantiers/commit_importants.log (cf v808).
File size: 5.0 KB
Line 
1#!/bin/bash
2
3# Directory where the NCL files are located  ** EDIT **
4
5export NCLDIR=~/LMDZ5/UTIL/NCL
6
7########################################################################
8# Mise en place des valeurs par defaut
9########################################################################
10
11inFile="nofile"
12var="liste"
13labelvar="dummy"
14var2="dummy"
15media="x11"
16outputname="quickview"
17typeplot="coupe"
18ntcutmin=0
19ntcutmax=$ntcutmin
20tavg=-999
21pcutmin=1.e4
22pcutmax=$pcutmin
23pavg=-999
24lcutmin=0
25lcutmax=$lcutmin
26lavg=-999
27Lcutmin=0
28Lcutmax=$Lcutmin
29Lavg=-999
30# (if limits are identical, average need to be made to trim the dimension)
31
32# -888 => default variable chosen by NCL
33valmin=-888
34valmax=-888
35step=-888
36valmin2=-888
37valmax2=-888
38step2=-888
39
40Lcenter=0
41lcenter=40
42grid="True"
43gLspc=30
44glspc=20
45gcol="white"
46
47########################################################################
48# lecture des options
49########################################################################
50
51while (($# > 0))
52  do
53  case $1 in
54      "-h") cat $NCLDIR/help.txt
55          exit;;
56      "-f")
57          inFile=$2 ; shift ; shift ;;
58      "-v")
59          var=${2%%,*}
60          tmp=${2#*,}
61          labelvar=${tmp%%,*}
62          if [ $labelvar != $var ]; then
63            labelvar=${tmp#*,}
64          fi
65          shift ; shift ;;
66      "-v2")
67          var2=$2 ; shift ; shift ;;
68      "-med")
69          media=$2 ; shift ; shift ;;
70      "-oname")
71          outputname=$2 ; shift ; shift ;;
72      "-type")
73          typeplot=$2
74          if [ $typeplot = "stereo" ]; then
75            lcenter=90
76          fi
77          shift ; shift ;;
78      "-tc")
79          ntcutmin=${2%%,*}
80          tmp=${2#*,}
81          ntcutmax=${tmp%%,*}
82          if [ $ntcutmin = $ntcutmax ]; then
83            tavg=-999
84          else
85            tavg=${tmp#*,}
86          fi
87          shift ; shift ;;
88      "-pc")
89          pcutmin=${2%%,*}
90          tmp=${2#*,}
91          pcutmax=${tmp%%,*}
92          if [ $pcutmin = $pcutmax ]; then
93            pavg=-999
94          else
95            pavg=${tmp#*,}
96          fi
97          shift ; shift ;;
98      "-lc")
99          lcutmin=${2%%,*}
100          tmp=${2#*,}
101          lcutmax=${tmp%%,*}
102          if [ $lcutmin = $lcutmax ]; then
103            lavg=-999
104          else
105            lavg=${tmp#*,}
106          fi
107          shift ; shift ;;
108      "-Lc")
109          Lcutmin=${2%%,*}
110          tmp=${2#*,}
111          Lcutmax=${tmp%%,*}
112          if [ $Lcutmin = $Lcutmax ]; then
113            Lavg=-999
114          else
115            Lavg=${tmp#*,}
116          fi
117          shift ; shift ;;
118       "-lim")
119          valmin=${2%%,*}
120          tmp=${2#*,}
121          valmax=${tmp%%,*}
122          if [ $valmin = $valmax ]; then
123          # si pas les deux, remises au default
124             valmin=-888
125             valmax=-888
126          else
127            step=${tmp#*,}
128          fi
129          if [ $step = $valmax ]; then
130          # si pas de step, il est remis au default
131             step=-888
132          fi
133          shift ; shift ;;
134       "-lim2")
135          valmin2=${2%%,*}
136          tmp=${2#*,}
137          valmax2=${tmp%%,*}
138          if [ $valmin2 = $valmax2 ]; then
139          # si pas les deux, remises au default
140             valmin2=-888
141             valmax2=-888
142          else
143            step2=${tmp#*,}
144          fi
145          if [ $step2 = $valmax2 ]; then
146          # si pas de step, il est remis au default
147             step2=-888
148          fi
149          shift ; shift ;;
150       "-center")
151          Lcenter=${2%%,*}
152          tmp=${2#*,}
153          lcenter=${tmp%%,*}
154          shift ; shift ;;
155       "-grid")
156          if [ $2 = "T" -o $2 = "y" -o $2 = "Y" ]; then
157            grid="True"
158          fi
159          if [ $2 = "F" -o $2 = "n" -o $2 = "N" ]; then
160            grid="False"
161          fi
162          shift ; shift ;;
163      "-gcol")
164          gcol=$2 ; shift ; shift ;;
165       "-gspc")
166          gLspc=${2%%,*}
167          tmp=${2#*,}
168          glspc=${tmp%%,*}
169          shift ; shift ;;
170  esac
171done
172
173###############################################################
174if [ $inFile = "nofile" ]; then
175        echo "You have to specify a file ! ( -f filename )"; exit
176fi 
177
178if [ -f $inFile ]; then
179        echo "File "$inFile" exists"
180else 
181        echo "File "$inFile" does not exists"; exit
182fi 
183
184if [ $typeplot = "ortho" ]; then
185    # lat and lon put to 2d, defaulting to full map in planeto.ncl
186        lcutmin=-90
187        lcutmax=90
188        lavg=0
189        Lcutmin=0
190        Lcutmax=180
191        Lavg=0
192fi 
193
194export SIMUDIR=`pwd`
195inFile=$SIMUDIR/$inFile
196cd $NCLDIR
197ncl inFile=\"$inFile\"  \
198    var=\"$var\" \
199    labelvar=\"$labelvar\" \
200    var2=\"$var2\" \
201    media=\"$media\" \
202    outputname=\"$outputname\" \
203    typeplot=\"$typeplot\" \
204    ntcutmin=$ntcutmin \
205    ntcutmax=$ntcutmax \
206    tavg=$tavg \
207    pcutmin=$pcutmin \
208    pcutmax=$pcutmax \
209    pavg=$pavg \
210    lcutmin=$lcutmin \
211    lcutmax=$lcutmax \
212    lavg=$lavg \
213    Lcutmin=$Lcutmin \
214    Lcutmax=$Lcutmax \
215    Lavg=$Lavg \
216    valmin=$valmin \
217    valmax=$valmax \
218    step=$step \
219    valmin2=$valmin2 \
220    valmax2=$valmax2 \
221    step2=$step2 \
222    Lcenter=$Lcenter \
223    lcenter=$lcenter \
224    grid=$grid \
225    gLspc=$gLspc \
226    glspc=$glspc \
227    gcol=\"$gcol\" \
228    planeto.ncl
229
230if [ -f *pdf ]; then 
231   mv *pdf $SIMUDIR
232fi
233
234if [ -f *ps ]; then 
235   mv *ps $SIMUDIR
236fi
237
Note: See TracBrowser for help on using the repository browser.