# Prove the earth is round DIY

• 20 Replies
• 2928 Views
?

#### Kasper Marstal

• 186
##### Prove the earth is round DIY
« on: August 15, 2008, 08:51:10 AM »
Hello.

Step 2: Run the code below and see that it predicts the correct current global mean temperature distribution (which is saved as the var G in the code).

Step 3: Rewrite the code to fit a flat earth model.

Step 4: Observe that it predicts severely incorrect temperatures.

Quote

% Two-Dimensional Simple Climate Model
% Version 1.0
% Matlab code by Kasper Korsholm Marstal (Copyright 2008)
%
% Based on an Energy Balance Model by North (1975).
% Enhaced with a latitude dependent albedo where no
% ice caps are present, a latitude dependent
% heat capacity, a temperature dependent meriodinal

%% INITIALIZE
clear all;
echo off;

%% Control (set to 0 to fall back to North parameters)

%% SETUP
sb = 5.67*10^(-8); % [sb] = W m^-2 K^-4
yr = 365*24*3600; % [yr] = s
whc = 4180*1000; % heat capacity for one cubmic metre of water, [whc] = J kg^-1 K^-1
A = 203.9; % outgoing radiation parameterization for North feedback
B = 2.09; % outgoing radiation parameterization for North feedback
C_0 = whc*50; % heatcapacity for 50 m^3 water, [C] = J K^-1
D = 0.445; % diffusion constant
dt = 20000; % time resolution, [dt] = s
k = 0.478; % value of climate sensitivity for advanced feedback
deltat_0 = 33.4; % [deltat_0] = K
numx = 181; % grid resolution
Q = 1367/4; % Solar constant [sc] = W m^-2
tmax = 93*yr; % [tmax] = yr

%% SIMULATION ALGORITHM
% Set arrays to speed up loops
% Set X-array
xsouth = -1;
xnorth = 1;
dx = (xnorth-xsouth)/(numx-1);
X = [];
for i = 0 : 1 : numx-1
X = [ X xsouth+dx*(i) ];
end

% Set solar distribution array
S = [];
for i = 1 : 1 : numx
S = [ S 1-0.241*(3*X(i)^2-1) ];
end

% Set latitude dependent albedo array if advalb = 1
ALB = [];
for i = 1 : 1 : numx
ALB = [ ALB 0.30+0.1*(3*X(i)^2-1) ];
end
end

% Set latitude dependant heat capacity array if advheatcap = 1
C = [];
for i = 1 : 1 : numx
C = [ C dt/((160-40*(3*X(i)^2-1))*whc) ]; %dt/((200-100*(3*X(i)^2-1))*whc)
end
end

% Set X1 array (used in evaluation of the diffusion)
X1 = [];
for i = 1 : 1 : numx
n = i-1;
if(n < 1)
n = n+numx;
end
X1 = [ X1 1-((1/2)*(abs(X(i))+abs(X(n))))^2 ];
end

% Set initial temperature conditions
% First guess
T = [];
for i = 1 : 1 : numx
T = [ T 303.16-45.*X(i)^2 ];
end

