- Timestamp:
- Mar 25, 2009 10:25:50 AM (16 years ago)
- Location:
- palm/trunk/SCRIPTS/NCL
- Files:
-
- 4 edited
-
UserGuide_for_NCLscripts.pdf (modified) (previous)
-
cross_sections.ncl (modified) (9 diffs)
-
ncl_preferences.ncl (modified) (7 diffs)
-
profiles.ncl (modified) (3 diffs)
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SCRIPTS/NCL/cross_sections.ncl
r218 r267 1143 1143 end if 1144 1144 if (xyc .EQ. 1)then 1145 cs_res@tiXAxisString = "x [ gridpoints]"1146 cs_res@tiYAxisString = "y [ gridpoints]"1145 cs_res@tiXAxisString = "x [m]" 1146 cs_res@tiYAxisString = "y [m]" 1147 1147 end if 1148 1148 if (xzc .EQ. 1)then 1149 cs_res@tiXAxisString = "x [ gridpoints]"1150 cs_res@tiYAxisString = "z [ gridpoints]"1149 cs_res@tiXAxisString = "x [m]" 1150 cs_res@tiYAxisString = "z [m]" 1151 1151 end if 1152 1152 if (yzc .EQ. 1)then 1153 cs_res@tiXAxisString = "y [ gridpoints]"1154 cs_res@tiYAxisString = "z [ gridpoints]"1153 cs_res@tiXAxisString = "y [m]" 1154 cs_res@tiYAxisString = "z [m]" 1155 1155 end if 1156 1156 end if … … 1492 1492 1493 1493 if ( sort .eq. "time" ) then 1494 if ( z_d( lo) .eq. -1)then1494 if ( z_d(zs+lo) .eq. -1)then 1495 1495 if (delta_z .EQ. -1) then 1496 1496 level = "-average" 1497 1497 else 1498 level = "=" + z_d( lo) + "m"1498 level = "=" + z_d(zs+lo) + "m" 1499 1499 end if 1500 1500 else 1501 level = "=" + z_d( lo) + "m"1501 level = "=" + z_d(zs+lo) + "m" 1502 1502 end if 1503 1503 cs_res@gsnCenterString = "t=" + decimalPlaces(t_all(li)/3600,2,True) +"h z"+level … … 1519 1519 1520 1520 if ( sort .eq. "layer" ) then 1521 if (z_d( li) .eq. -1) then1521 if (z_d(zs+li) .eq. -1) then 1522 1522 if (delta_z .EQ. -1) then 1523 1523 level = "-average" 1524 1524 else 1525 level = "=" + z_d( li) + "m"1525 level = "=" + z_d(zs+li) + "m" 1526 1526 end if 1527 1527 else 1528 level = "=" + z_d( li) + "m"1528 level = "=" + z_d(zs+li) + "m" 1529 1529 end if 1530 1530 cs_res@gsnCenterString = "t=" + decimalPlaces(t_all(lo)/3600,2,True) + "h z"+ level … … 1555 1555 1556 1556 if ( sort .eq. "time" ) then 1557 if ( y_d( lo) .eq. -1 ) then1557 if ( y_d(ys+lo) .eq. -1 ) then 1558 1558 level = "-average" 1559 1559 else 1560 level = "=" + y_d( lo) + "m"1560 level = "=" + y_d(ys+lo) + "m" 1561 1561 end if 1562 1562 cs_res@gsnCenterString = "t=" + decimalPlaces(t_all(li)/3600,2,True) + "h y"+ level … … 1578 1578 1579 1579 if ( sort .eq. "layer" ) then 1580 if ( y_d( li) .eq. -1 ) then1580 if ( y_d(ys+li) .eq. -1 ) then 1581 1581 level = "-average" 1582 1582 else 1583 level = "=" + y_d( li) + "m"1583 level = "=" + y_d(ys+li) + "m" 1584 1584 end if 1585 1585 cs_res@gsnCenterString = "t=" + decimalPlaces(t_all(lo)/3600,2,True) + "h y"+ level … … 1610 1610 1611 1611 if ( sort .eq. "time" ) then 1612 if ( x_d( lo) .eq. -1 ) then1612 if ( x_d(xs+lo) .eq. -1 ) then 1613 1613 level = "-average" 1614 1614 else 1615 level = "=" + x_d( lo) + "m"1615 level = "=" + x_d(xs+lo) + "m" 1616 1616 end if 1617 1617 cs_res@gsnCenterString = "t=" + decimalPlaces(t_all(li)/3600,2,True) + "h x"+ level … … 1633 1633 1634 1634 if ( sort .eq. "layer" ) then 1635 if ( x_d( li) .eq. -1 ) then1635 if ( x_d(xs+li) .eq. -1 ) then 1636 1636 level = "-average" 1637 1637 else 1638 level = "=" + x_d( li) + "m"1638 level = "=" + x_d(xs+li) + "m" 1639 1639 end if 1640 1640 cs_res@gsnCenterString = "t=" + decimalPlaces(t_all(lo)/3600,2,True) + "h x"+ level … … 1680 1680 end do 1681 1681 end if 1682 else 1682 else 1683 1683 if (format_out .EQ. "eps" .OR. format_out .EQ. "epsi") then 1684 1684 gsn_panel(wks_ps,plot(0:(no_time*no_layer*no_var)-1),(/no_var,no_layer*no_time/),cs_resP) … … 1687 1687 print(" ") 1688 1688 else 1689 do np = 0,no_layer*no_time*no_var-1,no_rows*no_columns 1689 do np = 0,no_layer*no_time*no_var-1,no_rows*no_columns 1690 1690 if ( np + no_rows*no_columns .gt. (no_layer*no_time*no_var)-1) then 1691 1691 gsn_panel(wks_ps, plot(np:(no_layer*no_time*no_var)-1),(/no_rows,no_columns/),cs_resP) -
palm/trunk/SCRIPTS/NCL/ncl_preferences.ncl
r236 r267 30 30 31 31 32 file_1 = "File in" 32 file_1 = "File in" 33 33 34 34 … … 677 677 678 678 679 file_1 = "File in" 679 file_1 = "File in" 680 680 681 681 … … 717 717 718 718 end if 719 ;***************************************************720 ; REQUIRED --prof3d-- 3D OR PROFILE DATA721 ;722 ; must be switched on [=1] if 3D data are used or switched off [=0] if profile data are used for profiles723 ;724 ; data type: integer725 ;726 ; default: -1727 ;***************************************************728 if (.not. isvar("prof3d"))then729 730 731 prof3d = -1732 733 734 end if735 719 ;*************************************************** 736 720 ; OPTIONAL --format_out-- OUTPUT FORMAT … … 1267 1251 1268 1252 1269 end_time_step = -1.d 1253 end_time_step = -1.d 1270 1254 1271 1255 … … 2038 2022 2039 2023 2040 file_1 = "File in" 2024 file_1 = "File in" 2041 2025 2042 2026 … … 2146 2130 2147 2131 2148 no_rows = 22132 no_rows = 1 2149 2133 2150 2134 … … 2204 2188 end_time_step = -1.d 2205 2189 2206 2190 2207 2191 end if 2208 2192 ;*************************************************** -
palm/trunk/SCRIPTS/NCL/profiles.ncl
r250 r267 73 73 start_f(0) = start_f_1 74 74 end_f(0) = end_f_1 75 76 if (prof3d .EQ. -1)then77 print(" ")78 print("Set 'prof3d' to 0 or 1")79 print(" ")80 exit81 else82 if (prof3d .NE. 0 .AND. prof3d .NE. 1)then83 print(" ")84 print("'prof3d'= "+prof3d+" is invalid; set 'prof3d' to 0 or 1")85 print(" ")86 exit87 end if88 end if89 75 90 76 if (no_files .GT. 1) then … … 372 358 373 359 vNam = getfilevarnames(f_att) 374 360 375 361 if (nof .EQ. 0)then 376 362 vNam0=vNam … … 407 393 print(" ") 408 394 end if 395 396 prof3d = 0 397 do varn = dim-1,0,1 398 if ( isStrSubset( vNam(varn), "zu_3d") .OR. isStrSubset( vNam(varn), "zw_3d")) then 399 prof3d = 1 400 break 401 end if 402 end do 409 403 410 404 ;-----------------------------------------------------below steps only for first file -> nof=0
Note: See TracChangeset
for help on using the changeset viewer.
