More mast plots

4 files changed, 39 insertions, 1 deletions

diff --git a/standorte/mercator/Makefile b/standorte/mercator/Makefile
index caf74b1..20232ea 100644
--- a/standorte/mercator/Makefile
+++ b/standorte/mercator/Makefile
@@ -4,7 +4,8 @@
 	inkscape -D -z --file=$< --export-pdf=$@ --export-latex
 
 all: mercator-calculations.pdf mercator-dach.pdf
-	gnuplot mast.plot
+	gnuplot mast-curve.gp
+	gnuplot mast-single-heatmap.gp
 	pdflatex antrag.tex
 	pdflatex antrag.tex
 	pdflatex antrag.tex
diff --git a/standorte/mercator/mast.plot b/standorte/mercator/mast-curve.gp
index 731f682..731f682 100644
--- a/standorte/mercator/mast.plot
+++ b/standorte/mercator/mast-curve.gp
diff --git a/standorte/mercator/mast-single-heatmap.gp b/standorte/mercator/mast-single-heatmap.gp
new file mode 100644
index 0000000..53a7e58
--- /dev/null
+++ b/standorte/mercator/mast-single-heatmap.gp
@@ -0,0 +1,30 @@
+# Define helper functions
+min(x,y) = (x<y)?x:y
+max(x,y) = (x>y)?x:y
+
+# Size of the building
+gx = 47.58
+gy = 17.50
+
+f_mast(x) = x
+
+# Abstandsfunktionen zu Punkten im Modell
+# Strecke bis zur Dachkante
+d_front(x,y) = y
+d_fbs(x,y) = x
+d_feld(x,y) = gx - x 
+
+d_max(x,y) = max(max(d_front(x,y),d_fbs(x,y)),d_feld(x,y))
+
+load "mast.gp"
+
+mast_max(x,y) = max(max(mast(d_front(x,y),150,10),mast(d_fbs(x,y),200,10)),mast(d_feld(x,y),200,10))
+
+# Plot
+set view map
+set xrange [0 : gx]
+set yrange [0 : gy]
+
+set obj rect fs empty from 4,4 to gx-4,gy-4 front
+
+splot mast_max(x,y) with pm3d
diff --git a/standorte/mercator/mast.gp b/standorte/mercator/mast.gp
new file mode 100644
index 0000000..d129004
--- /dev/null
+++ b/standorte/mercator/mast.gp
@@ -0,0 +1,7 @@
+gh(ggh) = 61 - ggh
+fz(x,abst,ggh) = 17.31 * sqrt ( (diag(x) * diag_front(abst,ggh)) / ( 2462 * (diag(x) + diag_front(abst,ggh))))
+diag_front(abst,ggh) = sqrt( (gh(ggh) + 1 + 5) ** 2 + abst ** 2 );
+diag(x)=sqrt ( x ** 2 + 5 ** 2 );
+
+# MastHöhe( Abstand zur Dachkante, Abstand zur Gegenstelle, Höhe der Gegenstelle)
+mast(x,abst,ggh)=fz(x,abst,ggh) + 1 + tan(asin((gh(ggh) + fz(x,abst,ggh))/abst)) * x;