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 5 years ago

Last modified 5 years ago

#2137 closed Bug

wrong formulae for Lupkes 2012 ? — at Version 1

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

Description (last modified by nicolasmartin)

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 * A**0

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

Commit History (0)

(No commits)

Change History (1)

comment:1 Changed 5 years ago by nicolasmartin

  • Description modified (diff)
Note: See TracTickets for help on using tickets.