#!/bin/csh
setenv OPENGGCMDIR $HOME/GGCM/head
$OPENGGCMDIR/bin/script.runme `pwd` `basename $0` 0; exit
#....... template
#--------------------------------------------------------------
#  This is the magic file that makes global MHD runs happen.
#  Don't change the first two lines -- they do the magic.
#  Everything below here is input to the model; obviously,
#  lines that start with a '#' are comments and blank lines
#  are ignored.
#  Unified input -- JR, IGPP/UCLA  12/28/2000
#--------------------------------------------------------------

#DESCRIPTION 1st test-run -avapirev

#------------- some top level options ----------------------------
RUN             default      #  run id, default sets to current dir
DOSMF           1            #  switch to create parameter database
DOINP           1            #  switch to create input file
DOGRD           1            #  switch to create the grid
KEEPTMP         1            #  switch to keep tmp files for debugging
TARGET          zaphod-eth   #  prepare run for which machine?
MHDCOORD        MHD

#------------- compiling, profiling and debugging ----------------
MPIF77_FLAGS    -g -fast     #  options for compilaton #-g -v -O -fast -Msave -Wl -Bstatic -Mextend #-Mextend -fast -v
LOCALF77        pgf77        #  the local f77 compiler

check_mhdvar    false            #  check mhd variables for sane values every step
CHECKSUMS       false            #  write checksums to compare code changes

#------------- directories ---------------------------------------
GGCMDIR         default                  # base directory for the code
OUTDIR          default                  # where the output goes

#------------- the code input files ------------------------------
GLOBU           default      #  code input (utilities), default: latest
GRIDA           default      #  ditto, grid/x-grid generation
GRIDB           default      #  ditto, grid/partitioning
GRIDC           default      #  ditto, grid/mapping coefficients
MHDSC           default      #  ditto, constants and definitions
MHDCa           default      #  ditto, MHD core code
MHDMI           default      #  ditto, MI coupling code/sigmas/CTIM interf.
IOSOL           default      #  ditto, ionosphere potential solver
IOPOS           default      #  ditto, ionosphere variables post-prpcessing
COTR1           default      #  ditto, coord. transf. and time routines
MACHS           default      #  ditto, machine dependent stuff
MINVA           default      #  ditto, minimum variance program
CTIM            none         #  1.0.10       #  ctim version (1.0.10) or "none" for sigaa2 sigma computation


#-------------- define the grid ------------------------------------

#---------------------------------------------------------------
#------  grid definition for 'tail1' 
#---------------------------------------------------------------
npx               7
npy               3
npz               1
nx                315
ny                90 
nz                100
dy0               0.30
dz0               0.30
xx1               -24.01
xx2               350
yy2               48.01
zz2               48.01
#.... how the x-grid is defined 
#     f = fak*(hm + (h0-hm)*tanhm(b1*(x-x1)) + (hn-hm)*tanhp(b2*(x-x5))
#    *  - (hm-hmm)*(#one-tanhm(b3*(x-(x3-dmm)))-tanhp(b3*(x-(x3+dmm)))) )
ggx1              $v="-n   $NX     -x0 $XX1     -xn $XX2 "
ggx2              $v="-x1  -26     -x3   -4     -dmm  8    -x5   50 "    
ggx3              $v="-h0  2.0     -hn   43    -hm 1.7    -hmm 0.70 "
ggx4              $v="-b1  0.15    -b2  0.025   -b3  0.3 "
#.... for tail:
#ggx1              $v="-n   $NX     -x0 $XX1     -xn $XX2 "
#ggx2              $v="-x1  -26     -x3   20     -dmm  16    -x5   80 "
#ggx3              $v="-h0  1.0     -hn   5      -hm 0.60   -hmm 0.60 "
#ggx4              $v="-b1  0.15    -b2  0.025   -b3  0.3 "

#-------------- resistivity ----------------------------------------
MAGDIFFU        nl1            # resistivity model res1/nl1
diffthreshold   0.75           # resistivity threshold
diffconstant    0.01           # resistivity constant  

#-------------- output control -------------------------------------
# ':' seperated lists, unit is seconds, +val is only that step, -val is modulus
outputmode      combined       # single/combined (use single only for jobs with many procs on SP2)
outtimerun      -1200000          # frequency for runtime graphics, 0 if not wanted
outtime2d       -60
outtimeio       -60
outtime3d       -1800
outtimeex       18000:1080000:1
outplanex       7:15:30:50:100:130           # output plane location in RE
outplaney       0
outplanez       0
ORBITFILE       orbits.txt
orbit_dt        10
modtty          1              # output line to log file every modtty

#-------------- time step control ----------------------------------
RRMIN           0.10           # minimum permitted number density in cm**-3
SPEEDLIMIT      1500.0         # maximum permitted Alfven speed in km/s
thx             0.40           # cfl time step factor
modnewstep      1              # recalculate cfl limit every 'modnewstep' step

#-------------- Times, Solar Wind / IMF ----------------------------
##.. UT times, use DIPOLTIME "1967:01:01:00:00:0.0" for generic runs with no dipole tilt
STARTTIME  2001:01:30:00:00:27.043    # UT when the run starts 2001  30  0  0 27.043
ENDTIME    2001:01:31:23:59:21.984    # 2001  31 23 59 21.984
DIPOLTIME  1967:01:01:00:00:00.000    #  time for dipole orientation 
alphafak   1.10                 # density factor to account for SW alphas
##
f107       220      # F10.7 solar flux
MOX        228.8    # solar wind monitor position (GSE)
MOY        40.0
MOZ        -22.0
swturb1    false    # solar wind turbulence (true/false)
swturb_rr  0.0      # solar wind trbulence random perturbation for density in 1/cc
swturb_vv  0.0      # solar wind trbulence random perturbation for velocity in km/sec

