[ncl-talk] Fwd: sip error

[Dawit Abebe]

Dear all,

Ignore the request for help below. I figured out the problem and fixed. 

Thanks

Dawit

> Begin forwarded message:
> 
> From: Dawit Abebe <dawit.kan at gmail.com>
> Subject: sip error
> Date: April 4, 2017 at 5:48:24 PM GMT+3
> To: ncl-talk at ucar.edu
> Cc: Mary Haley <haley at ucar.edu>
> 
> Dear nil-talk group,
> 
> I’m trying to calculate SPI but failing to produce the intended result. I pasted below the commands used and the respective returns from the interactive apprach below. The lines I need help with are those highlighted blue:
> 
> 1. on line 20) as per the examples provided in https://www.ncl.ucar.edu/Document/Functions/Built-in/dim_spi_n.shtml <https://www.ncl.ucar.edu/Document/Functions/Built-in/dim_spi_n.shtml> the use of the missing value of the variable (prc at _FillValue) returns an error. When I removed the prc at _FillValue, it  seems okay, which according to the NCL documentation on missing values (,https://www.ncl.ucar.edu/Document/Functions/Built-in/new.shtml <https://www.ncl.ucar.edu/Document/Functions/Built-in/new.shtml>) omitting the parameter from the command line leads to NCL to use the default missing value of the variable (i.e. pcr in this case). I’m just wondering if this is correct and appropriate path to the result I am expecting?    
> 
> 2. on line 26-28: the command gives me an error. I tried to figure out this but failed to fix it. 
> 
> I’m calling upon anyone to help me out on this two issues.
> 
> Thank you so much
> 
> Dawit
> 
> ——————————————————————————————
> (my_ncl) Dagemawes-MacBook-Pro:~ dawitabebe$ ncl -Q
> ncl 0> 
> ncl 1> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl" 
> ncl 2> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"   
> ncl 3> load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
> ncl 4> 
> ncl 5>    diri = "./"
> ncl 6>    fili = "pre_1901_2014.nc"
> ncl 7>    f    = addfile(diri+fili, "r")
> ncl 8> 
> ncl 9>    prc  = f->pre
> ncl 10>   printVarSummary(prc)
> 
> Variable: prc
> Type: double
> Total Size: 393984 bytes
>             49248 values
> Number of Dimensions: 3
> Dimensions and sizes:	[time | 1368] x [lat | 4] x [lon | 9]
> Coordinates: 
>             time: [ 380..41987]
>             lat: [3.75..5.25]
>             lon: [36.75..40.75]
> Number Of Attributes: 5
>   correlation_decay_distance :	 450
>   long_name :	precipitation
>   units :	mm/month
>   missing_value :	9.96920996839e+36
>   _FillValue :	9.96920996839e+36
> ncl 11> 
> ncl 12> 
> ncl 13> dimp = dimsizes(prc)
> ncl 14> ntim = dimp(0)
> ncl 15> nlat = dimp(1)
> ncl 16> mlon = dimp(2)
> ncl 17> 
> ncl 18> len  = (/3, 12, 17, 24, 29, 36/)
> ncl 19> klen = dimsizes(len)
> ncl 20> spi  = new((/klen, ntim, nlat, mlon/), float, prc at _FillValue)
> warning:New: Could not coerce missing value parameter into appropriate type, using default
> ncl 21> spi  = new((/klen, ntim, nlat, mlon/), float)                 
> ncl 22> 
> ncl 23> opt  = True
> ncl 24> opt at spi_typee = 3
> ncl 25> 
> ncl 26> do k=0,klen-1
> ncl 27>    spi(k,:,:,:) = dim_spi_n(prc, len(k), True, 0)
> ncl 28> end do
> fatal:Assignment type mismatch, right hand side can't be coerced to type of left hand side
> fatal:["Execute.c":8638]:Execute: Error occurred at or near line 27
> 
>  

-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20170405/7af653da/attachment.html