AMATH_383_Solution_6

# AMATH_383_Solution_6 - AMATH 383 Solution#6 Matlab program...

This preview shows pages 1–3. Sign up to view the full content.

AMATH 383 Solution #6: Matlab program codes 1. clear %parameters A=202 B=1.90 C=1.6*B Q=343 T_c=-10 %albedo: alpha1 is ice-free, alpha2 is ice-covered alpha1=0.32 alpha2=0.62 alpha0=(alpha1+alpha2)/2 axis ([250 550 0 1]) %y=sin(latitude) %s(y) is the latitudinal distribution of solar heating y=0:0.005:1 alpha_bar=alpha2-(alpha2-alpha1)*y.*(1-0.241*(y.^2-1)) s=1-0.482*(3*y.^2-1)/2 de=(1-alpha0)*s+C*(1-alpha_bar)/B Q=(B+C)*(T_c+A/B)./de plot(Q,y,'LineWidth',2) hold on plot(343,y) ylabel('location of iceline, y=sin(latitude)') xlabel('Q in watts per square meter') hold on plot([250:441],0.995,'.k','LineWidth',2) hold on plot([330:550],0.005,'.b','LineWidth',2) 2. clear %parameters A=202 B=1.90 C=1.6*B Q=343 T_c=-10 %ice line location y_s=0.95 s_0=1-0.482*(3*y_s^2-1)/2

This preview has intentionally blurred sections. Sign up to view the full version.

View Full Document
%albedo: alpha1 is ice-free, alpha2 is ice-covered alpha1=0.32 alpha2=0.62 alpha0=(alpha1+alpha2)/2 %y=sin(latitude) %s(y) is the latitudinal distribution of solar heating y=0:0.005:1 s=1-0.482*(3*y.^2-1)/2 %T_free is ice-free solution. T_ice is ice-covered. T_part is partly
This is the end of the preview. Sign up to access the rest of the document.

## This note was uploaded on 10/21/2011 for the course APPLIED MA AMATH 383 taught by Professor Hongqian during the Spring '10 term at University of Washington.

### Page1 / 3

AMATH_383_Solution_6 - AMATH 383 Solution#6 Matlab program...

This preview shows document pages 1 - 3. Sign up to view the full document.

View Full Document
Ask a homework question - tutors are online