Changeset 513 for palm/trunk/SCRIPTS/NCL
- Timestamp:
- Mar 16, 2010 1:37:48 PM (15 years ago)
- Location:
- palm/trunk/SCRIPTS/NCL
- Files:
-
- 5 edited
-
.ncl.config.default (modified) (42 diffs)
-
UserGuide_for_NCLscripts.odt (modified) (previous)
-
UserGuide_for_NCLscripts.pdf (modified) (previous)
-
cross_sections.ncl (modified) (5 diffs)
-
profiles.ncl (modified) (4 diffs)
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SCRIPTS/NCL/.ncl.config.default
r418 r513 52 52 end if 53 53 ;*************************************************** 54 ; REQUIRED if more than one file of a job chain54 ; OPTIONAL/REQUIRED if more than one file of a job chain 55 55 ; --start_f-- START OF CYCLIC NUMBER 56 56 ; … … 74 74 end if 75 75 ;*************************************************** 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 80 82 ; 81 83 ; data type: integer … … 94 96 end if 95 97 ;*************************************************** 96 ; REQUIRED--xyc--xzc--yzc-- SELECTION of XY or XZ or YZ CROSS SECTIONS98 ; OPTIONAL --xyc--xzc--yzc-- SELECTION of XY or XZ or YZ CROSS SECTIONS 97 99 ; 98 100 ; desired section has to be set to 1; … … 217 219 ; OPTIONAL --var-- OUTPUT VARIABLES 218 220 ; 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 220 223 ; each variable 221 224 ; … … 276 279 ; OPTIONAL --xs-- MINIMUM X-AXIS 277 280 ; 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 279 283 ; for x (as for XY or XZ 2D-DATA); 280 284 ; there is no need to specify a precise meter value from the input file, … … 297 301 ; OPTIONAL --xe-- MAXIMUM X-AXIS 298 302 ; 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 300 305 ; for x (as for XY or XZ 2D-DATA); 301 306 ; there is no need to specify a precise meter value from the input file, … … 318 323 ; OPTIONAL --ys-- MINIMUM Y-AXIS 319 324 ; 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 321 327 ; for y (as for XY or YZ 2D-DATA); 322 328 ; there is no need to specify a precise meter value from the input file, … … 360 366 ; OPTIONAL --zs-- MINIMUM Z-AXIS 361 367 ; 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 363 370 ; 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 365 373 ; 366 374 ; data type: integer … … 380 388 ; OPTIONAL --ze-- MAXIMUM Z-AXIS 381 389 ; 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 383 392 ; 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 385 395 ; 386 396 ; data type: integer … … 510 520 end if 511 521 ;*************************************************** 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 513 524 ; 514 525 ; to use several parameters to set up the axes, 'axes_explicit' … … 576 587 ; NORMALISE X-, Y- and Z-AXES 577 588 ; 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; 579 591 ; by default no normalising is used and units are [meter] 580 592 ; … … 680 692 ; example: plotvec = ",u,w," 681 693 ; 682 ; default: "plotvec" (no combining, the vector plot is created seperately) 694 ; default: "plotvec" (no combining, 695 ; the vector plot is created seperately) 683 696 ;*************************************************** 684 697 if(.not. isvar("plotvec"))then … … 743 756 end if 744 757 ;*************************************************** 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 746 760 ; 747 761 ; 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 749 764 ; 750 765 ; data type: integer … … 763 778 end if 764 779 ;*************************************************** 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 766 782 ; 767 783 ; 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 769 786 ; 770 787 ; data type: integer … … 926 943 ; for all timesteps of a job chain, 927 944 ; 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'); 929 947 ; otherwise declare path and full file name 930 948 ; … … 946 964 end if 947 965 ;*************************************************** 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 949 968 ; 950 969 ; 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 952 972 ; 953 973 ; data type: integer … … 965 985 end if 966 986 ;*************************************************** 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 968 989 ; 969 990 ; 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 971 993 ; 972 994 ; data type: integer … … 1005 1027 ; for all timesteps of a job chain, 1006 1028 ; 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'); 1008 1031 ; otherwise declare path and full file name 1009 1032 ; … … 1025 1048 end if 1026 1049 ;*************************************************** 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 1028 1052 ; 1029 1053 ; 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 1031 1056 ; 1032 1057 ; data type: integer … … 1044 1069 end if 1045 1070 ;*************************************************** 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 1049 1076 ; the last cyclicnumber which contains timesteps you'd like to use 1050 1077 ; … … 1084 1111 ; all timesteps of a job chain, 1085 1112 ; 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'); 1087 1115 ; otherwise declare path and full file name 1088 1116 ; … … 1104 1132 end if 1105 1133 ;*************************************************** 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 1107 1136 ; 1108 1137 ; 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 1110 1140 ; 1111 1141 ; data type: integer … … 1123 1153 end if 1124 1154 ;*************************************************** 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 1126 1157 ; 1127 1158 ; 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 1129 1161 ; 1130 1162 ; data type: integer … … 1163 1195 ; all timesteps of a job chain, 1164 1196 ; 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'); 1166 1199 ; otherwise declare path and full file name 1167 1200 ; … … 1183 1216 end if 1184 1217 ;*************************************************** 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 1186 1220 ; 1187 1221 ; 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 1189 1224 ; 1190 1225 ; data type: integer … … 1202 1237 end if 1203 1238 ;*************************************************** 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 1205 1241 ; 1206 1242 ; 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 1208 1245 ; 1209 1246 ; data type: integer … … 1242 1279 ; all timesteps of a job chain, 1243 1280 ; 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'); 1245 1283 ; otherwise declare path and full file name 1246 1284 ; … … 1262 1300 end if 1263 1301 ;*************************************************** 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 1265 1304 ; 1266 1305 ; 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 1268 1308 ; 1269 1309 ; data type: integer … … 1281 1321 end if 1282 1322 ;*************************************************** 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 1284 1325 ; 1285 1326 ; 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 1287 1329 ; 1288 1330 ; data type: integer … … 1564 1606 ; OPTIONAL --over-- OVERLAYING 1565 1607 ; 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]; 1567 1610 ; this feature cannot be used with statistic regions 1568 1611 ; … … 1725 1768 ; all timesteps of a job chain, 1726 1769 ; 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'); 1728 1772 ; otherwise declare path and full file name 1729 1773 ; … … 1745 1789 end if 1746 1790 ;*************************************************** 1747 ; REQUIRED if more than one file of a job chain1791 ; OPTIONAL/REQUIRED if more than one file of a job chain 1748 1792 ; --start_f-- START OF CYCLIC NUMBER 1749 1793 ; 1750 1794 ; 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 1752 1797 ; 1753 1798 ; data type: integer … … 1766 1811 end if 1767 1812 ;*************************************************** 1768 ; REQUIRED if more than one file of a job chain1813 ; OPTIONAL/REQUIRED if more than one file of a job chain 1769 1814 ; --end_f-- END OF CYCLIC NUMBER 1770 1815 ; 1771 1816 ; 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 1773 1819 ; 1774 1820 ; data type: integer … … 2036 2082 ; data type: float 2037 2083 ; 2038 ; example: norm_x = 1 2084 ; example: norm_x = 1. 2039 2085 ; 2040 2086 ; default: 1. (no normalising) … … 2087 2133 ; OPTIONAL --unit_x--unit_y-- UNITS OF X-, Y-AXES 2088 2134 ; 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 2090 2137 ; 2091 2138 ; data type: string, … … 2157 2204 ; all timesteps of a job chain, 2158 2205 ; 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'); 2160 2208 ; otherwise declare path and full file name 2161 2209 ; … … 2177 2225 end if 2178 2226 ;*************************************************** 2179 ; REQUIRED if more than one file of a job chain2227 ; OPTIONAL/REQUIRED if more than one file of a job chain 2180 2228 ; --start_f-- START OF CYCLIC NUMBER 2181 2229 ; 2182 2230 ; 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 2184 2233 ; 2185 2234 ; data type: integer … … 2198 2247 end if 2199 2248 ;*************************************************** 2200 ; REQUIRED if more than one file of a job chain2249 ; OPTIONAL/REQUIRED if more than one file of a job chain 2201 2250 ; --end_f-- END OF CYCLIC NUMBER 2202 2251 ; 2203 2252 ; 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 2205 2255 ; 2206 2256 ; data type: integer … … 2353 2403 ; OPTIONAL --over-- OVERLAYING 2354 2404 ; 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]; 2356 2407 ; this feature cannot be used with statistic regions 2357 2408 ; … … 2390 2441 ; OPTIONAL --norm_t-- NORMALISING X-AXIS OF TIMESERIES PLOTS (T-AXIS) 2391 2442 ; 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; 2393 2445 ; by default no normalising and unit is set to [s] 2394 2446 ; … … 2399 2451 ; default: 1. 2400 2452 ;*************************************************** 2401 if(.not. isvar("norm t"))then2453 if(.not. isvar("norm_t"))then 2402 2454 2403 2455 … … 2417 2469 ; default: " " 2418 2470 ;*************************************************** 2419 if(.not. isvar("unit t"))then2471 if(.not. isvar("unit_t"))then 2420 2472 2421 2473 -
palm/trunk/SCRIPTS/NCL/cross_sections.ncl
r418 r513 82 82 shape = 1 83 83 end if 84 84 85 85 if (xyc .NE. 0 .AND. xyc .NE. 1)then 86 86 print(" ") … … 1352 1352 data!4 = "x" 1353 1353 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))) 1356 1356 1357 1357 unit(varn) = data_att@units … … 1623 1623 1624 1624 if ( sort .eq. "layer" ) then 1625 if (z_d( zs+li) .eq. -1) then1625 if (z_d(li) .eq. -1) then 1626 1626 if (delta_z .EQ. -1) then 1627 1627 level = "-average" … … 1721 1721 1722 1722 if ( sort .eq. "time" ) then 1723 if ( x_d( xs+lo) .eq. -1 ) then1723 if ( x_d(lo) .eq. -1 ) then 1724 1724 level = "-average" 1725 1725 else … … 1744 1744 1745 1745 if ( sort .eq. "layer" ) then 1746 if ( x_d( xs+li) .eq. -1 ) then1746 if ( x_d(li) .eq. -1 ) then 1747 1747 level = "-average" 1748 1748 else -
palm/trunk/SCRIPTS/NCL/profiles.ncl
r418 r513 48 48 file_in_1 = new(no_files,logical) 49 49 start_f = new(no_files,integer) 50 end_f = new(no_files,integer)50 end_f = new(no_files,integer) 51 51 52 52 if (file_1 .EQ. "File in") then … … 433 433 end if 434 434 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 435 444 indices = new((/comma/),integer) 436 445 comma = 0 … … 1571 1580 do i=1,dimz-1 1572 1581 do j=0,np-1 1573 temp = f[:]->$vNam(varn)$1582 temp = f[:]->$vNam(varn)$ 1574 1583 temp_att = f_att->$vNam(varn)$ 1575 1584 data_temp = temp(ti_in(j),i,start_y:end_y,start_x:end_x) … … 1580 1589 do i=0,dimz-1 1581 1590 do j=0,np-1 1582 temp = f[:]->$vNam(varn)$1591 temp = f[:]->$vNam(varn)$ 1583 1592 temp_att = f_att->$vNam(varn)$ 1584 1593 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.
