Changeset 513 for palm/trunk/SCRIPTS


Ignore:
Timestamp:
Mar 16, 2010 1:37:48 PM (14 years ago)
Author:
heinze
Message:

Using the NCL scripts by means of the shell script palmplot

Location:
palm/trunk/SCRIPTS
Files:
1 added
5 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SCRIPTS/NCL/.ncl.config.default

    r418 r513  
    5252      end if
    5353      ;***************************************************
    54       ; REQUIRED if more than one file of a job chain
     54      ; OPTIONAL/REQUIRED if more than one file of a job chain
    5555      ;          --start_f-- START OF CYCLIC NUMBER
    5656      ;
     
    7474      end if
    7575      ;***************************************************
    76       ; REQUIRED if more than one file of a job chain --end_f-- END OF CYCLIC NUMBER
    77       ;
    78       ; if there is more than one file for all timesteps of a job chain, declare
    79       ; the last cyclicnumber which contains timesteps you'd like to use
     76      ; OPTIONAL/REQUIRED if more than one file of a job chain
     77      ;          --end_f-- END OF CYCLIC NUMBER
     78      ;
     79      ; if there is more than one file for all timesteps of a job chain,
     80      ; declare the last cyclicnumber which contains timesteps you'd like
     81      ; to use
    8082      ;
    8183      ; data type: integer
     
    9496      end if
    9597      ;***************************************************
    96       ; REQUIRED --xyc--xzc--yzc-- SELECTION of XY or XZ or YZ CROSS SECTIONS
     98      ; OPTIONAL --xyc--xzc--yzc-- SELECTION of XY or XZ or YZ CROSS SECTIONS
    9799      ;
    98100      ; desired section has to be set to 1;
     
    217219      ; OPTIONAL --var-- OUTPUT VARIABLES
    218220      ;
    219       ; names of output variables; inidcate them with one comma before and after
     221      ; names of output variables; inidcate them with one comma before
     222      ; and after
    220223      ; each variable
    221224      ;
     
    276279      ; OPTIONAL --xs-- MINIMUM X-AXIS
    277280      ;
    278       ; value for minimum of x-axis, selectable if there are NO preseted layers
     281      ; value for minimum of x-axis, selectable if there are
     282      ; NO preseted layers
    279283      ; for x (as for XY or XZ 2D-DATA);
    280284      ; there is no need to specify a precise meter value from the input file,
     
    297301      ; OPTIONAL --xe-- MAXIMUM X-AXIS
    298302      ;
    299       ; value for maximum of x-axis, selectable if there are NO preseted layers
     303      ; value for maximum of x-axis, selectable if there are
     304      ; NO preseted layers
    300305      ; for x (as for XY or XZ 2D-DATA);
    301306      ; there is no need to specify a precise meter value from the input file,
     
    318323      ; OPTIONAL --ys-- MINIMUM Y-AXIS
    319324      ;
    320       ; value for minimum of y-axis, selectable if there are NO preseted layers
     325      ; value for minimum of y-axis, selectable if there are
     326      ; NO preseted layers
    321327      ; for y (as for XY or YZ 2D-DATA);
    322328      ; there is no need to specify a precise meter value from the input file,
     
    360366      ; OPTIONAL --zs-- MINIMUM Z-AXIS
    361367      ;
    362       ; index for minimum of z-axis, selectable if there are NO preseted layers
     368      ; index for minimum of z-axis, selectable if there are
     369      ; NO preseted layers
    363370      ; for z (as for XZ or YZ 2D-DATA);
    364       ; you cannot specify a meter value from the input file due to grid stretching 
     371      ; you cannot specify a meter value from the input file due to
     372      ; grid stretching 
    365373      ;
    366374      ; data type: integer
     
    380388      ; OPTIONAL --ze-- MAXIMUM Z-AXIS
    381389      ;
    382       ; index for maximum of z-axis, selectable if there are NO preseted layers
     390      ; index for maximum of z-axis, selectable if there are
     391      ; NO preseted layers
    383392      ; for z (as for XZ or YZ 2D-DATA);
    384       ; you cannot specify a meter value from the input file due to grid stretching 
     393      ; you cannot specify a meter value from the input file due to
     394      ; grid stretching 
    385395      ;
    386396      ; data type: integer
     
    510520      end if
    511521      ;***************************************************
    512       ; OPTIONAL --axes_explicit-- EXPLICIT SETTINGS OF X- AND Y-AXES ARE SWITCHED ON
     522      ; OPTIONAL --axes_explicit--
     523      ;          EXPLICIT SETTINGS OF X- AND Y-AXES ARE SWITCHED ON
    513524      ;
    514525      ; to use several parameters to set up the axes, 'axes_explicit'
     
    576587      ;                              NORMALISE X-, Y- and Z-AXES
    577588      ;
    578       ; normalising the axes with the given value; units can explicitly be set;
     589      ; normalising the axes with the given value;
     590      ; units can explicitly be set;
    579591      ; by default no normalising is used and units are [meter]
    580592      ;
     
    680692      ; example: plotvec = ",u,w,"
    681693      ;
    682       ; default:   "plotvec" (no combining, the vector plot is created seperately)
     694      ; default:   "plotvec" (no combining,
     695      ;                       the vector plot is created seperately)
    683696      ;***************************************************
    684697      if(.not. isvar("plotvec"))then
     
    743756      end if
    744757      ;***************************************************
    745       ; REQUIRED --start_f_1-- START OF CYCLIC NUMBER OF THE 1ST INPUT FILE
     758      ; OPTIONAL/REQUIRED --start_f_1-- if more than one file of a job chain
     759      ;                    START OF CYCLIC NUMBER OF THE 1ST INPUT FILE
    746760      ;
    747761      ; if there is more than one file for all timesteps of a job chain,
    748       ; declare the first cyclicnumber which contains timesteps you'd like to use
     762      ; declare the first cyclicnumber which contains timesteps
     763      ; you'd like to use
    749764      ;
    750765      ; data type: integer
     
    763778      end if
    764779      ;***************************************************
    765       ; REQUIRED --end_f_1-- END OF CYCLIC NUMBER OF THE 1ST INPUT FILE
     780      ; OPTIONAL/REQUIRED --end_f_1--  if more than one file of a job chain
     781      ;            END OF CYCLIC NUMBER OF THE 1ST INPUT FILE
    766782      ;
    767783      ; if there is more than one file for all timesteps of a job chain,
    768       ; declare the last cyclicnumber which contains timesteps you'd like to use
     784      ; declare the last cyclicnumber which contains timesteps
     785      ; you'd like to use
    769786      ;
    770787      ; data type: integer
     
    926943      ; for all timesteps of a job chain,
    927944      ; declare path and file name without cycle number and ".nc" (first and
    928       ; last cycle numbers are given with the parameters 'start_f' and 'end_f');
     945      ; last cycle numbers are given with the parameters
     946      ; 'start_f' and 'end_f');
    929947      ; otherwise declare path and full file name
    930948      ;
     
    946964      end if
    947965      ;***************************************************
    948       ; REQUIRED --start_f_2-- START OF CYCLIC NUMBER OF THE 2ND INPUT FILE
     966      ; OPTIONAL/REQUIRED --start_f_2-- if more than one file of a job chain
     967      ;                       START OF CYCLIC NUMBER OF THE 2ND INPUT FILE
    949968      ;
    950969      ; if there is more than one file for all timesteps of a job chain,
    951       ; declare the first cyclicnumber which contains timesteps you'd like to use
     970      ; declare the first cyclicnumber which contains timesteps
     971      ; you'd like to use
    952972      ;
    953973      ; data type: integer
     
    965985      end if
    966986      ;***************************************************
    967       ; REQUIRED --end_f_2-- END OF CYCLIC NUMBER OF THE 2ND INPUT FILE
     987      ; REQUIRED --end_f_2-- if more than one file of a job chain
     988      ;                      END OF CYCLIC NUMBER OF THE 2ND INPUT FILE
    968989      ;
    969990      ; if there is more than one file for all timesteps of a job chain,
    970       ; declare the last cyclicnumber which contains timesteps you'd like to use
     991      ; declare the last cyclicnumber which contains timesteps
     992      ; you'd like to use
    971993      ;
    972994      ; data type: integer
     
    10051027      ; for all timesteps of a job chain,
    10061028      ; declare path and file name without cycle number and ".nc" (first and
    1007       ; last cycle numbers are given with the parameters 'start_f' and 'end_f');
     1029      ; last cycle numbers are given with the parameters
     1030      ; 'start_f' and 'end_f');
    10081031      ; otherwise declare path and full file name
    10091032      ;
     
    10251048      end if
    10261049      ;***************************************************
    1027       ; REQUIRED --start_f_3-- START OF CYCLIC NUMBER OF THE 3RD INPUT FILE
     1050      ; OPTIONAL/REQUIRED --start_f_3-- if more than one file of a job chain
     1051      ;                        START OF CYCLIC NUMBER OF THE 3RD INPUT FILE
    10281052      ;
    10291053      ; if there is more than one file for all timesteps of a job chain,
    1030       ; declare the first cyclicnumber which contains timesteps you'd like to use
     1054      ; declare the first cyclicnumber which contains timesteps
     1055      ; you'd like to use
    10311056      ;
    10321057      ; data type: integer
     
    10441069      end if
    10451070      ;***************************************************
    1046       ; REQUIRED --end_f_3-- END OF CYCLIC NUMBER OF THE 3RD INPUT FILE
    1047       ;
    1048       ; if there is more than one file for all timesteps of a job chain, declare
     1071      ; OPTIONAL/REQUIRED --end_f_3-- if more than one file of a job chain
     1072      ;                    END OF CYCLIC NUMBER OF THE 3RD INPUT FILE
     1073      ;
     1074      ; if there is more than one file for all timesteps of a job chain,
     1075      ; declare
    10491076      ; the last cyclicnumber which contains timesteps you'd like to use
    10501077      ;
     
    10841111      ; all timesteps of a job chain,
    10851112      ; declare path and file name without cycle number and ".nc" (first and
    1086       ; last cycle numbers are given with the parameters 'start_f' and 'end_f');
     1113      ; last cycle numbers are given with the parameters
     1114      ; 'start_f' and 'end_f');
    10871115      ; otherwise declare path and full file name
    10881116      ;
     
    11041132      end if
    11051133      ;***************************************************
    1106       ; REQUIRED --start_f_4-- START OF CYCLIC NUMBER OF THE 4TH INPUT FILE
     1134      ; OPTIONAL/REQUIRED --start_f_4-- if more than one file of a job chain
     1135      ;                          START OF CYCLIC NUMBER OF THE 4TH INPUT FILE
    11071136      ;
    11081137      ; if there is more than one file for all timesteps of a job chain,
    1109       ; declare the first cyclicnumber which contains timesteps you'd like to use
     1138      ; declare the first cyclicnumber which contains timesteps
     1139      ; you'd like to use
    11101140      ;
    11111141      ; data type: integer
     
    11231153      end if
    11241154      ;***************************************************
    1125       ; REQUIRED --end_f_4-- END OF CYCLIC NUMBER OF THE 4TH INPUT FILE
     1155      ; OPTIONAL/REQUIRED --end_f_4-- if more than one file of a job chain
     1156      ;                   END OF CYCLIC NUMBER OF THE 4TH INPUT FILE
    11261157      ;
    11271158      ; if there is more than one file for all timesteps of a job chain,
    1128       ; declare the last cyclicnumber which contains timesteps you'd like to use
     1159      ; declare the last cyclicnumber which contains timesteps
     1160      ; you'd like to use
    11291161      ;
    11301162      ; data type: integer
     
    11631195      ; all timesteps of a job chain,
    11641196      ; declare path and file name without cycle number and ".nc" (first and
    1165       ; last cycle numbers are given with the parameters 'start_f' and 'end_f');
     1197      ; last cycle numbers are given with the parameters
     1198      ; 'start_f' and 'end_f');
    11661199      ; otherwise declare path and full file name
    11671200      ;
     
    11831216      end if
    11841217      ;***************************************************
    1185       ; REQUIRED --start_f_5-- START OF CYCLIC NUMBER OF THE 5TH INPUT FILE
     1218      ; OPTIONAL/REQUIRED --start_f_5-- if more than one file of a job chain
     1219      ;                     START OF CYCLIC NUMBER OF THE 5TH INPUT FILE
    11861220      ;
    11871221      ; if there is more than one file for all timesteps of a job chain,
    1188       ;  declare the first cyclicnumber which contains timesteps you'd like to use
     1222      ; declare the first cyclicnumber which contains timesteps
     1223      ; you'd like to use
    11891224      ;
    11901225      ; data type: integer
     
    12021237      end if
    12031238      ;***************************************************
    1204       ; REQUIRED --end_f_5-- END OF CYCLIC NUMBER OF THE 5TH INPUT FILE
     1239      ; OPTIONAL/REQUIRED --end_f_5-- if more than one file of a job chain
     1240      ;END OF CYCLIC NUMBER OF THE 5TH INPUT FILE
    12051241      ;
    12061242      ; if there is more than one file for all timesteps of a job chain,
    1207       ; declare the last cyclicnumber which contains timesteps you'd like to use
     1243      ; declare the last cyclicnumber which contains timesteps
     1244      ; you'd like to use
    12081245      ;
    12091246      ; data type: integer
     
    12421279      ; all timesteps of a job chain,
    12431280      ; declare path and file name without cycle number and ".nc" (first and
    1244       ; last cycle numbers are given with the parameters 'start_f' and 'end_f');
     1281      ; last cycle numbers are given with the parameters
     1282      ; 'start_f' and 'end_f');
    12451283      ; otherwise declare path and full file name
    12461284      ;
     
    12621300      end if
    12631301      ;***************************************************
    1264       ; REQUIRED --start_f_6-- START OF CYCLIC NUMBER OF THE 6TH INPUT FILE
     1302      ; OPTIONAL/REQUIRED --start_f_6-- if more than one file of a job chain
     1303      ;                     START OF CYCLIC NUMBER OF THE 6TH INPUT FILE
    12651304      ;
    12661305      ; if there is more than one file for all timesteps of a job chain,
    1267       ; declare the first cyclicnumber which contains timesteps you'd like to use
     1306      ; declare the first cyclicnumber which contains timesteps
     1307      ; you'd like to use
    12681308      ;
    12691309      ; data type: integer
     
    12811321      end if
    12821322      ;***************************************************
    1283       ; REQUIRED --end_f_6-- END OF CYCLIC NUMBER OF THE 6TH INPUT FILE
     1323      ; OPTIONAL/REQUIRED --end_f_6-- if more than one file of a job chain
     1324      ; END OF CYCLIC NUMBER OF THE 6TH INPUT FILE
    12841325      ;
    12851326      ; if there is more than one file for all timesteps of a job chain,
    1286       ; declare the last cyclicnumber which contains timesteps you'd like to use
     1327      ; declare the last cyclicnumber which contains timesteps
     1328      ; you'd like to use
    12871329      ;
    12881330      ; data type: integer
     
    15641606      ; OPTIONAL --over-- OVERLAYING
    15651607      ;
    1566       ; predefined overlaying of standard variables is switched on [=1] or off [=0];
     1608      ; predefined overlaying of standard variables is switched
     1609      ; on [=1] or off [=0];
    15671610      ; this feature cannot be used with statistic regions
    15681611      ;
     
    17251768      ; all timesteps of a job chain,
    17261769      ; declare path and file name without cycle number and ".nc" (first and
    1727       ; last cycle numbers are given with the parameters 'start_f' and 'end_f');
     1770      ; last cycle numbers are given with the parameters
     1771      ; 'start_f' and 'end_f');
    17281772      ; otherwise declare path and full file name
    17291773      ;
     
    17451789      end if
    17461790      ;***************************************************
    1747       ; REQUIRED if more than one file of a job chain
     1791      ; OPTIONAL/REQUIRED if more than one file of a job chain
    17481792      ;           --start_f-- START OF CYCLIC NUMBER
    17491793      ;
    17501794      ; if there is more than one file for all timesteps of a job chain,
    1751       ; declare the first cyclicnumber which contains timesteps you'd like to use
     1795      ; declare the first cyclicnumber which contains timesteps
     1796      ; you'd like to use
    17521797      ;
    17531798      ; data type: integer
     
    17661811      end if
    17671812      ;***************************************************
    1768       ; REQUIRED if more than one file of a job chain
     1813      ; OPTIONAL/REQUIRED if more than one file of a job chain
    17691814      ;          --end_f-- END OF CYCLIC NUMBER
    17701815      ;
    17711816      ; if there is more than one file for all timesteps of a job chain,
    1772       ; declare the last cyclicnumber which contains timesteps you'd like to use
     1817      ; declare the last cyclicnumber which contains timesteps
     1818      ; you'd like to use
    17731819      ;
    17741820      ; data type: integer
     
    20362082      ; data type: float
    20372083      ;
    2038       ; example: norm_x = 1
     2084      ; example: norm_x = 1.
    20392085      ;
    20402086      ; default:   1. (no normalising)
     
    20872133      ; OPTIONAL --unit_x--unit_y-- UNITS OF X-, Y-AXES
    20882134      ;
    2089       ; unit of the x-axis can be set if norm_x is unequal 1, otherwise the unit is [1/m]; by default no unit is set for the y-axis
     2135      ; unit of the x-axis can be set if norm_x is unequal 1,
     2136      ; otherwise the unit is [1/m]; by default no unit is set for the y-axis
    20902137      ;
    20912138      ; data type: string,
     
    21572204      ; all timesteps of a job chain,
    21582205      ; declare path and file name without cycle number and ".nc" (first and
    2159       ; last cycle numbers are given with the parameters 'start_f' and 'end_f');
     2206      ; last cycle numbers are given with the parameters
     2207      ; 'start_f' and 'end_f');
    21602208      ; otherwise declare path and full file name
    21612209      ;
     
    21772225      end if
    21782226      ;***************************************************
    2179       ; REQUIRED if more than one file of a job chain
     2227      ; OPTIONAL/REQUIRED if more than one file of a job chain
    21802228      ;          --start_f-- START OF CYCLIC NUMBER
    21812229      ;
    21822230      ; if there is more than one file for all timesteps of a job chain,
    2183       ; declare the first cyclicnumber which contains timesteps you'd like to use
     2231      ; declare the first cyclicnumber which contains timesteps
     2232      ; you'd like to use
    21842233      ;
    21852234      ; data type: integer
     
    21982247      end if
    21992248      ;***************************************************
    2200       ; REQUIRED if more than one file of a job chain
     2249      ; OPTIONAL/REQUIRED if more than one file of a job chain
    22012250      ;           --end_f-- END OF CYCLIC NUMBER
    22022251      ;
    22032252      ; if there is more than one file for all timesteps of a job chain,
    2204       ; declare the last cyclicnumber which contains timesteps you'd like to use
     2253      ; declare the last cyclicnumber which contains timesteps
     2254      ; you'd like to use
    22052255      ;
    22062256      ; data type: integer
     
    23532403      ; OPTIONAL --over-- OVERLAYING
    23542404      ;
    2355       ; predefined overlaying of standard variables are switched on [=1] or off [=0];
     2405      ; predefined overlaying of standard variables are switched
     2406      ; on [=1] or off [=0];
    23562407      ; this feature cannot be used with statistic regions
    23572408      ;
     
    23902441      ; OPTIONAL --norm_t-- NORMALISING X-AXIS OF TIMESERIES PLOTS (T-AXIS)
    23912442      ;
    2392       ; normalising the x-axis with the given value; the unit can explicitly set;
     2443      ; normalising the x-axis with the given value;
     2444      ; the unit can explicitly set;
    23932445      ; by default no normalising and unit is set to [s]
    23942446      ;
     
    23992451      ; default:   1.
    24002452      ;***************************************************
    2401       if(.not. isvar("normt"))then
     2453      if(.not. isvar("norm_t"))then
    24022454     
    24032455         
     
    24172469      ; default:   " "
    24182470      ;***************************************************
    2419       if(.not. isvar("unitt"))then
     2471      if(.not. isvar("unit_t"))then
    24202472     
    24212473         
  • palm/trunk/SCRIPTS/NCL/cross_sections.ncl

    r418 r513  
    8282      shape = 1
    8383   end if
    84    
     84
    8585   if (xyc .NE. 0 .AND. xyc .NE. 1)then
    8686      print(" ")
     
    13521352         data!4 = "x" 
    13531353
    1354          MinVal(varn) = min(data(varn,:,:,:,:))
    1355          MaxVal(varn) = max(data(varn,:,:,:,:))
     1354         MinVal(varn) = min(data(varn,start_time_step:end_time_step,0:(ze-zs),0:(ye-ys),0:(xe-xs)))
     1355         MaxVal(varn) = max(data(varn,start_time_step:end_time_step,0:(ze-zs),0:(ye-ys),0:(xe-xs)))
    13561356         
    13571357         unit(varn) = data_att@units
     
    16231623                 
    16241624                  if ( sort .eq. "layer" ) then
    1625                      if (z_d(zs+li) .eq. -1) then
     1625                     if (z_d(li) .eq. -1) then
    16261626                        if (delta_z .EQ. -1) then
    16271627                           level = "-average"
     
    17211721
    17221722                  if ( sort .eq. "time" ) then
    1723                      if ( x_d(xs+lo) .eq. -1 ) then
     1723                     if ( x_d(lo) .eq. -1 ) then
    17241724                        level = "-average"
    17251725                     else
     
    17441744
    17451745                  if ( sort .eq. "layer" ) then
    1746                      if ( x_d(xs+li) .eq. -1 ) then
     1746                     if ( x_d(li) .eq. -1 ) then
    17471747                        level = "-average"
    17481748                     else
  • palm/trunk/SCRIPTS/NCL/profiles.ncl

    r418 r513  
    4848   file_in_1 = new(no_files,logical)
    4949   start_f   = new(no_files,integer)
    50    end_f = new(no_files,integer)
     50   end_f     = new(no_files,integer)
    5151   
    5252   if (file_1 .EQ. "File in") then
     
    433433      end if
    434434
     435      if(comma .gt. dim)
     436          print(" ")
     437          print("The variables 'var="+var+"'" )
     438          print("do not exist on your input file;")
     439          print("be sure to have one comma before and after each variable")
     440          print(" ")
     441          exit
     442      end if
     443
    435444      indices = new((/comma/),integer)
    436445      comma   = 0
     
    15711580                  do i=1,dimz-1
    15721581                     do j=0,np-1
    1573                         temp= f[:]->$vNam(varn)$
     1582                        temp = f[:]->$vNam(varn)$
    15741583                        temp_att = f_att->$vNam(varn)$
    15751584                        data_temp = temp(ti_in(j),i,start_y:end_y,start_x:end_x)
     
    15801589                  do i=0,dimz-1
    15811590                     do j=0,np-1
    1582                         temp= f[:]->$vNam(varn)$
     1591                        temp = f[:]->$vNam(varn)$
    15831592                        temp_att = f_att->$vNam(varn)$
    15841593                        data_temp = temp(ti_in(j),i,start_y:end_y,start_x:end_x)
Note: See TracChangeset for help on using the changeset viewer.