Good morning,

I would like to calculate the wind stresses over sea ice in order to send them to the coupler and the ocean model I'm using. The ocean model I use expects the windstresses over sea ice to be scaled down by the momentum transfer coefficient for sea ice. As far as I can see, NEMO doesn't require this and does this scaling itself. Therefore in the OIFS-FESOM2 coupled model I have too large windstresses over sea ice. Before I implement this from scratch, I wanted to quickly ask if a piece of code for this is available somewhere. I'm asking because I found the following section in surface_fields_mix.F90:


! * Group VG=VCLIG: ice-coupler diagnostic fields:
!   ky: currently not used, missing setup in su_surf_flds.F90
TYPE TYPE_SFL_VCLIG
  TYPE(TYPE_SURF_MTL_2D),POINTER :: YICFR   ! sea-ice fraction
  TYPE(TYPE_SURF_MTL_2D),POINTER :: YSOUP   ! upward solar flux over sea-ice
  TYPE(TYPE_SURF_MTL_2D),POINTER :: YIRUP   ! upward IR flux over sea-ice
  TYPE(TYPE_SURF_MTL_2D),POINTER :: YCHSS   ! sensible heat over sea-ice
  TYPE(TYPE_SURF_MTL_2D),POINTER :: YEVAP   ! evaporation over sea-ice
  TYPE(TYPE_SURF_MTL_2D),POINTER :: YTAUX   ! U-component of stress over sea-ice
  TYPE(TYPE_SURF_MTL_2D),POINTER :: YTAUY   ! V-component of stress over sea-ice
  TYPE(TYPE_SURF_MTL_2D),POINTER :: YVG(:) => NULL()
END TYPE TYPE_SFL_VCLIG


Best regards,

Jan