%%%Constants%%%%%%
Mdot_stream=4.5; %inches/month
deltaD_stream=-85;
delta18O_stream=-14;

Acache=10;
Runoff=0.15;

Z_lake_max=150; %inches
%%%%%Monthly Data%%%%%%%%%%%%
pan=[1.4 2.21 4.07 5.94 8.31 9.52 8.97 8.81 6.58 4.46 2.21 1.32];%inches
Ctemp=[47.8 53.3 57.4 63 71 78.2 84.1 82.6 76.8 67.8 55.8 47.5];%farenheit
Chum=[ 78 67 58 46 39 35 33 37 41 46 62 75];%relative humidity, percent
Cprecip=[ 0.86 1.06 1.04 .57 .2 .1 .01 .09 .17 .29 .7 .63];%inches


%%%%%%Initializations%%%%%%%%
Z_lake=125; %inches
delta_18O_lake=-2;
delta_D_lake=-10;

%%%%%%%Main loop%%%%%%%%%%%%
for itime=1:120
imonth=mod(itime,12)+1;
hum=Chum(imonth)/100;


DZ_lake=?????????????;
Z_lake=min(Z_lake+DZ_lake, Z_lake_max);
Z_lake_store(itime)=Z_lake;

end