[ferret_users] Re: Calculating stream function from known velocity field.

[FENG HE <fenghe@xxxxxxxx>]

Ferret could not support the exchange of integral and differential. 

To verify the streamfunction solution,  
Psi = Int{y0:y}u(x,y) dy - Int{x_0:x}v(x,y_0) dx  

it is neccessary to have d/dx(Int{y0:y}u(x,y) dy) =
Int{y0:y}(du(x,y)/dx) dy, but
Ferret could not support this.

Here is a test using coads_climatology

use coads_climatology
set region/x=340:350/y=40:50/l=1
let UWND_ddb = UWND[x=@ddb]
let UWND_ddb_iin = UWND_ddb[y=@iin]
let UWND_iin= UWND[y=@iin]
let UWND_iin_ddb = UWND_iin[x=@ddb]
list UWND_ddb_iin
             VARIABLE : UWND_DDB[Y=@IIN]
             FILENAME : coads_climatology.cdf
             FILEPATH : /contrib/ferret-5.22//data/
             SUBSET   : 5 by 5 points (LONGITUDE-LATITUDE)
             TIME     : 16-JAN      06:00
               19W    17W    15W    13W    11W   
              161    162    163    164    165
 49N   / 70:  0.231 -2.136 -0.327 -1.634 -2.852
 47N   / 69: -0.342 -1.077 -0.369 -1.615 -2.513
 45N   / 68: -0.641 -0.924  0.486 -2.642 -1.725
 43N   / 67: -0.236 -0.497 -0.051 -1.943 -1.921
 41N   / 66: -0.071 -0.255  0.198 -1.417 -1.029
list UWND_iin_ddb
             VARIABLE : UWND_IIN[X=@DDB]
             FILENAME : coads_climatology.cdf
             FILEPATH : /contrib/ferret-5.22//data/
             SUBSET   : 5 by 5 points (LONGITUDE-LATITUDE)
             TIME     : 16-JAN      06:00
               19W    17W    15W    13W    11W   
              161    162    163    164    165
 49N   / 70:  0.182 -2.242 -0.318 -1.921 -3.125
 47N   / 69: -0.376 -1.137 -0.346 -1.830 -2.680
 45N   / 68: -0.651 -0.950  0.491 -2.755 -1.825
 43N   / 67: -0.238 -0.505 -0.044 -1.988 -1.953
 41N   / 66: -0.071 -0.255  0.198 -1.417 -1.029


Best regards

Feng He 
Center for Climate Research
Atmospheric and Oceanic Sciences
University of Wisconsin - Madison