c	file name TEST_TROUGH.FOR
c	prog to look at the latitudinal variation of the ICED trough
	open (unit=1,file='data.out',type='old')
	pi = 3.141592654
	y1 = exp(-0.375)
	do 20 i = 1,30
	x = - float(i) / 10.
	z = exp((x - x*x) / 2.)
	t_iced = 1. - 0.4 * z
	yy = z
	if(abs(x).le.0.5) yy = y1 * sin(pi/2.*abs(x)/0.5)
	t_mod = 1. - 0.4 * yy
	write(6,100) i,x,z,yy,t_iced,t_mod
	write(1,100) i,x,z,yy,t_iced,t_mod
100	format(i5,7f10.3)
20	continue
c
c	now try X1 > 0
	y1 = exp(0.125)
	do 30 i = 1,30
	x = float(i) / 10.
	z = exp((x - x*x) / 2.)
	t_iced = 1. - 0.4 * z
	yy = z
	if(abs(x).le.0.5) yy = y1 * sin(pi/2.*abs(x)/0.5)
	t_mod = 1. - 0.4 * yy
	write(6,100) i,x,z,yy,t_iced,t_mod
	write(1,100) i,x,z,yy,t_iced,t_mod
30	continue
	end