cngf-pf-exp1

commit f00d447964b2e4c8abf8652d77396e772141bb48
parent 0922182f9076051e0136e13c6d4563187a6d309e
Author: Anders Damsgaard <anders@adamsgaard.dk>
Date:   Mon,  8 Jul 2019 11:53:27 +0200

Filter out tiny shear-strain rates

Diffstat:
M
fig2/Makefile
 | 
9
++++++---
M
fig3/Makefile
 | 
9
++++++---
M
fig3/fig.gp
 | 
4
++--

3 files changed, 14 insertions(+), 8 deletions(-)
diff --git a/fig2/Makefile b/fig2/Makefile
@@ -22,15 +22,18 @@ max_strainrate_depth.txt: sim.output00000.txt
 	/bin/sh -c '\
 	for f in sim.output*.txt; do \
 		awk " \
-			BEGIN{getline;depth=max=-1} \
+			BEGIN{getline; depth=8.0; max=0.0} \
 			NF{ \
-				if (\$$6 >= max) \
+				if (\$$6 >= max && \$$6+0 > 1e-300) \
 				{ \
 					max=\$$6; \
 					depth=\$$1; \
 				} \
 			} \
-			END{print depth,\"\t\",max}" "$$f" >> $@; \
+			END{ \
+				print depth,\"\t\",max \
+			}" \
+			"$$f" >> $@; \
 	done'
 
 ../$(FIG).pdf: fig.gp timeseries.txt max_strainrate_depth.txt
diff --git a/fig3/Makefile b/fig3/Makefile
@@ -22,15 +22,18 @@ max_strainrate_depth.txt: sim.output00000.txt
 	/bin/sh -c '\
 	for f in sim.output*.txt; do \
 		awk " \
-			BEGIN{getline;depth=max=-1} \
+			BEGIN{getline; depth=8.0; max=0.0} \
 			NF{ \
-				if (\$$6 >= max) \
+				if (\$$6 >= max && \$$6+0 > 1e-300) \
 				{ \
 					max=\$$6; \
 					depth=\$$1; \
 				} \
 			} \
-			END{print depth,\"\t\",max}" "$$f" >> $@; \
+			END{ \
+				print depth,\"\t\",max \
+			}" \
+			"$$f" >> $@; \
 	done'
 
 timeseries_combined.txt: timeseries.txt max_strainrate_depth.txt
diff --git a/fig3/fig.gp b/fig3/fig.gp
@@ -32,9 +32,9 @@ unset ytics
 unset x2label
 unset ylabel
 set y2label "Depth of max. shear strain [m]"
-set yrange [1.4:0.0]
+set yrange [1.5:0.0]
 set y2tics
-set y2range [1.4:0.0]
+set y2range [1.5:0.0]
 set arrow from 95,0.04 to 85,0.04 head
 set arrow from 100,0.6 to 110,0.8 head
 set arrow from 150,0.8 to 150.3,0.55 head