Re: Analysis over Irregular Region

From: Dennis Shea <shea_at_nyahnyahspammersnyahnyah>
Date: Wed, 10 Sep 2008 09:01:28 -0600

Also,

One approach to plotting the masks can be seen at:

http://www.ncl.ucar.edu/Applications/Scripts/eosdis_5.ncl

The data are completely different but the plot
code will be similar.

Another example approach can be seen for example 5 at

http://www.ncl.ucar.edu/Applications/mask.shtml

Dennis Shea wrote:
> 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
_______________________________________________
ncl-talk mailing list
ncl-talk_at_ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Wed Sep 10 2008 - 09:01:28 MDT

This archive was generated by hypermail 2.2.0 : Fri Sep 12 2008 - 21:42:04 MDT