c	prog to look at the latitudinal variation ofthe ICED trough
	open (unit=1,file='data.out',type='old')
	pi = 3.141592654
	y1 = exp(-0.375)
	do 20 i = 1,40
	x = - float(i) / 10.
	z = exp((x - x*x) / 2.)
	tz = 1. - 0.4 * z
	yy = z
	if(abs(x).le.0.5) yy = y1 * sin(pi/2.*abs(x)/0.5)
	tt = 1. - 0.4 * yy
	write(6,100) i,x,z,yy,tz,tt
	write(1,100) i,x,z,yy,tz,tt
100	format(i5,7f10.3)
20	continue
	end