help

Dear Ncl users, I want to plot pentad bases LH for the month of May. When I got the printVarSummary (u) it indicates only fist pentad but I don't know how to confirm the other pentad? Could anyone check my code where I did the mistake?
Thanks in advance,
kim
**********************
 ====================
===================
begin
  f = addfile ("lhtfl.sfc.gauss.2000.nc","r")
  u0 = short2flt(f->lhtfl(120:124,:,:))
  u1 = short2flt(f->lhtfl(125:129,:,:))
  u2 = short2flt(f->lhtfl(133:137,:,:))
  u3 = short2flt(f->lhtfl(138:142,:,:))
  u4 = short2flt(f->lhtfl(143:147,:,:))
  u5 = short2flt(f->lhtfl(148:152,:,:))
  u = u0
; u = u1
; u = u2
; u = u3
; u = u4
; u = u5
   u = (u0+u1+u2+u3+u4+u5)
   printVarSummary(u)
   time = u&time
   utc  = ut_calendar(time,0)    ; yr, mo, dy, hr, mn, sc
   mons = floattoint(utc(:, 1))  ; Values from 1 to 12
   days = floattoint(utc(:, 2))
   uRegion = wgt_areaave_Wrap(u(:,{0:40},{30:120}),1.0,1.0,0)
 
   wks   = gsn_open_wks ("x11","LH")
   res                   = True
   res@tmLabelAutoStride = True
   res@xyLineColors      =(/"red"/)
   res@xyLineThicknesses = 2.0
   res@tiYAxisString        = "W/m^2"
   plot  = gsn_csm_xy(wks,time,uRegion,res)
   end
**************
Variable: u
Type: float
Total Size: 360960 bytes
            90240 values
Number of Dimensions: 3
Dimensions and sizes:   [time | 5] x [lat | 94] x [lon | 192]
Coordinates:
            time: [17525784..17525880]
            lat: [88.542..-88.542]
            lon: [ 0..358.125]
Number Of Attributes: 15
  missing_value :       32766
  _FillValue :  32766
  parent_stat : Individual Obs
I
  statistic :   Mean
M
  level_desc :  Surface
0
  dataset :     NCEP Reanalysis Daily Averages
AJ
  var_desc :    Latent Heat Net Flux
CT
  GRIB_name :   LHTFL
  GRIB_id :     121
  least_significant_digit :     0
  precision :   1
  units :       W/m^2
  actual_range :        ( -218.1, 1140.3 )
  valid_range : ( -2200, 2200 )
  long_name :   mean Daily Latent Heat Net Flux at surface

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Sun Mar 7 21:50:59 2010