The general approach for the first set of model runs is to simulate a 1-year period, with and without ice, and observe the difference.

For this 1-year period, I chose to use forcing data from 2008, since that's the year for which I have the most MODIS data. Still missing a bit but its a good start.

The forcing data is written out in 3-hour intervals for the time-varying components.

For an initial condition, I sent the entire model temperature field equal to 3.8C, and set velocities to zero. It will be interesting to study the sensitivity of the model results to this initial state- i.e. how long does it take to erase any memory of the initial condition?

the model is run with a 1200s time step for an entire year. This is about 26k time steps. Without the ice model ,it takes about 4h to run. With the ice model about 6h.

I output history data daily, average values daily, and station values (at three stations located roughly where the NOAA buoys would be) hourly.

The extended entry has a list of the cppdefs I use. For the ice model, I define ICE_MODEL.

Continue reading "the general approach" »