glimmer ela model (EIS driver) from config file
0001 % glimmer ela model (EIS driver) 0002 % from config file 0003 ela_a = 14430.069930; 0004 ela_b = -371.765734; 0005 ela_c = 2.534965; 0006 bmax_mar = 4 0007 t=6; % time in model years 0008 lamb0=35; % latitude 0, 35deg here 0009 konstlat=60; % use 50deg N as input lat for timeseries 0010 % copy/pasted from mb2.data file 0011 %mb2=[-125000,350;-115000,350;-114000,350;-113000,400;-112000,225;-111000,250;-109000,225;-105000,350;-104000,100;-103000,75;-102000,50;-100000,100;-98000,250;-96000,250;-94000,300;-91000,150;-89000,125;-84000,50;-80000,350;-76000,500;-73000,600;-72000,70;-70000,400;-69000,125;-68000,300;-66000,-25;-65000,-250;-64000,-400;-62000,-400;-61000,-350;-59000,-50;-58000,0;-56000,200;-54000,300;-53000,200;-51000,250;-43001,250;-43000,300;-42000,300;-40000,200;-30000,200;-26000,150;-25000,-450;-20000,-500;-18000,-900;-17001,-900;-17000,-175;-15000,-50;-14999,200;-12000,200;-8000,200;-6000,300;-4999,300;-4000,700;-1000,700;0,300;]; 0012 0013 %% read in lat and grid 0014 [xcoord,ycoord,alp10k,alpheader]=asciigrid2mat('L:\glimmer_alpen\alpen10k.asc'); 0015 [lxcoord,lycoord,lat,latheader]=asciigrid2mat('L:\glimmer_alpen\alps10k_lat.asc'); 0016 0017 %% model: 0018 0019 ELA = ela_a + ela_b*lat+ela_c*lat^2+ mb2(t,2) ; 0020 0021 %% dump 0022 0023 mat2asciigrid(['L:\glimmer_alpen\ela_t',num2str(t),'.asc'],xcoord,ycoord,ELA); 0024 0025 0026 0027 %% 0028 %load('D:\Matlab\MCS_Scand\glimmer_mb'); 0029 lamb0=35; % 0030 lat=60; 0031 elaseries=zeros(size(mb2,1),2); 0032 for t=1:size(elaseries) 0033 elaseries(t,1)= mb2(t,1); 0034 elaseries(t,2) = ela_a + ela_b*konstlat+ela_c*konstlat^2 + mb2(t,2) ; 0035 end 0036 0037 % start -125000 at t=1 end at -80000 at 19 0038 figure, plot(1:41,elaseries(1:19,2))