From 7dfcc480ba1e19bd3232349fc733caef94034292 Mon Sep 17 00:00:00 2001 From: stainer_t Date: Mon, 8 Sep 2025 13:48:49 +0200 Subject: Initial commit from Polytechnique Montreal --- Dragon/src/SYBRII.f | 39 +++++++++++++++++++++++++++++++++++++++ 1 file changed, 39 insertions(+) create mode 100644 Dragon/src/SYBRII.f (limited to 'Dragon/src/SYBRII.f') diff --git a/Dragon/src/SYBRII.f b/Dragon/src/SYBRII.f new file mode 100644 index 0000000..2db117e --- /dev/null +++ b/Dragon/src/SYBRII.f @@ -0,0 +1,39 @@ +*DECK SYBRII + SUBROUTINE SYBRII(PP,SG,TAU0,TAUI) +* +*----------------------------------------------------------------------- +* +*Purpose: +* Evaluation of the $R_{ii}$ function in 1D slab geometry. +* +*Copyright: +* Copyright (C) 2005 Ecole Polytechnique de Montreal +* This library is free software; you can redistribute it and/or +* modify it under the terms of the GNU Lesser General Public +* License as published by the Free Software Foundation; either +* version 2.1 of the License, or (at your option) any later version +* +*Author(s): A. Hebert +* +*Parameters: input +* SG specular parity sign. +* TAU0 side-to-side optical path. +* TAUI optical path in volume i. +* +*Parameters: output +* PP value of the expression. +* +*----------------------------------------------------------------------- +* +*---- +* SUBROUTINE ARGUMENTS +*---- + REAL PP,SG,TAU0,TAUI +* + IF(TAUI.GT.1.0E-10) THEN + PP=SG*TABEN(2,TAU0)*TAUI-(TABEN(3,TAU0)-TABEN(3,TAU0+SG*TAUI)) + ELSE + PP=0.5*(TAUI**2)*TABEN(1,TAU0) + ENDIF + RETURN + END -- cgit v1.2.3