Changeset 15145 for NEMO/trunk/tests/ISOMIP+/MY_SRC/eosbn2.F90
- Timestamp:
- 2021-07-26T18:16:45+02:00 (3 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
NEMO/trunk/tests/ISOMIP+/MY_SRC/eosbn2.F90
r14995 r15145 1185 1185 z1_T0 = 1._wp/40._wp 1186 1186 ! 1187 DO_2D( 1, 1, 1, 1)1187 DO_2D( nn_hls, nn_hls, nn_hls, nn_hls ) 1188 1188 ! 1189 1189 zt = ctmp (ji,jj) * z1_T0 1190 zs = SQRT( ABS( psal(ji,jj) + zdeltaS ) * r1_S0 )1190 zs = SQRT( ABS( psal(ji,jj) + zdeltaS ) * z1_S0 ) 1191 1191 ztm = tmask(ji,jj,1) 1192 1192 ! … … 1249 1249 ! 1250 1250 z1_S0 = 1._wp / 35.16504_wp 1251 DO_2D( 1, 1, 1, 1)1251 DO_2D( nn_hls, nn_hls, nn_hls, nn_hls ) 1252 1252 zs= SQRT( ABS( psal(ji,jj) ) * z1_S0 ) ! square root salinity 1253 1253 ptf(ji,jj) = ((((1.46873e-03_wp*zs-9.64972e-03_wp)*zs+2.28348e-02_wp)*zs & … … 1356 1356 CASE( np_teos10, np_eos80 ) !== polynomial TEOS-10 / EOS-80 ==! 1357 1357 ! 1358 DO_3D( 1, 1, 1, 1, 1, jpkm1 )1358 DO_3D( nn_hls, nn_hls, nn_hls, nn_hls, 1, jpkm1 ) 1359 1359 ! 1360 1360 zh = gdept(ji,jj,jk,Kmm) * r1_Z0 ! depth … … 1415 1415 CASE( np_seos ) !== Vallis (2006) simplified EOS ==! 1416 1416 ! 1417 DO_3D( 1, 1, 1, 1, 1, jpkm1 )1417 DO_3D( nn_hls, nn_hls, nn_hls, nn_hls, 1, jpkm1 ) 1418 1418 zt = pts(ji,jj,jk,jp_tem) - 10._wp ! temperature anomaly (t-T0) 1419 1419 zs = pts (ji,jj,jk,jp_sal) - 35._wp ! abs. salinity anomaly (s-S0)
Note: See TracChangeset
for help on using the changeset viewer.