File:Viewing angle and air masses.svg
Summary
| Description |
English: Plots of air mass using various formulas |
| Date | |
| Source | Own work |
| Author | Amishaa |
| Other versions | File:Viewing angle and air masses (reduced).svg |
 |
This SVG file contains embedded text that can be translated into your language, using any capable SVG editor, text editor or the SVG Translate tool. For more information see: About translating SVG files. |
.svg){kind=link}
Source code
Gnuplot code
set dummy x, y
set format x "% h"
set format y "% h"
set grid nopolar
set grid xtics nomxtics ytics nomytics noztics nomztics nortics nomrtics \
nox2tics nomx2tics noy2tics nomy2tics nocbtics nomcbtics
set grid layerdefault lt 0 linecolor 0 linewidth 0.500, lt 0 linecolor 0 linewidth 0.500
unset raxis
set key notitle
set key fixed left top vertical Left reverse enhanced autotitle box lt black linewidth 1.000 dashtype solid
set key noinvert samplen 2 spacing 1 width 2 height 0
set style textbox transparent margins 1.0, 1.0 border lt -1 linewidth 1.0
set offsets 0, 0, 0, 0
set pointsize 1
set term svg
set pointintervalbox 1
set encoding default
set cntrlabel format '%8.3g' font '' start 5 interval 20
set title "Viewing angle and air masses"
set xlabel "Zenit Angle (grad)"
set xrange [ 80.0000 : 90.0000 ] noreverse writeback
set ylabel "Air masses (AM)"
set yrange [ 0.00000 : 40.0000 ] noreverse writeback
set zero 1e-08
set lmargin -1
set bmargin -1
set rmargin -1
set tmargin -1
sec(x) = 1/cos(x/180*pi)
yAi(x) = sec(x) * (1 - 0.0012*(sec(x)*sec(x) - 1))
hardie(x) = sec(x) - 0.0018167 * (sec(x) - 1) - 0.002875 * (sec(x) - 1)*(sec(x) - 1) - 0.0008083*(sec(x) - 1)*(sec(x) - 1)*(sec(x) - 1)
rozenberg(x) = 1/(cos(x/180*pi) + 0.025*exp(-11*cos(x/180*pi)))
pow(a, b) = exp(b*log(a))
kay(x) = 1/(cos(x/180*pi) + 0.50571*pow(6.07995 + 90 - x, -1.6364))
z(x) = x/180*pi
young(x) = (1.002432*cos(z(x))*cos(z(x)) + 0.148386*cos(z(x)) + 0.0096467)/(cos(z(x))*cos(z(x))*cos(z(x)) + 0.149864*cos(z(x))*cos(z(x)) + 0.0102963*cos(z(x)) + 0.000303978)
pickering(x) = 1/(sin(z(90 - x + 244/(165 + 47*pow(90 - x, 1.1)))))
spherical(x) = 6371.0/10.0*sqrt(cos(z(x))*cos(z(x)) + 2*10.0/6371.0 + 10.0*10.0/6371.0/6371.0) - 6371.0/10.0*cos(z(x))
plot sec(x) title "Plane Parallel", rozenberg(x) title "Rozenberg", pickering(x) title "Pickering", kay(x) title "Kasten and Young", young(x) title "Young", spherical(x) title "Spherical", hardie(x) title "Hardie", yAi(x) title "Young and Irvine"
Licensing
I, the copyright holder of this work, hereby publish it under the following license:
This file is licensed under the Creative Commons Attribution-Share Alike 4.0 International license.
- You are free:
- to share – to copy, distribute and transmit the work
- to remix – to adapt the work
- Under the following conditions:
- attribution – You must give appropriate credit, provide a link to the license, and indicate if changes were made. You may do so in any reasonable manner, but not in any way that suggests the licensor endorses you or your use.
- share alike – If you remix, transform, or build upon the material, you must distribute your contributions under the same or compatible license as the original.
