arr_a = rand(nx,ny)
arr_b = rand(nx,ny)
print('write arrays')
pprint(arr_a)
pprint(arr_b)
#------------------------------------------------------------------------------
# write a nc file
#------------------------------------------------------------------------------
ncf = Dataset('write_read.nc', 'w', format='NETCDF4')
# attributes
ncf.description = 'netCDF4 write/read test' # string
ncf.size = nx*ny # number
# dimensions
ncf.createDimension('nx', nx)
ncf.createDimension('ny', ny)
# variables
va = ncf.createVariable('va', 'f8', ('nx','ny'))
va.unit = 'm/s'
vb = ncf.createVariable('vb', 'f8', ('nx','ny'))
vb.unit = 'kg'
# write data to variables
va[:] = arr_a
vb[:] = arr_b
