Dear All,
I am trying to find the correlation between a station's monthly precipitation data and monthly ssts for a region that is in a lat,lon,time format. I read in my precipitation data using an readAsciiTable. However, when I tried plotting the correlation between the two I get the following with the error at the end. Attached is the script that I am using.
Can you help me please?
Regards,
Nkese.
Variable: time1
Type: double
Total Size: 13344 bytes
1668 values
Number of Dimensions: 1
Dimensions and sizes: [time | 1668]
Coordinates:
time: [16383360..17601072]
Number Of Attributes: 6
long_name : Time
bounds : time_bnds
standard_name : time
units : hours since 1-1-1 00:00:0.0
actual_range : ( 16383360, 17601072 )
delta_t : 0000-00-01 00:00:00
Variable: sstnew
Type: float
Total Size: 432345600 bytes
108086400 values
Number of Dimensions: 3
Dimensions and sizes: [lat | 180] x [lon | 360] x [time | 1668]
Coordinates:
lat: [89.5..-89.5]
lon: [-179.5..179.5]
time: [187001..200812]
Number Of Attributes: 10
description : HadISST 1.1 monthly average sea surface temperature
actual_range : ( -1.8, 34.76296 )
missing_value : -1e+30
_FillValue : -1e+30
scale_factor : 1
cell_methods : time: lat: lon: mean
add_offset : 0
units : degC
standard_name : sea_surface_temperature
long_name : Monthly 1 degree resolution SST
(0)
(0) Monthly 1 degree resolution SST: min=-1.8 max=34.763
Variable: sst_region
Type: float
Total Size: 45858336 bytes
11464584 values
Number of Dimensions: 3
Dimensions and sizes: [lat | 82] x [lon | 191] x [time | 732]
Coordinates:
lat: [40.5..-40.5]
lon: [-179.5..10.5]
time: [194601..200612]
Number Of Attributes: 10
long_name : Monthly 1 degree resolution SST
standard_name : sea_surface_temperature
units : degC
add_offset : 0
cell_methods : time: lat: lon: mean
scale_factor : 1
_FillValue : -1e+30
missing_value : -1e+30
actual_range : ( -1.8, 34.76296 )
description : HadISST 1.1 monthly average sea surface temperature
(0)
(0) Monthly 1 degree resolution SST: min=1.66987 max=32.8006
Variable: sst_anom
Type: float
Total Size: 45858336 bytes
11464584 values
Number of Dimensions: 3
Dimensions and sizes: [lat | 82] x [lon | 191] x [time | 732]
Coordinates:
lat: [40.5..-40.5]
lon: [-179.5..10.5]
time: [194601..200612]
Number Of Attributes: 12
anomaly_op_ncl : Annual Cycle Removed:function rmMonAnnCycLLT:contributed.ncl
description : HadISST 1.1 monthly average sea surface temperature
actual_range : ( -1.8, 34.76296 )
missing_value : -1e+30
_FillValue : -1e+30
scale_factor : 1
cell_methods : time: lat: lon: mean
add_offset : 0
units : degC
standard_name : sea_surface_temperature
long_name : Monthly 1 degree resolution SST
reference : function rmMonAnnCycLLT in contrib.ncl
(0)
(0) Monthly 1 degree resolution SST: min=-4.30175 max=6.04899
Variable: precipanom_data_ts
Type: float
Total Size: 2928 bytes
732 values
Number of Dimensions: 1
Dimensions and sizes: [732]
Coordinates:
Number Of Attributes: 1
_FillValue : 9.96921e+36
Variable: corr
Type: float
Total Size: 62648 bytes
15662 values
Number of Dimensions: 2
Dimensions and sizes: [82] x [191]
Coordinates:
Number Of Attributes: 1
_FillValue : -1e+30
fatal:No coordinate variable exists for dimension (lat) in variable (corr)
fatal:(lat) is not coordinate variable in variable(corr).
fatal:["Execute.c":7743]:Execute: Error occurred at or near line 82 in file hadsst1.ncl
CONFIDENTIALITY: This email (including any attachments) may contain confidential, proprietary and/or privileged information. Any duplication, copying, distribution, dissemination, transmission, disclosure or use in any manner of this email (including any attachments) without the authorisation of the sender is strictly prohibited. If you receive this email (including any attachments) in error, please notify the sender and delete this email (including any attachments) from your system. Thank you.
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
This archive was generated by hypermail 2.1.8 : Thu Feb 21 2013 - 11:26:43 MST