# Just calculate absorption coefficients for one continuum or full model
# and write results to a file.
#
# 2003-11-16 Stefan Buehler

Main {

# Set output format to ASCII:
output_file_formatSetAscii{}

# Create frequency grid:
VectorNLinSpace(f_grid){
  start =  58.4446e9
  stop  =  58.4486e9
  n     =  1000
#  start =  30.0000e9
#  stop  =  90.0000e9
#  n     =  120
#  start =  115.0000e9
#  stop  =  125.0000e9
#  n     =  10
}

# ... and write to file.
WriteXML(f_grid){""}

# 0=false 1=true
ZeemanO2Settings{
	         ZeemanO2OnOff = 1
	         ZeemanO2PressureLimit = 1.000e3
	        }

# Create a pressure vector [Pa]:
NumericSet(rte_pressure){1.200000e+00}

# Create a temperature vector [K]:
NumericSet(rte_temperature){1.741000e+02}

# Create a vmr vector [absolute numbers]:
#                                  O2        H2O
VectorSetExplicitly(rte_vmr_list){[0.2100e0, 1.000000e-06]}

ReadXML(species_index){"species_index.xml"}



# Write these to files:
#WriteXML(){""}
#WriteXML(abs_t){""}
#WriteXML(abs_vmr){""}

IndexSet(stokes_dim){
	value = 4
}

IndexSet(f_index){
	value = -1
}

#  AgendaSet(zeeman_prop_agenda)
# {
# 	StringSet(abs_model){"PWR98"}
# 	VectorSetExplicitly(abs_user_parameters){[]}
# 	absO2ZeemanModel{}
# } 
  	


#  AgendaSet(zeeman_prop_agenda)
#   {
#	StringSet(abs_model){"PWR98"}
#	VectorSetExplicitly(abs_user_parameters){[]}
#	absO2ZeemanModel{}
#   }
  
  AgendaSet( opt_prop_gas_agenda ) {
    StringSet(abs_model){"PWR98"}
    VectorSetExplicitly(abs_user_parameters){[]}
    absO2ZeemanModel{}
    ext_matInit{}
    abs_vecInit{}
    ext_matAddGasZeeman{}
    abs_vecAddGasZeeman{}
	
  }

ReadXML(geomag_los){"geomag_los.xml"}
WriteXML(geomag_los){""}

IndexSet(ppath_index){0}

 test_zeeman{}



# ------------------------------------------------------------------
# O2 Zeeman absorption model
#StringSet(abs_model){"PWR98"}
#VectorSetExplicitly(abs_user_parameters){[]}
#absO2Model{}
#WriteXML(abs){"absO2Model_file.xml"}
# ------------------------------------------------------------------
WriteXML(ext_mat_zee){""}
WriteXML(abs_vec_zee){""}
#WriteXML(abs_vec){""}
#WriteXML(ext_mat){""}
}