Changeset 5304
- Timestamp:
- 2015-05-28T12:05:28+02:00 (9 years ago)
- Location:
- branches/2015/dev_r5177_CNRS4_stopar/NEMOGCM/NEMO/OPA_SRC
- Files:
-
- 2 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 -
branches/2015/dev_r5177_CNRS4_stopar/NEMOGCM/NEMO/OPA_SRC/stopts.F90
r5296 r5304 20 20 21 21 PUBLIC sto_pts ! called by step.F90 22 PUBLIC sto_pts_init ! called by nemogcm.F90 22 PUBLIC sto_pts_init ! called by nemogcm.F90 23 23 24 24 ! Public array with random tracer fluctuations … … 37 37 !!---------------------------------------------------------------------- 38 38 !! *** ROUTINE sto_pts *** 39 !! 39 !! 40 40 !! ** Purpose : Compute current stochastic tracer fluctuations 41 41 !! … … 84 84 ! & zdtskp * MAX(sto2d(ji,jj,jsto_eosk),0._wp) + & 85 85 ! & zdtskm * MIN(sto2d(ji,jj,jsto_eosk),0._wp) 86 zdts = zdts * tmask(ji,jj,jk) *SIN( 2._wp *gphit(ji,jj) * rad )86 zdts = zdts * tmask(ji,jj,jk) *SIN( gphit(ji,jj) * rad ) 87 87 pts_ran(ji,jj,jk,jts,jdof) = zdts * 0.5_wp 88 88 ! … … 133 133 !!---------------------------------------------------------------------- 134 134 !! *** ROUTINE sto_pts_init *** 135 !! 135 !! 136 136 !! ** Purpose : Initialisation for stochastic tracer fluctuations 137 137 !!
Note: See TracChangeset
for help on using the changeset viewer.