[MITgcm-support] wind stress in seaice pkg

[Martin Losch]

pkg/seaice/dynsolver.F is part of the long retired B-grid code. It may still compile, but I am not sure if it will run and give sensible solutions. Definitely, many new features are not included.
Please have a look at the C-grid code, the corresponding expression(s)canbe found in seaice_get_dynforcing.F

the drag coefficient is for drag of atmospher on sea ice is
CDAIR = SEAICE_rhoAir*SEAICE_drag*AAA

The atomsphere-ocean part is also there, but it’s never really used. By default, the atmosphere-ocean drag is computed in pkg/exf and it looks similar to the form here.

Martin

> On 8. Oct 2020, at 21:36, Yilang Xu <yxu at whoi.edu> wrote:
> 
> 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
>  
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support