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

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