	% set parameter values
	
	Gcons=6.673e-11;
	rhom=3300.;
	rhow=1025.;
	rhoc=2800.;
	dseaf=4000.;
	dmoho=7000.;
	flwp=5000.;
	
	seafint=((rhoc-rhow)/(rhom-rhow))*exp(-dseaf/flwp);
	mohoint=((rhom-rhoc)/(rhom-rhow))*exp(-(dseaf+dmoho)/flwp);
	

	delgovw = 2*pi()*Gcons*(rhom-rhow)*(seafint+mohoint)*1.0e5
	
	