[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