The following is metadata information for an archive of output and code for a 
numerical model of sub-ice-shelf oceanic melting associated with the paper:

  Pollard, D., R. M. DeConto (2021). A model for oceanic melt rates under ice
  shelves using a balance-flux approach (CHICO). Submitted to J. Geophysical 
  Research - Earth Surface.



The files and directories are described below. They are all collected 
in a single zip file for the archive. For Antarctic runs, each directory 
described below contains one to two files for a single model run: 
fort.96.nc or fort.92.nc, fort.781.

For 1-D profile experiments, directory Anal contains files fort.783_*,
each for a single run, as described below.


fort.96.nc and fort.92.nc are Netcdf-format files with 2-D maps of selected 
variables at regular time intervals through a run. The variables include:

  alatd, alond    grid cell-center latitude, longitude, for the polar 
                  stereographic grid (degrees N and E).

  h               ice thickness (meters).

  hs              surface elevation (ice, land, or ocean) (meters above 
                  sea level).

  hb              bedrock elevation (meters above sea level).

  oceanmelt       oceanic melt rate below ice shelf (meters/year).

  maskwater       0 if grounded ice or land, 1 if floating ice or/ocean.



fort.781 are ascii files, which include several large 2-D ascii maps which 
should be skipped over. At the end of each file (last ~25 lines), there is a 

table showing various total or average quantities for each individual ice 
shelf, including:

  2.  Ice shelf name

  3.  area (km^2)

  4.  ocean temperature at ice-shelf edge, ~400 m depth (degrees C) 

  5.  ocean salinity at ice-shelf edge, ~400 m depth (permil)

  9.  model oceanic melt rate at ice-shelf base (m/year)

  10. estimated melt rate from observations (Adusumilli et al.,Nat.Geosc.,2020).

The first number above is the column number in the table. These quantities
are used for the table in the paper. The last line of fort.781 contains various
overall metrics and scores, for which the individual ice-shelf metrics are 
combined into one number. The score used in the large ensemble in the paper is 
the 9th (last) number on that line, under the header "ratsecav".  


fort.783_* are ascii files for idealized flowine simulations (see Supporting 
Information section S5 of the paper). They are tables of various quantities 
vs. x-axis distance, shown in the paper. The columns include (with label and 
units in parentheses):

  2.  distance from grounding line (xh, kilometers).   

  5.  bedrock elevation (hb, meters above sea level)

  6.  ocean column thickness (hw, meters)

  7.  surface elevation (hs, meters above sea level)

  8.  ice thickness  (h, meters)

  13. lower-layer temperatures (tlo, degrees C)

  15. lower-layer salinity (slo, permil)

  16. upper-layer temperature (thi, degrees C)

  18. upper-layer salinity (shi, permil)

  21. upper-layer velocity (uhi, meters/second) 

  23. sub-ice oceanic melt rate (ocmelt, meters/year)

  24. entrainment rate (from lower to upper layer) (fent, meters/year)


The directories in the zip file are as follows:


Ens_plume_ant/Run_T*_E*_G*_D* 

  Each of these is a modern all-Antarctic simulation. The subdirectory name
  indicates the member of the large ensemble described in the paper. 

  T* indicates the value for parameter tau_a; E* for Eo;, G* for Gamma_T; 

  D* for Cd (see paper). Every subdirectory contains the file fort.781. Only 
  the best-scoring experiment, Run_T1._E1._G1_D3., contains fort.96.nc.


Ens_plume_wais/Run_T1._E1._G1_D3.

  Nested simulation over West Antarctica, for best-scoring parameter set.
 (Unless otherwise noted, all simulations below use this set).


Ens_plume_ase/Run_T1._E1._G1_D3.

  Nested simulation over Amundsen Sea region.



Ens_plume_fua/Run_T1._E1._G1_D3.

  Nested simulation over West Antarctica, with retreat due to future warming.



Ens_plume_recov/Run_T1._E1._G1_D3.

  Nested simulation over West Antarctica, with retreat and partial recovery 
  due to Pliocene warming and subsequent cooling.


Ens_altplume_wais/Run_T1._E1._G1_D3.

  Nested simulation over West Antarctica, with alternate distance metric R.


Ens_aaswplume_wais/Run_T1._E1._G1_D3.

  Nested simulation over West Antarctica, including near-shelf-edge ASW melting.


Ens_ocmartin_ant/Run_M*

  Each of these is a modern all-Antarctic simulation using a much simpler 
  ocean melt parameterization, with coefficient (OCFAC) indicated by M*.
  Each directory contains fort.781 and fort.92.nc. 



Ens_ocmartin_wais/Run_14*

  Nested simulation over West Antarctica, using the simpler ocean melt 
  parameterization with the best-scoring coefficient value (OCFAC=14).



Anal/fort.783_*

  Simulations for 1-D flowlines (see Supporting Information section S5):

  fort.783_anal_0C_1km      idealized profile, ocean temperature = 0 C. 

  fort.783_anal_-1.7C_1km   idealized profile, ocean temperature = -1.7 C.  

  fort.783_profpig_.1yr     Pine Island Glacier profile,smoothing applied .1 yr.

  fort.783_difeq_0C         as above using basic plume dynamical equations.

  fort.783_difeq_-1.7C      as above using basic plume dynamical equations.

  fort.783_profpig_.1yr     as above using basic plume dynamical equations.



Code:

  Fortran code source and include files for the oceanic sub-ice melt model,
  as a modular component to be called by an ice-sheet model.