Re: Plotting after using escorc

From: Dennis Shea <shea_at_nyahnyahspammersnyahnyah>
Date: Mon Feb 18 2013 - 09:29:05 MST

;correlate between sst & precip
corr = escorc(sstnew,precipanom_data_ts)
printVarSummary(corr)
corr!0 = "lat"
corr!1 = "lon"
corr&lat = sst_region&lat
corr&lon = sst_region&lon
printVarSummary(corr)

On 2/18/13 7:30 AM, Nkese Mc Shine wrote:
> 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
>
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Mon Feb 18 09:29:12 2013

This archive was generated by hypermail 2.1.8 : Thu Feb 21 2013 - 11:26:43 MST