home *** CD-ROM | disk | FTP | other *** search
/ NetNews Usenet Archive 1992 #16 / NN_1992_16.iso / spool / comp / lang / fortran / 2779 < prev    next >
Encoding:
Text File  |  1992-07-21  |  2.1 KB  |  56 lines
  1. Newsgroups: comp.lang.fortran
  2. Path: sparky!uunet!cs.utexas.edu!torn!cunews!nrcnet0!cu23.crl.aecl.ca!wl.aecl.ca!harrisp
  3. From: harrisp@wl.aecl.ca
  4. Subject: Help! Partial Differential equation solvers
  5. Message-ID: <21JUL92.11592567@wl.aecl.ca>
  6. Sender: news@cu23.crl.aecl.ca (USENET News System)
  7. Nntp-Posting-Host: wc4.wl.aecl.ca
  8. Organization: AECL RESEARCH
  9. Date: Tue, 21 Jul 1992 17:59:25 GMT
  10. Lines: 44
  11.  
  12. To:   Greg Lampard
  13.  
  14. From: HARRISP@WL.AECL.CA   (Phil Harris)
  15. Date: 22-JUL-1992
  16. RE:   Solving partial differential equations numerically
  17.  
  18. Can anyone recommend a fortran callable  package for solving pde's and systems
  19. of pde's. I really need the source code, preferably written as standard 
  20. Fortran 77 as I want this to run on a PC (80486) using Lahey Fortran. Some 
  21. documentation would be nice also. I'm looking for something that is fast, 
  22. since the model I am developing has to run for several months of real time. 
  23. I wrote a simple program to solve the equations using explicit
  24. finite differences which works just fine, but I need to use very small time 
  25. steps for it to remain numerically stable, and takes almost a minute of 
  26. execution time on my '486 to simulate the results of a minute of real time.
  27. I am currently working on a code using implicit (Crank-Nicholson) finite 
  28. differencing, but am having some trouble getting it to execute correctly. 
  29. I'm hoping that once I get this working, I will be able to use larger time
  30. steps.
  31.  
  32. The equations I need to solve are fairly simple:
  33.  
  34. dQ(z,t)
  35. -------  = a1*( Const1(z) - Q(z,t) )*C(z,t) - a2*Q(z,t)
  36.   dt
  37.  
  38. dC(z,t)                dQ(z,t)    dC(z,t)
  39. -------- = -Const2 * ( ------- +  ------- )
  40.   dt                     dt         dz
  41.  
  42. a1, a2, Const2  are independant of z and t
  43. Const1 independant of t, but dependant on z.
  44.  
  45. Anybody have any suggestions. I've heard rumors of some NETLIB routines,
  46. but wouldn't know where to start in terms of which routines and where
  47. to get them from.
  48.  
  49.  
  50. Phil Harris                |  Imagination is more
  51. Whiteshell Laboratories            |  important than knowledge.
  52. Atomic Energy of Canada            |           Albert Einstein
  53. HARRISP@WL.AECL.CA            | 
  54.  
  55.  
  56.