------------------------------ CORHEL Parameters -----------------------------
 Carrington rotation    : 1913
 Resolution             : Low
 Case name              : corona_mas
 Simulation model       : Corona
 Corona model           : MAS
 ------------------------------------------------------------------------------


 ------------------------------ Corona Parameters -----------------------------
 Observatory            : Kitt Peak
 Polefitting - Th 0     : 0.25 radian
 Polefitting - Th 1     : 0.45 radian
 Polefitting - Th R     : 0.35 radian
 Maximum mode           : 9
 Number of filter passes: 3
 Minimum image scale    : -50.00 Gauss
 Maximum image scale    : 55.00 Gauss
 Palette                : Seismic
 List of output plots   : br,bt,bp,jr,jt,jp,rho,t,p,vr,vt,vp
 Maximum time           : 120 code units
 Output interval        : 120 code units
 Minimum time step      : 0.005 code units
 Maximum time step      : 0.050 code units
 Lundquist number       : 1.00e3
 Viscosity              : 0.01
 Density                : 2.00 x 10^8 cm^-3
 Temperature            : 1.80e6 K

 * 1 code units in time = 24 minutes.
 ------------------------------------------------------------------------------

--14:21:55--  ftp://nsokp.nso.edu/kpvt/synoptic/mag/m1913f.fits
           => `m1913f.fits'
Resolving nsokp.nso.edu... 140.252.1.61
Connecting to nsokp.nso.edu|140.252.1.61|:21... connected.
Logging in as anonymous ... Logged in!
==> SYST ... done.    ==> PWD ... done.
==> TYPE I ... done.  ==> CWD /kpvt/synoptic/mag ... done.
==> PASV ... done.    ==> RETR m1913f.fits ... done.
Length: 138,240 (135K) (unauthoritative)

    0K .......... .......... .......... .......... .......... 37%  148.69 KB/s
   50K .......... .......... .......... .......... .......... 74%  412.61 KB/s
  100K .......... .......... .......... .....                100%  388.07 KB/s

14:21:56 (246.51 KB/s) - `m1913f.fits' saved [138240]

 
 ### KPSYN, version 1.07, 11/02/2007 ###
 
 Reading file: kittpeak1913.hdf
 
 ### Kitt Peak synoptic file dimensions:
 
 Number of points in longitude =          360
 Number of points in latitude  =          180
 
 ############################# SELECT OUTPUT MESH #############################
 
 Specify the type of output latitude mesh:
 
 y ==> use the original "Wilcox data" latitude mesh
 n ==> use a uniform latitude mesh
 
 Use Wilcox latitude mesh? [y|n]:  
 ######################### OUTPUT RAW KITT PEAK DATA ##########################
 
 ### Kitt Peak (raw) Br ###
 Enter HDF output file name [CR to skip]:  
 Writing Kitt Peak (raw) Br to file: br_raw.hdf
 
 ############################ POLAR FIELD FITTING #############################
 
 Would you like to fit the polar field? [y|n]:  
 The data in the theta range [TH0,TH1] is used to fit the lowest-m
 components of Br with a two-term or three-term geometric Taylor
 expansion near the poles.
 
 Br is reset to the fitted field for theta in the range [0,THR] 
 for each pole.
 
 ### Fitting region ###
 NORTH POLE: Enter TH0, TH1 [radians]:  
 ### Fitting region ###
 SOUTH POLE: Enter TH0, TH1 [radians]:  
 ### Region to reset ###
 NORTH POLE: Enter THR [radians]:  
 ### Region to reset ###
 SOUTH POLE: Enter THR [radians]:  
 ### Type of fit ###
 
 y ==> use a two-term Taylor series fit
 n ==> use a three-term Taylor series fit
 
 Use a two-term fit? [y|n]:  
 ### Maximum m-mode ###
 Maximum m-mode to fit:  
 ### Fitting the polar fields ...
 
 ### From Kitt Peak synoptic data:
 
 Br (monopole) =  -0.1759853    
 
 ####################### DIAGNOSTICS AFTER POLE FITTING #######################
 
 ### Kitt Peak (corrected) Br ###
 Enter HDF output file name [CR to skip]:  
 Writing Kitt Peak (corrected) Br to file: br_polefit.hdf
 
 ############################### INTERPOLATION ################################
 
 ### Interpolating Br ...
 
 ### After interpolation:
 
 Br (monopole) =   1.4003520E-04
 
 ### Computing the Fourier transform of Br ...
 
 ### RMS value of Br (before filter) =    7.603417    
 
 ### Br_rms vs. m (before filter) ###
 Enter text output file name [CR to skip]:  
 Writing Br_rms to file: br_rms_before_filt.dat
 
 ### Sine part of Fourier modes of Br ###
 Enter HDF output file name [CR to skip]:  
 Writing Br_sin to file: brs.hdf
 
 ### Cosine part of Fourier modes of Br ###
 Enter HDF output file name [CR to skip]:  
 Writing Br_cos to file: brc.hdf
 
 ############################ LONGITUDINAL FILTER #############################
 
 ### Longitudinal filter parameters:
 
 Filter type: exp[-(0.7*m/MMAX)**8)]
 
 MMAX specifies the mode that is attenuated to 95% of its original amplitude.
 Modes with M smaller than MMAX are not attenuated appreciably.
 The mode with M = 1.64*MMAX is attenuated to 5% of its original amplitude.
 If filtering is not desired, set MMAX=0.
 
 Enter MMAX:  
 ### Applying the longitudinal filter ...
 
 ### RMS value of Br (after filter) =    5.397908    
 
 ### Br_rms vs. m (after filter) ###
 Enter text output file name [CR to skip]:  
 Writing Br_rms to file: br_rms_after_filt.dat
 
 ################################ FINAL OUTPUT ################################
 
 File name for interpolated Br:  
 Writing interpolated Br to file: br_m9.hdf
 
 #################################### DONE ####################################
 
  
 Code: MAS
 Version: 0.66d_ccmc_v01
 Updated on: 04/30/2003
 Source file: mas47d_ccmc_v01.F
 Run ID: mas
 Run started on: 2009/04/07
 Run started at: 14:22:00
 Ran on machine: zmac.predsci.com
 Machine type: Darwin
  
 ### Parameters:
  
 NR =     61
 NT =     71
 NP =     64
 NHISTMAX =   2002
  
 ### Input file contents:
  
  &invars
  option='streamer'
  fldtype='potential' bingauss=.true. bnfile='br_m9_filt3.dat'
  eqtype='parker' np1d=130 onedfile='parker1.8.pw' rhor0=2.00
  bcr0type='1dchar' bcr1type='1dchar'
  tmax=120 ntmax=40000 dtmax=0.050 dtmin=0.005
  ifideal=0 slund=1.00e3 visc=0.01
  rsifile=' '
  rl=29. g0=.823 ifrho=1 iftemp=1 ifvdgv=1 ifpc=1
  rfrac=.4,1. drratio=60.,4. nfrmesh=5
  tfrac=.5,1. dtratio=.5,2. nftmesh=5
  mmodes=0,1,2,3,4,5,6,7,8,9
  tpltxint=120. ihistint=1 ifprec=1 trsdump=20.
  upwindv=1. cfl=.4 isitype=1
  dformat='hdf'
  plotlist='br','bt','bp','jr','jt','jp','rho','t','p',
   'vr','vt','vp'
  tnode=0.,10.0 vnode=0.,0. ishearprof=2
  arotate=1. brotate=0. crotate=0.
  ihst=23 jhst=33 khst=1
  parchar=.false. ubzero=.true. nfiltub=3
  he_frac=0. radloss=0. tcond=0. ifaw=0
  tnode_ch=0.,500. q0phys_ch=0.,0.
  tbc0=1.80e6 tbc1=0.
  upwinda=0. upwindar=1. upwindat=1. upwindap=1.
  upwindphi=1. strmax=.4
  fdphi=.false.
  si_aggressive=.true.
  /
 
 
 >running
  
 Code: MAS
 Version: 0.66d_ccmc_v01
 Updated on: 04/30/2003
 Source file: mas47d_ccmc_v01.F
 Run ID: mas
 Run ended on: 2009/04/07
 Run ended at: 16:06:18
 Ran on machine: zmac.predsci.com
 Machine type: Darwin
  
 ### End of run ...
  
 NTIME =   3109
 TIME =  1.2001248E+02
  
 ### CPU time used =    6243.33 seconds.