tic
[hmbc_spec hmbc_axes] = read_nmrp('hmbc.ft2',1); 
[tocsy_spec tocsy_axes] = read_nmrp('tocsy.ft2',1);
     % as in previous scripts, read in files for indirect covariance
[Z_spec Z_axes] = ...
     covar(tocsy_spec,hmbc_spec,'axes1',tocsy_axes,'axes2',hmbc_axes,'power',1,'Z',1);  
     % Z-matrix transform of HMBC-TOCSY unsymmetric covariance spectrum
     % value of 1 for 'Z' flag produces Z matrix (2 produces 'mrs' result)
write_covar2pipe(Z_spec, Z_axes,'hmbctocsy.zmat'); 
toc