% Equilibrium with present day parameters
G = [244.2380,246.5491,248.1093,249.5257,250.8844,252.2142,253.5281,254.8332,256.1338,257.4330,258.7326,260.0344,261.3395,262.6489,263.9635,265.2408,266.4844,267.6973,268.8821,270.0406,271.1745,272.2853,273.3741,274.4418,275.4894,276.5175,277.5267,278.5177,279.4907,280.4463,281.3848,282.3064,283.2115,284.1001,284.9726,285.8291,286.6697,287.4946,288.3038,289.0975,289.8756,290.6384,291.3858,292.1179,292.8347,293.5363,294.2226,294.8937,295.5496,296.1903,296.8158,297.4262,298.0213,298.6012,299.1660,299.7155,300.2498,300.7688,301.2726,301.7611,302.2344,302.6923,303.1350,303.5622,303.9742,304.3708,304.7519,305.1177,305.4681,305.8030,306.1224,306.4264,306.7149,306.9879,307.2454,307.4874,307.7138,307.9246,308.1199,308.2996,308.4637,308.6122,308.7451,308.8624,308.9640,309.0501,309.1205,309.1752,309.2143,309.2378,309.2456,309.2378,309.2143,309.1752,309.1205,309.0501,308.9640,308.8624,308.7451,308.6122,308.4637,308.2996,308.1199,307.9246,307.7138,307.4874,307.2454,306.9879,306.7149,306.4264,306.1224,305.8030,305.4681,305.1177,304.7519,304.3708,303.9742,303.5622,303.1350,302.6923,302.2344,301.7611,301.2726,300.7688,300.2498,299.7155,299.1660,298.6012,298.0213,297.4262,296.8158,296.1903,295.5496,294.8937,294.2226,293.5363,292.8347,292.1179,291.3858,290.6384,289.8756,289.0975,288.3038,287.4946,286.6697,285.8291,284.9726,284.100,283.2115,282.3064,281.3848,280.4463,279.4907,278.5177,277.5267,276.5175,275.4894,274.4418,273.3741,272.2853,271.1745,270.0406,268.8821,267.6973,266.4844,265.2408,263.9635,262.6489,261.3395,260.0344,258.7326,257.4330,256.1338,254.8332,253.5281,252.2142,250.8844,249.5257,248.1093,246.5491,244.2380];

tmod = [ G ; G ]; % temperature distribution arrays for evalution
Z = G ; % temperature distribution matrix for save

% Initiate vars
c = 0; % time counter for temperature distribution save
p1 = 1; % pointer 1
p2 = 2; % pointer 2

% Further speed up loops
dtC = dt/C_0;
end
dx2 = dx^2;
yrdt = yr/dt;
Tmean = mean(T);

% Find temperatures
for t = 1 : dt : tmax % time-dependent loop
dc = D*(1+0.03*(mean(tmod(p1,:))-288)); % Temperature dependent diffusion cofficient
else
dc = D;
end
for i = 1 : 1 : numx % latitude-dependent loop
in = i-1;
ip = i+1;
if (in < 1);  % these two ifs ensures symmetry around poles
in = in+numx;
end
if (ip > numx);
ip = ip-numx;
end
if (tmod(p1,i) < 263.16)
a = 0.6;
else
a = ALB(i);
end
else
if (tmod(p1,i) < 263.16)
a = 0.6;
else
a = 0.3;
end
end
dtC = C(i);
end
d = (dc/dx2)*(X1(ip)*(tmod(p1,ip)-tmod(p1,i))-X1(i)*(tmod(p1,i)-tmod(p1,in)));
f = -sb*(tmod(p1,i)-(deltat_0+k*(tmod(p1,i)-Tmean)))^4; % advanced feedback
else
f = -A-B*(tmod(p1,i)-273.16); % North feedback
end
tmod(p2,i) = tmod(p1,i)+dtC*(Q*S(i)*(1-a)+f+d); % evaluate and save T(x,t)
end
% Save temperature distribution once a year
c = c+1;
if (c > yrdt)
Z = [ Z ; tmod(p2,:) ];
c = 0;
end
% Switch pointers
psave = p1;
p1 = p2;
p2 = psave;
end

%% OUTPUT

% Surfaceplot of temperature distribution matrix
surf(Z)
xlabel('Latitude');
ylabel('Year');
zlabel('Temperature');

% Global mean temperature array
M = [];
for i = 1 : 1 : (tmax/yr)-1;
M = [ M mean(Z(i,:)) ];
end

% Global mean temperature change
delta_T = M((tmax/yr)-1)-M(1);
[/size]

?

#### semperround

• 2903
##### Re: Prove the earth is round DIY
« Reply #1 on: August 15, 2008, 08:52:42 AM »
let me see here, where's my damn calculator.
an vir

?

#### semperround

• 2903
##### Re: Prove the earth is round DIY
« Reply #2 on: August 15, 2008, 09:09:12 AM »
i think mr. common sense put it simpler in this thread http://theflatearthsociety.org/forum/index.php?topic=22241.0. this is easier for us lazy folk to understand.
an vir

?

#### Kasper Marstal

• 186
##### Re: Prove the earth is round DIY
« Reply #3 on: August 15, 2008, 11:37:54 AM »
Bump.

