!copyright (C) 2001 MSC-RPN COMM %%%RPNPHY%%% ***S/R INCTPHY - INITIALISE LE COMMON 'CTESPHY'. *SUBROUTINE INCTPHY 2 * #include "impnone.cdk"
* *Author * N. Brunet (Jan91) * *Revision * 001 B. Bilodeau (Nov 1995) - Change KARMAN to 0.40 * *Object * to initialize the variables in common block CTESPHY by * reading the file "CONSTANTES". * *Arguments * None. * *PARAMETRES * NBRE - NOMBRE DE CONSTANTES DANS LE FICHIER INTEGER NBRE PARAMETER(NBRE=31) * *IMPLICITES #include "consphy.cdk"
*MODULES EXTERNAL CONSTNT * ** *---------------------------------------------------------------------- INTEGER FLAG, I REAL TEMP1(NBRE) * EQUIVALENCE (TEMP1(1),CPD) * DATA INIT/.FALSE./ * CHARACTER *10 NOM1(NBRE) * DATA NOM1/ 'CPD', 'CPV', 'RGASD', 'RGASV', 'TRPL', $ 'TCDK', 'RAUW', 'EPS1', 'EPS2', 'DELTA', $ 'CAPPA', 'TGL', 'CONSOL', 'GRAV', 'RAYT', $ 'STEFAN', 'PI', 'OMEGA', $ 'KNAMS', 'STLO', 'KARMAN', 'RIC', 'CHLC', 'CHLF', $ 'T1S', 'T2S', 'AW', 'BW', 'AI', 'BI', 'SLP'/ * *---------------------------------------------------------------------- * DO 10 I=1,NBRE CALL CONSTNT(TEMP1(I),FLAG,NOM1(I),0) IF(FLAG.EQ.0)THEN WRITE(6,600)NOM1(I) 600 FORMAT(/,5X,'LA CONSTANTE',2X,A10,1X,'N EXISTE PAS',/) CALL QQEXIT(1) ELSE ** WRITE(6,610)NOM1(I),TEMP1(I) 610 FORMAT(1X,'VALEUR DE',1X,A10,2X,'=',1X,E15.7) END IF 10 CONTINUE * * * DONNER A LA CONSTANTE "KARMAN" LA VALEUR 0.40 * --------------------------------------------- * CALL CONSTNT(0.40 ,FLAG,'KARMAN',3) CALL CONSTNT(KARMAN,FLAG,'KARMAN',0) * * WRITE(6,1000) WRITE(6,1010) 'THE VALUE OF THE CONSTANT KARMAN HAS BEEN *' WRITE(6,1010) ' ------ *' WRITE(6,1010) 'CHANGED TO ', KARMAN,' IN S/R INCTPHY' WRITE(6,1000) PRINT *,' ' * INIT=.TRUE. * * 1000 FORMAT ( ' ', + / ' ****************************************************', + / ' ') * 1010 FORMAT ( ' * ',A,F4.2,A,' *') * * RETURN END