dgoldberg/streamice/slope_limiter.F

#include "CPP_OPTIONS.h"

      _RL FUNCTION SLOPE_LIMITER (num, denom)
 
       _RL num, denom
!       _RL SLOPE_LIMITER
       _RL r

! used for flux limiting in advective subroutines
! Van Leer limiter (source: Wikipedia!)

       if (denom .eq. 0.0) then
        SLOPE_LIMITER = 0.0
       elseif (num*denom .le. 0.0) then
        SLOPE_LIMITER = 0.0
       else
        r = num/denom
        SLOPE_LIMITER = (r+abs(r))/(1+abs(r))
       endif

       RETURN
       END 
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 = (r+abs(r))/(1+abs(r))
19	endif
20
21	RETURN
22	END