search

NOAA-PMEL / Ferret

The Ferret program from NOAA/PMEL
https://ferret.pmel.noaa.gov/Ferret/
The Unlicense
55 stars 20 forks source link

wrong contour lines with Curvilinear data when units of X coord not Longitude #960

Closed karlmsmith closed 6 years ago

karlmsmith commented 6 years ago

Reported by @AnsleyManke on 18 Sep 2009 17:22 UTC Reported by Patrick McGuire 

let xb = x[x=0:290:10]
let yb = y[y=0:290:10]
let xb2 = xb + 0*yb
let zb = I[x=0:290:10] * J[y=0:290:10]

contour/hlimits=-100:400:50 xb2+zb,xb2,zb

The plot draws some of the contours over to the left past 0.

Migrated-From: http://dunkel.pmel.noaa.gov/trac/ferret/ticket/1688

karlmsmith commented 6 years ago

Comment by @AnsleyManke on 22 Sep 2009 16:18 UTC Another example with FILL


let xb = x[x=0:290:10]
let yb = y[y=0:290:10]
let xb2 = xb + 0*yb
let zb = I[x=0:290:10] * J[y=0:290:10]
fill/hlimits=-100:400:50 0.0000001*(xb2+zb)^3,xb2,zb

This is fixed in routines fer/plt/load_curvilienar.F, fer/plt/curv_coord_sub.F, fer/plt/curv_coord_add_sub.F, taking care to separate the cases of curvilinear coordinates in units of longitude and modulo operations on longitudes, from the case where the curvilienar X coordinate is not modulo longitude.

karlmsmith commented 6 years ago

Attachment from @AnsleyManke on 18 Sep 2009 17:22 UTC the plot from the example script contourbug