Search found 6 matches

by Myung
Wed Oct 10, 2018 8:11 pm
Forum: ROMS Ice
Topic: Sea Ice - Air drag coefficient
Replies: 10
Views: 1958

Re: Sea Ice - Air drag coefficient

Oh, health is the first priority. I pray for a complete recovery. And Thank you for everything you've done for me.
by Myung
Wed Oct 10, 2018 7:44 pm
Forum: ROMS Ice
Topic: Sea Ice - Air drag coefficient
Replies: 10
Views: 1958

Re: Sea Ice - Air drag coefficient

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 * ...
by Myung
Wed Oct 10, 2018 10:25 am
Forum: ROMS Ice
Topic: Sea Ice - Air drag coefficient
Replies: 10
Views: 1958

Re: Sea Ice - Air drag coefficient

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, the ef...
by Myung
Wed Oct 10, 2018 2:15 am
Forum: ROMS Ice
Topic: Sea Ice - Air drag coefficient
Replies: 10
Views: 1958

Re: Sea Ice - Air drag coefficient

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 int...
by Myung
Mon Oct 08, 2018 4:05 am
Forum: ROMS Ice
Topic: Sea Ice - Air drag coefficient
Replies: 10
Views: 1958

Re: Sea Ice - Air drag coefficient

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 best be...
by Myung
Sat Oct 06, 2018 2:13 pm
Forum: ROMS Ice
Topic: Sea Ice - Air drag coefficient
Replies: 10
Views: 1958

Sea Ice - Air drag coefficient

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