T_in_file = 'WOCE_PTEMP_MODIFIED_gx3v6_climatology.nc'
S_in_file = 'WOCE_SALT_MODIFIED_gx3v6_climatology.nc'
out_file = 'WOCE_MODIFIED_ic_gx3v6.ieeed'

vart = 'TEMPERATURE_MODIFIED'
vars = 'SALINITY_MODIFIED'

ncid = ncdf_open(T_in_file)
varid = ncdf_varid(ncid,vart)
ncdf_varget,ncid,varid,temp
ncdf_close,ncid
km = (size(temp))(3)

openw,1,out_file
for i=0,km-1 do begin
  writeu,1,double(reform(temp[*,*,i]))
endfor

ncid = ncdf_open(S_in_file)
varid = ncdf_varid(ncid,vars)
ncdf_varget,ncid,varid,salt
ncdf_close,ncid

for i=0,km-1 do begin
  writeu,1,double(reform(salt[*,*,i]))
endfor

close,1

end