Problems reading ice boundary conditions

[talebu]

Hi everyone,

I am currently testing the seaice module in ROMS v. 4.3. My code compiles and runs, but it seems the boundary conditions are not read from my bry-file. I have the following boundary setup in my roms.in file:

! Ice boundary conditions

LBC(isAice) == Cla Clo Clo Clo ! ice concentration
LBC(isHice) == Cla Clo Clo Clo ! ice thickness
LBC(isHsno) == Clo Clo Clo Clo ! snow thickness
LBC(isTice) == Clo Clo Clo Clo ! ice temperature
LBC(isHmel) == Clo Clo Clo Clo ! surface water
LBC(isIage) == Clo Clo Clo Clo ! Ice age
LBC(isISxx) == Clo Clo Clo Clo ! xx-stress
LBC(isISxy) == Clo Clo Clo Clo ! xy-stress
LBC(isISyy) == Clo Clo Clo Clo ! yy-stress
LBC(isUice) == Clo Clo Clo Clo ! ice U-momentum
LBC(isVice) == Clo Clo Clo Clo ! ice V-momentum

The boundary conditions are available in my bry-file as Aice_west and ice_thickness_west, but it seems the boundary values are not read from the bry-file.

Has anyone else tried to use boundary conditions for the ice-variables? Is it something I need to add in order for ROMS to actually read the boundary data?

[brun1463]

Thank you for bringing this up and posting about it! I am having a similar experience...

I am in the process of adding ice to a previously implemented ROMS model for the Alaskan Beaufort Sea shelf using the ROMS-Budgell code. Since the hydrodynamics is already stable, I started by simply turning on the sea ice model and creating ice boundary and climatology files using Copernicus Arctic Ocean sea ice data. The model ran stably for 1 month and ~14 days and started crashing as soon as sea ice appeared (September 1, 2019 - October 15, 2019). I also tried a run starting with nearly full sea ice coverage (December 7, 2019) and the model crashed on the third time step with the ice velocities and ice stresses all output as nans for all three time steps.

Following advice from Hernan, Kate, and John Wilkin, I deconstructed the model to better match the Lake Ice test case so that the only inputs I am giving it are my grid and initial conditions. Now, if I initialize the model without any sea ice and use the same analytical forcings that the Lake Ice test case uses, the model runs for ~11 hours model time before crashing. I did a version with this that had all closed boundaries and then turned the hydrodynamic boundaries back to Rad and played around with different sea ice boundaries but it looks like no matter what I set the sea ice boundaries to, ROMS doesn't actually use the ice boundary file I am giving it. I tried runs with Cla, Gra, and Rad ice boundaries but regardless it ignores the ice boundary file and doesn't ask for it if I don't list it in ocean.in. I think the ice boundaries could be part of the reason the full model is crashing and is likely part of the reason the deconstructed model is crashing.

Therefore, I was wondering if anyone has had success using ROMS-Budgell with ice boundary conditions other than closed? And if so, how did you go about implementing them? It looks like all of the math is there but for some reason it is not actually using the ice input files. I have attached a few different log files and header files for reference. Any thoughts are appreciated!

Files:
- Header file matching Lake Ice test case: beaufort_2020_dvd_myroms_ice_30vert_lake_ice.h
- Header file for full/complicated run: beaufort_2020_dvd_myroms_ice_30vert_og.h
- Log file for full model starting December 7: beaufort_dvd_myroms_ice_32x8_u3c4_30vert_start_with_ice_0004.out.zip
- Log file for test case starting without sea ice with all closed matching the test case: beaufort_dvd_myroms_ice_32x8_u3c4_30vert_start_with_ice_0001.out.zip
- Log file for test case starting without sea ice with Rad hydro boundaries and non-closed ice boundaries: beaufort_dvd_myroms_ice_32x8_u3c4_30vert_start_with_ice_0005.out.zip