; Observational parameters
lambda    =    2.196       ;; wavelength of image (micron)

; Model parameters
Lstar     =  5.5d6         ;; Luminosity of star (in Lsun)
Rstar     =     60.        ;; radius of central star (in Rsun) 59.7=Rcore John
Rin       =      0.1       ;; smallest inpact parameter (in Rstar)
Rout      =  59.3d3        ;; outer radius of model (in Rsun) john:59.3d3
distance  =   2300.        ;; distance to object (in parsec)
tau_core  =    0.667       ;; 'core radius' defined as tau_ff=tau_core
include_ff =   x           ;; include ff opacity
include_es =   x           ;; include e- scattering opacity

; Wind parameters:
Mdot      =    0.0016      ;; Mass Loss Rate [Msun/yr]
Twind     =   10d3         ;; Temperature of the wind [K]
v_inf     =    500d0       ;; Terminal velocity of wind
v0        =     50d0       ;; Initial velocity of wind at stellar surface
Beta      =    1.0         ;; Beta of velocity law
f         =    0.256       ;; Constant clumping factor
use_f_law =                ;; Use velocity dependant clumping fac. (like John)
fmin      =    0.30625     ;; minimum clumping
v_cl      =    100.        ;; 'clumping velocity' (in km/s)
X         =    0.56        ;; mass fraction of H
Y         =    0.44        ;; mass fraction of He
qH        =    1.0         ;; Ionisation fraction of Hydrogen
qHe       =    0.0         ;; Ionisation fraction of Helium
; Grid parameters
np_inRstar=     16         ;; number of impact parameters for p < Rstar
nR        =    256         ;; number of points in Radial direction
Rfine_fac =     16         ;; for determining best R grid, initial grid is
;; taken logarithmic with Rfine_fac*nR points

; General constants
mH        =  1.67d-24      ;; mass of H atom [g]
Rsun      =  7d10          ;; Radius of sun [cm]
sig_thom  =  6.65d-25      ;; Thomson cross section of free electron [cm2]
Msun      =  1.98892d33    ;; Solar mass [g]
Lsun      =  3.846d33      ;; Solar luminosity [erg/s]
year      =  3.15576d7     ;; lenght of 1 year [s]
sig_SB    =  5.6705d-5     ;; Stefan-Boltzmann constant [ers/s/cm2/K4]
c         =  3d10          ;; speed of light [cm/s]
h         =  6.626176d-27  ;; Planck constan [erg s]
k         =  1.3807d-16    ;; Boltzmann constant [erg/K]
au        =  1.496d13      ;; Astronomical Unit [cm]
pc        =  3.086d18      ;; parsec [cm]

gIII      =  1.2           ;; ff gaunt factor to be used
gII       =  0.94          ;; bf gaunt factor to be used

; PARAMETERS for IMAGE and VISIBILITY CURVE COMPUTATION
make_im_uv =  x            ;; compute a synthetic image and visibilty curve?
npix       = 1024          ;; image of wind will have nxn pixels
pixelsize  = 0.25          ;; pixelsize in mas
npix_ft    = 1024          ;; number of pixels to be used in fourier transform
phi        = 0.            ;; baseline angle (degrees), at the moment obsolete
resolved_f = 0.43          ;; fraction of resolved flux
vis_ps     =               ;; make a postscript in POSTSCRIPST?
vis_dataf  =  x            ;; make a data file of vis curve?



;john_file  =  IP_2P2_1min3 ;; file where John's intensities are
;john_scale = 1.268         ;; spatial scale factor for john's model
;john_file  =  IP_2P2_1_25min3 ;; file where John's intensities are
;john_scale = 1.146         ;; spatial scale factor for john's model
 john_file  =  ip_2p2_mum1.5e-3.dat ;;
;john_scale = 1.053         ;; spatial scale factor for john's model
 john_scale = 1.0365        ;; spatial scale factor for john's model
plot_john   =     x         ;; plot also john's model, 1.25e-3 (scaled)
