Philippe 02/05/2016: moved all LIBTOOLS files in LIBTOOLS directory
[MNH-git_open_source-lfs.git] / LIBTOOLS / readme / LATEX / lfi2grb.tex
1 \section{Conversion to GRIB or Vis5D files}
2
3 \subsection{Presentation}
4 FM synchronous file can be convert into \underline{GRIB} 
5 or \underline{Vis5D} format.
6 This section aims at describ how the converter works and how use it.
7
8 The GRIB (GRId in Binary) format is a standard meteorological one, defined 
9 by the WMO. GRIB files can be plotted with METVIEW
10 %\footnote{available on {\tt xdata} workstation in CNRM}
11 graphic interface (developped at ECMWF), or 
12 R2\footnote{used in the GMME/MICADO team at CNRM} software.
13
14 The Vis5D format is specified for using Vis5D\footnote{home page 
15 {\tt http://www.ssec.wisc.edu/\~ billh/vis5d.html}}
16 software (following the GNU General Public License): 3 spatial
17 dimensions, time dimension, 5$^{th}$ dimension for enumeration of variables.
18 It is rather designed for animation of 3D plotting.
19
20 Choice was made to put together the two file formats in a same conversion
21 program because in both cases specificities of Meso-NH grids have to be
22 treated in the same way (horizontally: Arakawa C-grid, vertically: Gal-Chen
23 coordinate $\hat z$ following terrain). However, the user has to choose one 
24 of the two formats available when running the tool
25 (see section \ref{s:execution}).
26
27
28 \subsection{Usage} \label{s:execution}
29 The interactive tool is
30 called {\tt lfi2grb} or {\tt lfi2v5d} according the wanted output
31 file format, but it runs the same program. Some questions are to be
32 answered to indicate the number and type of vertical levels, the type of
33 horizontal domain,
34 and the name of the variables to write into the output file.
35 All that is typed on keyboard is saved in {\tt dirconv.grb} or {\tt dirconv.v5d}
36 file, it can be appended and used as input (after renaming it) for the next call
37 of the tool (e.g. {\tt mv dirconv.grb dirgrb ; lfi2grb < dirgrb}).
38
39 For historical reasons, a program with the same goal of conversion to GRIB or 
40 Vis5d has been first developped as a main program
41 of MesoNH, as DIAG program is. This program called {\bf CONVLFI} runs with
42 the MesoNH procedure {\bf prepmodel} and
43 a namelist file {\tt CONVLFI1.nam} (see \ref{ss:convlfi}). 
44
45 To use the converter after a {\bf DIAG prepmodel} job, the Meso-NH file must
46 remain a synchronous file, not transformed onto a diachronic file:
47 in {\bf prepmodelrc} specify {\tt OUTFILE\_TOOLS='fm'} 
48 (default is 'conv2dia' to convert with {\tt conv2dia}).
49
50
51 \subsubsection{{\tt lfi2grb} tool}
52 When {\tt lfi2grb} tool is invoked, you must indicate, 
53 after the name of the input file, first 
54 the horizontal grid (type, eventually type of interpolation and domain),
55 the vertical grid (type and levels), 
56 then the list of the 3-dimensional fields to convert, 
57 and the list of the 2-dimensional ones.
58
59 For the \underline{horizontal grid}, you can either keep the one of MesoNH file
60 (cartesien or conformal projection) or interpolate onto a lat-lon regular grid.
61 In the first case, you can replace all the fields on mass points (A-grid)
62 or keep the native grid (C-grid).
63 In the second case, you have to indicate
64 the bounds of the domain with north and south latitudes and west and east
65 longitudes, as well as the type of horizontal interpolation:
66 nearest-neighbour value or bilinear interpolation with the 4 surrounding values.
67 The resolution of the lat.-lon. grid is automatically initialized 
68 with the equivalent value of the grid-mesh where the map scale is minimum. 
69 The program also indicates the number of grid points of the Meso-NH domain 
70 inside the prescribed lat-lon domain. If there are points of lat-lon domain
71 outside Meso-NH one, the value of the interpolated fields at these points
72 will be a missing one.
73
74 The \underline{vertical grid} can be either the native K levels or pressure
75 levels. 
76 In the first case ({\tt K}), all levels are kept and no interpolation is done:
77 the height specified in the GRIB header is the one of the grid without orography.
78 In the second case ({\tt P}), the list of pressure levels is either specified
79 manually or computed using a linear function from user-specified
80 minimum, maximum and increment values. If a prescribed level is below the lower
81 Meso-NH level or above the upper MesoNH level, the value of the field at this 
82 level will be a missing one. Otherwise, the value is computed from
83 a linear interpolation in log(P).
84
85 The \underline{3-dimensional fields} to convert are specified as follows: 
86 one field per line with first the name of the record in the input file
87 following by its grib code (tabular character is allowed). Note that no test
88 is done on the value of grib code (GRIB header {\sf ISEC1(6)}): you choose it
89 to easily identify the field with the software used after the conversion.
90 The end of the list is indicated by the keyword {\tt END}.
91
92 The \underline{2-dimensional fields} to convert are specified as follows:
93 one field per line with first the name of the record in the input file
94 (it can be a K-level of a 3-dimensional field too),
95 following by its grib code and possibly level indicator and level value
96 (tabular character is allowed).
97 Note that the value of the level indicator ({\sf ISEC1(7)}) is optional 
98 (the default value is 105: {\sf 'specified height above ground'}).
99 So is the level value ({\sf ISEC1(8)}), the default value is the altitude of 
100 the first mass point of the K-levels.
101 The end of the list is indicated by the keyword {\tt END}.
102
103 \subsubsection{Example of {\tt lfi2grb} use}
104 \begin{itemize}
105 \item to convert onto a GRIB file with horizontal and vertical interpolations in P levels:\\
106 (all that is typed on keyboard (in {\it italic} in the example below)
107 is saved in {\tt dirconv.grb})
108 \end{itemize}
109 \small
110 {\tt - ENTER FM synchronous FILE NAME (without .lfi) ?}  \\
111 {\tt\it CEXP.1.CSEG.001d } \hspace{3.5cm} $<$- the input file must be splitted in .des and .lfi \\
112 {\tt - Horizontal interpolation to lat-lon regular grid? (Y/y/O/o/N/n)}\\
113 {\tt\it y }  \\
114 {\tt - Type of interpolation? NEARest-neighbour (default) or BILInear }\\
115 {\tt\it NEAR }  \\
116 {\tt  - NSWE target domain bounds (in degrees)? }\\
117 {\tt\it 55. 35. -20. 10. }  \\
118 {\tt - Vertical grid: type K or P ? }\\
119 {\tt\it P }  \\
120 {\tt - Type of vertical grid: given by linear FUNCTN (default) or MANUALly ?}\\
121 {\tt\it FUNCTN }  \\
122 {\tt - Enter number of P levels ?} \\
123 {\tt\it 5 }  \\
124 {\tt - Values of the  5  P levels (hPa, from bottom to top):} \\
125 {\tt\it 1000. 850. 700. 500. 300. }  \\
126 {\tt - Enter 3D variables to CONVERT (1/1 line, end by END): }\\
127 {\tt  MesoNH field name, grib parameter indicator }\\
128 {\tt\it UM  33 }\\
129 {\tt - next 3D field or END ? }\\
130 {\tt\it VM  34 }\\
131 {\tt - next 3D field or END ? }\\
132 {\tt\it END }\\
133 {\tt - Enter 2D variables to CONVERT (1/1 line, end by END): }\\
134 {\tt  MesoNH field name, grib parameter indicator, eventually level indicator and level value}\\
135 {\tt\it T2M  13  105  2}\\
136 {\tt - next 2D field or END ? }\\
137 {\tt\it THM\_K\_2  13}\\
138 {\tt - next 2D field or END ? }\\
139 {\tt\it END}\\
140 {\tt 2 fields (3D), and   2 fields (2D) written in CEXP.1.CSEG.001d.GRB }\\
141  
142 \normalsize
143 \subsubsection{{\tt lfi2v5d} tool}
144 When {\tt lfi2v5d} tool is invoked, you must indicate, 
145 after the name of the input file, first 
146 the vertical grid (type and levels), 
147 then the list of the 3-dimensional fields to convert, 
148 and the list of the 2-dimensional ones.
149
150 No horizontal interpolation is available for the Vis5D format output: all the
151 converted fields are replaced on mass points (A-grid) of the MesoNH grid
152 (cartesien or conformal projection).
153
154 The \underline{vertical grid} can be either the native K levels, altitude
155 levels or pressure levels. 
156 In the first case ({\tt K}), all levels are kept and the fields are interpolated
157 on the levels of the lowest point of the domain.
158 In the second and third cases ({\tt Z} and {\tt P}), the list of levels is
159 either specified
160 manually or computed using a linear function from user-specified
161 minimum, maximum and increment values. The value of the field is computed from
162 a linear interpolation in Z or in log(P).
163
164 The \underline{3-dimensional fields} to convert are specified with 
165 one record name per line.
166 The end of the list is indicated by the keyword {\tt END}.
167
168 Then the \underline{2-dimensional fields},
169 or a K-level of 3-dimensional fields,
170 to convert are specified in the same way.
171
172 \subsubsection{Example of {\tt lfi2v5d} use}
173 \begin{itemize}
174 \item to convert onto a Vis5D file with vertical interpolation in Z levels:\\
175 (all that is typed on keyboard (in {\it italic} in the example below)
176 is saved in {\tt dirconv.v5d})
177 \end{itemize}
178 \small
179 {\tt - ENTER FM synchronous FILE NAME (without .lfi) ?}  \\
180 {\tt\it CEXP.1.CSEG.001 } \hspace{3.5cm} $<$- the input file must be splitted in .des and .lfi \\
181 {\tt - Verbosity level ?}  \\
182 {\tt\it 5 }  \\
183 {\tt - File 2D (xz): L2D=T or F ?}  \\
184 {\tt\it F }  \\
185 {\tt - Vertical grid: type K,Z or P ?}  \\
186 {\tt\it Z }  \\
187 {\tt - Type of vertical grid: given by linear FUNCTN (default) or MANUALly ?} \\
188 {\tt\it FUNCTN }  \\
189 {\tt - Vertical grid: min, max, int (m for Z, hPa for P)?} \\
190 {\tt\it 1500 9000 3000 }  \\
191 {\tt - Enter 3D variables to CONVERT (1/1 line, end by END): }\\
192 {\tt\it THM }  \\
193 {\tt - next 3D field or END ? }\\
194 {\tt\it POVOM }  \\
195 {\tt - next 3D field or END ? }\\
196 {\tt\it END }\\
197 {\tt - Enter 2D variables to CONVERT (1/1 line, end by END): }\\
198 {\tt\it ZS }  \\
199 {\tt - next 2D field or END ? }\\
200 {\tt\it END }\\
201 {\tt 2 fields (3D), and   1 fields (2D) written in CEXP.1.CSEG.001d.V5D }\\
202
203 \subsubsection{{\bf CONVLFI} program} \label{ss:convlfi}
204 The MesoNH program {\bf CONVLFI} allows conversion onto GRIB 
205 (the horizontal grid is either the native 
206 MesoNH grid (Arakawa C-grid) of the field, the MesoNH mass grid
207 (Arakawa A-grid),
208 the vertical grid is either the native K levels or pressure levels), or
209 conversion onto Vis5D (the horizontal grid is the MesoNH mass grid
210 (A-grid), the vertical grid is either the native K levels without orography,
211 altitude or pressure levels).
212
213 The conversion is done with the Meso--NH procedure {\bf prepmodel} used with 
214 the {\bf CONVLFI} program and the {\tt CONVLFI1.nam} namelist file.
215 Up to 24 FM files can be treated identically in a single prepmodel job.
216 \\
217
218 A) In the file \underline{\bf prepmodelrc}, the input and output host, directories
219 and login control variables refer to the input and output files as usual.
220 The other control variables to initialize specifically in this file are:
221 \begin{itemize}
222 \item MAINPROG=CONVLFI
223 \item LOAD\_OPT='location\_of\_v5d\_library'
224 \item OUTHOST=name\_workstation  (for example) \\
225 this allows future use of {\tt vis5d} or {\tt metview} on your local host.
226 \end{itemize}
227
228 B) In the \underline{\tt CONVLFI1.nam} namelist file, the user must indicate
229 the format type wanted, the number and type of vertical levels, 
230 the type of horizontal interpolation on a lat/lon domain 
231 as well as the name of the variables to write into the output file:
232 \begin{enumerate}
233 \item\underline{Namelist NAM\_OUTFILE}: 
234
235 \begin{center}
236 \begin{tabular} {|l|l|l|}
237 \hline
238 Fortran name & Fortran type & default value\\
239 \hline
240 \hline
241 CMNHFILE     & array of character (len=28)  & none   \\
242 COUTFILETYPE & character (len=3)   & none   \\
243 NVERB        & integer    & 5  \\
244 LAGRID       & logical    & .TRUE.  \\
245 CLEVTYPE     & character (len=1)   & 'P' if COUTFILETYPE='GRB'   \\
246              &                     & 'K' if COUTFILETYPE='V5D'   \\
247 CLEVLIST     & character (len=6)   & 'FUNCTN'   \\
248 XVLMIN       & real    & 10000.  if COUTFILETYPE='GRB'  \\
249 XVLMAX       & real    & 100000. if COUTFILETYPE='GRB'  \\
250 XVLINT       & real    & 10000.  if COUTFILETYPE='GRB'  \\
251 LLMULTI      & logical    & .TRUE. \\
252 \hline
253 \end{tabular}
254 \end{center}
255
256 \begin{itemize}
257 \item CMNHFILE: name of the input FM file (from an initialization sequence, or
258 a model simulation, or after diagnostics computation).
259 \index{CMNHFILE!\innam{NAM\_OUTFILE}}
260 \item COUTFILETYPE: type of the output file, appended
261 to CMNHFILE to generate the name of the output file.
262 \begin{itemize}
263 \item 'V5D'
264 \item 'GRB'
265 \end{itemize}
266 \index{COUTFILETYPE!\innam{NAM\_OUTFILE}}
267 \item NVERB: verbosity level 
268 \begin{itemize}
269 \item  0 for minimum of prints
270 \item  5 for intermediate level of prints
271 \item  10 for maximum of prints.
272 \end{itemize}
273 \index{NVERB!\innam{NAM\_OUTFILE}}
274 \item LAGRID: switch to interpolate fields on an Arakawa A-grid (mass grid),
275 \subitem forced to .TRUE. if Vis5D file or horizontal interpolation.
276 \index{LAGRID!\innam{NAM\_OUTFILE}}
277 \item CLEVTYPE: type of vertical levels in output file,
278 \index{CLEVTYPE!\innam{NAM\_OUTVER}}
279 \begin{itemize}
280 \item 'P' pressure levels
281 \item 'Z' z levels (only used for COUTFILETYPE='V5D')
282 \item 'K' 
283 \subitem if COUTFILETYPE='GRB': native vertical grid of Meso-NH (no
284 interpolation, height specified in GRIB message is the one of the grid 
285 without orography),
286 \subitem if COUTFILETYPE='V5D': native vertical grid of Meso-NH (fields are
287 interpolated on the levels of the lowest point of the domain).
288 \end{itemize}
289 \item CLEVLIST: how vertical levels are specified
290 \begin{itemize}
291 \item 'MANUAL' number and list of levels specified in the 1$^{st}$ free-format
292 part,
293 \item 'FUNCTN' using a linear function, with the next 3 parameters. 
294 \end{itemize}
295 \index{CLEVLIST!\innam{NAM\_OUTVER}}
296 \item XVLMIN: minimum value for the vertical grid 
297 \subitem (in m for CLEVTYPE = 'Z', in Pa for CLEVTYPE = 'P'),
298 \item XVLMAX: maximum value for the vertical grid (`'),
299 \item XVLINT: increment value for the vertical grid (`').
300 \item LLMULTI: switch to produce a multigrib file (.T.) or monogrib files (.F.),
301 only used for COUTFILETYPE='GRB' (each monogrib file name is composed with the
302 date, the variable name and the level).
303 \index{LLMULTI!\innam{NAM\_OUTFILE}}
304
305 \end{itemize}
306
307 \item\underline{Free-format part}: (number and list of vertical levels) \\
308 This part is only used if CLEVLIST='MANUAL':
309 \begin{enumerate}
310 \item first the number of vertical levels,
311 \item then the list of levels, by increasing values in m if CLEVTYPE = 'Z', or decreasing
312 values in Pa if CLEVTYPE = 'P'
313 \end{enumerate}
314
315 \item\underline{Free-format part}: (variable names)
316 This part indicates the record name of the variables of the input file to
317 write in the output file. It is specified in two parts:
318 \begin{enumerate}
319 \item between the keywords BEGIN\_3D and END\_3D: the name of the 3D fields,
320 following by their grib code if COUTFILETYPE='GRB' (separed by tabular 
321 character).
322 \item between the keywords BEGIN\_2D and END\_2D: the name of the 2D fields,
323 following by their grib code, and possibly level indicator and level value
324 if COUTFILETYPE='GRB' (separed by tabular character).
325 \end{enumerate}
326 {\bf N.B.:} do not forget the comment line after the keyword BEGIN\_3D
327 and BEGIN\_2D.
328
329
330 \end{enumerate}
331
332 \underline{C) Example of namelist file CONVLFI1.nam}
333 \begin{itemize}
334 \item
335 to convert into a Vis5d file:
336 \end{itemize}
337
338 \begin{verbatim}
339 &NAM_OUTFILE  CMNHFILE(1)='T1E20.2.09B24.002',
340               CMNHFILE(2)='T1E20.2.09B24.003',
341               COUTFILETYPE='V5D',
342               CLEVTYPE='Z', CLEVLIST='MANUAL',
343               LAGRID=T, NVERB=10 /
344 15
345 30.
346 100.
347 250.
348 500.
349 1000.
350 1500.
351 2000.
352 2500.
353 3000.
354 3500.
355 4000.
356 4500.
357 5000.
358 6000.
359 8000.
360
361 BEGIN_3D
362 #variables 3D (MesoNH field name)
363 UM
364 VM
365 WM
366 THM
367 END_3D
368 BEGIN_2D
369 #variables 2D (MesoNH field name)
370 ZS
371 END_2D
372 \end{verbatim}
373
374 \begin{itemize}
375 \item
376 to convert into a GRIB file:
377 \end{itemize}
378 \begin{verbatim}
379 &NAM_OUTFILE  CMNHFILE(1)='T1E20.2.09B24.002',
380               CMNHFILE(2)='T1E20.2.09B24.003',
381               COUTFILETYPE='GRB', 
382               CLEVTYPE='P', CLEVLIST='FUNCTN',
383               XVLMAX=100000., XVLMIN=10000., XVLINT=10000.,
384               LAGRID=T, NVERB=5 /
385
386 BEGIN_3D
387 #variables 3D (MesoNH field name, grib parameter indicator)
388 UM  33
389 VM  34
390 WM  40
391 THM 13
392 END_3D
393 BEGIN_2D
394 #variables 2D (MesoNH field name, grib parameter indicator)
395 ZS 8
396 END_2D
397 next lines are ignored
398 codes example:
399 MSLP    1  
400 ACPRR   61 
401 INPRR   59 
402 PABSM   1
403 ALT 6
404 TEMP    11
405 REHU    52
406 RVM 53
407 RCM 153
408 RRM 170
409 RIM 178
410 RSM 171
411 RGM 179
412 RHM 226
413 RARE    230
414 HHRE    231
415 VVRE    232
416 VDOP    233
417 POVOM   234
418 \end{verbatim}
419
420
421 \normalsize
422 \subsection{Short description of the program}
423 Two main tasks are performed by the program: 
424 \begin{enumerate}
425  \item \subitem After the specification of the name of the input file, a `light'
426 initialization subroutine {\tt init\_for\_convlfi.f90 } is called to initialize
427 the I/O interface, the geometry, dimensions, grids, metric coefficients, times,
428 and to read pressure field.
429  \subitem According the output grids choosen, extra arrays are allocated for 
430 interpolations. 
431 \ignore{
432 If horizontal interpolation is required, the equivalent 
433 resolution and the number of usefull points are computed by the subroutine 
434 {\tt ini2lalo.f90}.
435 }%ignore
436  \item Then fields are treated one after another: first 3D fields, then 
437 2D fields. 
438  \subitem In the case of GRIB conversion, fields are interpolated and written
439 one after another (subroutine {\tt code\_and\_write\_grib.f90 } called for each
440 horizontal level of each field). 
441  \subitem For Vis5D conversion, fields are interpolated and written
442 all together (subroutine \newline {\tt code\_and\_write\_vis5d.f90 } called at the end).
443 \end{enumerate}
444 Using a `light' initialization routine and reading fields name from standard
445 input allows the conversion program not to be dependant of a MesoNH version 
446 or program.
447
448
449 \subsection{Some tips to use Vis5D}
450 See the complete guide for using Vis5D: file README.ps in the Vis5D package.
451
452 \subsubsection{Utilities} (section 5 of README.ps)
453 \begin{itemize}
454 \item
455 {\tt v5dinfo filename}: shows summary of the v5d file: number and name of
456 the variables, size of the 3-D grid, number of time steps, vertical
457 grid definition and projection definition.
458 \item
459 {\tt v5dstats filename}: shows statistics of the v5d file:
460 minimum value, maximum value, mean value, standard deviation of
461  each variable.
462 \item
463 {\tt v5dedit filename}: edits the header of the v5d file and allows to change
464 it: variables names, variables units, times and dates, projection, vertical
465 coordinate system, low levels. \\
466 {\it Useful to set the variable's units since they are not set by the program
467  CONVLFI.}
468 \item
469 {\tt v5dappend [-var] filename1 ... targetfile}: joins v5d files together: 
470 {\it useful since the {\bf prepmodel} job generates a separate v5d file for each
471  timestep}, {\tt var} indicates list of variables to omit in the target file,
472 the dimensions of 3-D grids must be the same in each input file.
473 \end{itemize}
474
475 \subsubsection{Options} \label{ss:opt} (section 6.1 of README.ps) \\
476
477 To call Vis5D: {\tt vis5d file1 [options] file2 [options] ...} \\
478 Options can be be specified here when calling, or by pressing the {\sf DISPLAY}
479 button of the main control panel and then the 'Options' menu.
480
481 Options useful to set when calling: \\
482 {\tt [-date]} use 'dd month yy' instead of julian 'yyddd' date, \\
483 {\tt [-box x y z]} specify the aspect ratio of the 3-D box (default is 2 2 1), \\
484 {\tt [-mbs n]} override the assumed system memory size of 32 megabytes (Vis5D
485 tells you value to specify if not enough), \\
486 {\tt [-topo file]} use a topography file other than the default EARTH.TOPO
487
488
489 \subsubsection{Control panel} (section 6.2 of README.ps) \\
490 The top buttons control primary functions of Vis5D (see section
491 \ref{sss:funct}). \\
492 The middle ones control the viewing modes (see section \ref{sss:viewing}).\\
493 The bottom 2-D matrix of buttons contains physical variables on the rows, and
494 types of graphic representation on the columns. To control any type of graphic,
495 click on the button with the left mouse button. 
496 A pop-up window appears when clicking with the middle mouse button, and
497 one window to modify colors with the right button  
498 (see section \ref{sss:graph}).
499 \\
500
501 \underline{\bf Primary functions} \label{sss:funct}(section 6.3 of README.ps)
502 \begin{itemize}
503 \item{\sf SAVE PIC} to save the image in a file: first toggle the {\sf REVERSE}
504 button to reverse black and white, then toggle the {\sf SAVE PIC} button and
505 choose {\tt xwd} (X Window Dump) format. The file can be visualised with
506  {\tt xv} utility and transformed into {\tt postscript} format.
507
508 \item{\sf GRID\#s} to display the grid indices instead of latitude, longitude and
509 vertical units along the edges of the box.
510
511 \item{\sf CONT\#s, LEGENDS} to toggle on or off the isoline values, the colorbar
512 legends.
513
514 \item{\sf BOX, CLOCK} to toggle on or off the display of the box and the clock.
515
516 \item{\sf TOP, SOUTH, WEST} to set a top (or bottom), a south (or north), a west
517 (or east) view.
518 {\it Select} {\sf SOUTH} {\it to visualise 2D file.}
519
520 \item{\sf SAVE, RESTORE, SCRIPT} to save and restore isolines, colors, labels,
521 view (write and read a Tcl script).
522
523 \item{\sf UVW VARS} to specify the names of the variables to use to display wind
524 slices and trajectories, several triplets of variables can be used.
525
526 \item{\sf NEW VAR..} to duplicate variables or create new ones by specifying
527 mathematical expressions (formulas use names of existing variables, numbers,
528 arithmetic operations, functions such as $SQRT,EXP,LOG,SIN,COS,TAN,ABS,MIN,MAX$,
529 ex: horizontal wind speed, $spd=SQRT(UM*UM+VM*VM)$
530 see section 6.13 of README.ps).
531
532 \item{\sf ANIMATE} when several time steps: left mouse button: forward,
533 right button: backward, S key: slower, F key: faster.
534
535 \item{\sf STEP} when several time steps: left mouse button: one step ahead,
536  middle button: first step, right button: one step back.
537
538 \item{\sf DISPLAY} to change the number of displays, the display options
539 (see section \ref{ss:opt}), the display parameters (as with the {\tt v5dedit}
540 utility).
541
542 \end{itemize}
543
544 \underline{\bf Viewing modes} \label{sss:viewing}(section 6.4 of README.ps) \\
545 The underlined modes are the most useful (the others are much better displayed
546 with {\tt diaprog} Meso-NH graphics).
547 \begin{itemize}
548 \item\underline{\sf Normal} 
549  to rotate, zoom and translate the graphics in the 3D window.
550
551 %\item{\sf Trajectory}
552 % to create and display wind trajectories.
553 %
554 \item\underline{\sf Slice}
555  to reposition horizontal and vertical slices.
556
557 \item\underline{\sf Label}
558  to create and edit text labels in the 3D window.
559
560 \item{\sf Probe}
561  to inspect individual grid values with a cursor moving through the 3D grid.
562
563 \item{\sf Sounding}
564  to display a vertical sounding at the location of the moveable cursor.
565
566 \item{\sf Clipping}
567  to reposition the six bounding planes of the 3-D box. Select one plane (top, bottom,
568  north, south, west or east) with the middle mouse button, and reposition it
569  with the right mouse button.
570
571 \end{itemize}
572
573 \underline{\bf Types of graphic representations} \label{sss:graph}(sections 6.5 to 6.9 of README.ps) \\
574 The underlined types are the most useful (the others are much better displayed
575 with {\tt diaprog} Meso-NH graphics).
576 \begin{itemize}
577 \item\underline{\sf Isosurfaces}: 
578  A 3-D contour surface showing the volume bounding by a particular value of the
579 field (set with the left mouse button). The isosurface is either monocolor
580 or colored according to the values of another variable (right mouse button). 
581
582 \item\underline{\sf Slices}: 
583 Planar cross section (horizontally or vertically) can be moved in this mode.
584 To replace geographic coordinates by grid
585 coordinates, press the {\sf "GRID \#s"} button on the control panel.
586
587 \subitem contour line: interval can be changed
588 and min/max values specified in the pop-up window. {\tt -10 (-30,20)} will
589 plot values between -30 and 20 at intervals 10 with negative values dashed.
590 Color can be changed with the right mouse button.
591
592 \subitem colored slice: colors can be changed in the pop-up window 
593 (with the mouse buttons or arrow keys). Color table is displayed in the
594 3-D window if the {\sf "LEGEND \#s"} button is selected. 
595 %Transparency can be changed by pressing the SHIFT key while using mouse.
596 To change limits of plotted values, use the keyboard array buttons when in 
597 the variable control panel (left and right for limits in the extend of the 
598 variable values, up and down for colors inside it).
599
600 \subitem wind vector slice: (buttons {\sf Hwind1, Vwind1, Hwind2, Vwind2})
601 the scale parameter multiplies the length of vectors drawn
602 (double: 2, half: 0.5), the density parameter controls the number of vectors
603 (between zero and one, 0.5 for one vector of two, 0.25 for one of four).
604
605 \subitem wind stream slice: (buttons {\sf HStream, VStream})
606 the density parameter controls the number of streamlines
607 (between zero and two).
608
609 \item\underline{\sf Volume rendering}: {\it for powerful workstations..}
610  
611 \end{itemize}
612
613
614 \subsubsection{Advanced use}
615
616 \begin{itemize}
617 \item generate your own topography file, with the {\tt maketopo.c} program
618 in the {\tt util} directory (see 5 of README.ps).
619
620 \item Tcl language, to write script (button {\sf SCRIPT}) or
621 interactively (button {\sf INTERP..}) (see 6.16 of README.ps).
622
623 \item external analysis functions written in Fortran,
624 in {\tt userfuncs} directory (see 6.13.3 of README.ps).
625
626 \end{itemize}
627
628 \subsection{State of art}
629 The converter only runs on Linux and VPP.
630 In HP, right compilation options have to be found to use the external library...