for (ia,alt) in enumerate(altval): # Altitudes.
# Pull out value.
[locval, ialt] = gitmLib.locValue(mObj, t, args.var, alt, \
args.lat*np.pi/180.0, \
args.lon*np.pi/180.0)
# Fill array
arr[it,ia] = locval
# Write output to a file.
try:
root = Dataset(outfile, 'w', clobber=True)
except:
sys.stderr.write("FATAL ERROR: tseriesPrfLoc: error opening file "+outfile+". Terminating\n")
sys.exit(1)
root.description = "Produced by tseriesPrfLoc. Vertical profiles at a location versus time"
root.latitude = args.lat
root.longitude = args.lon
root.physicalvar = args.var
root.createDimension('time',ntimes)
root.createDimension('altitude',nalt)
tvar = root.createVariable('time','f8',('time',))
avar = root.createVariable('altitude','f8',('altitude',))
var = root.createVariable('variable','f8',('time','altitude'))
tvar[:] = times.copy()
avar[:] = altval.copy()
var[:] = arr.copy()
root.close()
# fout.write("#Variable {:s} at {:6.2f} Longitude and {:6.2f} Latitude and Altitude {:6.2f}\n".format(args.var, args.lon, args.lat, alt))
# fout.write("#TimeJ2000 TimeStr LocVal\n")
