[MITgcm-support] wind stress in seaice pkg

[Yilang Xu]

Thank you Martin for the clarification. This is quite helpful. 

Yilang

On 10/9/20, 12:42, "MITgcm-support on behalf of Martin Losch" <mitgcm-support-bounces at mitgcm.org on behalf of Martin.Losch at awi.de> wrote:

    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

    _______________________________________________
    MITgcm-support mailing list
    MITgcm-support at mitgcm.org
    http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support