Opened 8 years ago

Closed 8 years ago

Last modified 3 years ago

#1041 closed Bug (fixed)

Corrected diagnostic for the ice shear rate

Reported by: bouillon Owned by: gm
Priority: low Milestone:
Component: LIM3 Version: trunk
Severity: Keywords: LIM*
Cc: Branch review:
MP ready?: Task progress:

Description

The ice shear rate invariant computed in limrhg.F90 is not correct.

It is only a diagnostic and it is not used in the rest of the model, so there will be no change in the results.

Line 695
shear_i(ji,jj) = zds (ji,jj)
should be replaced by

zdst = ( e2u( ji , jj ) * v_ice1(ji,jj) &
& - e2u( ji-1, jj ) * v_ice1(ji-1,jj) &
& + e1v( ji , jj ) * u_ice2(ji,jj) &
& - e1v( ji , jj-1 ) * u_ice2(ji,jj-1) &
& ) &
& / area(ji,jj)
shear_i(ji,jj) = SQRT( zdt(ji,jj)*zdt(ji,jj) + zdst*zdst )

and line 702
CALL lbc_lnk( shear_i(:,:), 'F', 1. )
by
CALL lbc_lnk( shear_i(:,:), 'T', 1. )

Commit History (2)

ChangesetAuthorTimeChangeLog
3791gm2013-02-10T14:16:41+01:00

dev_MERGE_2012: #1041 : corrected diagnostic for the ice shear rate

3789gm2013-02-10T14:11:44+01:00

dev_v3_4_STABLE_2012: #1041 : corrected diagnostic for the ice shear rate

Attachments (1)

new_limrhg.F90 (38.6 KB) - added by bouillon 8 years ago.
limrhg.F90 corrected

Download all attachments as: .zip

Change History (7)

Changed 8 years ago by bouillon

limrhg.F90 corrected

comment:1 Changed 8 years ago by clevy

  • Owner changed from NEMO team to gm

comment:2 Changed 8 years ago by gm

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

Thanks Sylvain for this correction

The changes have been made. There is a minor difference compare to what you propose (introduction of a 2D local array zdst) in order to avoid duplicate calculation.

see Changeset 3789 for branches/2012/dev_v3_4_STABLE_2012 (i.e. v3.4.1 stable)

and 

see Changeset 3791 for branches/2012/dev_MERGE_2012 (i.e. v3.5 alpha)

gurvan

comment:3 Changed 5 years ago by nicolasmartin

  • Keywords invariant removed

comment:4 Changed 5 years ago by nicolasmartin

  • Keywords rate removed

comment:5 Changed 5 years ago by nicolasmartin

  • Keywords shear removed

comment:6 Changed 3 years ago by nemo

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