Dear Ferret users,
I was surprised that a construct like temp[x=120:150] refers to strange longitudes when 120 or 150 is out of the actual defined range. A self-contained example is pasted below.
This is problematic when you save a subset of data because the longitude axis of the new dataset will include these strange longitude points unless you carefully specify the longitude range when saving it.
Is this an expected behavior or is it something that should be fixed in a future release?
Ryo
----
! PyFerret v7.63 (optimized)
! Darwin 19.6.0 - 12/10/21
yes? set data coads_climatology
yes? save/clobber/file=tmp.nc SST[x=130:140,y=20,L=1]
LISTing to file tmp.nc
yes? set data tmp.nc
yes? list sst[x=120:135]
[ . . . ]
45W / 0: ....
131E / 1: 25.68
133E / 2: 25.90
135E / 3: 26.13
yes? list sst[x=135:140]
[ . . . ]
135E / 3: 26.13
137E / 4: 26.31
139E / 5: 26.39
yes? list sst[x=135:145]
[ . . .]
135E / 3: 26.13
137E / 4: 26.31
139E / 5: 26.39
45W / 6: ....
yes? save/file=tmp-new.nc sst[x=120:145] ! acquires strange longitudes.
yes? plot sst