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
[covar_spec covar_axes] = ...
     covar(tocsy_spec,hmbc_spec,'axes1',tocsy_axes,'axes2',hmbc_axes,'power',0.5);  
     % GIC performed is [TOCSY*HMBC]^1/2
write_covar2pipe(covar_spec, covar_axes,'hmbctocsy.gic'); 
toc

% example of plotting a region from a GIC spectrum

Plot2DSpec([3.0 4.5],[125.5 179.5],hmbc_spec,hmbc_axes,5*(10^7),1.3,30);
print -dpdf gic.pdf