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

Tue Dec 8 13:10:18 MST 2015

```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/"
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)

;------------------------------------------------------------------
;------------------------------------------------------------------

diri = "./"
fili = "Ecearth320x160x1961-2005.nc"
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