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