Parent Directory
| 
Revision Log
| 
Revision Graph
Initial check-in of Dan Goldberg's streamice package
| 1 | #include "CPP_OPTIONS.h" |
| 2 | |
| 3 | _RL FUNCTION SLOPE_LIMITER (num, denom) |
| 4 | |
| 5 | _RL num, denom |
| 6 | ! _RL SLOPE_LIMITER |
| 7 | _RL r |
| 8 | |
| 9 | ! used for flux limiting in advective subroutines |
| 10 | ! Van Leer limiter (source: Wikipedia!) |
| 11 | |
| 12 | if (denom .eq. 0.0) then |
| 13 | SLOPE_LIMITER = 0.0 |
| 14 | elseif (num*denom .le. 0.0) then |
| 15 | SLOPE_LIMITER = 0.0 |
| 16 | else |
| 17 | r = num/denom |
| 18 | SLOPE_LIMITER = 0.0*(r+abs(r))/(1+abs(r)) |
| 19 | endif |
| 20 | |
| 21 | RETURN |
| 22 | END |
| ViewVC Help | |
| Powered by ViewVC 1.1.22 |