In MATLAB I'm trying to write a function that determines the design height and width of the snowdrift (see attached image) given Hc, Lu and Llower. Also calculate the maximum pressure due to the snow. The function should have Three inputs, three outputs. [out1, out2, out3] = snow_drift (in1, in2, in3).

Given:

- Drift snow loads (the triangular portion) shall be superimposed on the balanced snow load (the rectangular portion).

- If Hc/Hb

- For leeward drifts ((left to right in picture) : Hd = 0.43 * (Lu^1/3) * ((pg+10)^1/4 )-15

- where Hd is the height of the drift (ft),

- Lu is the length of the upper roof (ft),

- pg is the ground snow load (psf). EXCEPT: if lu

- IF Hd Hc, the drift height is equal to Hd and the drift width, w, equals 4hd. If Hd > Hc, the drift height equals Hc and: W = (4*Hd^2)/Hc

- For windward drifts (right to left) use the same equations for Hd and W, but use the length of the lower roof for Lu and multiply the resulting value of Hd by 0.75.

- The design drift shall be the larger of the windward and leeward drifts

- pg = 90 psf

- Knowing the uniform load, we can calculate the height if we know the density. Fortunately, ASCE 7 gives us density, y = 0.13pg + 25 but no more than 30 pcf. The balanced snow load is 0.7pg; given y, you can calculate Hb.

So Far I have coded for this:

function [Hd, W, Hb] = snow_drift(Hc, Lu, Llower) % this function determines the design height and width of the now drifft % when given Hc, Lu, and Llower pg = 90; Hd = 0.43 * (Lu^1/3) * ((pg+10)^1/4 )-15; if Lu

if Hd

if Hd>Hc Hd = Hd; W = (4*Hd^2)/Hc; end

Hb =

if Hd/Hb

end

Thanks for the help!

Surcharge Load Due to Drifting Balanced Snow Load 1 1 1 1 1 Surcharge Load Due to Drifting Balanced Snow Load 1 1 1 1 1