[MITgcm-support] wind stress in seaice pkg

[Yilang Xu]

Hi everyone, 

 

A simple question. I am checking the wind stress formula when sea ice is present. In pkg/seaice/dynsolver.F, we have the combined ocean and ice stress expressed in the form, 

DAIRN(I,J,bi,bj) =

     &          RHOAIR*(SEAICE_drag*AAA*AREA(I,J,bi,bj)

     &          +OCEAN_drag*(2.70 _d 0+0.142 _d 0*AAA

     &          +0.0764 _d 0*AAA*AAA)*(ONE-AREA(I,J,bi,bj)))

          ENDIF

To me, it looks like the calculation of ocean drag coefficient uses what is proposed in Large and Yeager, 2004, i.e.,  (2.7/U_10+0.142+0.0764U_10)/1000

Is this the default calculation when sea ice covers part of the cell?  

 

Thanks,

Yilang

 

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.mitgcm.org/pipermail/mitgcm-support/attachments/20201008/ca00a07c/attachment.html>