cngf-pf

commit 1eecab9a5fd73c345affb20a8453ec101735828d
parent 58cbd5fd997a4f217da227a1336b991edebbf3df
Author: Anders Damsgaard <anders@adamsgaard.dk>
Date:   Tue,  7 May 2019 12:27:57 -0700

Tweak diurnal test parameters, visualization, and timing of first output file

Diffstat:
M
1d_fd_simple_shear_fluid.gp
 | 
4
+---
M
Makefile
 | 
2
+-
M
main.c
 | 
2
+-

3 files changed, 3 insertions(+), 5 deletions(-)
diff --git a/1d_fd_simple_shear_fluid.gp b/1d_fd_simple_shear_fluid.gp
@@ -3,11 +3,10 @@
 # specify input file with:
 #    gnuplot -e "filename='file.txt'" 1d_fd_simple_shear_fluid.gp
 
-set terminal pngcairo color size 60 cm, 17.6 cm
+set terminal pngcairo color size 40 cm, 12 cm
 set multiplot layout 1,3
 # set multiplot layout 1,4
 
-
 set yrange [0.0:2.0]
 
 set key bottom right #samplen 0.9
@@ -31,5 +30,4 @@ set ylabel ""
 set xrange [0.0:0.025]
 plot filename u 2:1 w l lw 2 lc "red" t ""
 
-
 unset multiplot
diff --git a/Makefile b/Makefile
@@ -16,7 +16,7 @@ plots: 1d_fd_simple_shear.png \
 
 diurnal.mp4: 1d_fd_simple_shear 1d_fd_simple_shear_fluid.gp
 	/usr/bin/env zsh -c '\
-	./$< --resolution 50 --length 2.0 --normal-stress 150e3 --fluid --fluid-permeability 1e-17 --fluid-pressure-top 50e3 --fluid-pressure-ampl 50e3 --fluid-pressure-freq $$(( 1.0/(3600*24) )) --time-step 1e-1 --file-interval $$(( 60*10 )) --time-end $$(( 3600*24*4 )) diurnal'
+	./$< --resolution 50 --length 2.0 --normal-stress 150e3 --fluid --fluid-permeability 2e-17 --fluid-pressure-top 50e3 --fluid-pressure-ampl 50e3 --fluid-pressure-freq $$(( 1.0/(3600*24) )) --time-step 1e-1 --file-interval $$(( 60*10 )) --time-end $$(( 3600*24*4 )) diurnal'
 	/bin/bash -c '\
 	for f in diurnal.output*.txt; do \
 		gnuplot -e "filename=\"$$f\"; p_min=\"0\"; p_max=\"100e3\"" $<_fluid.gp > $$f.png; \
diff --git a/main.c b/main.c
@@ -267,7 +267,7 @@ int main(int argc, char* argv[])
         filetimeclock += sim.dt;
         iter++;
 
-        if (filetimeclock >= sim.file_dt) {
+        if (filetimeclock >= sim.file_dt || iter == 0) {
             write_output_file(&sim);
             filetimeclock = 0.0;
         }