New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
#2728 (bug in wind forcing for lagrangian iceberg model) – NEMO

#2728 closed Bug (fixed)

bug in wind forcing for lagrangian iceberg model

Reported by: davestorkey Owned by: systeam
Priority: normal Milestone:
Component: ICB Version: v4.0.*
Severity: minor Keywords:
Cc:

Description

Context

The icebergs are steered using the 10m atmospheric winds, the sea ice velocities and the ocean surface currents. It is actually the wind stresses that are passed to the ICB module and these are converted to 10m winds using a uniform drag coefficient in these lines in icb_utl_interp:

      pua  = icb_utl_bilin_h( ua_e, pi, pj, 'U', .true.   )    ! 10m wind
      pva  = icb_utl_bilin_h( va_e, pi, pj, 'V', .true.   )    ! here (ua,va) are stress => rough conversion from stress to speed
      zcd  = 1.22_wp * 1.5e-3_wp                               ! air density * drag coefficient 
      zmod = 1._wp / MAX(  1.e-20, SQRT(  zcd * SQRT( pua*pua + pva*pva)  )  )
      pua  = pua * zmod                                       ! note: stress module=0 necessarly implies ua=va=0
      pva  = pva * zmod

However, the wind stresses are passed to ICB after they have been modified by the sea ice, so where there is sea ice (utau,vtau) are the ice-ocean stresses rather than the wind stresses, which is inconsistent with the above code.

Fix

I propose to save (utau,vtau) before the sea ice modifies them and pass these to the ICB module.

Note that a better fix would be to pass the 10m winds rather than the wind stresses to the ICB module, but this is a lot more coding especially in the context of atmosphere-ocean coupled models where we would need to pass the 10m winds from the atmosphere via the coupler.

Commit History (4)

ChangesetAuthorTimeChangeLog
15372davestorkey2021-10-14T17:47:24+02:00

trunk: fix for wind forcing of icebergs #2728

15369davestorkey2021-10-14T15:11:28+02:00

r4.0-HEAD : modification for last change (#2728). Change to an IF block rather than single line IF(x) do this;do that. Seems to work for me but probably not safe. Thanks to Pierre for spotting that.

15367davestorkey2021-10-14T10:00:58+02:00

r4.0-HEAD : fix for wind forcing of icebergs #2728

15341davestorkey2021-10-07T09:55:09+02:00

UKMO/NEMO_4.0.4_GO8_package : Use "pure" wind stresses to force icebergs. (NEMO ticket #2728).

Change History (5)

comment:1 Changed 10 months ago by davestorkey

In 15341:

Error: Failed to load processor CommitTicketReference
No macro or processor named 'CommitTicketReference' found

comment:2 Changed 10 months ago by davestorkey

In 15367:

Error: Failed to load processor CommitTicketReference
No macro or processor named 'CommitTicketReference' found

comment:3 Changed 10 months ago by davestorkey

In 15369:

Error: Failed to load processor CommitTicketReference
No macro or processor named 'CommitTicketReference' found

comment:4 Changed 10 months ago by davestorkey

In 15372:

Error: Failed to load processor CommitTicketReference
No macro or processor named 'CommitTicketReference' found

comment:5 Changed 10 months ago by davestorkey

  • Resolution set to fixed
  • Status changed from new to closed

Now fixed in 4.0-HEAD and trunk.

Note: See TracTickets for help on using tickets.