Dear Ncl users, I want to plot pentad bases LH for the two months May and June over one plot. When I got the printVarSummary(u) it indicates only fist five days but how I confirm the other pentad all together? Could anyone check the code given below?.
Thanks in advance,
kim
========================================
========================================
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,:,:))
u6 = short2flt(f->lhtfl(153:157,:,:))
u7 = short2flt(f->lhtfl(158:162,:,:))
u8 = short2flt(f->lhtfl(163:167,:,:))
u9 = short2flt(f->lhtfl(168:172,:,:))
u10 = short2flt(f->lhtfl(173:177,:,:))
u11 = short2flt(f->lhtfl(178:182,:,:))
u = u0
u = u1
u = u2
u = u3
u = u4
u = u5
u = u6
u = u7
u = u8
u = u9
u = u10
u = u11
u = (/u0+u1+u2+u3+u4+u5+u6+u7+u8+u9+u10+u11/)
printVarSummary(u)
time = u&time
printVarSummary(time)
**************************
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: [17527176..17527272]
lat: [88.542..-88.542]
lon: [ 0..358.125]
Number Of Attributes: 15
missing_value : 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
_FillValue : 32766
=============
Variable: time
Type: double
Total Size: 40 bytes
5 values
Number of Dimensions: 1
Dimensions and sizes: [time | 5]
Coordinates:
Number Of Attributes: 5
avg_period : 0000-00-01 00:00:00
delta_t : 0000-00-01 00:00:00
actual_range : ( 17522904, 17531664 )
long_name : Time
units : hours since 1-1-1 00:00:0.0
(0) 17527176
(1) 17527200
(2) 17527224
(3) 17527248
(4) 17527272
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Fri Mar 5 07:27:01 2010
This archive was generated by hypermail 2.1.8 : Fri Mar 12 2010 - 09:11:56 MST