Re: Analysis over Irregular Region

I guess there was no online response to the original question.
That is unusual. Usually every question gets answered.

---
The mask must be created manually.
Assuming that the grids have 2-dimensional lat/lon
arrays (say XLAT and XLON, eg lambert-conformal)
I think you will have todo something like:
     regMask = new ((/nlat,mlon/), integer)
Then set each point to a mask value:
0=NW, 1=N, 2=A, 3=NO, 4=SO
Since A looks like the largest region, I'd first set
      regMask = 2
Then use something like the "where" function
http://www.ncl.ucar.edu/Document/Functions/Built-in/where.shtml
to set large groups of grid points to the other regions,
       regmask = where ( xlat.ge.47.25 .and. xlon.le.8.5, 0, regmask)
       regmask = where ( xlat.ge.48.   .and.
                        (xlon.ge.9.0 .and. xlon.le.14) ,  1, regmask)
       etc
Without some external infomation, I doubt there is anyway
that you will avoid this type of approach. Probably even
setting individual grid points.
I would save to a netCDF file.
      regMask_at_long_name = "Region Mask"
      regMask_at_units     = "0=NW, 1=N, 2=A, 3=NO, 4=SO"
      regMask!0 = " ... "    ; use same name as xlat
      regMask!1 = " ... "
      system("/bin/rm -f REGMASK.nc")   ; remove any previous instance
      fout = addfile ("REGMASK.nc", "c")
      fout->rmask = regMask
IMRAN NADEEM wrote:
> Dear All,
> 
>     I want to do some analysis over irregular regions similar to shown 
> in the attached picture.How can I create mask for such regions on an RCM 
> grid. In Talk archives, I found a similar question asked in the past, 
> but I couldn't find reply to that question.
> 
> http://www.ncl.ucar.edu/Support/talk_archives/2005/0360.html
> 
> Any suggestion is welcome
> Best Regards
> Nadeem
> 
> ----------------
> Imran Nadeem
> PhD Student
> Institute of Meteorology
> Univ. of Natural Resources and Applied Life Sciences (BOKU)
> 
> Peter-Jordan Strasse 82
> 1190 Vienna, Austria
> 
> 
> ------------------------------------------------------------------------
> 
> 
> ------------------------------------------------------------------------
> 
> _______________________________________________
> ncl-talk mailing list
> ncl-talk_at_ucar.edu
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
_______________________________________________
ncl-talk mailing list
ncl-talk_at_ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk