/[necpp]/necpp/FORTRAN/somnec_readme.txt
ViewVC logotype

Contents of /necpp/FORTRAN/somnec_readme.txt

Parent Directory Parent Directory | Revision Log Revision Log


Revision 1.1 - (hide annotations) (download)
Sat Sep 17 02:35:02 2005 UTC (7 years, 9 months ago) by tcam-guest
Branch: MAIN
CVS Tags: HEAD
File MIME type: text/plain
Added Sommerfeld Norton Code
1 tcam-guest 1.1 The accompanying file SOMNEC.FOR contains the FORTRAN source code for
2     the SOMNEC program. The SOMNEC program generates grids used by the NEC2
3     program for the Sommerfeld-Norton approximation to realistic grounds.
4     I hand-coded the source code from the listing given in NOSC TD-116,
5     Volume I, Part II, pages 390 -- 433. I changed the code given in the listing
6     to enable interactive user input of relative dielectric constant, ground
7     conductivity, and frequency in MHz. The source code as given relied on the use
8     of a tape drive to accomplish communication of the results of the SOMNEC
9     program to the NEC2 program (old-timers, takes you back, doesn't it?).
10     I elected to accomplish inter-program communication on the DEC VAX and Unix
11     workstations available to me by using a binary unformatted file. Hence the
12     user is also prompted for a file name, and the grid arrays are written to the
13     file. I also had to make changes in the NEC2 code, to tell NEC2 the name of
14     the SOMNEC file, to open the file, and to read in the file contents. Users
15     wishing to employ the SOMNEC code with implementations of NEC2 or NEC81
16     expecting SOMNEC file input will have to determine the file type and format
17     expected and recode accordingly. I eventually incorporated the SOMNEC
18     program in NEC2 as a set of subroutines, and that is the version I am now
19     using.
20     The file SOMNEC_OUTPUT.TXT contains the results of a run of SOMNEC
21     using the inputs in example 10 in Volume II of NOSC TD 116. Note the
22     agreement with the results given on pages 154 -- 156 of Volume II. The
23     time for execution of this case was 257 seconds on a DEC VAX 8600.
24     On a Silicon Graphics 4D/440 workstation, I obtained around 100 seconds
25     execution time for the same case, and I would expect tens of seconds when I
26     port the code to a DEC AXP, as I will do soon.
27     I hope this code is useful to you -- 73!
28    
29     George B. Christianson
30     Princeton University Plasma Physics Laboratory
31    
32     gchristianson@pppl.gov
33    
34     (Amateur Radio NJ2P)
35    
36     October 5, 1994

  ViewVC Help
Powered by ViewVC 1.1.5