Changeset 5366 for trunk/NEMOGCM/NEMO
- Timestamp:
- 2015-06-05T13:21:13+02:00 (9 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMOGCM/NEMO/OPA_SRC/storng.F90
r5329 r5366 1 1 MODULE storng 2 !$AGRIF_DO_NOT_TREAT 2 3 !!====================================================================== 3 4 !! *** MODULE storng *** … … 89 90 y = m( m( m(y,13_8), -17_8 ), 43_8 ) 90 91 z = 6906969069_8 * z + 1234567_8 91 92 92 93 kiss = x + y + z 93 94 … … 266 267 267 268 uran = half * ( one + REAL(kiss(),8) / huge64 ) 268 269 269 270 END SUBROUTINE kiss_uniform 270 271 … … 334 335 CALL kiss_uniform(u2) 335 336 zz = u1 * u1 * u2 336 337 337 338 accepted = rr .GE. (zz*p1-p2) 338 339 IF (.NOT.accepted) accepted = rr .GE. LOG(zz) … … 345 346 ! (see Devroye, Non-Uniform Random Variate Generation, p. 415) 346 347 c = one/k ; d = (one-k) * EXP( (k/(one-k)) * LOG(k) ) 347 348 348 349 accepted=.FALSE. 349 350 DO WHILE (.NOT.accepted) … … 396 397 397 398 END SUBROUTINE kiss_sample 398 399 !$AGRIF_END_DO_NOT_TREAT 399 400 END MODULE storng
Note: See TracChangeset
for help on using the changeset viewer.