Geology Reference
In-Depth Information
%alphanoise.m for 1/f^alpha noise models, alpha =
1 and 2
% Note: outputs alpha1 and alpha2 are time series
realizations;
% mag1 and mag2 are the theoretical noise
spectra
npts=2048;
fmag=0:1/npts:0.5-1/npts;
fmag=fmag';
[alpha1,mag1]=alphanoise(1,npts);
[alpha2,mag2]=alphanoise(2,npts);
mag1=mag1';
mag2=mag2';
mag12=mag1.^2;
mag22=mag2.^2;
mag12a=(length(mag12)/sum(mag12)) * mag12;
mag22a=(length(mag22)/sum(mag22)) * mag22;
figure; plot(fmag,mag12a);hold all; plot(fmag,mag22a)
FigureĀ  4.19 Hypothesis testing of the signal + red noise time series in
FigureĀ 4.17 .
dt=1;num=1000;nw=4;
[Mspecred,specred,po,fd,fr,tabMC,tabchi,tabtchi,t
heored,rho]=redconf(signalmarkov5,tt,dt,num,nw);
figure;plot(fd,po);hold all; plot(fd,theored);hold
all;plot(fd,tabtchi);
A.4 Computation of the Obliquity
and Precession Index
For the La2010a (or b, c, or d) solution calculations for the obliquity and
precession index can be made using FORTRAN programs as follows:
(1) The Laskar et al. (1993) FORTRAN code and input parameters files may
be downloaded from the online VizieR Catalogue Service:
http://vizier.cfa.harvard.edu/ftp/cats/VI/63/prepa.f
http://vizier.cfa.harvard.edu/ftp/cats/VI/63/prepa.par
http://vizier.cfa.harvard.edu/ftp/cats/VI/63/integ.f
http://vizier.cfa.harvard.edu/ftp/cats/VI/63/integsub.f
http://vizier.cfa.harvard.edu/ftp/cats/VI/63/integ.par
http://vizier.cfa.harvard.edu/ftp/cats/VI/63/climavar.f
http://vizier.cfa.harvard.edu/ftp/cats/VI/63/ climavar.par
Mac users must apply mac2unix ( one time only ) to convert all files after
downloading, which may be obtained here:
http://sourceforge.net/projects/cs-cmdtools/files/
Search WWH ::




Custom Search