density_model_Goldstein.m

%% Density model of Goldstein et al., 2016, JGR

% Cosine Coefficients vs L: cm(L) for L=[4:19]
Cm = [...
    25.7 21.8 12.1 5.66   3.09  1.98 0.712 0.559 0.328 0.207 0.142 0.131 0.123 0.109 0.0894 0.0796;
    4.54 5.27 4.57 0.131 -0.295 0.166 0.0337 -0.00525 -0.0844 -0.00714 0.0225 0.0151 0.0215 0.0309 0.0246 0.0229;
    1.8 1.96 -0.456 0.106 0.0714 -0.319 -0.0662 -0.0993 0.0328 -0.0292 0.0164 -0.00806 -0.00806 -0.00129 0.00632 0.00112;
    4.48 1.3 -0.924 -0.00976 0.201 -0.14 -0.047 -0.0305 -0.0231 0.0233 -0.00718 0.00244 0.00463 -0.000285 0.00273 0.00118;
    -4.59 -1.4 -0.526 0.0815 -0.0384 -0.0474 0.127 0.0508 -0.0148 -0.0123 0.00593 0.00342 0.000359 0.00319 0.00516 0.0038;
    3.28 1.46 0.117 0.115 0.0155 -0.0313 -0.0191 0.0129 0.0177 0.00446 -0.000558 0.0037 0.00256 0.000964 -0.000815 0.0106];

% Sine Coefficients vs L: sm(L) for L=[4:19]
Sm = [...
    0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0;
    -2.26 0.667 0.0393 -0.932 -0.0963 0.196 -0.128 -0.16  -0.0587 -0.0305 -0.0122 -0.0355 -0.0419 -0.0243 -0.00542 -0.00561;
    5.61 0.96 -2.06 -0.85 0.24 -0.0342 0.0501 0.0134 0.0505 -0.017 0.00144 0.00879 -0.00178 -0.00784 0.00851 0.00132;
    -3.25 -0.744 0.446 0.138 0.0429 0.0159 -0.0375 0.0327 -0.0269 0.0247 0.00145 0.00222 -0.00568 -0.000499 0.00107 0.0032;
    -4.16 0.191 0.722 0.267 -0.109 -0.204 -0.0235 0.0134 -0.0133 -0.0181 -0.00731 -0.00483 -0.00489 0.003 0.0045 0.00507;
    2.6 1.11 0.579 0.0242 0.035 0.13 -0.0334 0.000942 -0.00108 0.00536 -0.00391 -0.00418 -0.00426 -0.00453 -0.00526 -0.01];

save('density_goldstein2016.mat', 'Cm', 'Sm')