[ncl-talk] Questions for dsgrid3
Dennis Shea
shea at ucar.edu
Mon Nov 9 21:57:26 MST 2015
You should look at:
http://www.ncl.ucar.edu/Applications/ESMF.shtml
Examples 10, 14, 17 & 18
There is no need to use degrid3.ncl.
On Mon, Nov 9, 2015 at 9:40 PM, Sujeong Lim <sjlim1202 at gmail.com> wrote:
> 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
>
> --
>
>
>
> _______________________________________________
> ncl-talk mailing list
> ncl-talk at ucar.edu
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>
More information about the ncl-talk
mailing list