VAX/PC TEST USING QP PROFILE	12-mar092
  1  0.                 OF DUPLICATE W CARDS, THE LAST ONE DOMINATES            
  1  1.                 ordinary ray
  1  -1.                EXTRAORDINARY RAY                                       
 24  78.6        1      LAT. OF NORTH MAGC POLE, DEG NORTH (to agree with GGM)
 25  290.2       1      LONG. OF NORTH MAGC POLE, DEG EAST (to agree with GGM)
  1  0.                 no magnetic field
 24  90.0        1      put the magnetic pole at the geog pole for no-field case
 25  0.          1      ditto
  3  0.                 TRANSMITTER HEIGHT, KM                                  
  4   44.65      1      TRANSMITTER LATITUDE, DEG NORTH                         
  5   296.4      1      TRANSMITTER LONGITUDE, DEG EAST                         
  4   82.60      1      TRANSMITTER LATITUDE, DEG NORTH                         
  5   297.4      1      TRANSMITTER LONGITUDE, DEG EAST                         
  4   32.90      1      TRANSMITTER LATITUDE, DEG NORTH                         
  5   263.5      1      TRANSMITTER LONGITUDE, DEG EAST                         
  7   20.00             INITIAL FREQUENCY, MC/S                                 
  8   29.00             final frequency
  9   9.00                 STEP FREQUENCY                                     
 11   321.2      1      INITIAL AZIMUTH ANGLE, DEGS CLOCKWISE FROM NORTH POLE   
 12   321.2      1      final azimuth angle
 13   0.         1       STEP AZIMUTH ANGLE                                 
 15    0.0       1      INITIAL ELEVATION ANGLE, DEG                            
 16   12.0       1      FINAL ELEVATION ANGLE, DEG                              
 17    4.0       1      STEP IN ELEVATION ANGLE, DEG                            
 20   110.              RECEIVER HEIGHT ABOVE THE EARTH, KM                     
 20   0.                rx height above the earth
 22  6.0
 22  1.                 NUMBER OF HOPS  (+1 FOR N > 1)        
 23  5000.              Maximum no of steps per hop
 41  1.                 Runge Kutta integration
 41  3.                 Adams-Moulton with relative error checking
 42  0.001               max allowable single step error
 42  0.0001              max allowable single step error
 45  10.                max step length
 46  0.0001             minimum step length in km
 46  0.000001           minimum step length in km
 46  0.001              minimum step length in km
 57  2.                 INTEGRATE AND PRINT PHASE PATH                          
 57  0.                  do not do phase path
 58  2.                 INTEGRATE AND PRINT ABSORPTION                          
 58  0.                  do not do absorption
 71  500.0                NUMBER OF STEPS FOR EACH PRINTING                   
 71   50.0                NUMBER OF STEPS FOR EACH PRINTING                   
 71   10.0                NUMBER OF STEPS FOR EACH PRINTING                   
 72  1.                 PUNCH RAYSETS                                           
 81  1.                 PLOT PROJECTION OF RAY PATH ON A VERTICAL PLANE         
 83  40.0        1      LEFT LATITUDE OF PLOT, DEG                              
 84  -105.       1      LEFT LONGITUDE OF PLOT, DEG                             
 85  52.12       1      RIGHT LATITUDE OF PLOT, DEG                             
 86  -81.8       1      RIGHT LONGITUDE OF PLOT, DEG                            
 87  100.0        1     DISTANCE BETWEEN TIC MARKS, KM                          
101  6.5                CRITICAL FREQUENCY, MC/S                                
101  8.977750275        critical freq MHz (1.e12 per cubic metre)
102  300.0              HMAX, KM                                                
103  60.                SCALE HEIGHT, KM                                        
103  100.0              ym for qp layer
104  0.5                     ALPHA CHAPMAN LAYER                                
104  0.                 not a quasi-parabolic layer
104  1.                 quasi-parabolic layer
150  1.                 CALL PERTURBATION SUBROUTINE                            
150  0.                 no perturbation
151  250.               Z0, KM                                                  
152 100.                SH, SCALE HEIGHT, KM                                    
153  0.1                DELTA                                                   
155  100.               LAMBDAX, HORIZONTAL WAVELENGTH, KM                      
156  100.               LAMBDAZ, VERTICAL WAVELENGTH, KM                        
201  0.01               GYROFREQ ON THE GROUND AT THE EQUATOR - no field !!!!
201  0.8                GYROFREQUENCY ON THE GROUND AT THE EQUATOR, MHZ         
201  0.9                GYROFREQUENCY ON THE GROUND AT THE EQUATOR, MHZ (75W)         
251  3.65      E4       COLLISION FREQUENCY AT H1, /SEC                         
252  100.0              H1, REFERENCE HEIGHT, KM                                
253  .148               A1, EXPONENTIAL DECREASE OF NU WITH HEIGHT, /KM         
254  30.                COLLISION FREQUENCY AT H2, /SEC                         
255  140.               H2, REFERENCE HEIGHT, KM                                
256  .0183              A2, EXPONENTIAL DECREASE OF NU WITH HEIGHT, /KM         
300  0.                 data for q-p layer with latitudinal gradients
300  0.5                gradient of FC with latitude (A)
300  0.0                gradient of FC with latitude (A)
301  0.5                gradient of HMAX with latitude (B)
301  0.0                gradient of HMAX with latitude (B)
302  -34.7       1      starting latitude for gradients (TH00)
303  10.                starting FC for gradients (FC0)
304  250.               starting HMAX for gradients (HM0)
305  100.               semi-thickness of q-p (YM)
  0                     (A BLANK IN COL. 1-3 ENDS THE CURRENT  W ARRAY)