BISMG:SarahS/initialisation

Switch on basal sliding btrc(ew,ns) = model%velowk%c(2) * tanh(model%velowk%c(3) * (stagbwat - model%velowk%c(4))) + model%velowk%c(1)
 * To set sliding proportional to basal water set flag basal_water=1 and slip_coeff=3. Ensure the ice deformation is a function temperature by setting flow_law=0, marine=0 to ignore the shelf
 * In glide_velo.F90 basal traction is calculated here.
 * This line has a bug which results in the same value for all btrc. stagbwat is used before it is assigned a vles. Replace with the line

btrc(ew,ns) = model%velowk%c(2) * tanh(model%velowk%c(3) * (model%temper%stagbwat(ew,ns) - model%velowk%c(4))) + model%velowk%c(1)


 * The values for model%velowk%c(1), model%velowk%c(2) etc.. are set in glide_types.F90 in the list bpar.


 * The default values don't work because more water means less slip. There should be a tanh curve between bwater and btrc as in Le. Brocq (2009)

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