#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: |
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)
Changeset | Author | Time | ChangeLog |
---|---|---|---|
3791 | gm | 2013-02-10T14:16:41+01:00 | dev_MERGE_2012: #1041 : corrected diagnostic for the ice shear rate |
3789 | gm | 2013-02-10T14:11:44+01:00 | dev_v3_4_STABLE_2012: #1041 : corrected diagnostic for the ice shear rate |
Attachments (1)
Change History (7)
Changed 11 years ago by bouillon
comment:1 Changed 11 years ago by clevy
- Owner changed from NEMO team to gm
comment:2 Changed 11 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 8 years ago by nicolasmartin
- Keywords invariant removed
comment:4 Changed 8 years ago by nicolasmartin
- Keywords rate removed
comment:5 Changed 8 years ago by nicolasmartin
- Keywords shear removed
comment:6 Changed 6 years ago by nemo
- Keywords LIM* added
limrhg.F90 corrected