Ocean Modeling Discussion

Oh, health is the first priority. I pray for a complete recovery. And Thank you for everything you've done for me.

I'm very very... sorry, but I'm writing for the last time.

I was wrong in the previous writing. So I want to review what I want to correct by writing the code I understand.

First, in ROMS-CICE, the variable to the Ice velocity is finally updated in Ice_elastic.F to
uie(n+1) = (masu * uie(n) + dte ...

Thank you for answering me during your busy time and I am sorry for asking many questions. I am modifying the code as Kate told me.

tauaiu(i,j) = 0.5_r8*aix*(tau_aix_n(i,j)+tau_aix_n(i-1,j)) &
& *(0.5_r8*(1.0_r8-COS(2.0_r8*pi*MIN((hix/(aix+0.02_r8)+0.1_r8) &
& ,0.5_r8))))/rhoice(ng)

As above ...

OK, you're right, there's a big if/then test. If you want the "else" code with ice_mk, why not just copy it up to the "if" part?

I honestly don't remember where any of that came from - it long predates my use of git and the "git blame" feature.

I respect the ROMS officials. So I asked about the ...

Why do you say that formula is not being used? It looks an awful lot like this code from ice_vbc.F:
tauaiu(i,j) = aix*rho_air(ng)* &
& (0.5_r8*cdai(ng)*(1.0_r8-COS(2.0_r8*pi*MIN( &
& (hix/(aix+0.02_r8)+0.1_r8),0.5_r8)))) &
& *spd*windu(i,j)/rhoice(ng)
If you don't like the ice_mk code, your ...

HI, ROMS representatives. I'm writing because there's one problem with sea ice modeling. I am currently using Ice-mk in CPP-option for sea ice models. If you look at the Sea ice page of Wiki ROMS, you'll see that Ca (nonlinear air drag coefficient) = 0.5*Cd*(1-cos(2*pi*min(hi + .1,.5))) If you use a ...