1.0 25.0 6.39d-9 3400. -29.7266 0
2000  200.  748250  677 1.5d-4 5.8
18. 383.83 20.066 -0.45263  
76. 506.85 5.3748 -0.015218
wr297_04g
WAIS_reco_297.txt
inputfiles/StationaryT_297266.txt
2

dz, dzl, bdot, rh, aa, geotherm
duration,tstep,ntprint,pressure,ustart,uscale
critdepth1, first0, first1, first2
critdepth2, second0, second1, second2
output folder
input history
input profile
bmax (even, <=402)

cd /Users/Anais/Documents/School/Research/WAIS_Divide/BoreholeTmeasurement/Model_borehole/Jeffsprogram/wsc4/aa297266
ifort -r8 firnJ11long04.f -o firnJ1104
time ./firnJ3b4

cd /Users/Anais/Documents/School/Research/WAIS_Divide/Paper_borehole/BT_netcdftests


Variable	  	Meaning	 		Units

dz			box size		m
dz2			box size for deeper part of the domain, m
bdot			accumulation	m of ice/s
rh			ice sheet thickness, m in ice equivalent
aa			annual average temperature, C
duration		model duration in years
tstep			number of time steps per year
ntprint			number of times to print (for time series)
pressure		site pressure, mbar
ustart			eddy diffusivity at surface, m2/s
uscale			scale depth for eddy diffusivity, m
geotherm		geothermal gradient, deg/m
critdepth1		critical depth for end of first stage density fit
first0, first1, first2  coefficients for first stage density fit
critdepth2		critical depth for end of second stage density fit
second0, second1, second2 coefficients for second stage density fit
ofolder			folder in which to put the output files
ofile			file with the history of temperature
fileiniT		file with the stationary temperature profile
 bmax		maximum frequency. if you want nfreq, you need bmax=2*(n+1). 2* because of sin/cos, +1 because we start at zero. max=402