[ncl-talk] Sub: Interpolation
Adv
advita6 at gmail.com
Fri Dec 18 00:24:45 MST 2015
Hi,
I have been trying to do objective analysis only over the region of
longitude from -116 to -90.35 and latitude of 40.01 to 49. While plotting
the output, the values are spreaded over the global region. I don't
understand where is wrong. It's been long I am kinda stuck with this
objective analysis. I have attached here the output figure and code too.
Could someone please help me to fix this?
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)
mnlat=(min(lat))
mxlat=(max(lat))
mnlon=(min(lon))
mxlon=(max(lon))
;print(lon)
print(mnlon)
print(mxlon)
print(mnlat)
print(mxlat)
mlon = 390
dlon = -90.35 / mlon
lonn = fspan (-116, (mlon - 1) * dlon, mlon)
print(lonn)
mlat = 390
dlat = 49.0 / mlat
latt = fspan (40.01, (mlat - 1) * dlat, mlat)
print(latt)
;gcdist = gc_latlon(mnlat,mnlon, mxlat,mxlon, 300,2)
;print(gcdist)
;print (gcdist at gclat+" "+gcdist at gclon )
;return
;latt=gcdist at gclat
;lonn=gcdist at gclon
;print(min(latt))
;print(max(latt))
;print(max(lonn))
;print(min(lonn))
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/)
;390 rows by 50 columns
a0 = onedtond(z1,(/390,50/))
;print(a0)
x=dim_avg(a0)
print(x)
set_default_fillvalue("string", "")
zi = natgrid(lon, lat, x, lonn, latt)
print(zi)
zi!0="lat"
zi!1="lon"
zi&lat=latt
zi&lon=lonn
printVarSummary(zi)
print(zi)
opt = True
opt at timing = True
rscan = (/20,15,10, 5,3 /)
grid = obj_anal_ic_deprecated(lon, lat, x, lonn, latt, rscan, opt)
printVarSummary(grid)
print(grid)
print(min(grid))
print(max(grid))
grid!0="lt"
grid!1="ln"
grid<=latt
grid&ln=lonn
printVarSummary(grid)
;return
;------------------------------------------------------------------
; plot
;------------------------------------------------------------------
wks = gsn_open_wks("pdf","rdm2grid")
gsn_define_colormap(wks,"amwg")
plot = new(2,graphic)
res = True
res at gsnDraw = False ; don't draw
res at gsnFrame = False ; don't advance frame
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 gsnSpreadColors = True ; spread out color table
res at lbLabelBarOn = False ; turn off individual cb's
res at gsnCenterString = "Original grid"
plot(0) = gsn_csm_contour_map_ce(wks,zi,res)
plot(1) = gsn_csm_contour_map_ce(wks,grid,res)
;************************************************
; create panel
;************************************************
resP = True ; modify the panel plot
;resP at txString = title
resP at gsnMaximize = True ; make large
resP at gsnPanelLabelBar = True ; add common colorbar
resP at lbLabelAutoStride = True ; Let NCL decide spacing
;resP at lbLabelStride = 2 ; force every other label
;resP at lbLabelFontHeightF = 0.0125 ; make labels smaller [0.2
default]
gsn_panel(wks,plot,(/2,1/),resP) ; now draw as one plot
end
Warning message:
(0) check_for_y_lat_coord: Warning: Data either does not contain a valid
latitude coordinate array or doesn't contain one at all.
(0) A valid latitude coordinate array should have a 'units' attribute
equal to one of the following values:
(0) 'degrees_north' 'degrees-north' 'degree_north' 'degrees north'
'degrees_N' 'Degrees_north' 'degree_N' 'degreeN' 'degreesN' 'deg north'
(0) check_for_lon_coord: Warning: Data either does not contain a valid
longitude coordinate array or doesn't contain one at all.
(0) A valid longitude coordinate array should have a 'units' attribute
equal to one of the following values:
(0) 'degrees_east' 'degrees-east' 'degree_east' 'degrees east'
'degrees_E' 'Degrees_east' 'degree_E' 'degreeE' 'degreesE' 'deg east'
(0) check_for_y_lat_coord: Warning: Data either does not contain a valid
latitude coordinate array or doesn't contain one at all.
(0) A valid latitude coordinate array should have a 'units' attribute
equal to one of the following values:
(0) 'degrees_north' 'degrees-north' 'degree_north' 'degrees north'
'degrees_N' 'Degrees_north' 'degree_N' 'degreeN' 'degreesN' 'deg north'
(0) check_for_lon_coord: Warning: Data either does not contain a valid
longitude coordinate array or doesn't contain one at all.
(0) A valid longitude coordinate array should have a 'units' attribute
equal to one of the following values:
(0) 'degrees_east' 'degrees-east' 'degree_east' 'degrees east'
'degrees_E' 'Degrees_east' 'degree_E' 'degreeE' 'degreesE' 'deg east'
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20151218/b750ff44/attachment-0001.html
-------------- next part --------------
A non-text attachment was scrubbed...
Name: rdm2grid.pdf
Type: application/pdf
Size: 697435 bytes
Desc: not available
Url : http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20151218/b750ff44/attachment-0001.pdf
More information about the ncl-talk
mailing list