[ferret_users] ASCII to netCDF issue

Hi Ferret Users,

I’m working on converting an ASCII file to netCDF for use in the LAS. I’ve included the file in this message (EV000-87.dat), and the netCDF file this script creates.

Here is the script:

def axis/y=90.0:-90.0/npoints=37/units="degrees_north" yax

def axis/x=0:355/npoints=72/units="degrees_east" xax

define axis/t="23-oct-2003":"23-oct-2003":1/units=hours tax1

def grid/x=xax/y=yax/t=tax1 my_grid

file/skip=2/columns=73/order=yx/var="dummy1,effective_dose"/grid=my_grid "/home/doostra/NAIRAS/EV000-87.dat"

set variable/title="Effective Dose (NAIRAS)" effective_dose

set variable/units="uSv/hr"/bad=-999.0 effective_dose

save/file=nairas_hal-87.nc/clobber effective_dose

The issue is that the netcdf file only has 18 of the 37 latitude points, so in the LAS it only plots half of the grid. I’ve tried defining y with modulo and this gets me all the lat points, but the data is missing for half again. I’ve also tried defining y with modulo and x=180:-180, as well as some other things—basically grasping at straws.

Does anyone see anything I’m doing wrong? Or, have a better way to do this? I’m open to any and all feedback.

Thank you!

Dan

Daniel Oostra (SSAI)

NASA/LaRC Atmospheric Sciences Data Center

2 South Wright Street

Hampton, VA 23681-2199

Ph: 757.864.6157

daniel.h.oostra@xxxxxxxx