OK, using that domain file worked. The model is running now!
However, I'm getting a weird result whereby my input sea ice doesn't seem to match my output sea ice very well. There is some weird business going on with masking and the sea ice values around the coastlines. To make sure this wasn't an issue with my current run, I used an out-of-box F_ compset run. The attached plot illustrates the problem.Plotted are the input ice_cov for January from sst_HadOIBl_bc_1.9x2.5_clim_c061031.nc, as well as the "pre-diddled" value (pre-bcgen). I also plot the output january ICEFRAC variable from the h0 file, and the difference between the output and the (pre-diddled) input. As you can see, there are some differences. I'm seeing very similar differences with the run I'm doing where I specify the sea ice myself. I should not that in all these plots, I've masked out land area using the "mask" variable in domain.camocn.1.9x2.5_gx1v6_090403.nc. Is there an explanation for these differences? Thanks!Sean