[ncl-talk] sub: css grid and objective analysis...

[Adv]

Hi,
I have been trying to compute objective analysis of station data sets.

begin
 strll  = asciiread("stationmac.txt", -1, "string") ;
  print(strll)   ; STATION","STATION_NAME", "LAT", "LON"
  print("=====")

sloc = str_get_field(strll , 2, ",")
  lat  = tofloat( str_get_field(strll , 3, ",")  )
  lon  = tofloat( str_get_field(strll , 4, ",")  )
  print(sloc +" : "+lat+" "+lon)

 NOBS  = 65
diri="/home/dailysimms/prec-obs/01/"
z1 = asciiread(diri+"PrecJan1965-2014",-1,"float")
print(z1)
printVarSummary(z1)
;a = (/1,2,3,4,5,6,7,8/)
;391 rows by 50 columns
a0 = onedtond(z1,(/391,50/))
;print(a0)
x=dim_avg(a0)

;------------------------------------------------------------------
; Read data
;------------------------------------------------------------------

  diri = "./"
  fili = "Ecearth320x160x1961-2005.nc"
  f    = addfile(diri+fili, "r")
  Z    = f->t2(0,{40:49},{242:264})
  Z1    = f->t2(0,{40:49},{242:264})*86400
printVarSummary(Z)

;------------------------------------------------------------------
; Generate N random values from the above
;------------------------------------------------------------------
 dimz  = dimsizes( Z )
nlat  = dimz(0)
  mlon  = dimz(1)
  nx = 8
  ny = 19
  plat =  90.*fspan(-1.,1.,nx)
  plon = 180.*fspan(-1.,1.,ny)
  N     = nlat*mlon
  iZ    = generate_unique_indices( N )
 trip  = cssgrid(lat, lon, x,nlat,mlon)
dlon=0.1
  reps  = random_uniform (-dlon, dlon, N)  ; minor location perturbations

  IOBS  = iZ(0:NOBS)   ; convenience
                       ; clarity only ... create explicit variables
 rlon  = trip(0,IOBS) + reps(IOBS)
 rlat  = trip(1,IOBS) + reps(IOBS(::-1))
 rZ    = trip(2,IOBS)
;******************************************************************
; call Wrap version of obj_anal_ic
;******************************************************************

  opt    = True
  opt at timing = True

 dcrit  = (/45,35,25, 15,10 /)
  grid   = obj_anal_ic_Wrap(lon, lat, x, rlon, rlat, dcrit, opt)
printVarSummary(grid)
print(grid)
return

ERROR:

CSSGRID - error  -319, coordinates   319 and M coincide for some M  >   319
>= 1 (coordinate numbering starting at 1).
Error number -319 returned from c_cssgridd
fatal:cssgrid: Error number -319.
fatal:["Execute.c":8575]:Execute: Error occurred at or near line 45 in file
oban.ncl

Could someone help me to fix this bug?
Thank you
Adv
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20151208/feb5ad86/attachment.html