This is the evidence you've been asking for; anyone disprove flat earth theory in their own livingroom. This topic made it to the 2nd page without a single FE comment.

The earth is round.

?

#### Kasper Marstal

• 186
##### Re: Prove the earth is round DIY
« Reply #4 on: August 16, 2008, 05:04:37 AM »
Bump.

Boy, it feels good. I'm able to say the earth is round without being contradicted! This topic has made it to the second page 2 times without a single FE comment.

The earth is round.

#### Parsifal

• Official Member
• 36115
• Bendy Light specialist
##### Re: Prove the earth is round DIY
« Reply #5 on: August 16, 2008, 05:06:55 AM »

Any advice on how to do this?
I'm going to side with the white supremacists.

?

#### Kasper Marstal

• 186
##### Re: Prove the earth is round DIY
« Reply #6 on: August 16, 2008, 05:25:30 AM »

Any advice on how to do this?

You can buy it although it is expensive. If you study you can get it for almost no money through your university. You can download it on mininova.org.

#### Parsifal

• Official Member
• 36115
• Bendy Light specialist
##### Re: Prove the earth is round DIY
« Reply #7 on: August 16, 2008, 05:34:06 AM »

Only for Windows, it appears.
I'm going to side with the white supremacists.

?

#### Kasper Marstal

• 186
##### Re: Prove the earth is round DIY
« Reply #8 on: August 16, 2008, 06:10:46 AM »
I'm sorry, you'll have to figure out another way to get it or you won't be able do disprove the earth is flat yourself.

#### Parsifal

• Official Member
• 36115
• Bendy Light specialist
##### Re: Prove the earth is round DIY
« Reply #9 on: August 16, 2008, 06:29:47 AM »
I'm sorry, you'll have to figure out another way to get it or you won't be able do disprove the earth is flat yourself.

Demonoid has it for Linux, it seems. Unfortunately, the torrent for the 2008 version doesn't seem to be connecting to the server (it's externally tracked), and the 2007 version only has two seeders.
I'm going to side with the white supremacists.

?

#### Kasper Marstal

• 186
##### Re: Prove the earth is round DIY
« Reply #10 on: August 16, 2008, 08:05:26 AM »
isohunt.com has a matlab torrent.

#### sokarul

• 17785
• Discount Chemist
##### Re: Prove the earth is round DIY
« Reply #11 on: August 16, 2008, 08:35:38 AM »
Can you recode that for mathmatica?
ANNIHILATOR OF  SHIFTER

It's no slur if it's fact.

?

#### Kasper Marstal

• 186
##### Re: Prove the earth is round DIY
« Reply #12 on: August 16, 2008, 09:04:52 AM »
Can you recode that for mathmatica?

I'm not familiar with mathematica so that would be rather timeconsuming. I guess the answer is no

?

#### cbarnett97

• 2746
##### Re: Prove the earth is round DIY
« Reply #13 on: August 16, 2008, 05:55:05 PM »
if you do like to program just get a telescope any will do (the better the more accurate the results) and chart the orbital paths of the planets for a couple of years and then show that the earth is not rotating, the planets rotate about us, and the distances are not what we thought. If you do not want to take the time to do that then just get the observastions of others and use those to calculate the above, I would suggest the observations made by Brahe
Only 2 things are infinite the universe and human stupidity, but I am not sure about the former.

?

#### Kasper Marstal

• 186
##### Re: Prove the earth is round DIY
« Reply #14 on: August 16, 2008, 06:21:47 PM »
if you do like to program just get a telescope any will do (the better the more accurate the results) and chart the orbital paths of the planets for a couple of years and then show that the earth is not rotating, the planets rotate about us, and the distances are not what we thought. If you do not want to take the time to do that then just get the observations of others and use those to calculate the above, I would suggest the observations made by Brahe

I wouldn't suggest using observations from Brahe since they were conducted in the late 15th century. Do you really think this hasn't been done before? By astronomers and hobbyists alike?

?

#### cbarnett97

• 2746
##### Re: Prove the earth is round DIY
« Reply #15 on: August 16, 2008, 06:33:41 PM »
if you do like to program just get a telescope any will do (the better the more accurate the results) and chart the orbital paths of the planets for a couple of years and then show that the earth is not rotating, the planets rotate about us, and the distances are not what we thought. If you do not want to take the time to do that then just get the observations of others and use those to calculate the above, I would suggest the observations made by Brahe

