NCL 缺测值填平

在画位势高度的时候发现等压截取会出现缺测，参考官网案例，修正了下，再次记录下代码。

  geop_plane   = wrf_user_interp_level( gmph,p,850.,False)/9.8

  guess     = 1                ; use zonal means
  is_cyclic = True             ; cyclic [global]
  nscan     = 1500             ; usually much less than this
  eps       = 1.e-2            ; variable dependent
  relc      = 0.6              ; relaxation coefficient
  opt       = 0                ; not used

  poisson_grid_fill( geop_plane, is_cyclic, guess, nscan, eps, relc, opt)

原结果：

填平修正后：