Changeset 267 for palm/trunk/SCRIPTS/NCL/cross_sections.ncl
- Timestamp:
- Mar 25, 2009 10:25:50 AM (16 years ago)
- File:
-
- 1 edited
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)
Note: See TracChangeset
for help on using the changeset viewer.