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