{\rtf1\mac\ansicpg10000\cocoartf102
{\fonttbl\f0\fswiss\fcharset77 Helvetica;}
{\colortbl;\red255\green255\blue255;}
\margl1440\margr1440\vieww9000\viewh9000\viewkind0
\pard\tx720\tx1440\tx2160\tx2880\tx3600\tx4320\tx5040\tx5760\tx6480\tx7200\tx7920\tx8640\ql\qnatural

\f0\fs24 \cf0 % runs a one-dimensional thermal/plume cloud model, based on\
% Rogers&Yau, stepping forward in time in 60second increments\
% and for different initial values of \
%beta, the product of thermal volume*buoyancy*gravity\
% february 21, 2006, MSC409\
\
beta=0.1 ;\
a= 6.6   ;\
%a is another constant, empirically based,relating time,beta, and height\
b=0.5 ;\
%b is a constant relating radius of the thermal to height\
%initialize arrays\
z=fltarr(7)\
r=z\
FOR itime = 0:420:60     %time in seconds\
	z_squared=a*sqrt(beta)*itime ;\
	z(itime)=sqrt(z_squared) ;\
	radius(itime)=b*z(itime) ;\
END\
time=(0,60,120,180,240,300,360,420) ;\
% plot,z and r as a function of time, include value of beta in the title\
figure; plot(z,time) ;\
figure;plot(r,time) ;\
\
%save values to an ascii file\
filename='thermal_model_out.txt' ;\
SAVE(filename,'time','z','r', '-ascii','-tabs') ;