jstults
/ ice-code.R
Created
February 4, 2012 16:26
Comparison of Northern Hemisphere Perennial and Seasonal Sea Ice
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| ssq = function (x) {sum(x ^ 2)} | |
| plotTrend = function (x, st = 1978, en = c(2012, 12), y.pos = NA, x.pos = NA, main.t = "Untitled") { | |
| ### Get trend | |
| trend = lm(window(x, st, en) ~ I(time(window(x, st, en)))) | |
| ### Initialize variables | |
| N = length(window(x, st, en)) | |
| I = seq(1:N) / frequency(x) |
jstults
/ falknerskan.mac
Created
January 29, 2012 20:50
Maxima definition of Falkner-Skan governing equations
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| /* Falkner-Skan ODE | |
| Zhang, J., Chen, B., 'An iterative method for solving the | |
| Falkner-Skan equation' | |
| Nachtsheim, P., Swigert, P., 'Satisfaction of Asymptotic Boundary | |
| Conditions in Numerical Solution of Systems of Nonlinear Equations | |
| of Boundary-Layer Type,' NASA TN D-3004, Lewis Research Center, | |
| Cleveland, Oct, 1965. |
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| import scipy as sp | |
| from scipy.integrate import ode | |
| from scipy.optimize import fsolve | |
| from matplotlib import rc | |
| #rc('text', usetex=True) | |
| tick_size = 'large' | |
| label_size = 'xx-large' | |
| rc('xtick', labelsize=8) | |
| rc('ytick', labelsize=8) | |
| rc('legend', fontsize=12) |
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| module fs | |
| implicit none | |
| contains | |
| subroutine rates(F, t, y, arg) | |
| double complex, intent(in) :: t | |
| double complex, intent(in), dimension(3) :: y, arg | |
| double complex, intent(out), dimension(3,1) :: F | |
| F(1,1) = y(2)*arg(3) |
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| <?xml version="1.0" encoding="UTF-8" ?> | |
| <!DOCTYPE html> | |
| <html b:version='2' class='v2' expr:dir='data:blog.languageDirection' xmlns='http://www.w3.org/1999/xhtml' xmlns:b='http://www.google.com/2005/gml/b' xmlns:data='http://www.google.com/2005/gml/data' xmlns:expr='http://www.google.com/2005/gml/expr'> | |
| <head> | |
| <meta content='IE=EmulateIE7' http-equiv='X-UA-Compatible'/> | |
| <b:if cond='data:blog.isMobile'> | |
| <meta content='width=device-width,initial-scale=1.0,minimum-scale=1.0,maximum-scale=1.0' name='viewport'/> | |
| <b:else/> | |
| <meta content='width=1100' name='viewport'/> | |
| </b:if> |
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
| import scipy as sp | |
| from scipy import fft, ifft | |
| from scipy.integrate import cumtrapz | |
| from scipy.integrate import simps | |
| import csv | |
| from matplotlib import rc | |
| fig_width_pt = 469.75 | |
| inch_per_pt = 1.0 / 72.27 | |
| figw = fig_width_pt * inch_per_pt | |
| figh = (620.43 / 2.1) * inch_per_pt |
NewerOlder