[Thread Prev][Thread Next][Index]

[ferret_users] Sverdrup Streamfunction - indefinite integral with irregular boarders



I am trying to calculate and plot with FERRET the Sverdrup Streamfunction from the wind stress curl curl_tau

acording to the formel:

Psi = 1/beta Int_xlim(y) ^x curl(tau)dx, xlim (y) is the longitude of the boundary at each latitude.

I use the following definition (for the pacific ocean):

LET PSI = IF curl_tau[x=160E:0@iin] THEN (1/(beta*rho_ocean))*( (-1)*curl_tau[x=160E:0@iin] + curl_tau[x=160E:0@din] )/1E06 ELSE 0

fill PSI ; go land

the result is as in the attached GIF-image, but it is not sutisfing. I have the following problem:

- the islands(Hawaii, New Zeland, ...) and the Continent (Australia) disdurb the calculation (no streamfunction behind the land-boarders).
The problem is computing an indefinite integral in an area with irregular boundaries.
How can I solve this? Can someone help?

best regards
Theodor



--
--------------------------------------
Dr. Theodoros Yemenis
Center for Marine and Atmospheric Research
University of Hamburg
Institute of Oceanography
Bundesstrasse 53
D-20146 Hamburg
Germany
---------------------------------------

phone: +49 40 42838 6104
fax: +49 40 42838 7063
+49 40 36 801 801
email: yemenis@ifm.zmaw.de ---------------------------------------


GIF image


[Thread Prev][Thread Next][Index]

Dept of Commerce / NOAA / OAR / PMEL / TMAP

Contact Us | Privacy Policy | Disclaimer | Accessibility Statement