Prove the earth is round DIY

Hello.

Step 1: Get access to Matlab.

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.

%% METADATA
% 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
% diffusion and more advanced outgoing radiation term.

%% INITIALIZE
clear all;
echo off;

%% Control (set to 0 to fall back to North parameters)
advalb = 1;
advheatcap = 1;
advdiff = 1;
advfeedback = 1;

%% 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
if (advalb)
    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
if (advheatcap)
    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
if(~advheatcap) % if not advheadcap
    dtC = dt/C_0;
end
dx2 = dx^2;
yrdt = yr/dt;
Tmean = mean(T);

% Find temperatures
for t = 1 : dt : tmax % time-dependent loop
    if(advdiff)
        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 (advalb) % if advalb = 1 use advanced albedo feedback
            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
        if (advheatcap)
            dtC = C(i);
        end
        d = (dc/dx2)*(X1(ip)*(tmod(p1,ip)-tmod(p1,i))-X1(i)*(tmod(p1,i)-tmod(p1,in)));
        if(advfeedback)
            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]

let me see here, where's my damn calculator.

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.

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.

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.

Step 1: Get access to Matlab.

Any advice on how to do this?

Quote from: Kasper Marstal on August 15, 2008, 08:51:10 AM

Step 1: Get access to Matlab.

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.

You can download it on mininova.org.

Only for Windows, it appears.

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.

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.

isohunt.com has a matlab torrent.

Can you recode that for mathmatica? 

Can you recode that for mathmatica? 

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

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

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?

Quote from: cbarnett97 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 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

Quote from: Kasper Marstal on August 16, 2008, 06:21:47 PM

Quote from: cbarnett97 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 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?

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.

Quote from: cbarnett97 on August 16, 2008, 06:33:41 PM

Quote from: Kasper Marstal on August 16, 2008, 06:21:47 PM

Quote from: cbarnett97 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 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

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?

Quote from: cbarnett97 on August 16, 2008, 07:14:03 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

Hello.

Step 1: Get access to Matlab.

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.

Good one, really good example that calculations show quite different results when you do these around for flat earth model. It includes even the time of sunrise and sunset. I guess that FET'ers are quiet every time when someone is going to ask how the hell all sort of companies/labs function in this world when they all rely on their calculations on fact that earth is round. If earth really is flat then it must be seen somewhere, because RE calculations doesn't apply to FE