#-------------- MI coupling and ionosphere -------------------------
ISPHERE       3.5       # inner boundary radius (RE)
io_shield_lat 50.0       # ionosphere shielding latitude
BND_IO_RR     2.0       # density at inner boundary in cc
BND_IO_TT     200       # temperature at inner bndry in eV
NTHETA        361       # ionosphere grid latitude points
NPHI          121       # ionosphere grid longitude points
NIOX          40        # no of panels in theta for one hemisphere
NIOY          8         # no of modes in phi
NIOGO         4         # order of gauss legendre integration
njslide       2         # fac sliding time averages
modiono       20        # recalc ionosphere every soso step
sigmatrans    3000      # after which timestep begin transition to s-c sigma
sig_const     5.00        
sigh_const    5.00        
sigmafak      1.0       # factor to multiply all ionosphere conductances   
mape1         1.0       # MI coupling factor for mapped E
fac_filt      2         # how often to spatially filter FAC
sig_floor     0.5       # minimum sigmas
fak_fac_dyn   1.0       # factor for dynamo currents (only w/ctim)
sig_smooth    2         # spatially smoothing sigmas
sig_h_back    0.3
sig_p_back    0.5       # nightside H/P background conductances
sigfac_fak    1.00      # fudge factor for discrete e-flux
sigdif_fak    1.00      # fudge factor for diffuse  e-flux

BIGDIM        30000000   # global memory

#-------------- CTIM inputs ----------------------------------------
ctim_potfak       0.6   # factor for driving ctim potential
CTIM_IN           ./
ctim_f_prsxxxx    {"$CTIM_IN/inp_ctim.fields.20000921000.asc"}
ctim_f_mmmxxxx    {"$CTIM_IN/inp_ctim.mmm20000921000"}
ctim_f_velxxxx    {"$CTIM_IN/inp_ctim.vel20000921000"}
ctim_f_bmagfld    {"$CTIM_IN/inp_ctim.bmagfld"}
ctim_f_substorm   {"$CTIM_IN/inp_ctim.substorm"}
ctim_f_cormag     {"$CTIM_IN/inp_ctim.cormag"}
ctim_f_datam      {"$CTIM_IN/inp_ctim.datam"}
ctim_f_levpi      {"$CTIM_IN/inp_ctim.levpi"}
ctim_f_gwatts     {"$CTIM_IN/inp_ctim.gwatts"}
ctim_f_cpcpd      {"$CTIM_IN/inp_ctim.cpcpd"}
ctim_f_ionprof    {"$CTIM_IN/inp_ctim.ionprof"}
ctim_f_prof2      {"$CTIM_IN/inp_ctim.prof2"}
ctim_f_dmsp       {"$CTIM_IN/inp_ctim.dmsp"}
ctim_f_holt       {"$CTIM_IN/inp_ctim.holt"}
ctim_f_dmspmod    {"$CTIM_IN/inp_ctim.dmspmod"}
ctim_f_rice91     {"$CTIM_IN/inp_ctim.rice91"}
ctim_f_riceprof   {"$CTIM_IN/inp_ctim.riceprof"}
ctim_f_hprof      {"$CTIM_IN/inp_ctim.hprof"}
ctim_f_hough      {"$CTIM_IN/inp_ctim.hough"}
ctim_f_inshaun    {"$CTIM_IN/inp_ctim.d1d2"}
ctim_i10          5       # I10 is the smothing frequency.
ctim_i5           1       #  i5 is mwt calling frequency.
ctim_i8           6       #  i8 is shaun frequency.
ctim_windmx       1000.0  # maximum wind
ctim_tmpmin       130.0   # minimum temperature
ctim_ampl22       608.0   # amplitudes of the tidal model
ctim_ampl23       500.0
ctim_ampl24       138.0
ctim_ampl25       138.0
ctim_ampl11       200.0
ctim_phase22      9.6     #  phases of tidal model.
ctim_phase23      4.6
ctim_phase24      1.2
ctim_phase25      2.2
ctim_phase11      12.0
ctim_kp           3   #  kp index
ctim_xx           91  #  global setting???
ctim_amp          1.3 #  ratio of E foster field to real E field.
ctim_switch_1    default    #  ctim boolean switches
ctim_switch_2    default
ctim_switch_3    default
ctim_switch_4    default
ctim_switch_5    default
ctim_switch_6    default
ctim_switch_7    default
ctim_switch_8    default
ctim_switch_9    default
ctim_switch_10   default
ctim_switch_11   default
ctim_switch_12   default
ctim_switch_13   default
ctim_switch_14   default
ctim_switch_15   default
ctim_switch_16   default
ctim_switch_17   default
ctim_switch_18   default
ctim_switch_19   default
ctim_switch_20   default
ctim_switch_21   default
ctim_switch_22   default
ctim_switch_23   default
ctim_switch_24   default
ctim_switch_25   default
ctim_switch_26   default
ctim_switch_27   default
ctim_switch_28   default
ctim_switch_30   default
ctim_switch_31   default
ctim_switch_32   default
ctim_switch_33   default
ctim_switch_34   default
ctim_switch_35   default
ctim_switch_37   default
#-------------------------- The End ------------------------------------


