C
C Inputs for thunder:
C
C For potential:
C   HISTVOL = 24-char mss path to tiegcm dynamo history volume
C   DDHHMM = model time (3 ints: day,hr,min) of desired history from histvol
C
C     Z1 IS THE HEIGHT OF THE SIGMA 1 SURFACE
C     SIGSL IS THE SEA LEVEL CONDUCTIVITY
C     SIGO IS THE SEA LEVEL CONDUCTIVITY THAT DEFINES THE UPPER
C            EXPONENTIAL CONDUCTIVITY IN REGIONS 1 AND 2
C     ZM IS THE HEIGHT OF THE MAGNETOSPHERE IN METERS
C     ALP DEFINES THE EXPONENTIAL FALL OFF OF CONDUCTIVITY IN
C            REGIONS 1 AND 2
C     DSIGO GIVES THE LATITUDINAL VARIATION OF CONDUCTIVITY ALONG
C            THE SIGMA1 SURFACE
C     IUBCD IS THE PARAMETER THAT SPECIFIES WHETHER THE IONOSPHERIC
C            POTENTIAL IS COUPLED INTO MODEL CALCULATIONS
C            IUBCD = 0 FOR NO IONOSPHERIC POTENTIAL
C            IUBCD = 1 INCLUDES COUPLING WITH IONOSPHERIC POTENTIAL
C     ITILT IS THE PARAMETER THAT GOVERNS GEOMAGNETIC AND GEOGRAPHIC
C            COORDINATES
C            ITILT = 0 GEOGRAPHIC POLE IS COINCIDENT WITH
C                      GEOMAGNETIC POLE
C            ITILT = 1 GEOGRAPHIC AND GEOMAGNETIC POLES ARE
C                      SEPARATED
C     IDELT IS THE PARAMETER CONTROLLING THE PLOTS BELOW THE
C            SIGMA1 SURFACE
C            IDELT = 0 PLOT FULL SOLUTION WITH NEGATIVE REGIONS
C            IDELT = 1 PLOT ONLY POSITIVE REGIONS GIVING SURFACE DETAILS
C     ICLDS IS THE PARAMETER CONTROLLING THE PRESENCE OF CLOUDS
C            ICLDS = 0 NO CLOUDS IN THE CALCULATION
C            ICLDS = 1 CLOUDS INCLUDED GCM OUTPUT
C     ILWRG IS THE PARAMETER CONTROLLING THE SOLUTION TECHNIQUE IN
C            REGION 0
C            ILWRG = 0 ANALYTIC SOLUTION USING EXPONENTIAL FORM
C            ILWRG = 1 TWO POINT BOUNDARY VALUE NUMERICAL SOLUTION
C
 HISTVOL = '/ROBLE/RGR93/TESW01     '
 DDHHMM = 1,0,0
 Z1 = 8.E+3
 SIGSL = 1.14E-14
 SIGO = 1.14E-13
 ZM = 105.E+3
CALP = 1./12.E+3
 ALP = 8.3333E-5
 DSIGO = 0.4
 IUBCD = 1
 ITILT = 1
 IDELT = 1
 ILWRG = 0
 ICLDS = 0
 ENDOFREAD