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.
Changeset 6505 – NEMO

Changeset 6505


Ignore:
Timestamp:
2016-05-01T18:11:06+02:00 (8 years ago)
Author:
gm
Message:

#1722 - trunk: TEOS-10 freezing point bug correction

File:
1 edited

Legend:

Unmodified
Added
Removed
  • trunk/NEMOGCM/NEMO/OPA_SRC/TRA/eosbn2.F90

    r6504 r6505  
    10151015      REAL(wp), DIMENSION(jpi,jpj), INTENT(out  )           ::   ptf    ! freezing temperature [Celcius] 
    10161016      ! 
    1017       INTEGER  ::   ji, jj   ! dummy loop indices 
    1018       REAL(wp) ::   zt, zs   ! local scalars 
     1017      INTEGER  ::   ji, jj          ! dummy loop indices 
     1018      REAL(wp) ::   zt, zs, z1_S0   ! local scalars 
    10191019      !!---------------------------------------------------------------------- 
    10201020      ! 
     
    10231023      CASE ( np_teos10, np_seos )      !==  CT,SA (TEOS-10 and S-EOS formulations) ==! 
    10241024         ! 
     1025         z1_S0 = 1._wp / 35.16504_wp 
    10251026         DO jj = 1, jpj 
    10261027            DO ji = 1, jpi 
    1027                zs= SQRT( ABS( psal(ji,jj) ) * r1_S0 )           ! square root salinity 
     1028               zs= SQRT( ABS( psal(ji,jj) ) * z1_S0 )           ! square root salinity 
    10281029               ptf(ji,jj) = ((((1.46873e-03_wp*zs-9.64972e-03_wp)*zs+2.28348e-02_wp)*zs & 
    10291030                  &          - 3.12775e-02_wp)*zs+2.07679e-02_wp)*zs-5.87701e-02_wp 
     
    10741075      CASE ( np_teos10, np_seos )      !==  CT,SA (TEOS-10 and S-EOS formulations) ==! 
    10751076         ! 
    1076          zs  = SQRT( ABS( psal ) * r1_S0 )           ! square root salinity 
     1077         zs  = SQRT( ABS( psal ) / 35.16504_wp )           ! square root salinity 
    10771078         ptf = ((((1.46873e-03_wp*zs-9.64972e-03_wp)*zs+2.28348e-02_wp)*zs & 
    10781079                  &          - 3.12775e-02_wp)*zs+2.07679e-02_wp)*zs-5.87701e-02_wp 
Note: See TracChangeset for help on using the changeset viewer.