[ncl-talk] unable to plot correlation

[Dennis Shea]

Interpolate so that both grids have the same dimensions.

On Tue, Feb 23, 2016 at 3:27 PM, Ipshita Majhi <ipmajhi at alaska.edu> wrote:

>
> Dear NCL,
>
> I used escorc to calculate correlation. The lat and lon dimension of the
> two variables are different so the correlation output had four dimension.
> How can i resolve it. Thank you, listed below is the code:-
>
>
> ;********************************************
> ;Reading in monthly mean and original weekly data
> ;for sce
> ;********************************************
>
> a=addfile("~/Documents/NCL_files/Snow_cover_IMS/NCL_codes/sce_mon_mean.nc","r")
> ; open output netCDF file
> b=addfile("~/Documents/NCL_files/Sea_ice/nhsce_v01r01_19661004_20140602.nc
> ","r")
>
> ;Extracting time and snow cover extent (sce) from respective data files
>
> sce=byte2flt(a->sce)
>
> sce at _FillValue = 0
>
> sce!1="lat"
> sce!2="lon"
>
> lat=b->latitude
> lon=b->longitude
>
> sce at lat2d=lat
> sce at lon2d=lon
>
> ;*********************************************
>
> ;*********************************************
> ;Reading in Sea ice data
> ;*********************************************
>
> c=addfile("~/Documents/NCL_files/Sea_ice/HadISST_ice.nc","r")
> lat_sic=c->latitude
> lon_sic=c->longitude
>
> sic=c->sic
>
> ;Extracting 1967-2012
>
> sic_1967_2012=sic(1164:1715,:,:)
>
> copy_VarAtts(sic,sic_1967_2012)
> copy_VarCoords_1(sic,sic_1967_2012)
>
> ;*******************************************
>
> ;*******************************************
> ;Detrending and reording both the datasets
> ;*******************************************
>
> jan_sce = new((/46,88,88/),float)
> jan_sic=new((/46,180,360/),float)
>
> jan_sce at _FillValue = -0.99
> jan_sic at _FillValue= -1e+30
>
> jan_sce at lat2d=lat
> jan_sce at lon2d=lon
> ;*****************************************
>
> ;*****************************************
> ;Extracting monthly data for sce
> ;*****************************************
>
> do nyr=0,551,12
>
> jan_sce(nyr/12,:,:) =sce(nyr,:,:)
> jan_sic(nyr/12,:,:)=sic_1967_2012(nyr,:,:)
>
> end do
>
> ;****************************************
>
> ;****************************************
> ;Reordering
>
> jan_reorder_sce=jan_sce(lat|:,lon|:,time|:); Dimension is (88,88,:)
> jan_reorder_sic=jan_sic(latitude|:,longitude|:,time|:); Dimension is
> (180,360,:)
>
> ;****************************************
>
> ;****************************************
> ;Detrending now
> ;****************************************
> jan_sce_dt=dtrend(jan_reorder_sce,False)
> jan_sic_dt=dtrend(jan_reorder_sic,False)
> ;****************************************
>
> ;****************************************
> ;Correlation coefficient
>
> jan_ccr = escorc(jan_sce_dt,jan_sic_dt)
> print(jan_ccr)
>
>
> _______________________________________________
> ncl-talk mailing list
> ncl-talk at ucar.edu
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20160223/b0f6ecce/attachment.html