# Re: how to determine whether a gridpoint is in Europe base on lat & lon information?

From: Dennis Shea <shea_at_nyahnyahspammersnyahnyah>
Date: Wed, 16 Sep 2009 11:58:38 -0600

True, there is no function the will explicitly return a yes/no if you
are over Europe.
What about a country like Turkey ... part Europe/part Asia?

The question is not clear:
Do you mean the entire area [land+ocean+Mediterranean] that
encompasses Europe?
Do you mean the just land? Iceland/Svalbard/Novaya Zemlya [ISNZ]
are officially part
of Europe. Should these
be included?

Let's say that a *crude* regional approximation that excludes ISNZ but
includes
water area is sufficient. Your data array is x(time,lat,lon).

The region could be ....
latS = 36
latN = 70
lonW = -10
lonE = 60

work = x(0,:,:)
LAT = conform( work, work&lat, 0)
LON = conform( work, work&lon, 1)
eMask = new ( (/nlat,mlon/), "logical", "No_FillValue")
eMask = where(LAT.ge.latS .and. LAT.le.latN .and. \
LON.ge.lonE .and. LON.le.lonW, True,
False)
delete(work)
delete(LAT)
delete(LON)

========
If you want to exclude water [ocean/Sea] areas ... Do you have a
If so, this would have to be included in the above 'where' atatement.

If you do not have a land/sea mask, you could use NCL's

See, Example 5 at
==========

Wei Huang wrote:
> Kai,
>
> Unfortunately, there is no such function available yet.
>
> Wei Huang
> huangwei_at_ucar.edu
> VETS/CISL
> National Center for Atmospheric Research
> P.O. Box 3000 (1850 Table Mesa Dr.)
> Boulder, CO 80307-3000 USA
> (303) 497-8924
>
>
>
>
>
> On Sep 16, 2009, at 6:29 AM, Kai Zhang wrote:
>
>
>> Dear all,
>>
>> I want to calculate the total amount of emission over Europe based on
>> gridded data. Is there a way to determine whether a grid point is in
>> Europe base on lat & lon information?
>>
>> Thanks and best regards,
>>
>> Kai
>> _______________________________________________
>> ncl-talk mailing list
>> List instructions, subscriber options, unsubscribe:
>> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>>
>
> _______________________________________________
> ncl-talk mailing list
> List instructions, subscriber options, unsubscribe:
> http://mailman.ucar.edu/mailman/listinfo/ncl-talk
>

--
======================================================
Dennis J. Shea                  tel: 303-497-1361    |
P.O. Box 3000                   fax: 303-497-1333    |
Climate Analysis Section                             |
Climate & Global Dynamics Div.                       |
National Center for Atmospheric Research             |
Boulder, CO  80307                                   |
USA                        email: shea 'at' ucar.edu |
======================================================
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Wed Sep 16 2009 - 11:58:38 MDT

This archive was generated by hypermail 2.2.0 : Wed Sep 16 2009 - 12:58:13 MDT