Changeset 6505
- Timestamp:
- 2016-05-01T18:11:06+02:00 (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMOGCM/NEMO/OPA_SRC/TRA/eosbn2.F90
r6504 r6505 1015 1015 REAL(wp), DIMENSION(jpi,jpj), INTENT(out ) :: ptf ! freezing temperature [Celcius] 1016 1016 ! 1017 INTEGER :: ji, jj ! dummy loop indices1018 REAL(wp) :: zt, zs ! local scalars1017 INTEGER :: ji, jj ! dummy loop indices 1018 REAL(wp) :: zt, zs, z1_S0 ! local scalars 1019 1019 !!---------------------------------------------------------------------- 1020 1020 ! … … 1023 1023 CASE ( np_teos10, np_seos ) !== CT,SA (TEOS-10 and S-EOS formulations) ==! 1024 1024 ! 1025 z1_S0 = 1._wp / 35.16504_wp 1025 1026 DO jj = 1, jpj 1026 1027 DO ji = 1, jpi 1027 zs= SQRT( ABS( psal(ji,jj) ) * r1_S0 ) ! square root salinity1028 zs= SQRT( ABS( psal(ji,jj) ) * z1_S0 ) ! square root salinity 1028 1029 ptf(ji,jj) = ((((1.46873e-03_wp*zs-9.64972e-03_wp)*zs+2.28348e-02_wp)*zs & 1029 1030 & - 3.12775e-02_wp)*zs+2.07679e-02_wp)*zs-5.87701e-02_wp … … 1074 1075 CASE ( np_teos10, np_seos ) !== CT,SA (TEOS-10 and S-EOS formulations) ==! 1075 1076 ! 1076 zs = SQRT( ABS( psal ) * r1_S0) ! square root salinity1077 zs = SQRT( ABS( psal ) / 35.16504_wp ) ! square root salinity 1077 1078 ptf = ((((1.46873e-03_wp*zs-9.64972e-03_wp)*zs+2.28348e-02_wp)*zs & 1078 1079 & - 3.12775e-02_wp)*zs+2.07679e-02_wp)*zs-5.87701e-02_wp
Note: See TracChangeset
for help on using the changeset viewer.