Inicio Explorar Axuda
Iniciar sesión
AndDT
/
slackbuilds
Seguir 1
Destacar 0
Fork 0
Ficheiros Incidencias 0 Pull Requests 0 Wiki
Árbore: 9303d7f91a
Ramas Etiquetas
slackbuilds
/
academic
/
xfoil
/
xoper.f.patch

xoper.f.patch 2.3 KB
Histórico Raw

1234567891011121314151617181920212223242526272829303132333435363738394041424344454647484950515253545556575859606162 
  1. --- src/xoper.f	2013-12-08 13:21:22.000000000 -0500
    2. 

  2. +++ src.xoper.f.new	2016-09-15 20:50:03.000000000 -0400
    3. 

  3. @@ -3064,32 +3064,32 @@
    4. 

  4.          ione = int( (fnum-float(10*iten)) / 0.99999 )
    5. 

  5.          idec = int( (fnum-float(10*iten)-float(ione)) / 0.09999 )
    6. 

  6.  
    7. 

  7. -        fname = char(iten+izero) 
    8. 

  8. -     &       // char(ione+izero) 
    9. 

  9. -     &       // char(idec+izero) // '.bl'
    10. 

  10. -        lu = 44
    11. 

  11. -        open(lu,file=fname,status='unknown')
    12. 

  12. -        rewind(lu)
    13. 

  13. -        write(lu,'(a,a)') 
    14. 

  14. -     &'#       s         ue          H          P         K ',
    15. 

  15. -     &'        x    -m du/dx'
    16. 

  16. -c       1234567890 1234567890 1234567890 1234567890 1234567890 1234567890
    17. 

  17. -        do ibl = 2, iblte(is)
    18. 

  18. -          iblm = max( ibl-1 , 2 )
    19. 

  19. -          iblp = min( ibl+1 , iblte(is) )
    20. 

  20. -          i  = ipan(ibl ,is)
    21. 

  21. -          hk = dstr(ibl,is) / thet(ibl,is)
    22. 

  22. -          ddef = dstr(ibl,is)*uedg(ibl,is)
    23. 

  23. -          pdef = thet(ibl,is)*uedg(ibl,is)**2
    24. 

  24. -          edef = tstr(ibl,is)*uedg(ibl,is)**3 * 0.5
    25. 

  25. -          duds = (uedg(iblp,is)-uedg(iblm,is))
    26. 

  26. -     &         / (xssi(iblp,is)-xssi(iblm,is))
    27. 

  27. -          dpds = -ddef*duds
    28. 

  28. -          write(lu,9977) 
    29. 

  29. -     &       xssi(ibl,is), uedg(ibl,is), hk, pdef, edef, x(i), dpds
    30. 

  30. - 9977     format(1x, 3f11.4, 2f11.6, f11.3, e14.6 )
    31. 

  31. -        enddo
    32. 

  32. -        close(lu)
    33. 

  33. +c        fname = char(iten+izero) 
    34. 

  34. +c     &       // char(ione+izero) 
    35. 

  35. +c     &       // char(idec+izero) // '.bl'
    36. 

  36. +c        lu = 44
    37. 

  37. +c        open(lu,file=fname,status='unknown')
    38. 

  38. +c        rewind(lu)
    39. 

  39. +c        write(lu,'(a,a)') 
    40. 

  40. +c     &'#       s         ue          H          P         K ',
    41. 

  41. +c     &'        x    -m du/dx'
    42. 

  42. +cc       1234567890 1234567890 1234567890 1234567890 1234567890 1234567890
    43. 

  43. +c        do ibl = 2, iblte(is)
    44. 

  44. +c          iblm = max( ibl-1 , 2 )
    45. 

  45. +c          iblp = min( ibl+1 , iblte(is) )
    46. 

  46. +c          i  = ipan(ibl ,is)
    47. 

  47. +c          hk = dstr(ibl,is) / thet(ibl,is)
    48. 

  48. +c          ddef = dstr(ibl,is)*uedg(ibl,is)
    49. 

  49. +c          pdef = thet(ibl,is)*uedg(ibl,is)**2
    50. 

  50. +c          edef = tstr(ibl,is)*uedg(ibl,is)**3 * 0.5
    51. 

  51. +c          duds = (uedg(iblp,is)-uedg(iblm,is))
    52. 

  52. +c     &         / (xssi(iblp,is)-xssi(iblm,is))
    53. 

  53. +c          dpds = -ddef*duds
    54. 

  54. +c          write(lu,9977) 
    55. 

  55. +c     &       xssi(ibl,is), uedg(ibl,is), hk, pdef, edef, x(i), dpds
    56. 

  56. +c 9977     format(1x, 3f11.4, 2f11.6, f11.3, e14.6 )
    57. 

  57. +c        enddo
    58. 

  58. +c        close(lu)
    59. 

  59.  
    60. 

  60.  
    61. 

  61.  
    62. 

  62.