00001 c artrsmrascala_.F
00002 c*********************************************************************
00003 c* (c) Copyright 2000-2002 Brown Deer Technology, LLC.
00004 c* All rights reserved.
00005 c*********************************************************************
00006
00007 subroutine artrsmrascala(
00008 & Asize,Msize
00009 & ,amax,mmax
00010 & ,mrastep,mracoef
00011 & ,mracoefrs
00012 & ,m,scala
00013 & )
00014
00015 implicit none
00016
00017 c ****************************************
00018 c ***** COMMON BLOCKS AND PARAMETERS *****
00019 c ****************************************
00020
00021 #include "stormdef.h"
00022 real INVSQRT2
00023 parameter( INVSQRT2 = STORMDEF_MATHCONST_INVSQRT2 )
00024
00025 c *******************************
00026 c ***** VARIABLES EXCHANGED *****
00027 c *******************************
00028
00029 integer Asize
00030 integer Msize
00031 integer amax
00032 integer mmax
00033 integer mrastep
00034 real mracoef(0:1,0:Msize-1,0:Asize-1)
00035 real mracoefrs(0:1,0:Msize-1,0:Asize-1)
00036 integer m
00037 real scala(0:Asize-1)
00038
00039 c ****************************
00040 c ***** LOCAL VARIABLES ******
00041 c ****************************
00042
00043 integer a
00044
00045 c ******************************
00046 c ***** EXTERNAL FUNCTIONS *****
00047 c ******************************
00048
00049 c if ((m.ge.0).and.(m.le.mmax)) then
00050 do a=0,amax
00051 scala(a)=mracoef(0,m,a)*INVSQRT2**(m)
00052 enddo
00053 c else
00054 c do a=0,amax
00055 c scala(a)=1.0/0.0
00056 c enddo
00057 c endif
00058
00059 return
00060 end
00061 c********************************************************
