Changeset 5065


Ignore:
Timestamp:
2015-02-06T12:43:59+01:00 (6 years ago)
Author:
davestorkey
Message:

Fix bug in cd_neutral_10m in sbcblk_core.F90. See ticket #1465.

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/NEMOGCM/NEMO/OPA_SRC/SBC/sbcblk_core.F90

    r4990 r5065  
    848848      rgt33 = 0.5_wp + SIGN( 0.5_wp, (zw10 - 33._wp) )   ! If zw10 < 33. => 0, else => 1   
    849849      cd_neutral_10m = 1.e-3 * ( & 
    850          &       (rgt33 + 1._wp)*( 2.7_wp/zw10 + 0.142_wp + zw10/13.09_wp - 3.14807E-10*zw10**6) & ! zw10< 33. 
     850         &       (1._wp - rgt33)*( 2.7_wp/zw10 + 0.142_wp + zw10/13.09_wp - 3.14807E-10*zw10**6) & ! zw10< 33. 
    851851         &      + rgt33         *      2.34   )                                                    ! zw10 >= 33. 
    852852      ! 
Note: See TracChangeset for help on using the changeset viewer.