[ncl-talk] Questions for dsgrid3
Sujeong Lim
sjlim1202 at gmail.com
Mon Nov 9 21:40:44 MST 2015
Hello, I have a question for dsgrid3.
I want to convert my unstructured data to 1deg data. (ex. T(time, lev,
ncol) -> T(time, lev, lat, lon)
As I know, I can use ESMF_unstruct_to_0.25deg.ncl and degrid3.ncl. Is it
right?
As a practice, I tested degrid3 example.
Though I copied the example in the ncl homepage, it showed error message.
(Red line)
fatal:Number of dimensions on right hand side do not match number of
dimension in left hand side
fatal:Execute: Error occurred at or near line 56 in file degrid3_org.ncl
How can I do?
Sincerely,
Sujeong
--
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
begin
NUM = 1000
NX = 21
NY = 21
NZ = 21
RAND_MAX = 32767.0
xi = new((/NUM/), float)
yi = new((/NUM/), float)
zi = new((/NUM/), float)
u = new((/NUM/), float)
xo = new((/NX/), float)
yo = new((/NY/), float)
zo = new((/NZ/), float)
xmin = -2.0
ymin = -2.0
zmin = -2.0
xmax = 2.0
ymax = 2.0
zmax = 2.0
;
; Create random data in three space and define a function.
;
rand1 = new((/NUM/), float)
rand2 = new((/NUM/), float)
rand3 = new((/NUM/), float)
srand(1)
do i = 0,NUM - 1
rand1(i) = rand
rand2(i) = rand
rand3(i) = rand
end do
xi = xmin + (xmax - xmin) * (rand1 / RAND_MAX)
yi = ymin + (ymax - ymin) * (rand2 / RAND_MAX)
zi = zmin + (zmax - zmin) * (rand3 / RAND_MAX)
u = (xi * xi) + (yi * yi) + (zi * zi)
;
; Create the output grid.
;
ii = fspan(0, 20.0, 21)
xo = xmin + (ii / (NX - 1)) * (xmax - xmin)
yo = ymin + (ii / (NY - 1)) * (ymax - ymin)
zo = zmin + (ii / (NZ - 1)) * (zmax - zmin)
;
; Interpolate.
;
* u = dsgrid3s(xi, yi, zi, u, xo, yo, zo)*
end
--
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20151110/014f55c7/attachment.html
More information about the ncl-talk
mailing list