* 32, 1.826, 0.879, 0.169, 0.0,
* 36, 2.090, 0.818, 0.221, 0.0 /
C-------------------------------------------------------------
c ... Set the wound corrections of Green and Clothier (1988)
C-------------------------------------------------------------
DATA WC2/ 0, 0, 1, 0, 0.0,
* 16, -0.171, 1.299, 0.0194, -0.000093,
* 20, -0.159, 1.318, 0.0270, -0.000140,
* 24, -0.135, 1.326, 0.0367, -0.000194,
* 28, -0.143, 1.306, 0.0488, -0.000267,
* 32, -0.067, 1.355, 0.0571, -0.000203,
* 36, -0.013, 1.379, 0.0670, -0.000105/
PI=3.14159
TWOPI=2.*PI
C------------------------------------------------------------
c ... Set the graphics mode
C------------------------------------------------------------
CTITLE(1:40) = ' HEAT-PULSE analysis program '
CALL GRAPHICSMODE(CTITLE)
1 CALL CLEARSCREEN($GCLEARSCREEN)
CALL SETTEXTPOSITION(5,1,CURPOS)
NXPIX = WINC.NUMXPIXELS !-20
NYPIX = WINC.NUMYPIXELS !- 100
NCOL(1)=4
NCOL(2)=10
NCOL(3)=1
NCOL(4)=15
BEGYR = -9999
C------------------------------------------------------------
C ... Open input/output files
C------------------------------------------------------------
WRITE(*,'(a\)') 'Enter the Input filename (.dat): '
READ(*,'(A)') FILNAM
CALL NCSTRING(FILNAM,NC)
OPEN(UNIT=30, FILE=FILNAM(1:NC)//'.DAT', STATUS='UNKNOWN')
OPEN(UNIT=31, FILE=FILNAM(1:NC)//'.VEL', STATUS='UNKNOWN')
OPEN(UNIT=32, FILE=FILNAM(1:NC)//'.FLX', STATUS='UNKNOWN')
OPEN(UNIT=33, FILE=FILNAM(1:NC)//'.SUM', STATUS='UNKNOWN')
C------------------------------------------------------------
c ... Read input data
C------------------------------------------------------------
WRITE(*,*)
WRITE(*,*) 'Enter the number of heat pulse probes (1..4)'
READ(30,1099) ANS, NHPV
1099 FORMAT(A16,4I8)
WRITE(*,1099) ANS, NHPV
WRITE(*,*)
WRITE(*,*) 'Enter the number of sensors in each probe'
READ(30,1099) ANS, (NTC(J), J=1,NHPV)
WRITE(*,*) 'Enter the wound width from 1.6 to 3.6 (mm)'
READ(30,1098) ANS, (WW(I), I=1,NHPV)
1098 FORMAT(A16,4F8.2)
WRITE(*,1098) ANS, (WW(I), I=1,NHPV)
SWANCOR = .FALSE.
READ(30,1099) ANS, ICORR
IF(ICORR.GT.(0.0)) THEN
WRITE(*,*) 'Use Swansons corrections: YES'
SWANCOR=.TRUE.
ELSE
WRITE(*,*) 'Use Swansons corrections: NO'
SWANCOR=.FALSE.
ENDIF
WRITE(*,*)
DO 12 J=1,NHPV
IF(WW(J).EQ.(0.0)) IW(J) = 1
IF(WW(J).GE.(1.6)) IW(J) = INT((WW(J)-1.6)/0.4+2.5)
IF(IW(J).GT.7) IW(J) = 7
12 CONTINUE
C------------------------------------------------------------
(4 pages)
Manymanuals.com
Manymanuals.de
Manymanuals.fr
Manymanuals.it
Manymanuals.pl
Manymanuals.cz
Manymanuals.es
Manymanuals-pt.com
Comments to this Manuals