newflxio_nc = '/ptmp/yeager/sfwf_tx1v1.nc'
newsalx_nc = '/ptmp/yeager/tx1v1/PHC2_salx_tx1v1.nc'
outbin = '/ptmp/yeager/tx1v1/sfwf_mm_PHC2_salx_flxio_20080502.ieeer8'
outnc = '/ptmp/yeager/tx1v1/sfwf_mm_PHC2_salx_flxio_20080502.nc'

nx = 360
ny = 240

ncid = ncdf_open(newsalx_nc)
varid = ncdf_varid(ncid,'SALT')
ncdf_varget,ncid,varid,newsss
newsss = reform(newsss)
ncdf_close,ncid

ncid = ncdf_open(newflxio_nc)
varid = ncdf_varid(ncid,'flxio')
ncdf_varget,ncid,varid,newflxio
ncdf_close,ncid

ncid = ncdf_create(outnc,/clobber)
dim1 = ncdf_dimdef(ncid,'lat',ny)
dim2 = ncdf_dimdef(ncid,'lon',nx)
dim3 = ncdf_dimdef(ncid,'month',12)
var1 = ncdf_vardef(ncid,'sss',[dim2,dim1,dim3],/double)
var2 = ncdf_vardef(ncid,'flxio',[dim2,dim1,dim3],/double)
ncdf_attput,ncid,var1,'missing_value',-99.
ncdf_attput,ncid,var2,'missing_value',-99.
ncdf_control,ncid,/ENDEF
ncdf_varput,ncid,var1,newsss
ncdf_varput,ncid,var2,newflxio
ncdf_close,ncid


openw,1,outbin,/SWAP_END
writeu,1,double(newsss)
writeu,1,double(newflxio)
close,1

end