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	heimbach	1.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	dgoldberg	1.2	SLOPE_LIMITER = (r+abs(r))/(1+abs(r))
19	heimbach	1.1	endif
20
21			RETURN
22			END