Dear all,
Could anyone give me help how to solve the problem for calculating RH from
JRA-25 data.
thanks and regards,
wei
***********************************************
p = a->pres
printVarSummary(p)
t = a1->tmp
printVarSummary(t)
q = a2->spf
printVarSummary(q)
w = q/(1-q)
w!0="time"
w!1="lev"
w!2="lat"
w!3="lon"
printVarSummary(w)
rh = relhum (t, w, conform(t,p,1))
printVarSummary(rh)
end
***********************************************
print Var Summary:
Variable: p
Type: float
Total Size: 4008960 bytes
1002240 values
Number of Dimensions: 3
Dimensions and sizes: [time | 24] x [lat | 145] x [lon | 288]
Coordinates:
time: [0.5..143.5]
lat: [90..-90]
lon: [ 0..358.75]
Number Of Attributes: 11
center : Japanese Meteorological Agency - Tokyo (RSMC)
long_name : Pressure reduced to MSL
units : Pa
_FillValue : 1e+20
level_indicator : 102
gds_grid_type : 0
parameter_table_version : 3
parameter_number : 2
forecast_time : 0
forecast_time_units : hours
initial_time : 07/26/2010 (00:00)
Variable: t
Type: float
Total Size: 92206080 bytes
23051520 values
Number of Dimensions: 4
Dimensions and sizes: [time | 24] x [lev | 23] x [lat | 145] x [lon |
288]
Coordinates:
time: [0.5..143.5]
lev: [0..1000]
lat: [90..-90]
lon: [ 0..358.75]
Number Of Attributes: 11
center : Japanese Meteorological Agency - Tokyo (RSMC)
long_name : Temperature
units : K
_FillValue : 1e+20
level_indicator : 100
gds_grid_type : 0
parameter_table_version : 3
parameter_number : 11
forecast_time : 0
forecast_time_units : hours
initial_time : 07/26/2010 (00:00)
Variable: q
Type: float
Total Size: 48107520 bytes
12026880 values
Number of Dimensions: 4
Dimensions and sizes: [time | 24] x [lev | 12] x [lat | 145] x [lon |
288]
Coordinates:
time: [0.5..143.5]
lev: [100..1000]
lat: [90..-90]
lon: [ 0..358.75]
Number Of Attributes: 11
center : Japanese Meteorological Agency - Tokyo (RSMC)
long_name : Specific humidity
units : kg kg^-1
_FillValue : 1e+20
level_indicator : 100
gds_grid_type : 0
parameter_table_version : 3
parameter_number : 51
forecast_time : 0
forecast_time_units : hours
initial_time : 07/26/2010 (00:00)
Variable: w
Type: float
Total Size: 48107520 bytes
12026880 values
Number of Dimensions: 4
Dimensions and sizes: [time | 24] x [lev | 12] x [lat | 145] x [lon |
288]
Coordinates:
Number Of Attributes: 1
_FillValue : 1e+20
fatal:conform: The array to be conformed must have the same number of
dimensions as indicated by the length of the last argument
fatal:Execute: Error occurred at or near line 24 in file rh.ncl
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Fri Jun 15 07:49:58 2012
This archive was generated by hypermail 2.1.8 : Fri Jun 15 2012 - 14:51:31 MDT