!  Convert watts/cm^2 to volts/meter for a material of
!  known constitutive parameters.

!  Program inputs are:
!    Material complex permittivity Ke' and Ke"
!    Material complex permeability Ku' and Ku"
!    Irradiance in watts/cm2

!  Program outputs are:
!    Program input values
!    Material complex index of refraction
!    Volts/meter in material for specified irradiance

!  For my purposes I found it easier to just key the input values 
!  directly into the macro rather than use input boxes.  Feel free 
!  to modify this macro as you see fit.

!  Author:  Mike I. Jones  3/24/95  updated 8/28/99

OUTPUT SCREEN

PI = 3.141592653589793

! VELOCITY OF LIGHT IN A VACUUM IN METERS/SEC.

VC = 2.99792458E8

! PERMEABILITY OF FREE SPACE IN HENRYS/METER

U0 = 4 * PI * 1E-7

! PERMITTIVITY OF FREE SPACE IN FARADS/METER

E0 = 1E7 / (4 * PI * VC * VC)

FREESPACEIMPEDANCE = SQRT( U0 / E0 )

!  POWERCM2 = IRRADIANCE IN WATTS/CM^2

POWERCM2 = 500

!  POWERM2 = IRRADIANCE IN WATTS/M^2

POWERM2 = POWERCM2 * 10000

!  CONSTITUTIVE PARAMETERS

KE1 = 9.00
KE2 = 0.001

KU1 = 1
KU2 = 0

FORMAT 10.6
PRINT
PRINT "CONSTITUTIVE PARAMETERS:"
PRINT
PRINT "ke'   = ",KE1,"    ke'' = ",KE2
PRINT
PRINT "ku'   = ",KU1,"    ku'' = ",KU2
PRINT

!  CONVERT ke', ke", ku', ku" TO COMPLEX INDEX

Q1 = SQRT(  ( KE1*KE1 + KE2*KE2 ) * ( KU1*KU1 + KU2*KU2 ) )

Q2 = KE1 * KU1 - KE2 * KU2

INDEX = SQRT ( Q1 + Q2 ) / SQRT(2)

Q3 = KE1*KU2 + KE2*KU1

IF Q3==0 
    ABSORPTION_INDEX = 0
  ELSE
    ABSORPTION_INDEX = ABSO(Q3)/Q3 * SQRT ( Q1 - Q2 ) / SQRT(2)
ENDIF

PRINT "COMPLEX INDEX OF REFRACTION:"
PRINT
PRINT "INDEX = ",INDEX,"    K    = ",ABSORPTION_INDEX
PRINT

FORMAT 10.3
PRINT "IRRADIANCE IN WATTS/CM^2 = ",POWERCM2
PRINT

VPM = SQRT( 2 * FREESPACEIMPEDANCE ) * SQRT( POWERM2 / INDEX )

PRINT "VOLTS / METER = ",VPM

END