I wouldn't suggest using observations from Brahe since they were conducted in the late 15th century. Do you really think this hasn't been done before? By astronomers and hobbyists alike?
That is my point, in all that time not one person has noticed anything wrong with them? not to mention he was setting out to prove a geocentric universe while kepler was trying to prove circular obits and yet they were both wrong
Only 2 things are infinite the universe and human stupidity, but I am not sure about the former.

?

#### Kasper Marstal

• 186
##### Re: Prove the earth is round DIY
« Reply #16 on: August 16, 2008, 06:45:21 PM »
if you do like to program just get a telescope any will do (the better the more accurate the results) and chart the orbital paths of the planets for a couple of years and then show that the earth is not rotating, the planets rotate about us, and the distances are not what we thought. If you do not want to take the time to do that then just get the observations of others and use those to calculate the above, I would suggest the observations made by Brahe

I wouldn't suggest using observations from Brahe since they were conducted in the late 15th century. Do you really think this hasn't been done before? By astronomers and hobbyists alike?
That is my point, in all that time not one person has noticed anything wrong with them? not to mention he was setting out to prove a geocentric universe

Jesus fucking christ. You proposed I should use them. And of course they are not as accurate as today. Of course "someone" noticed. You think you can find a person today who will claim otherwise?

Quote
while kepler was trying to prove circular obits and yet they were both wrong

Kepler was not trying to prove circular orbits, he proved that planets orbit in elliptic shapes and he was right. Jesus fucking christ.

And for the record, FE is more geocentric than RE since it violates the cosmological principle.

?

#### cbarnett97

• 2746
##### Re: Prove the earth is round DIY
« Reply #17 on: August 16, 2008, 07:14:03 PM »
if you do like to program just get a telescope any will do (the better the more accurate the results) and chart the orbital paths of the planets for a couple of years and then show that the earth is not rotating, the planets rotate about us, and the distances are not what we thought. If you do not want to take the time to do that then just get the observations of others and use those to calculate the above, I would suggest the observations made by Brahe

I wouldn't suggest using observations from Brahe since they were conducted in the late 15th century. Do you really think this hasn't been done before? By astronomers and hobbyists alike?

That is my point, in all that time not one person has noticed anything wrong with them? not to mention he was setting out to prove a geocentric universe

Jesus fucking christ. You proposed I should use them. And of course they are not as accurate as today. Of course "someone" noticed. You think you can find a person today who will claim otherwise?

Quote
while kepler was trying to prove circular obits and yet they were both wrong

Kepler was not trying to prove circular orbits, he proved that planets orbit in elliptic shapes and he was right. Jesus fucking christ.

And for the record, FE is more geocentric than RE since it violates the cosmological principle.
He was tyring to prove circular orbits but he discovered that they had elliptical orbits instead while Brahe was trying to prove a geocentric universe and eded up helping to prove otherwise. Which is my point, that somtimes the obsevations do not fit the theory
Only 2 things are infinite the universe and human stupidity, but I am not sure about the former.

?

#### Kasper Marstal

• 186
##### Re: Prove the earth is round DIY
« Reply #18 on: August 16, 2008, 07:21:45 PM »
Which is my point, that somtimes the obsevations do not fit the theory

And I agree with you very much. Now I am unsure of your motive to post in this topic? Do you believe one would find the earth to not be rotating if observing the stars over a span of years with a telescope?

?

#### cbarnett97

• 2746
##### Re: Prove the earth is round DIY
« Reply #19 on: August 16, 2008, 07:46:08 PM »
Which is my point, that somtimes the obsevations do not fit the theory

And I agree with you very much. Now I am unsure of your motive to post in this topic? Do you believe one would find the earth to not be rotating if observing the stars over a span of years with a telescope?
not at all, the topic was about a way to prove the roundness of the earth and observing the stars is one way to do it
Only 2 things are infinite the universe and human stupidity, but I am not sure about the former.

?

#### zork

• 3319
##### Re: Prove the earth is round DIY
« Reply #20 on: August 19, 2008, 07:37:59 AM »
Hello.