#504 closed Bug (fixed)
bug in the coriolis term in LIM2
Reported by: | bouillon | Owned by: | nemo |
---|---|---|---|
Priority: | low | Milestone: | |
Component: | LIM2 | Version: | trunk |
Severity: | Keywords: | LIM* | |
Cc: |
Description
There is a bug in line 175 of limrhg.f90 in the trunk version of LIM2.
!!! zcorl(ji,jj) = zmass(ji,jj) * fcor(ji,jj)
The coriolis paramater fcor is located at center of grid cells, not on the corners. The expression below seems the correct one to me (mamm).
zcorl(ji,jj) = zmass(ji,jj) &
& *( fcor(ji,jj) * wght(ji,jj,2,2) + fcor(ji-1,jj)*wght(ji,jj,1,2) &
& + fcor(ji,jj-1) * wght(ji,jj,2,1) + fcor(ji-1,jj-1)*wght(ji,jj,1,1)) * zusw
It is ok in LIM3.
Sylvain
Commit History (1)
Changeset | Author | Time | ChangeLog |
---|---|---|---|
1569 | ctlod | 2009-07-31T18:37:36+02:00 | correction of a bug in the Coriolis term in the rheology of LIM2, see ticket: #504 |
Change History (5)
comment:1 Changed 15 years ago by gm
- Component changed from OPA to LIM2
comment:2 Changed 15 years ago by gm
comment:3 Changed 15 years ago by ctlod
- Resolution set to fixed
- Status changed from new to closed
comment:4 Changed 8 years ago by nicolasmartin
- Keywords coriolis removed
comment:5 Changed 7 years ago by nemo
- Keywords LIM* added
Note: See
TracTickets for help on using
tickets.
I agree, this is a bug.
A less time consuming solution should be to change the definition of fcor to set it to the coriolis value at I-point.
fcor(ji,jj) = ff(ji-1,jj-1)
Be carefull, when applied the lbc_blk to fcor ....
gurvan