%% make an estimate of precipitation efficiency
%% according to Smith P 175 this is 
%% obs rainfall over computed condensed water
%% run munich version to set the variables

%% columns 0-250 are the z values up to the 3 km freezing level
%% plot(z(1:250))

     estpe=sum(qr(1:250))/sum(qc(1:250))
     
% table
% w=2 m/s maxz=5 km, estpe=0.1456
% w=2 m/s maxz=7 km, estpe=0.1280
% w=5 m/s maxz=5 km, estpe=0.3971
% w=5 m/s maxz=7 km, estpe=0.3590
% w=7 m/s maxz=5 km, estpe=0.4165
% w=7 m/s maxz=7 km, estpe=0.4188
% w=8 m/s maxz=5 km, estpe=0.3719
% w=8 m/s maxz=7 km, estpe=0.4072
% w=9 m/s maxz=5 km, estpe=0.3132
% w=9 m/s maxz=7 km, estpe=0.3717
% w=10 m/s maxz=5 km, estpe= 0.2558
% w=10 m/s maxz=7 km, estpe=0.3239
