Hi, I'm trying to solve a little oceanographic problem. I wanna plot the stream lines of a system composed by the Sverdrup equation (equation regulating water motion) and the continuity equation (conservation of mass); and all this in spherical coordinates, so in a shell DR around the Earth (but keeping only the radial and longitudinal coordinates, which in the end is just a slice of the spherical shell). I attach also the code. 2 omega sin (phi) is the Coriolis parameter. w is the radial water velocity and v is the longitudinal one.

Any idea?

thanks a lot in advance for any suggestion! Lorenzo

Attachments: