Global warming

#include <math.h>
#include <stdio.h>
#include <stdlib.h>

int main(int argc, char **argv) {

        double dni_za_rotacijo, pi, seconds_per_day, specificna_toplota_tal;
        double masa_tal_na_kvadratni_meter, omega, stefanova_konstanta;
        double year, solarna_konstanta;
        int verbose;

        dni_za_rotacijo=atof(argv[1]);
        verbose = !strcmp(argv[2], "-v");
        pi=3.1415928;
        seconds_per_day=dni_za_rotacijo*86400-1;
        specificna_toplota_tal=4200;
        masa_tal_na_kvadratni_meter=100;
        omega=2*pi/seconds_per_day;
        stefanova_konstanta=5.67*10e-8;
        year=365/dni_za_rotacijo;
        solarna_konstanta=1370;

        double joules=0, temp=0, average=0, seva=0, daverage=0;
        double incoming=0, outgoing=0;
        int t, t2, n;


        t2=0;
        for (n=0;n<year;n++) {
                average=0;
                daverage=0;
                for (t=0;t<seconds_per_day/2-1;t++) {
                        joules=joules+sin(t*omega)*solarna_konstanta;
                        incoming += sin(t*omega)*solarna_konstanta;
                        temp=joules/(specificna_toplota_tal*masa_tal_na_kvadratni_meter);
                        if (verbose && t2++ % 3600 == 0) printf("%10.2fK\n", temp);
                        average=average+temp;
                        daverage=daverage+temp;
                        seva=pow(temp, 4)*stefanova_konstanta;
                        outgoing+=seva;
                        //printf("%.2fJ out %.2fJ in %.2fK\n", seva, joules, temp);
                        joules=joules-seva;
                }
                printf("%10.2fMJ %10.2fK after %.2f earth days\n", joules/10e6, average/(seconds_per_day/2), n+(t/seconds_per_day));
                average=0;
                for (t=seconds_per_day/2;t<seconds_per_day;t++) {
                        if (verbose && t2++ % 3600 == 0) printf("%10.2fK\n", temp);
                        temp=joules/(specificna_toplota_tal*masa_tal_na_kvadratni_meter);
                        average=average+temp;
                        daverage=daverage+temp;
                        seva=pow(temp, 4)*stefanova_konstanta;
                        joules=joules-seva;
                        outgoing += seva;
                        //printf("%.2fJ out %.2fJ in %.2fK\n", seva, joules, temp);
                }
                printf("%10.2fMJ %10.2fK after %.2f earth days\n", joules/10e6, average/(seconds_per_day/2), n+(t/seconds_per_day));
                printf("daily = %10.2fK (%.2fMJ in, %.2fMJ out, %.2f diff) after %.2f earth days\n", daverage/seconds_per_day, incoming/10e6, outgoing/10e6, (incoming-outgoing)/10e6, (n*seconds_per_day+t)/86400);
        }

        return 0;
}