source: palm/trunk/SCRIPTS/NCL/ncl_preferences.ncl @ 190

Last change on this file since 190 was 190, checked in by letzel, 14 years ago
  • NCL scripts in trunk/SCRIPTS/NCL updated to version 2.0 (.ncl_preferences replaced by ncl_preferences.ncl)
File size: 43.1 KB
Line 
1begin
2
3;*********************************************************
4; Please specify the used script by setting it to 1:
5;*********************************************************
6 
7   cross_sections = 0
8   profiles       = 0
9   timeseries     = 0
10   spectra        = 0
11   
12;*********************************************************
13
14   ;******************************************************
15   ; parameters used by CROSS_SECTIONS
16   ;******************************************************
17   
18   if (cross_sections .EQ. 1)then
19   
20      ;***************************************************
21      ; REQUIRED --file_1-- INPUT FILE
22      ;
23      ; input file produced by PALM; please declare path and name
24      ;
25      ; data type: string
26      ;
27      ; default:   "File in"
28      ;***************************************************
29      if(.not. isvar("file_1"))then
30     
31     
32         file_1 = "File in"
33         
34       
35      end if
36      ;***************************************************
37      ; REQUIRED --xyc, xzc, yzc-- SECTION of XY or XZ or YZ
38      ;
39      ; desired section has to be set to 1; e.g.: xyc=1, xzc=0 and yzc=0 for xy-section
40      ;
41      ; data type: integer
42      ;
43      ; default:   0
44      ;***************************************************
45      if(.not. isvar("xyc"))then     
46       
47         xyc = 0
48       
49      end if
50      if (.not. isvar("xzc"))then
51     
52         xzc = 0
53       
54      end if
55      if (.not. isvar("yzc"))then
56     
57         yzc = 0
58                 
59      end if             
60      ;***************************************************
61      ; OPTIONAL --format_out-- OUTPUT FORMAT
62      ;
63      ; supported file formats: x11, pdf, ps, eps, epsi or ncgm
64      ;
65      ; data type: string
66      ;
67      ; default:   "x11"
68      ;***************************************************
69      if(.not. isvar("format_out"))then 
70     
71     
72         format_out = "x11"
73         
74         
75      end if
76      ;***************************************************
77      ; OPTIONAL --file_out-- OUTPUT FILE
78      ;
79      ; please declare path and name
80      ;
81      ; data type: string
82      ;
83      ; default:   "~/test"
84      ;***************************************************
85      if(.not. isvar("file_out"))then   
86     
87     
88         file_out = "~/test_cs"
89         
90         
91      end if     
92      ;***************************************************
93      ; OPTIONAL --no_columns-- NUMBER OF ROWS
94      ;
95      ; number of plots in one row
96      ;
97      ; data type: integer
98      ;
99      ; default:   1
100      ;***************************************************
101      if(.not. isvar("no_columns"))then 
102     
103     
104         no_columns = 1
105         
106         
107      end if     
108      ;***************************************************
109      ; OPTIONAL --no_lines-- NUMBER OF COLUMNS
110      ;
111      ; number of plots in one column
112      ;
113      ; data type: integer
114      ;
115      ; default:   2
116      ;***************************************************
117      if(.not. isvar("no_lines"))then   
118     
119     
120         no_lines = 2
121         
122         
123      end if     
124      ;***************************************************
125      ; OPTIONAL --sort-- TYPE OF SORTING
126      ;
127      ; sequence of plots; sorting either by time step [="time"] or by layer/height [="layer"]
128      ;
129      ; data type: string
130      ;
131      ; default:   "layer"
132      ;***************************************************
133      if(.not. isvar("sort"))then       
134     
135     
136         sort = "layer"
137         
138         
139      end if     
140      ;***************************************************
141      ; OPTIONAL --var-- OUTPUT VARIABLES
142      ;
143      ; names of output variables; please inidcate them with one comma before and after each variable;
144      ; e.g.: var=",u,v,w," (the windcomponents u, v and w will be plotted)
145      ;
146      ; data type: string
147      ;
148      ; default:   "all" (all available variables will be plotted)
149      ;***************************************************
150      if(.not. isvar("var"))then
151     
152         
153         var = "all"
154         
155         
156      end if     
157      ;***************************************************
158      ; OPTIONAL --start_time_step-- FIRST TIME STEP
159      ;
160      ; first time step of plot; you don't have to specify a precise time value from the input file,
161      ; it will be rounded to the next existent time step
162      ;
163      ; data type: double; dimension: [hour]
164      ;
165      ; default:   -1.d (first available time step)
166      ;***************************************************
167      if(.not. isvar("start_time_step"))then
168     
169         
170         start_time_step = -1.d
171         
172         
173      end if     
174      ;***************************************************
175      ; OPTIONAL --end_time_step-- LAST TIME STEP
176      ;
177      ; last time step of plot; you don't have to specify a precise time value from the input file,
178      ; it will be rounded to the next existent time step
179      ;
180      ; data type: double; dimension: [hour]
181      ;
182      ; default:   -1.d (last available time step)
183      ;***************************************************
184      if(.not. isvar("end_time_step"))then       
185     
186     
187         end_time_step = -1.d
188         
189         
190      end if     
191      ;***************************************************
192      ; OPTIONAL --xs-- MINIMUM X-AXIS
193      ;
194      ; value for minimum of x-axis, selectable if there are NO preseted layers for x (as for XY or XZ 2D-DATA);
195      ; you don't have to specify a precise meter value from the input file, it will be rounded to the next
196      ; existent value
197      ;
198      ; data type: double; dimension: [meter]
199      ;
200      ; default:   -1.d (minimum of range)
201      ;***************************************************
202      if(.not. isvar("xs"))then     
203     
204     
205         xs = -1.d
206         
207         
208      end if     
209      ;***************************************************
210      ; OPTIONAL --xe-- MAXIMUM X-AXIS
211      ;
212      ; value for maximum of x-axis, selectable if there are NO preseted layers for x (as for XY or XZ 2D-DATA);
213      ; you don't have to specify a precise meter value from the input file, it will be rounded to the next
214      ; existent value
215      ;
216      ; data type: double; dimension: [meter]
217      ;
218      ; default:   -1.d (maximum of range)
219      ;***************************************************
220      if(.not. isvar("xe"))then     
221     
222     
223         xe = -1.d
224         
225         
226      end if     
227      ;***************************************************
228      ; OPTIONAL --ys-- MINIMUM Y-AXIS
229      ;
230      ; value for minimum of y-axis, selectable if there are NO preseted layers for y (as for XY or YZ 2D-DATA);
231      ; you don't have to specify a precise meter value from the input file, it will be rounded to the next
232      ; existent value
233      ;
234      ; data type: double; dimension: [meter]
235      ;
236      ; default:   -1.d (minimum of range)
237      ;***************************************************
238      if(.not. isvar("ys"))then     
239     
240     
241         ys = -1.d
242         
243         
244      end if     
245      ;***************************************************
246      ; OPTIONAL --ye-- MAXIMUM Y-AXIS
247      ;
248      ; value for maximum of y-axis, selectable if there are NO preseted layers for y (as for XY or YZ 2D-DATA);
249      ; you don't have to specify a precise meter value from the input file, it will be rounded to the next
250      ; existent value
251      ;
252      ; data type: double; dimension: [meter]
253      ;
254      ; default:   -1.d (maximum of range)
255      ;***************************************************
256      if(.not. isvar("ye"))then   
257     
258         
259         ye = -1.d
260         
261         
262      end if     
263      ;***************************************************
264      ; OPTIONAL --zs-- MINIMUM Z-AXIS
265      ;
266      ; value for minimum of z-axis, selectable if there are NO preseted layers for z (as for XZ or YZ 2D-DATA);
267      ; you don't have to specify a precise meter value from the input file, it will be rounded to the next
268      ; existent value
269      ;
270      ; data type: integer
271      ;
272      ; default:   -1 (minimum of range)
273      ;***************************************************
274      if(.not. isvar("zs"))then     
275     
276     
277         zs = -1
278         
279         
280      end if     
281      ;***************************************************
282      ; OPTIONAL --ze-- MAXIMUM Z-AXIS
283      ;
284      ; value for maximum of z-axis, selectable if there are NO preseted layers for z (as for XZ or YZ 2D-DATA);
285      ; you don't have to specify a precise meter value from the input file, it will be rounded to the next
286      ; existent value
287      ;
288      ; data type: integer
289      ;
290      ; default:   -1 (maximum of range)
291      ;***************************************************
292      if(.not. isvar("ze"))then
293     
294           
295         ze = -1
296         
297         
298      end if     
299      ;***************************************************
300      ; OPTIONAL --mode-- TYPE OF PLOT
301      ;
302      ; isoline plots [="Line"], contour plots [="Fill"] or both [="Both"] can be choosen
303      ;
304      ; data type: string
305      ;
306      ; default:   "Fill"
307      ;***************************************************
308      if(.not. isvar("mode"))then       
309     
310     
311         mode = "Fill"
312         
313         
314      end if     
315      ;***************************************************
316      ; OPTIONAL --fill_mode-- TYPE OF FILLING
317      ;
318      ; area fill [="AreaFill"], raster fill [="RasterFill"] or cell fill [="CellFill"] can be choosen
319      ;
320      ; data type: string
321      ;
322      ; default:   "AreaFill"
323      ;***************************************************
324      if(.not. isvar("fill_mode"))then   
325     
326     
327         fill_mode = "AreaFill"
328         
329         
330      end if     
331      ;***************************************************
332      ; OPTIONAL --shape-- ASPECT RATIO
333      ;
334      ; aspect ratio of axis is kept [=1] or  not[=0]
335      ;
336      ; data type: integer
337      ;
338      ; default:   1
339      ;***************************************************
340      if(.not. isvar("shape"))then
341     
342         
343         shape = 1
344         
345         
346      end if     
347      ;***************************************************
348      ; OPTIONAL --vector-- VECTOR PLOT
349      ;
350      ; a vector plot for one defined vector will be switched on [=1] or off [=0]
351      ;
352      ; data type: integer
353      ;
354      ; default:   0
355      ;***************************************************
356      if(.not. isvar("vector"))then
357     
358         
359         vector = 0
360         
361         
362      end if     
363      ;***************************************************
364      ; REQUIRED IF vector=1 --vec1, vec2-- VECTOR FOR VECTOR PLOT
365      ;
366      ; variables for the vector, vec1 hast to be set to component 1 and vec2 to component 2;
367      ; please inidcate them with one comma before and after each variable;
368      ; e.g.: vec1=",u," and vec2=",v," for plotting the horizontal windvector
369      ;
370      ; data type: string
371      ;
372      ; default:   vec1 = "component1", vec2 = "component2"
373      ;***************************************************
374      if(.not. isvar("vec1"))then     
375         
376         vec1 = "component1"
377         
378      end if
379      if (.not. isvar("vec2"))then
380     
381         vec2 = "component2"     
382         
383      end if     
384      ;***************************************************
385      ; OPTIONAL --plotvec-- VARIABLE FOR EXTRA VECTOR PLOT
386      ;
387      ; the vector plot can be combined with other variables; if desired inidcate the variables
388      ; with one comma before and after it; e.g. plotvec=",u,v,w," (the vector plot
389      ; will combined with the contour plots of u, v and w)
390      ;
391      ; data type: string
392      ;
393      ; default:   plotvec (no combining, the vector plot will be created seperately)
394      ;***************************************************
395      if(.not. isvar("plotvec"))then
396     
397         
398         plotvec = "plotvec"
399         
400         
401      end if     
402      ;***************************************************
403      ; OPTIONAL --ref_mag-- REFERENCE VECTOR
404      ;
405      ; length of the reference vector in a vector plot
406      ;
407      ; data type: float
408      ;
409      ; default:   0.05
410      ;***************************************************
411      if(.not. isvar("ref_mag"))then
412     
413         
414         ref_mag = 0.05
415         
416         
417      end if             
418      ;***************************************************
419     
420   end if
421
422   ;******************************************************
423   ; parameters used by PROFILES
424   ;******************************************************
425   
426   if (profiles .EQ. 1)then
427   
428      ;***************************************************
429      ; REQUIRED --file_1-- INPUT FILE
430      ;
431      ; input file produced by PALM; please declare path and name
432      ;
433      ; data type: string
434      ;
435      ; default:   "File in"
436      ;***************************************************
437      if (.not. isvar("file_1"))then
438     
439     
440         file_1 = "File in"
441         
442         
443      end if
444      ;***************************************************
445      ; REQUIRED --prof3d-- 3D OR PROFILE DATA
446      ;
447      ; must be switched on [=1] if 3D data are used or switched off [=0] if profile data are used for profiles
448      ;
449      ; data type: integer
450      ;
451      ; default:   -1
452      ;***************************************************
453      if (.not. isvar("prof3d"))then 
454     
455         
456         prof3d = -1
457         
458         
459      end if     
460      ;***************************************************
461      ; OPTIONAL --format_out-- OUTPUT FORMAT
462      ;
463      ; supported file formats: x11, pdf, ps, eps, epsi or ncgm
464      ;
465      ; data type: string
466      ;
467      ; default:   "x11"
468      ;***************************************************
469      if (.not. isvar("format_out"))then         
470     
471     
472         format_out = "x11"
473         
474         
475      end if     
476      ;***************************************************
477      ; OPTIONAL --file_out-- OUTPUT FILE
478      ;
479      ; please declare path and name
480      ;
481      ; data type: string
482      ;
483      ; default:   "~/test"
484      ;***************************************************
485      if (.not. isvar("file_out"))then 
486     
487       
488         file_out = "~/test_pr"
489         
490         
491      end if     
492      ;***************************************************
493      ; OPTIONAL --no_columns-- NUMBER OF ROWS
494      ;
495      ; number of plots in one row
496      ;
497      ; data type: integer
498      ;
499      ; default:   1
500      ;***************************************************
501      if (.not. isvar("no_columns"))then         
502     
503     
504         no_columns = 1
505         
506         
507      end if     
508      ;***************************************************
509      ; OPTIONAL --no_lines-- NUMBER OF COLUMNS
510      ;
511      ; number of plots in one column
512      ;
513      ; data type: integer
514      ;
515      ; default:   2
516      ;***************************************************
517      if (.not. isvar("no_lines"))then 
518     
519       
520         no_lines = 2
521         
522         
523      end if     
524      ;***************************************************
525      ; OPTIONAL --var-- OUTPUT VARIABLES
526      ;
527      ; names of output variables; please inidcate them with one comma before and after each variable;
528      ; e.g.: var=",u,v,w," (the windcomponents u, v and w will be plotted)
529      ;
530      ; data type: string
531      ;
532      ; default:   "all" (all available variables will be plotted)
533      ;***************************************************
534      if (.not. isvar("var"))then
535     
536         
537         var = "all"
538         
539         
540      end if   
541      ;***************************************************
542      ; OPTIONAL ---no_files-- NUMBER OF INPUT FILES
543      ;
544      ; up to 12 input files can be used; all input files must contain the same variables and dimensionalties;
545      ; identical variables will be merged in one plot;
546      ; combining or overlaying of any variables is not feasible;
547      ; an extra legend (besides the time step legend) will be drawn where you can explicitily name the items
548      ;
549      ; data type: integer
550      ;
551      ; default:   1
552      ;***************************************************
553      if (.not. isvar("no_files"))then
554     
555         
556         no_files = 1
557         
558         
559      end if     
560      ;***************************************************
561      ; OPTIONAL --name_legend_1-- 1st LEGEND ITEM
562      ;
563      ; if no_files > 1: name of 1st item in extra legend
564      ;
565      ; data type: string
566      ;
567      ; default:   "File 1"
568      ;***************************************************
569      if (.not. isvar("name_legend_1"))then
570     
571           
572         name_legend_1 = "File 1"
573         
574         
575      end if     
576      ;***************************************************
577      ; REQUIRED IF no_files > 1 --file_2-- 2nd INPUT FILE
578      ;
579      ; 2nd input file produced by PALM; please declare path and name
580      ;
581      ; data type: string
582      ;
583      ; default:   "File in"
584      ;***************************************************
585      if (.not. isvar("file_2"))then 
586     
587         
588         file_2 = "File in"
589         
590         
591      end if     
592      ;***************************************************
593      ; OPTIONAL --name_legend_2-- 2nd LEGEND ITEM
594      ;
595      ; if no_files > 1: name of 2nd item in extra legend
596      ;
597      ; data type: string
598      ;
599      ; default:   "File 2"
600      ;***************************************************
601      if (.not. isvar("name_legend_2"))then   
602     
603         
604         name_legend_2 = "File 2"
605         
606         
607      end if     
608      ;***************************************************
609      ; REQUIRED IF no_files > 2 --file_3-- 3rd INPUT FILE
610      ;
611      ; 3rd input file produced by PALM; please declare path and name
612      ;
613      ; data type: string
614      ;
615      ; default:   "File in"
616      ;***************************************************
617      if (.not. isvar("file_3"))then   
618     
619       
620         file_3 = "File in"
621         
622         
623      end if     
624      ;***************************************************
625      ; OPTIONAL --name_legend_3-- 3rd LEGEND ITEM
626      ;
627      ; if no_files > 2: name of 3rd item in extra legend
628      ;
629      ; data type: string
630      ;
631      ; default:   "File 3"
632      ;***************************************************
633      if (.not. isvar("name_legend_3"))then
634     
635           
636         name_legend_3 = "File 3"
637         
638         
639      end if     
640      ;***************************************************
641      ; REQUIRED IF no_files > 3 --file_4-- 4th INPUT FILE
642      ;
643      ; 4th input file produced by PALM; please declare path and name
644      ;
645      ; data type: string
646      ;
647      ; default:   "File in"
648      ;***************************************************
649      if (.not. isvar("file_4"))then
650     
651           
652         file_4 = "File in"
653         
654         
655      end if     
656      ;***************************************************
657      ; OPTIONAL --name_legend_4-- 4th LEGEND ITEM
658      ;
659      ; if no_files > 3: name of 4th item in extra legend
660      ;
661      ; data type: string
662      ;
663      ; default:   "File 4"
664      ;***************************************************
665      if (.not. isvar("name_legend_4"))then   
666     
667       
668         name_legend_4 = "File 4"
669         
670         
671      end if     
672      ;***************************************************
673      ; REQUIRED IF no_files > 4 --file_5-- 5th INPUT FILE
674      ;
675      ; 5th input file produced by PALM; please declare path and name
676      ;
677      ; data type: string
678      ;
679      ; default:   "File in"
680      ;***************************************************
681      if (.not. isvar("file_5"))then
682     
683           
684         file_5 = "File in"
685         
686         
687      end if     
688      ;***************************************************
689      ; OPTIONAL --name_legend_5-- 5th LEGEND ITEM
690      ;
691      ; if no_files > 4: name of 5th item in extra legend
692      ;
693      ; data type: string
694      ;
695      ; default:   "File 5"
696      ;***************************************************
697      if (.not. isvar("name_legend_5"))then
698     
699           
700         name_legend_5 = "File 5"
701         
702         
703      end if     
704      ;***************************************************
705      ; REQUIRED IF no_files > 5 --file_6-- 6th INPUT FILE
706      ;
707      ; 6th input file produced by PALM; please declare path and name
708      ;
709      ; data type: string
710      ;
711      ; default:   "File in"
712      ;***************************************************
713      if (.not. isvar("file_6"))then   
714     
715         
716         file_6 = "File in"
717         
718         
719      end if     
720      ;***************************************************
721      ; OPTIONAL --name_legend_6-- 6th LEGEND ITEM
722      ;
723      ; if no_files > 5: name of 6th item in extra legend
724      ;
725      ; data type: string
726      ;
727      ; default:   "File 6"
728      ;***************************************************
729      if (.not. isvar("name_legend_6"))then
730     
731           
732         name_legend_6 = "File 6"
733         
734         
735      end if     
736      ;***************************************************
737      ; REQUIRED IF no_files > 6 --file_7-- 7th INPUT FILE
738      ;
739      ; 7th input file produced by PALM; please declare path and name
740      ;
741      ; data type: string
742      ;
743      ; default:   "File in"
744      ;***************************************************
745      if (.not. isvar("file_7"))then
746     
747           
748         file_7 = "File in"
749         
750         
751      end if     
752      ;***************************************************
753      ; OPTIONAL --name_legend_7-- 7th LEGEND ITEM
754      ;
755      ; if no_files > 6: name of 7th item in extra legend
756      ;
757      ; data type: string
758      ;
759      ; default:   "File 7"
760      ;***************************************************
761      if (.not. isvar("name_legend_7"))then
762     
763           
764         name_legend_7 = "File 7"
765         
766         
767      end if     
768      ;***************************************************
769      ; REQUIRED IF no_files > 7 --file_8-- 8th INPUT FILE
770      ;
771      ; 8th input file produced by PALM; please declare path and name
772      ;
773      ; data type: string
774      ;
775      ; default:   "File in"
776      ;***************************************************
777      if (.not. isvar("file_8"))then
778     
779           
780         file_8 = "File in"
781         
782         
783      end if     
784      ;***************************************************
785      ; OPTIONAL --name_legend_8-- 8th LEGEND ITEM
786      ;
787      ; if no_files > 1: name of 8th item in extra legend
788      ;
789      ; data type: string
790      ;
791      ; default:   "File 8"
792      ;***************************************************
793      if (.not. isvar("name_legend_8"))then   
794     
795         
796         name_legend_8 = "File 8"
797         
798         
799      end if     
800      ;***************************************************
801      ; REQUIRED IF no_files > 8 --file_9-- 9th INPUT FILE
802      ;
803      ; 9th input file produced by PALM; please declare path and name
804      ;
805      ; data type: string
806      ;
807      ; default:   "File in"
808      ;***************************************************
809      if (.not. isvar("file_9"))then
810     
811           
812         file_9 = "File in"
813         
814         
815      end if     
816      ;***************************************************
817      ; OPTIONAL --name_legend_9-- 9th LEGEND ITEM
818      ;
819      ; if no_files > 8: name of 9th item in extra legend
820      ;
821      ; data type: string
822      ;
823      ; default:   "File 9"
824      ;***************************************************
825      if (.not. isvar("name_legend_9"))then
826     
827           
828         name_legend_9 = "File 9"
829         
830         
831      end if     
832      ;***************************************************
833      ; REQUIRED IF no_files > 9 --file_10-- 10th INPUT FILE
834      ;
835      ; 10th input file produced by PALM; please declare path and name
836      ;
837      ; data type: string
838      ;
839      ; default:   "File in"
840      ;***************************************************
841      if (.not. isvar("file_10"))then
842     
843           
844         file_10 = "File in"
845         
846         
847      end if     
848      ;***************************************************
849      ; OPTIONAL --name_legend_10-- 10th LEGEND ITEM
850      ;
851      ; if no_files > 9: name of 10th item in extra legend
852      ;
853      ; data type: string
854      ;
855      ; default:   "File 10"
856      ;***************************************************
857      if (.not. isvar("name_legend_10"))then 
858     
859         
860         name_legend_10 = "File 10"
861         
862         
863      end if     
864      ;***************************************************
865      ; REQUIRED IF no_files > 10 --file_11-- 11th INPUT FILE
866      ;
867      ; 11th input file produced by PALM; please declare path and name
868      ;
869      ; data type: string
870      ;
871      ; default:   "File in"
872      ;***************************************************
873      if (.not. isvar("file_11"))then
874     
875           
876         file_11 = "File in"
877         
878         
879      end if     
880      ;***************************************************
881      ; OPTIONAL --name_legend_11-- 11th LEGEND ITEM
882      ;
883      ; if no_files > 1: name of 11th item in extra legend
884      ;
885      ; data type: string
886      ;
887      ; default:   "File 11"
888      ;***************************************************
889      if (.not. isvar("name_legend_11"))then
890     
891           
892         name_legend_11 = "File 11"
893         
894         
895      end if     
896      ;***************************************************
897      ; REQUIRED IF no_files > 11 --file_12-- 12th INPUT FILE
898      ;
899      ; 12th input file produced by PALM; please declare path and name
900      ;
901      ; data type: string
902      ;
903      ; default:   "File in"
904      ;***************************************************
905      if (.not. isvar("file_12"))then
906     
907           
908         file_12 = "File in"
909         
910         
911      end if     
912      ;***************************************************
913      ; OPTIONAL --name_legend_12-- 12th LEGEND ITEM
914      ;
915      ; if no_files > 11: name of 12th item in extra legend
916      ;
917      ; data type: string
918      ;
919      ; default:   "File 12"
920      ;***************************************************
921      if (.not. isvar("name_legend_12"))then
922     
923           
924         name_legend_12 = "File 12"
925         
926         
927      end if     
928      ;***************************************************
929      ; OPTIONAL --start_time_step-- FIRST TIME STEP
930      ;
931      ; first time step of plot; you don't have to specify a precise time value from the input file,
932      ; it will be rounded to the next existent time step
933      ;
934      ; data type: double; dimension: [hour]
935      ;
936      ; default:   -1.d (first available time step)
937      ;***************************************************
938      if (.not. isvar("start_time_step"))then   
939     
940       
941         start_time_step = -1.d
942         
943         
944      end if     
945      ;***************************************************
946      ; OPTIONAL --end_time_step-- LAST TIME STEP
947      ;
948      ; last time step of plot; you don't have to specify a precise time value from the input file,
949      ; it will be rounded to the next existent time step
950      ;
951      ; data type: double; dimension: [hour]
952      ;
953      ; default:   -1.d (last available time step)
954      ;***************************************************
955      if (.not. isvar("end_time_step"))then
956     
957         
958         end_time_step = -1.d     
959         
960         
961      end if     
962      ;***************************************************
963      ; OPTIONAL --xs-- MINIMUM X-AXIS
964      ;
965      ; value for minimum of x-axis
966      ;
967      ; data type: float; dimension of variable
968      ;
969      ; default:   -1. (minimum of range)
970      ;***************************************************
971      if (.not. isvar("xs"))then   
972     
973         
974         xs = -1.
975         
976         
977      end if     
978      ;***************************************************
979      ; OPTIONAL --xe-- MAXIMUM X-AXIS
980      ;
981      ; value for maximum of x-axis
982      ;
983      ; data type: float; dimension of variable
984      ;
985      ; default:   -1. (maximum of range)
986      ;***************************************************
987      if (.not. isvar("xe"))then 
988     
989         
990         xe = -1.
991         
992         
993      end if     
994      ;***************************************************
995      ; OPTIONAL --min_z-- MINIMUM Z-AXIS
996      ;
997      ; value for minimum height of z-axis
998      ;
999      ; data type: double; dimension: [meter]
1000      ;
1001      ; default:   -1.d (related to the minimum of z-axis in input file)
1002      ;***************************************************
1003      if (.not. isvar("min_z"))then   
1004     
1005         
1006         min_z = -1.d
1007         
1008         
1009      end if     
1010      ;***************************************************
1011      ; OPTIONAL --max_z-- MAXIMUM Z-AXIS
1012      ;
1013      ; value for maximum height of z-axis
1014      ;
1015      ; data type: double; dimension: [meter]
1016      ;
1017      ; default:   -1.d (related to the maximum of z-axis in input file)
1018      ;***************************************************
1019      if (.not. isvar("max_z"))then   
1020     
1021         
1022         max_z = -1.d
1023         
1024         
1025      end if     
1026      ;***************************************************
1027      ; OPTIONAL --logy-- LOGARITHMIC SCALE
1028      ;
1029      ; logarithmic scale for y-axis will be switched on [=1] or off [=0]
1030      ;
1031      ; data type: integer
1032      ;
1033      ; default:   0
1034      ;***************************************************
1035      if (.not. isvar("logy"))then
1036     
1037           
1038         logy = 0
1039         
1040         
1041      end if     
1042      ;***************************************************
1043      ; OPTIONAL --norm-- NORMALISE Y-AXIS
1044      ;
1045      ; value for normalising the y-axis
1046      ;
1047      ; data type: float
1048      ;
1049      ; default:   1 (no normalising)
1050      ;***************************************************
1051      if (.not. isvar("norm"))then
1052     
1053           
1054         norm = 1
1055         
1056         
1057      end if     
1058      ;***************************************************
1059      ; OPTIONAL --over-- OVERLAYING
1060      ;
1061      ; predefined overlaying of standard variables will be switched on [=1] or off [=0];
1062      ; this feature cannot be used with statistic regions
1063      ;
1064      ; data type: integer
1065      ;
1066      ; default:   0
1067      ;***************************************************
1068      if (.not. isvar("over"))then
1069     
1070         
1071         over = 0
1072         
1073         
1074      end if     
1075      ;***************************************************
1076      ; OPTIONAL --combine-- COMBINING
1077      ;
1078      ; combining of two or three individual variables is possible and switched on [=1] or off[=0]
1079      ;
1080      ; data type: integer
1081      ;
1082      ; default:   0
1083      ;***************************************************
1084      if (.not. isvar("combine"))then 
1085     
1086         
1087         combine = 0
1088         
1089         
1090      end if     
1091      ;***************************************************
1092      ; REQUIRED IF combine=1 --number_comb-- NUMBER FOR COMBINING
1093      ;
1094      ; number of variables for combining; 2 or 3 is valid
1095      ;
1096      ; data type: integer
1097      ;
1098      ; default:   -1
1099      ;***************************************************
1100      if (.not. isvar("number_comb"))then
1101     
1102           
1103         number_comb = -1
1104         
1105         
1106      end if     
1107      ;***************************************************
1108      ; REQUIRED IF combine=1 --c_var-- VARIABLES FOR COMBINING
1109      ;
1110      ; names of variables; please inidcate them with one comma before and after each variable;
1111      ; e.g.: var=",u,v,w," (the windcomponents u, v and w will be combined)
1112      ;
1113      ; data type: string
1114      ;
1115      ; default:   "c_variables"
1116      ;***************************************************
1117      if (.not. isvar("c_var"))then
1118     
1119           
1120         c_var = "c_variables"
1121         
1122         
1123      end if     
1124      ;***************************************************
1125      ; OPTIONAL --black-- BLACK OR COLOR LINES
1126      ;
1127      ; colored [=0] or black and white [=1] plots
1128      ;
1129      ; data type: integer
1130      ;
1131      ; default:   0
1132      ;***************************************************
1133      if (.not. isvar("black"))then
1134     
1135           
1136         black = 0
1137         
1138         
1139      end if     
1140      ;***************************************************
1141      ; OPTIONAL --dash-- DASHED LINES
1142      ;
1143      ; use of different line patterns [=1] or continous lines [=0] for different time steps
1144      ;
1145      ; data type: integer
1146      ;
1147      ; default:   0
1148      ;***************************************************
1149      if (.not. isvar("dash"))then     
1150     
1151       
1152         dash = 0     
1153         
1154         
1155      end if     
1156      ;***************************************************
1157       
1158   end if
1159   
1160   ;******************************************************
1161   ; parameters used by SPECTRA
1162   ;******************************************************
1163   
1164   if (spectra .EQ. 1) then
1165   
1166      ;***************************************************
1167      ; REQUIRED --file_1-- INPUT FILE
1168      ;
1169      ; input file produced by PALM; please declare path and name
1170      ;
1171      ; data type: string
1172      ;
1173      ; default:   "File in"
1174      ;***************************************************
1175      if (.not. isvar("file_1"))then 
1176     
1177         
1178         file_1 = "File in"
1179         
1180         
1181      end if     
1182      ;***************************************************
1183      ; OPTIONAL --format_out-- OUTPUT FORMAT
1184      ;
1185      ; supported file formats: x11, pdf, ps, eps, epsi or ncgm
1186      ;
1187      ; data type: string
1188      ;
1189      ; default:   "x11"
1190      ;***************************************************
1191      if (.not. isvar("format_out"))then       
1192     
1193       
1194         format_out = "x11"
1195         
1196         
1197      end if     
1198      ;***************************************************
1199      ; OPTIONAL --file_out-- OUTPUT FILE
1200      ;
1201      ; please declare path and name
1202      ;
1203      ; data type: string
1204      ;
1205      ; default:   "~/test"
1206      ;***************************************************
1207      if (.not. isvar("file_out"))then 
1208     
1209       
1210         file_out = "~/test_sp"
1211         
1212         
1213      end if     
1214      ;***************************************************
1215      ; OPTIONAL --no_columns-- NUMBER OF ROWS
1216      ;
1217      ; number of plots in one row
1218      ;
1219      ; data type: integer
1220      ;
1221      ; default:   1
1222      ;***************************************************
1223      if (.not. isvar("no_columns"))then         
1224     
1225     
1226         no_columns = 1
1227         
1228         
1229      end if     
1230      ;***************************************************
1231      ; OPTIONAL --no_lines-- NUMBER OF COLUMNS
1232      ;
1233      ; number of plots in one column
1234      ;
1235      ; data type: integer
1236      ;
1237      ; default:   2
1238      ;***************************************************
1239      if (.not. isvar("no_lines"))then
1240     
1241         
1242         no_lines = 2
1243         
1244         
1245      end if     
1246      ;***************************************************
1247      ; OPTIONAL --var-- OUTPUT VARIABLES
1248      ;
1249      ; names of output variables; please inidcate them with one comma before and after each variable;
1250      ; e.g.: var=",u,v,w," (the windcomponents u, v and w will be plotted)
1251      ;
1252      ; data type: string
1253      ;
1254      ; default:   "all" (all available variables will be plotted)
1255      ;***************************************************
1256      if (.not. isvar("var"))then
1257     
1258         
1259         var = "all"
1260         
1261         
1262      end if     
1263      ;***************************************************
1264      ; OPTIONAL --height_level-- HEIGHT LEVELS
1265      ;
1266      ; defines if all heights or selected heights will be output; a 1D array with up to 10 elements
1267      ; from 0 to 9 must be indicated, e.g.: height_level=(/0,3,7/) (level 0, 3 and 7 as in input file
1268      ; will be output)
1269      ;
1270      ; data type: integer array
1271      ;
1272      ; default:   (/-1/) (all levels will be output)
1273      ;***************************************************
1274      if (.not. isvar("height_level"))then
1275     
1276         
1277         height_level = (/-1/)
1278         
1279         
1280      end if     
1281      ;***************************************************
1282      ; OPTIONAL --sort-- TYPE OF SORTING
1283      ;
1284      ; defines if all heights [="height"] or all timesteps [="time"] will be merged in one plot
1285      ;
1286      ; data type: string
1287      ;
1288      ; default:   "height"
1289      ;***************************************************
1290      if (.not. isvar("sort"))then
1291     
1292         
1293         sort = "height"
1294         
1295         
1296      end if     
1297      ;***************************************************
1298      ; OPTIONAL --start_time_step-- FIRST TIME STEP
1299      ;
1300      ; first time step of plot; you don't have to specify a precise time value from the input file,
1301      ; it will be rounded to the next existent time step
1302      ;
1303      ; data type: double; dimension: [hour]
1304      ;
1305      ; default:   -1.d (first available time step)
1306      ;***************************************************
1307      if (.not. isvar("start_time_step"))then
1308     
1309         
1310         start_time_step = -1.d
1311         
1312         
1313      end if     
1314      ;***************************************************
1315      ; OPTIONAL --end_time_step-- LAST TIME STEP
1316      ;
1317      ; last time step of plot; you don't have to specify a precise time value from the input file,
1318      ; it will be rounded to the next existent time step
1319      ;
1320      ; data type: double; dimension: [hour]
1321      ;
1322      ; default:   -1.d (last available time step)
1323      ;***************************************************
1324      if (.not. isvar("end_time_step"))then
1325     
1326         
1327         end_time_step = -1.d
1328         
1329         
1330      end if   
1331      ;***************************************************
1332      ; OPTIONAL --black-- BLACK OR COLOR LINES
1333      ;
1334      ; colored [=0] or black and white [=1] plots
1335      ;
1336      ; data type: integer
1337      ;
1338      ; default:   0
1339      ;***************************************************
1340      if (.not. isvar("black"))then
1341     
1342           
1343         black = 0
1344         
1345         
1346      end if     
1347      ;***************************************************
1348      ; OPTIONAL --dash-- DASHED LINES
1349      ;
1350      ; use of different line patterns [=1] or continous lines [=0]
1351      ;
1352      ; data type: integer
1353      ;
1354      ; default:   0
1355      ;***************************************************
1356      if (.not. isvar("dash"))then   
1357     
1358         
1359         dash = 0
1360         
1361         
1362      end if
1363      ;***************************************************
1364      ; OPTIONAL --logx-- LOGARITHMIC SCALE X-AXIS
1365      ;
1366      ; logarithmic scale for x-axis will be switched on [=1] or off [=0]
1367      ;
1368      ; data type: integer
1369      ;
1370      ; default:   1
1371      ;***************************************************
1372      if (.not. isvar("logx"))then   
1373     
1374         
1375         logx = 1
1376         
1377         
1378      end if     
1379      ;***************************************************
1380      ; OPTIONAL --logy-- LOGARITHMIC SCALE Y-AXIS
1381      ;
1382      ; logarithmic scale for y-axis will be switched on [=1] or off [=0]
1383      ;
1384      ; data type: integer
1385      ;
1386      ; default:   1
1387      ;***************************************************
1388      if (.not. isvar("logy"))then   
1389     
1390         
1391         logy = 1
1392         
1393         
1394      end if     
1395      ;***************************************************
1396      ; OPTIONAL --normx-- NORMALISE X-AXIS
1397      ;
1398      ; value for normalising the x-axis
1399      ;
1400      ; data type: float
1401      ;
1402      ; default:   1. (no normalising)
1403      ;***************************************************
1404      if (.not. isvar("normx"))then   
1405     
1406         
1407         normx = 1.
1408         
1409         
1410      end if     
1411      ;***************************************************
1412      ; OPTIONAL --normy-- NORMALISE Y-AXIS
1413      ;
1414      ; value for normalising the y-axis
1415      ;
1416      ; data type: float
1417      ;
1418      ; default:   1. (no normalising)
1419      ;***************************************************
1420      if (.not. isvar("normy"))then     
1421     
1422       
1423         normy = 1.
1424         
1425         
1426      end if     
1427      ;***************************************************
1428     
1429   end if
1430   
1431   
1432   ;******************************************************
1433   ; parameters used by TIMESERIES
1434   ;******************************************************
1435
1436   if (timeseries .EQ. 1) then
1437   
1438      ;***************************************************
1439      ; REQUIRED --file_1-- INPUT FILE
1440      ;
1441      ; input file produced by PALM; please declare path and name
1442      ;
1443      ; data type: string
1444      ;
1445      ; default:   "File in"
1446      ;***************************************************
1447      if (.not. isvar("file_1"))then 
1448     
1449         
1450         file_1 = "File in"
1451         
1452         
1453      end if     
1454      ;***************************************************
1455      ; OPTIONAL --format_out-- OUTPUT FORMAT
1456      ;
1457      ; supported file formats: x11, pdf, ps, eps, epsi or ncgm
1458      ;
1459      ; data type: string
1460      ;
1461      ; default:   "x11"
1462      ;***************************************************
1463      if (.not. isvar("format_out"))then         
1464     
1465     
1466         format_out = "x11"
1467         
1468         
1469      end if     
1470      ;***************************************************
1471      ; OPTIONAL --file_out-- OUTPUT FILE
1472      ;
1473      ; please declare path and name
1474      ;
1475      ; data type: string
1476      ;
1477      ; default:   "~/test"
1478      ;***************************************************
1479      if (.not. isvar("file_out"))then 
1480     
1481       
1482         file_out = "~/test_ts"
1483         
1484         
1485      end if     
1486      ;***************************************************
1487      ; OPTIONAL --no_columns-- NUMBER OF ROWS
1488      ;
1489      ; number of plots in one row
1490      ;
1491      ; data type: integer
1492      ;
1493      ; default:   1
1494      ;***************************************************
1495      if (.not. isvar("no_columns"))then       
1496     
1497       
1498         no_columns = 1
1499         
1500         
1501      end if     
1502      ;***************************************************
1503      ; OPTIONAL --no_lines-- NUMBER OF COLUMNS
1504      ;
1505      ; number of plots in one column
1506      ;
1507      ; data type: integer
1508      ;
1509      ; default:   2
1510      ;***************************************************
1511      if (.not. isvar("no_lines"))then
1512     
1513         
1514         no_lines = 2
1515         
1516         
1517      end if     
1518      ;***************************************************
1519      ; OPTIONAL --var-- OUTPUT VARIABLES
1520      ;
1521      ; names of output variables; please inidcate them with one comma before and after each variable;
1522      ; e.g.: var=",u,v,w," (the windcomponents u, v and w will be plotted)
1523      ;
1524      ; data type: string
1525      ;
1526      ; default:   "all" (all available variables will be plotted)
1527      ;***************************************************
1528      if (.not. isvar("var"))then
1529     
1530         
1531         var = "all"
1532         
1533         
1534      end if     
1535      ;***************************************************
1536      ; OPTIONAL --start_time_step-- FIRST TIME STEP
1537      ;
1538      ; first time step of plot; you don't have to specify a precise time value from the input file,
1539      ; it will be rounded to the next existent time step
1540      ;
1541      ; data type: double; dimension: [hour]
1542      ;
1543      ; default:   -1.d (first available time step)
1544      ;***************************************************
1545      if (.not. isvar("start_time_step"))then
1546     
1547         
1548         start_time_step = -1.d
1549         
1550         
1551      end if     
1552      ;***************************************************
1553      ; OPTIONAL --end_time_step-- LAST TIME STEP
1554      ;
1555      ; last time step of plot; you don't have to specify a precise time value from the input file,
1556      ; it will be rounded to the next existent time step
1557      ;
1558      ; data type: double; dimension: [hour]
1559      ;
1560      ; default:   -1.d (last available time step)
1561      ;***************************************************
1562      if (.not. isvar("end_time_step"))then
1563     
1564         
1565         end_time_step = -1.d
1566         
1567         
1568      end if     
1569      ;***************************************************
1570      ; OPTIONAL --over-- OVERLAYING
1571      ;
1572      ; predefined overlaying of standard variables will be switched on [=1] or off [=0];
1573      ; this feature cannot be used with statistic regions
1574      ;
1575      ; data type: integer
1576      ;
1577      ; default:   0
1578      ;***************************************************
1579      if (.not. isvar("over"))then
1580     
1581         
1582         over = 0
1583         
1584         
1585      end if
1586      ;***************************************************
1587
1588   end if
1589
1590         
1591;*********************************************************       
1592   
1593end
Note: See TracBrowser for help on using the repository browser.