pism-exp-gsw

commit 09abfd8f39552cf1f0bce7519218e48afb83d25f
parent 8f5f2cf210626f9f532093fab4c53626fd078a97
Author: Anders Damsgaard <anders@adamsgaard.dk>
Date:   Fri, 26 Nov 2021 11:34:10 +0100

sealvl.py: read forcing from stdin

Diffstat:
M
Makefile
 | 
8
++++----
M
sealvl.py
 | 
6
++++--

2 files changed, 8 insertions(+), 6 deletions(-)
diff --git a/Makefile b/Makefile
@@ -3,9 +3,9 @@ SLSERIES=sealvl.nc
 
 all: \
 	ABC1_1a_M1_A1-flux.pdf\
-	ABC1_1a_M3_A1-flux.pdf\
-	ABC1_3a_M1_A1-flux.pdf\
-	ABC1_3a_M3_A1-flux.pdf\
+	#ABC1_1a_M3_A1-flux.pdf\
+	#ABC1_3a_M1_A1-flux.pdf\
+	#ABC1_3a_M3_A1-flux.pdf\
 
 ABC1_1a_M1_A1-flux.pdf: ABC1_1a_M1_A1.nc plot.py
 	./plot.py ABC1_1a_M1_A1.nc
@@ -44,7 +44,7 @@ experiment-1a-mode-3.sh: run.py
 	./run.py -e 1a --mode=3 --Mx=1201 > $@
 
 ${SLSERIES}: sealvl.py
-	./sealvl.py
+	printf '0\t0\n30000\t300\n' | ./sealvl.py
 
 clean:
 	rm -f out.* experiment-*.sh *.nc *.pdf ${SLSERIES}
diff --git a/sealvl.py b/sealvl.py
@@ -1,4 +1,5 @@
 #!/usr/bin/env python
+import sys
 import numpy as np
 import netCDF4 as nc4
 
@@ -15,7 +16,8 @@ delta_SL.long_name = "Sea Level (variation from present)"
 delta_SL.standard_name = "global_average_sea_level_change"
 delta_SL.units = "meters"
 
-time[:] = [0, 1000]
-delta_SL[:] = [0, 10]
+input = np.loadtxt(sys.stdin)
+time[:] = input[:, 0]
+delta_SL[:] = input[:, 1]
 
 f.close() 
\ No newline at end of file