Ignore:
Timestamp:
Mar 25, 2009 10:25:50 AM (16 years ago)
Author:
weinreis
Message:

NCL scripts profiles cross_sections ncl_preference and UserGuide? modified

File:
1 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SCRIPTS/NCL/cross_sections.ncl

    r218 r267  
    11431143      end if
    11441144      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]"
    11471147      end if
    11481148      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]"
    11511151      end if
    11521152      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]"
    11551155      end if
    11561156   end if
     
    14921492                 
    14931493                  if ( sort .eq. "time" ) then
    1494                      if ( z_d(lo) .eq. -1)then
     1494                     if ( z_d(zs+lo) .eq. -1)then
    14951495                        if (delta_z .EQ. -1) then
    14961496                           level = "-average"
    14971497                        else
    1498                            level = "=" + z_d(lo) + "m"
     1498                           level = "=" + z_d(zs+lo) + "m"
    14991499                        end if
    15001500                     else
    1501                         level = "=" + z_d(lo) + "m"
     1501                        level = "=" + z_d(zs+lo) + "m"
    15021502                     end if
    15031503                     cs_res@gsnCenterString = "t=" + decimalPlaces(t_all(li)/3600,2,True) +"h  z"+level             
     
    15191519                 
    15201520                  if ( sort .eq. "layer" ) then
    1521                      if (z_d(li) .eq. -1) then
     1521                     if (z_d(zs+li) .eq. -1) then
    15221522                        if (delta_z .EQ. -1) then
    15231523                           level = "-average"
    15241524                        else
    1525                            level = "=" + z_d(li) + "m"
     1525                           level = "=" + z_d(zs+li) + "m"
    15261526                        end if
    15271527                     else
    1528                         level = "=" + z_d(li) + "m"
     1528                        level = "=" + z_d(zs+li) + "m"
    15291529                     end if
    15301530                     cs_res@gsnCenterString = "t=" + decimalPlaces(t_all(lo)/3600,2,True) + "h  z"+ level               
     
    15551555               
    15561556                  if ( sort .eq. "time" ) then
    1557                      if ( y_d(lo) .eq. -1 ) then
     1557                     if ( y_d(ys+lo) .eq. -1 ) then
    15581558                        level = "-average"
    15591559                     else
    1560                         level = "=" + y_d(lo) + "m"
     1560                        level = "=" + y_d(ys+lo) + "m"
    15611561                     end if
    15621562                     cs_res@gsnCenterString = "t=" + decimalPlaces(t_all(li)/3600,2,True) + "h  y"+ level
     
    15781578
    15791579                  if ( sort .eq. "layer" ) then
    1580                      if ( y_d(li) .eq. -1 ) then
     1580                     if ( y_d(ys+li) .eq. -1 ) then
    15811581                        level = "-average"
    15821582                     else
    1583                         level = "=" + y_d(li) + "m"
     1583                        level = "=" + y_d(ys+li) + "m"
    15841584                     end if
    15851585                     cs_res@gsnCenterString = "t=" + decimalPlaces(t_all(lo)/3600,2,True) + "h  y"+ level
     
    16101610
    16111611                  if ( sort .eq. "time" ) then
    1612                      if ( x_d(lo) .eq. -1 ) then
     1612                     if ( x_d(xs+lo) .eq. -1 ) then
    16131613                        level = "-average"
    16141614                     else
    1615                         level = "=" + x_d(lo) + "m"
     1615                        level = "=" + x_d(xs+lo) + "m"
    16161616                     end if
    16171617                     cs_res@gsnCenterString = "t=" + decimalPlaces(t_all(li)/3600,2,True) + "h  x"+ level
     
    16331633
    16341634                  if ( sort .eq. "layer" ) then
    1635                      if ( x_d(li) .eq. -1 ) then
     1635                     if ( x_d(xs+li) .eq. -1 ) then
    16361636                        level = "-average"
    16371637                     else
    1638                         level = "=" + x_d(li) + "m"
     1638                        level = "=" + x_d(xs+li) + "m"
    16391639                     end if
    16401640                     cs_res@gsnCenterString = "t=" + decimalPlaces(t_all(lo)/3600,2,True) + "h  x"+ level
     
    16801680         end do
    16811681      end if
    1682    else        
     1682   else       
    16831683      if (format_out .EQ. "eps" .OR. format_out .EQ. "epsi") then
    16841684         gsn_panel(wks_ps,plot(0:(no_time*no_layer*no_var)-1),(/no_var,no_layer*no_time/),cs_resP)
     
    16871687         print(" ")
    16881688      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 
    16901690            if ( np + no_rows*no_columns .gt. (no_layer*no_time*no_var)-1) then
    16911691               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.