last update: 9-28-92
1. THE APPLICATION

This application simulates scattering of electromagnetic waves  from
loaded multiple apertures. 

There are two slots on an infinitive extension metal plane. The TE plane 
wave  is incident to the structures from left. In the metal plane (to the 
right), there are two waveguides terminated by one microwave network. 
We intend to calculate  the scattering fields in the space (left the metal 
plane). The coordinates are set as follows: 

The point of origin is in the middle of the two slots. The x-axis of 
abscissas is perpendicular to the slots. The z-axis is in the slot axis 
direction.

The results give the pattern of the scattering magnetic field. The field point 
are allocated in a semicircle in x-y coordinate plane. The diameter of the 
circle is R and the points are equally divided. 

For details see SCCS-327 (in the directory benchmark/doc)
============================================================================

2. ORGANIZATION

Four versions of the code are provided here (each stored in its own 
subdirectory)

(1) sequential code named 'emf77.f' is written in standard Fortran-77.

(2) CM code named 'emfCM-2.fcm' is written in CM Fortran and use the CM 
    Scientific Subroutine Library (CMSSL). When you compile that, use the 
    command "cmf emfCM-2.fcm -lcmssl -o emfCM". There is also a version
    without calls to CMSSL ("pure Fortran-90", actually mpfortran code
    adapted for the Connection Machine). It is named efmF90.fcm and can
    be compiled on CM-2/CM-5 issuing the following command:
    "cmf emfF90.fcm -o emfF90.

(3) mpfortran  code is named 'emfF90.f'. To run it, use 
     "mpfortran -pesize=8 -o emfF90 emfF90.f" on MasPar/DECmpp 120000. 
   
(4) F77 with explicit message passing codes named 'emphost.f', 'empnode.f', 
    'empnode1.f', 'emphosts.f', 'empnodes.f', 'input', and 'makefile'.
    in which, 'emphost.f', 'empnode.f', 'empnode1.f' are double precision
    version and 'emphosts.f', 'empnodes.f' are single precision version.
    First, you use 'make rxnode' on resi.cs.rice.edu. Then, you can use
    'make rxhost' on io.cs.rice.edu. After you get 'emphost*' and 'empnode*'
    you can use 'getcube -td2' and then 'bphost < input' to run the program.
    In 'empnode.f' and 'empnode1.f', we use different distribution methods.
    Do not forget using 'relcube' to release the cube after you finish.


3. EXAMPLE INPUT

When you run the code. you will be asked to input the original data. The 
suggestion is:
I=5 J=5 M=20 N=20 N2=901  
a=0.25 b=0.25 d=0.5 l1=0.3 l2=0.3  
(1) space permittivity= 1.     (2) space permeability= 1. 
(3) waveguide permittivity= 1. (4) waveguide permeability= 1.  
wave incident angle = 0.5236                  wave length= 1. 
the distance from the measurement points= 10. the beginning angle = 0.  
Y11,Y12,Y21,Y22 are complex. The values can be:
Y11= (0.,0.) Y12= (1.,0.) Y21= (1.,0.) Y22= (0.,0.)  

4. EXAMPLE OUTPUT

The results generated by this application, that is the scattered field are 
stored in em9outs1.dat and em9outs2.dat created during the run time. For
the example input described above, the program generated the following results:

em9outs1.dat:
the current values are
  (-.8750663, .2970908) ( 8.7245658E-002, 6.7622304E-002)
 ( 9.3538225E-002,-6.3795924E-002) ( 3.9160445E-002, 7.7010617E-002)
 ( 9.4122812E-002,-4.0824637E-002) (-.8750663, .2970908)
 ( 8.7245658E-002, 6.7622304E-002) ( 9.3538225E-002,-6.3795924E-002)
 ( 3.9160445E-002, 7.7010617E-002) ( 9.4122812E-002,-4.0824637E-002)
and the voltages values are 
  (-109.4576, 51.48106) ( 21.57812,-13.64320) (-33.15636,-13.35294)
 ( 46.15070, 8.904148) (-23.28301,-33.85236) (-109.4576, 51.48106)
 ( 21.57812,-13.64319) (-33.15636,-13.35293) ( 46.15070, 8.904151)
 (-23.28302,-33.85237)

In em9outs2.dat:
the first 10 values are
0.127711	0.127714	0.127724	0.127740	0.127763
0.127793	0.127828	0.127871	0.127920	0.127975
You can use these data to check the success of your work.


Acknowledgment:  
This work is sponsored by DARPA under contract #DABT63-91-k-0005.
The content of the information does not necessary
reflect the position or the policy of the Government and no official
endorsement should be inferred.


                                     Good luck!