! err64_packed_data.jnl
! Test writing packed data. Prior to v6.4 this example shows the
! crash described in bug 1715.
! 4/15/19: Output changes with fix to issue 334: coordinates are not 
! regular but have a point_spacing="even" attribute. Ferret now checks
! and overrides this.

! ncpdq from NCO operators 
! sp ncpdq fine_x_axis.nc err64_packed_data.nc
use err64_packed_data.nc
sh att/all sst

! check that the coordinates packed by ncpdq are unpacked ok.
list/i=1:15 sst

! we can write the data out - saving the scaling to re-pack.
set att/output sst.scale_factor
set att/output sst.add_offset

! With double-precision Ferret, need to write as float,
! or the FillValue cant be represnted. 

set var/outtype=float sst

save/clobber/file=mypack.nc/i=1:15 sst
can data/all; use mypack.nc
sh att/all sst

list/i=1:15 sst

