use axis_lat.nc

let a=y[gy=lat]

save/clobber/file=axis_lat_new.nc a

! compare precision with
! ncdump axis_lat.nc
! and
! ncdump axis_lat_new.nc