- Timestamp:
- 2012-05-17T18:33:12+02:00 (12 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2012/dev_r3385_NOCS04_HAMF/NEMOGCM/NEMO/OPA_SRC/TRA/eosbn2.F90
r3294 r3396 121 121 REAL(wp) :: zd , zc , zaw, za ! - - 122 122 REAL(wp) :: zb1, za1, zkw, zk0 ! - - 123 REAL(wp) :: zrau0r ! - -124 123 REAL(wp), POINTER, DIMENSION(:,:,:) :: zws 125 124 !!---------------------------------------------------------------------- … … 133 132 ! 134 133 CASE( 0 ) !== Jackett and McDougall (1994) formulation ==! 135 zrau0r = 1.e0 / rau0136 134 !CDIR NOVERRCHK 137 135 zws(:,:,:) = SQRT( ABS( pts(:,:,:,jp_sal) ) ) … … 174 172 ! masked in situ density anomaly 175 173 prd(ji,jj,jk) = ( zrhop / ( 1.0_wp - zh / ( zk0 - zh * ( za - zh * zb ) ) ) & 176 & - rau0 ) * zrau0r* tmask(ji,jj,jk)174 & - rau0 ) * r1_rau0 * tmask(ji,jj,jk) 177 175 END DO 178 176 END DO … … 254 252 INTEGER :: ji, jj, jk ! dummy loop indices 255 253 REAL(wp) :: zt, zs, zh, zsr, zr1, zr2, zr3, zr4, zrhop, ze, zbw ! local scalars 256 REAL(wp) :: zb, zd, zc, zaw, za, zb1, za1, zkw, zk0 , zrau0r! - -254 REAL(wp) :: zb, zd, zc, zaw, za, zb1, za1, zkw, zk0 ! - - 257 255 REAL(wp), POINTER, DIMENSION(:,:,:) :: zws 258 256 !!---------------------------------------------------------------------- … … 265 263 ! 266 264 CASE( 0 ) !== Jackett and McDougall (1994) formulation ==! 267 zrau0r = 1.e0 / rau0268 265 !CDIR NOVERRCHK 269 266 zws(:,:,:) = SQRT( ABS( pts(:,:,:,jp_sal) ) ) … … 309 306 ! masked in situ density anomaly 310 307 prd(ji,jj,jk) = ( zrhop / ( 1.0_wp - zh / ( zk0 - zh * ( za - zh * zb ) ) ) & 311 & - rau0 ) * zrau0r* tmask(ji,jj,jk)308 & - rau0 ) * r1_rau0 * tmask(ji,jj,jk) 312 309 END DO 313 310 END DO
Note: See TracChangeset
for help on using the changeset viewer.