!copyright (C) 2001 MSC-RPN COMM %%%RPNPHY%%%SUBROUTINE INIRAS 1 * #include "impnone.cdk"
* *Author * B. Bilodeau (Apr 97) * *Revisions * *Object * initialize common block RASDATA for RAS convection scheme * *Arguments * *Notes * The code of INIRAS has been extracted from * s/r MRACRITN in order to make it reentrant. * ** * INTEGER L REAL TEM * #include "rasdata.cdk"
* DATA PH/150.0, 200.0, 250.0, 300.0, 350.0, 400.0, 450.0, 500.0 *, 550.0, 600.0, 650.0, 700.0, 750.0, 800.0, 850.0/ DATA A/ 1.6851, 1.1686, 0.7663, 0.5255, 0.4100, 0.3677 *, 0.3151, 0.2216, 0.1521, 0.1082, 0.0750, 0.0664 *, 0.0553, 0.0445, 0.0633/ DATA AD/15*0.0/, AC/15*0.0/, ACTOP/0.0/ * ACTOP = ACTP*FACM DO 5 L=1,15 A(L) = A(L)*FACM 5 CONTINUE DO 10 L=2,15 TEM = PH(L) - PH(L-1) AU(L) = A(L-1) / TEM AD(L) = A(L) / TEM AC(L) = PH(L)*AU(L) - PH(L-1)*AD(L) AD(L) = AD(L) - AU(L) 10 CONTINUE * RETURN END