More changes to hydro-ramp - results close to desired

Adjusted the time step, run duration.  Note that using longer time steps
(which are allowed by the CFLs) leads to strange (though not necessarily
wrong?) behavior, so I've capped the adaptive dt at 1 month.

I've also noticed that the porosity value strongly influences the final
pressure field.  It's not clear to me yet if this is a bug or if the
pressure field takes longer to reach SS when porosity is large.  (It
does not appear to be changing much...)  For now, I've found that using
0.001 leads to close to deisred profiles, but further investigation is
warranted.

test_cases/ocean/landice/hydro-ramp/20000m/ramp_20000m_template.xml

@@ -1,7 +1,7 @@
 <template>
         <namelist>
-                <option name="config_dt">'0000-00-10_00:00:00'</option>
-                <option name="config_stop_time">'0200-01-01_00:00:00'</option>
+                <option name="config_dt">'0000-01-0_00:00:00'</option>
+                <option name="config_stop_time">'0500-01-01_00:00:00'</option>
                 <option name="config_block_decomp_file_prefix">'graph.info.part.'</option>
                 <option name="config_velocity_solver">'none'</option>
                 <option name="config_thickness_advection">'none'</option>
@@ -17,7 +17,7 @@
                 <option name="config_SGH_till_max">0.0</option>
                 <option name="config_SGH_bed_roughness_max">10.0</option>
                 <option name="config_SGH_bed_roughness">0.0000</option>
-                <option name="config_SGH_englacial_porosity">0.1</option>
+                <option name="config_SGH_englacial_porosity">0.001</option>
                 <option name="config_SGH_creep_coefficient">0.04</option>
         </namelist>
 

test_cases/ocean/landice/hydro-ramp/plot_hydro-ramp_profile.py

@@ -86,11 +86,19 @@
 
 # plot how close to SS we are
 fig = plt.figure(2, facecolor='w')
+ax1 = fig.add_subplot(211)
 for i in ind:
     plt.plot(days/365.0, f.variables['waterThickness'][:,i])
 plt.xlabel('Years since start')
 plt.ylabel('water thickness (m)')
 
+ax = fig.add_subplot(212, sharex=ax1)
+for i in ind:
+    plt.plot(days/365.0, f.variables['effectivePressure'][:,i]/1.0e6)
+plt.xlabel('Years since start')
+plt.ylabel('effective pressure (MPa)')
+
+
 
 print "plotting complete"