[ncl-talk] Dimension size mismatch

Rashed Mahmood rashidcomsis at gmail.com
Fri Aug 9 04:25:20 MDT 2019


Hi Jalal,
I think the error message is likely related to the lines:

  poly_lat = clat(nc,nr,:)
  poly_lon = clon(nc,nr,:)

Changed those to:
  poly_lat * :=* clat(nc,nr,:)
  poly_lon* :=* clon(nc,nr,:)

hope that helps.

Rashed


On Thu, Aug 8, 2019 at 6:39 PM Md. Jalal Uddin via ncl-talk <
ncl-talk at ucar.edu> wrote:

>  do nr=0,nRadius-1
>     do nc=0,nCenter-1
>        poly_lat = clat(nc,nr,:)
>        poly_lon = clon(nc,nr,:)
>
>        min_lat  = min(poly_lat)
>        max_lat  = max(poly_lat)
>        min_lon  = min(poly_lon)
>        max_lon  = max(poly_lon)
>
> ;---Extract the desired rectangle of data
>
>        PRC  := jancont({min_lat:max_lat},{min_lon:max_lon})
>
> ;---Set points that are outside of the circle of data to missing
>
>        lat2d := conform(PRC, PRC&latitude, 0)
>        lon2d := conform(PRC, PRC&longitude, 1)
>
>        latlon_circle := gc_inout(lat2d,lon2d, poly_lat,poly_lon)
>        PRC    = where(latlon_circle,  PRC, PRC at _FillValue)
> end do
> end do
>
>
> Here, when center lat changed (nc), the dimension size of PRC also changed
> that means PRC has different size. Similar to
> http://www.ncl.ucar.edu/Applications/Scripts/polyg_29.ncl
>
> I tried to mask overland for PRC data but failed because of different
> dimension sizes.
>
> How can I make PRC as the same size? Any suggestions would be appreciated.
>
> Jalal
> --
> *Md. Jalal Uddin*
> MSc in Applied Meteorology
> Nanjing University of Information, Science and Technology, China
> B.Sc. in Disaster Management (Hons.)
> Patuakhali Science and Technology University, Bangladesh.
> Cell: +8613260859092, +8801792052662
> Web: www.dmjalal90.weebly.com
> Facebook: jalal.hossen.39
> LinkedIn: https://bd.linkedin.com/in/md-jalal-uddin-80a026b0
> _______________________________________________
> ncl-talk mailing list
> ncl-talk at ucar.edu
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20190809/6fc3c7ca/attachment.html>


More information about the ncl-talk mailing list