[ncl-talk] trend_query

[Dennis Shea]

I would recommend:

http://www.ncl.ucar.edu/Document/Functions/Built-in/regCoef_n.shtml

This is the most recent version. As noted in the documentation:
Note: with NCL V6.2.1 or later, you can use *regCoef_n*
<http://www.ncl.ucar.edu/Document/Functions/Built-in/regCoef_n.shtml> to
avoid having to reorder the arrays first:

   rc   = *regCoef_n*
<http://www.ncl.ucar.edu/Document/Functions/Built-in/regCoef_n.shtml>(x,
y, 0, 0)

===

****PLease See: Example 2***

rc   = *regCoef_n*(time, ts, 0, 0) ; *rc(nlat,mlon)*

copy_VarCoords(ts(0,:,),rc)
printVarSummary(rc)

=====
Using the original regCoef rather than regCoef_n may require the arrays to
be reordered.

http://www.ncl.ucar.edu/Document/Functions/Built-in/regCoef-1.shtml

****See: Examples 2 and 3***

   rc   = *regCoef*(x, y(lat|:,lon|:,time|:) )

If *y* has coordinate variables
<http://www.ncl.ucar.edu/Document/Language/cv.shtml> these may readily be
assigned via NCL syntax:

   rc!0   = "lat"    ; name dimensions
   rc!1   = "lon"
   rc&lat = y&lat    ; assign coordinate values to named dimensions
   rc&lon = y&lon




On Mon, Jan 2, 2017 at 3:01 AM, Kunal Bali <kunal.bali9 at gmail.com> wrote:

> Dear NCL users,
>
> I have the data with time steps > 1000
>
> I want to plot the spatial map of trend analysis of the given dataset.
> So I tried the script given below for getting the desired plot but
> couldn't successed.
>
> the given error is
>
>
> *fatal:regCoef: The rightmost dimension of x must be at least
> 2fatal:["Execute.c":8567]:Execute: Error occurred at or near line 32 in
> file trend.ncl*
>
> Could anyone please help me out
>
>
>
>    tmp  = a->SWTNTCLR
>    ts   = tmp(time|:,lat|:,lon|:)               ; reorder variable
>
>    delete(tmp)                                  ; no longer needed
>
> ;************************************************
> ; create x and calculate the regression coefficients (slopes, trends)
> ;************************************************
>    time         = ts&time                       ; days since 1850-01-01
>    rc           = regCoef(time,ts)
>
>    rc at long_name = "regression coefficient (trend)"
>    rc at units     = ts at units+"/day"
> ;   copy_VarCoords(ts(:,:,0), rc)                ; copy lat,lon coords
>
>
> ;************************************************
> ; plotting parameters
> ;************************************************
>    wks  = gsn_open_wks("ps" ,"regress")
>    gsn_define_colormap(wks,"ViBlGrWhYeOrRe")    ; choose colormap
>
>    res                       = True
>    res at gsnMaximize           = True             ; make large
>
>    res at cnFillOn              = True             ; turn on color
>    res at cnLinesOn             = False            ; turn off contour lines
>    res at cnLineLabelsOn        = False            ; turn off contour line
> labels
>  ;;res at cnFillMode            = "RasterFill"
>
> ;   res at cnLevelSelectionMode  = "ManualLevels"   ; set manual contour
> levels
> ;   res at cnMinLevelValF        =  -1.00           ; set min contour level
> ;   res at cnMaxLevelValF        =   1.00           ; set max contour level
>    res at cnLevelSpacingF       =   0.10           ; set contour interval
> ;
>    res at mpFillOn              = False            ; turn off default
> background gray
>  ;  res at mpCenterLonF          = 210
>
>    res at tiMainString          = fili
>    plot = gsn_csm_contour_map_ce(wks,rc,res)
>   end
>
>
>
>
>
>
>
>
> Regards
> Kunal Bali
>
>
>
>
>
> _______________________________________________
> 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/20170102/1b1e9508/attachment.html