mirror of
https://github.com/Ponce/slackbuilds
synced 2024-11-13 21:56:09 +01:00
62 lines
2.3 KiB
Diff
62 lines
2.3 KiB
Diff
|
--- 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)
|
||
|
|
|
||
|
|
|
||
|
|
|