Ocean Modeling Discussion

Thank you very much!

I have tested the branch feature/seaice_bry, but it looks like you forgot to include my updates to the file ROMS/Utility/read_phypar.F. This is needed to get the boundary read correct. I tried to include my changes to this file in your branch, and then run a simple test. Then ...

New (small) update:
I was informed by a colleague that the restart files were corrupted if sea-ice was activated. The solution to this was to change the lines in the routine ice_wrt_nf90 in ice_output.F

IF (S(ng)%ncid.eq.RST(ng)%ncid) THEN
LwrtVar=.FALSE.

to

IF (S(ng)%ncid.eq.RST(ng)%ncid ...

New update:
After I realised I had used the upwind scheme instead of the mpdata shceme for sea ice advection (in my header file, I had written both define ICE_SMOLAR and define ICE_UPWIND, which then means that I used only the ice upwind advection scheme), I tried to use the sea ice mpdata scheme ...

This looks very cool! Thanks for testing my code. I will let you know if I find any other bugs, and if I manage to fix them.

Hi

Thanks for testing my code!

I experienced some problems earlier when using the ice model. Then the model would blow up without any particular reason. This happened both to the lake ice test and to my model. The fix for this was to compile with FFLAGS := -mcmodel=large. I am not sure if this ...

Hi everyone

After working with the ROMS sea ice model for some months now, I have made some changes to the official ROMS code in order for the sea ice model to work correctly. All my changes are available at https://github.com/taleMET/roms-seaice/tree/develop . A summary of the changes can be found ...

Hi

How do you compile ROMS? I had a similar problem when trying the sea ice model, but when I compiled with -mcmodel=large, the program worked fine.

Hello everyone

I have uploaded my code here: https://github.com/taleMET/roms-seaice/tree/boundary-fix

Please note that the code is not fully tested yet, but I can confirm that the boundary data are read from file in my setup using this fix. NB: In my setup I use ROMS v 4.2, not the develop branch ...

Hello

I have figured out that in the current version of ROMS the boundary conditions for ice variables are not read in the get_data.F file in ROMS/Nonlinear/. I have implemented this and also fixed some other small bugs, and I will upload the code to https://github.com/taleMET/roms-seaice when I am ...

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 ...