!-------------------------------------- LICENCE BEGIN ------------------------------------
!Environment Canada - Atmospheric Science and Technology License/Disclaimer,
! version 3; Last Modified: May 7, 2008.
!This is free but copyrighted software; you can use/redistribute/modify it under the terms
!of the Environment Canada - Atmospheric Science and Technology License/Disclaimer
!version 3 or (at your option) any later version that should be found at:
!http://collaboration.cmc.ec.gc.ca/science/rpn.comm/license.html
!
!This software is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY;
!without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.
!See the above mentioned License/Disclaimer for more details.
!You should have received a copy of the License/Disclaimer along with this software;
!if not, you can write to: EC-RPN COMM Group, 2121 TransCanada, suite 500, Dorval (Quebec),
!CANADA, H9P 1J3; or send e-mail to service.rpn@ec.gc.ca
!-------------------------------------- LICENCE END --------------------------------------
***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
* 002 B. Dugas (Aug 2005) - INIT defined in INCTPHY_DATA
*
*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,INCTPHY_DATA
*
**
*----------------------------------------------------------------------
INTEGER FLAG, I
REAL TEMP1(NBRE)
*
EQUIVALENCE (TEMP1(1),CPD)
*
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
BLOCK DATA INCTPHY_DATA
*
#include "consphy.cdk"
*
DATA INIT/.FALSE./
*
END BLOCK DATA INCTPHY_DATA