%% ***************************************** % GCH 2535 - Hiver 2019 - Devoir No2 % Equipe: % % % %% ***************************************** clear all %% ***************************************** % Prétraitement % ***************************************** %Constantes physiques k_eau=; %Conductivité thermique de l'eau (W/(m.K)) rho_eau=; %Densité de l'eau (kg/m3) c_eau=4186; %Capacité thermique massique de l'eau (J/(kg.K)) alpha_eau=; %Diffusivité thermique de l'eau (m2/s) L_eau=; %Entalpie de fusion de l'eau (J/kg) h_air=; %Coefficient de convection thermique (W/((m2.K)) enso =; %Ensoleillement (W/m2) albd_eau =; %Albédo des lacs (-) %Conditions initiales T_glace=; delta_T=; T_air=; %Température de l'air (K) %Discrétisation du domaine H=1; %Profondeur du domaine (m) h0=; %Profondeur d'eau (m) h=; %Profondeur d'eau (m) N=; %Nombre de noeuds du domaine dx=; %Pas en espace (m) %Déclaration des vecteurs et matrices %Initialisation des vecteurs %% ***************************************** % Solveur % ***************************************** %Paramètres temporels tfin=3600*24*200; % Temps final (s) dt=; % pas de temps (s) numero_iteration=; frequence_sortie=; t=; %Boucle en temps while (t