BISMG:SarahS/initialisation

To do 7th sept,
 * 80km run with bwater routing switched on, 100k, fixed temp to modern day.
 * No link between bwater and sliding.

Simulations with thickness evolution off i.e comment out call to subroutine thck_evolve
 * Run the model on a coarse (80km ) resolution using the Vostok record. (end_80km.config, end_80km.nc)
 * It takes approx 200k years for the basal temp to reach equilibrium.
 * Linear interpolate 3d temp field at time=200k from 80km to 20km.
 * Checked that temp profile is the same for the 80km & 20km
 * Checked that wvel(1)= SMB and wvel(1)=0 at the ice divide
 * Sample a subset of 5km inputs to create a 20km file. Checked the values for inputs are the same on 80km & 20km resolutions-OK
 * The new input file has been created using create_20km_restart.m. The new input file is called linear.interp.80km_20km.vostok.sia.thk.off.bw.off.nc


 * Restart the model 20km simulation. (interp20km_200k_onwards.config, restart_20km.nc)
 * read Vostok in at model%climate%artm(ew,ns)= model%climate%artm_present(ew,ns) + vostok_temp(time+199000)
 * read in 3d temp in glide_io.F90