- Timestamp:
- 2015-05-28T12:05:28+02:00 (9 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2015/dev_r5177_CNRS4_stopar/NEMOGCM/NEMO/OPA_SRC/TRA/eosbn2.F90
r5296 r5304 315 315 REAL(wp), DIMENSION(jpi,jpj,jpk ), INTENT(in ) :: pdep ! depth [m] 316 316 ! 317 INTEGER :: ji, jj, jk, jsmp, jdof ! dummy loop indices 317 INTEGER :: ji, jj, jk, jsmp ! dummy loop indices 318 INTEGER :: jdof 318 319 REAL(wp) :: zt , zh , zstemp, zs , ztm ! local scalars 319 320 REAL(wp) :: zn , zn0, zn1, zn2, zn3 ! - - … … 345 346 ! (2) for t-dt, s-ds (with the opposite fluctuation) 346 347 DO jsmp = 1, nn_sto_eos*2 348 jdof = (jsmp + 1) / 2 347 349 zh = pdep(ji,jj,jk) * r1_Z0 ! depth 348 zt = pts (ji,jj,jk,jp_tem) * r1_T0 + pts_ran(ji,jj,jk,jp_tem,jdof) * zsign(jsmp)! temperature349 zstemp = pts 350 zt = (pts (ji,jj,jk,jp_tem) + pts_ran(ji,jj,jk,jp_tem,jdof) * zsign(jsmp)) * r1_T0 ! temperature 351 zstemp = pts (ji,jj,jk,jp_sal) + pts_ran(ji,jj,jk,jp_sal,jdof) * zsign(jsmp) 350 352 zs = SQRT( ABS( zstemp + rdeltaS ) * r1_S0 ) ! square root salinity 351 353 ztm = tmask(ji,jj,jk) ! tmask … … 372 374 & + (((((EOS600*zs+EOS500)*zs+EOS400)*zs+EOS300)*zs+EOS200)*zs+EOS100)*zs+EOS000 373 375 ! 374 zn_sto(jsmp) = ( ( zn3 * zh + zn2 ) * zh + zn1 ) * zh + zn0 376 zn_sto(jsmp) = ( ( zn3 * zh + zn2 ) * zh + zn1 ) * zh + zn0_sto(jsmp) 375 377 END DO 376 !377 378 ! 378 379 ! compute stochastic density as the mean of the (2*nn_sto_eos) densities … … 380 381 DO jsmp = 1, nn_sto_eos*2 381 382 prhop(ji,jj,jk) = prhop(ji,jj,jk) + zn0_sto(jsmp) ! potential density referenced at the surface 382 ! 383 ! 383 384 prd(ji,jj,jk) = prd(ji,jj,jk) + ( zn_sto(jsmp) * r1_rau0 - 1._wp ) ! density anomaly (masked) 384 385 END DO
Note: See TracChangeset
for help on using the changeset viewer.