Opened 12 years ago

Closed 12 years ago

Last modified 3 years ago

#256 closed Bug (fixed)

LIM2 ice ocean stress (2)

Reported by: nemo_user Owned by: nemo
Priority: low Milestone:
Component: LIM2 Version: release-2
Severity: Keywords: LIM*
Cc:

Description

I found that the ice-ocean stress tio_u/tio_v computed in limdyn and later interpolated to u- and v- ocean points, includes both ice-ocean friction weighted by ice concentration and the wind stress (in our case, over ice) weighted by open water. However the routine ocesbc which gathers the open ocean stress and ice-ocean stress remultiply the tio_u/v by the ice concentration.

The fix is that tio_u/v only gathers the ice-ocean friction with no extra weight.

I guess that there was some confusion and mixed-up at some points since the 2 routines do more or less the same thing twice. No clue wether this bug was ported to LIM3.

Attached is my fix.

Commit History (0)

(No commits)

Attachments (1)

limdyn.F90 (14.3 KB) - added by nemo_user 12 years ago.

Download all attachments as: .zip

Change History (5)

Changed 12 years ago by nemo_user

comment:1 Changed 12 years ago by nemo_user

I want to add as well that the bug is minor in the sense that in case of complete ice cover or open water the resulting stress over ocean would OK.

Fred.

comment:2 in reply to: ↑ description ; follow-up: Changed 12 years ago by ctlod

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

Replying to nemo_user:

I found that the ice-ocean stress tio_u/tio_v computed in limdyn and later interpolated to u- and v- ocean points, includes both ice-ocean friction weighted by ice concentration and the wind stress (in our case, over ice) weighted by open water. However the routine ocesbc which gathers the open ocean stress and ice-ocean stress remultiply the tio_u/v by the ice concentration.


Thank you for this contribution, it has been corrected in a clean way in the next relase nemo_v3 version.

The fix is that tio_u/v only gathers the ice-ocean friction with no extra weight.

I guess that there was some confusion and mixed-up at some points since the 2 routines do more or less the same thing twice. No clue wether this bug was ported to LIM3.

Attached is my fix.

comment:3 in reply to: ↑ 2 Changed 12 years ago by ctlod

Replying to ctlod:

Replying to nemo_user:

I found that the ice-ocean stress tio_u/tio_v computed in limdyn and later interpolated to u- and v- ocean points, includes both ice-ocean friction weighted by ice concentration and the wind stress (in our case, over ice) weighted by open water. However the routine ocesbc which gathers the open ocean stress and ice-ocean stress remultiply the tio_u/v by the ice concentration.


Thank you for this contribution, it has been corrected in a clean way in the next relase nemo_v3 version.



Now, wind stress fields has been entirely reviewed and is computed into limsbc_2.F90 (if using LIM 2.0) or limsbc.F90 (if using LIM 3.0) modules.

The fix is that tio_u/v only gathers the ice-ocean friction with no extra weight.

I guess that there was some confusion and mixed-up at some points since the 2 routines do more or less the same thing twice. No clue wether this bug was ported to LIM3.

Attached is my fix.

comment:4 Changed 3 years ago by nemo

  • Keywords LIM* added
Note: See TracTickets for help on using tickets.