[ferret_users] bug in scat2gridlaplace_xy

[Lev Tarasov <lev@xxxxxxxxxxxxxxxxxxxxxxxxxx>]

Hi;
 	when I use scat2gridlaplace_xy on a file with datapoints over 
europe, it offset points east of the 0 latitude by about 35 degrees in an 
eastward direction. Is there a newer version of this function that works
properly? Or can somebody at least give me the exact offset value and I
could stick in a brute force adjustment.

thanks,
Lev

The relevant script is:

FILE/SKIP=0/VAR="xlnb,ylt,wgt,id"  RSLweightmap10gNBb.EU
let xln = if xlnb lt 0. then xlnb+360. else xlnb
LET wgtl = scat2gridlaplace_xy(xln,ylt, wgt, x[gx=xlonef], y[gy=ylatef], 2., 1)


yes? show axis xlonef
 name       axis              # pts   start                end
 XLONEF    LONGITUDE          133mr   12.5W                119.5E(479.5)
   Axis span (to cell edges) = 133 (modulo length = 360)
yes? show axis ylatef
 name       axis              # pts   start                end
 YLATEF    LATITUDE            73 r   47.25N               83.25N
   Axis span (to cell edges) = 36.5


        NOAA/PMEL TMAP
        FERRET v5.81
        Linux(g77) 2.4.21-32 - 07/14/05
        10-May-06 14:19

snippet of RSLweightmap10gNBb.EU
---------------------------------
  -4.06 52.47 0.00 6044
  -5.07 51.66 0.00 6045
  -3.77 51.58 0.00 6046
  -3.02 51.35 0.00 6047
   0.83 51.70 0.00 6048
   0.31 51.46 0.00 6049
   1.02 51.07 0.00 6050
   0.33 50.84 0.00 6051
  -1.35 50.80 0.00 6052
  -2.54 50.62 0.00 6053
  -3.75 50.39 0.00 6054
  -5.48 50.13 0.00 6055
  -5.25 56.07 0.03 6056
   2.98 55.02 0.02 6057
   3.50 52.50 0.00 6058
  36.15 69.10 1.34 7001
  29.40 70.09 1.17 7002
  19.45 70.15 0.22 7003
  19.00 69.80 0.49 7004
---------------------------------


Lev Tarasov -   Dept of Physics, University of Toronto,
                60 St. George St., Toronto, Ontario, CANADA, M5S 1A7
                Tel (519)-821-3555  Fax (416)-978-8905
                email: lev.tarasoff@utoronto.ca
		http://www.atmosp.physics.utoronto.ca/people/lev/lev.html