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.
#2137 (wrong formulae for Lupkes 2012 ?) – NEMO

Opened 6 years ago

Last modified 5 years ago

#2137 closed Bug

wrong formulae for Lupkes 2012 ? — at Initial Version

Reported by: cbricaud Owned by: systeam
Priority: low Milestone:
Component: SBC Version: trunk
Severity: minor Keywords:
Cc:

Description

BE CAREFUL !!! Due to dynamic behaviour of this ticket creation page, it is highly recommend to set first all other fields before writing the ticket description below.
If you have lost your draft after an unwanted reload, you can click on the link 'Restore Form' in the contextual menu upper right to recover it.

Context

REAL(wp), DIMENSION(:,:), INTENT(inout) :: Cd

REAL(wp), PARAMETER
zCe = 2.23e-03_wp
REAL(wp), PARAMETER
znu = 1._wp
REAL(wp), PARAMETER
zmu = 1._wp
REAL(wp), PARAMETER
zbeta = 1._wp
REAL(wp)
zcoef !!---------------------------------------------------------------------- zcoef = znu + 1._wp / ( 10._wp * zbeta )


! ice-atm drag
Cd(:,:) = Cd_ice + & ! pure ice drag

& zCe * ( 1._wp - at_i_b(:,:) )zcoef * at_i_b(:,:)(zmu-1._wp) ! change due to sea-ice morphology

so
Cd(:,:) = Cd_ice + 2.23e-03 * ( 1._wp - A )1.1 * A0

but in the paper ( https://agupubs.onlinelibrary.wiley.com/doi/full/10.1029/2012JD017630 ) whe have:

Cd(:,:) = 2.23e-03 * ( 1._wp - A )1.1 * A

so whe shoud have in the code:

Cd(:,:) = Cd_ice + & ! pure ice drag

& zCe * ( 1._wp - at_i_b(:,:) )zcoef * at_i_b(:,:)zmu ! change due to sea-ice morphology


Analysis

...

Fix

...

Commit History (0)

(No commits)

Change History (0)

Note: See TracTickets for help on using tickets.