[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