summaryrefslogtreecommitdiff
path: root/academic/xfoil/xoper.f.patch
blob: cff6f6fd30297bcb7dfd57946af99b94603f3060 (plain)
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
--- src/xoper.f	2013-12-08 13:21:22.000000000 -0500
+++ src.xoper.f.new	2016-09-15 20:50:03.000000000 -0400
@@ -3064,32 +3064,32 @@
         ione = int( (fnum-float(10*iten)) / 0.99999 )
         idec = int( (fnum-float(10*iten)-float(ione)) / 0.09999 )
 
-        fname = char(iten+izero) 
-     &       // char(ione+izero) 
-     &       // char(idec+izero) // '.bl'
-        lu = 44
-        open(lu,file=fname,status='unknown')
-        rewind(lu)
-        write(lu,'(a,a)') 
-     &'#       s         ue          H          P         K ',
-     &'        x    -m du/dx'
-c       1234567890 1234567890 1234567890 1234567890 1234567890 1234567890
-        do ibl = 2, iblte(is)
-          iblm = max( ibl-1 , 2 )
-          iblp = min( ibl+1 , iblte(is) )
-          i  = ipan(ibl ,is)
-          hk = dstr(ibl,is) / thet(ibl,is)
-          ddef = dstr(ibl,is)*uedg(ibl,is)
-          pdef = thet(ibl,is)*uedg(ibl,is)**2
-          edef = tstr(ibl,is)*uedg(ibl,is)**3 * 0.5
-          duds = (uedg(iblp,is)-uedg(iblm,is))
-     &         / (xssi(iblp,is)-xssi(iblm,is))
-          dpds = -ddef*duds
-          write(lu,9977) 
-     &       xssi(ibl,is), uedg(ibl,is), hk, pdef, edef, x(i), dpds
- 9977     format(1x, 3f11.4, 2f11.6, f11.3, e14.6 )
-        enddo
-        close(lu)
+c        fname = char(iten+izero) 
+c     &       // char(ione+izero) 
+c     &       // char(idec+izero) // '.bl'
+c        lu = 44
+c        open(lu,file=fname,status='unknown')
+c        rewind(lu)
+c        write(lu,'(a,a)') 
+c     &'#       s         ue          H          P         K ',
+c     &'        x    -m du/dx'
+cc       1234567890 1234567890 1234567890 1234567890 1234567890 1234567890
+c        do ibl = 2, iblte(is)
+c          iblm = max( ibl-1 , 2 )
+c          iblp = min( ibl+1 , iblte(is) )
+c          i  = ipan(ibl ,is)
+c          hk = dstr(ibl,is) / thet(ibl,is)
+c          ddef = dstr(ibl,is)*uedg(ibl,is)
+c          pdef = thet(ibl,is)*uedg(ibl,is)**2
+c          edef = tstr(ibl,is)*uedg(ibl,is)**3 * 0.5
+c          duds = (uedg(iblp,is)-uedg(iblm,is))
+c     &         / (xssi(iblp,is)-xssi(iblm,is))
+c          dpds = -ddef*duds
+c          write(lu,9977) 
+c     &       xssi(ibl,is), uedg(ibl,is), hk, pdef, edef, x(i), dpds
+c 9977     format(1x, 3f11.4, 2f11.6, f11.3, e14.6 )
+c        enddo
+c        close(lu)