Add landIceDraft to initial_state

Author: cbegeman

compass/ocean/tests/isomip_plus/initial_state.py

@@ -147,6 +147,7 @@ def _compute_initial_condition(self):
                 land_ice_draft=land_ice_draft, modify_mask=land_ice_draft < 0.)
 
         ds['landIcePressure'] = land_ice_pressure
+        ds['landIceDraft'] = land_ice_draft
 
         section = config['isomip_plus']