#DEFINITIONS:  -*-sh-*-
#This file was generated by arts_file_from_command_list 
#in the PyARTS python module

Main{

output_file_formatSetAscii{
}

#### LOAD DATA: these files were created with MCDataPrepare.arts ######

ReadXML( f_grid ){
filename="data/f_grid.xml.generated"
} 

IndexSet(f_index){
value=0
}

ReadXML(p_grid){
filename="data/p_grid.xml"
}

AtmosphereSet3D{
}

ReadXML(lat_grid){
filename="data/lat_grid.xml"
}

ReadXML(lon_grid){
filename="data/lon_grid.xml"
}

ReadXML( t_field ){
filename="data/t_field.xml.generated"
} 

ReadXML( z_field ){
filename="data/z_field.xml.generated"
} 

ReadXML( vmr_field ){
filename="data/vmr_field.xml.generated"
} 

r_geoidSpherical{
r=-1
}

ReadXML(z_surface){
filename="data/z_surface.xml.generated"
}

ReadXML(abs_lookup){
filename="data/gas_abs_lookup_adapted.xml.generated"
}

abs_speciesSet{
species=["O2","N2","H2O"]
}

IndexSet (abs_lookup_is_adapted ){
value=1
}

ReadXML(cloudbox_limits){
filename="data/cloudbox_limits.xml.generated"
}

IndexSet(cloudbox_on){
value=1
}

ReadXML(pnd_field){
filename="data/pnd_field.xml.generated"
}

ReadXML(scat_data_mono){
filename="data/scat_data_mono.xml.generated"
}

#### Define Agendas #################################################

AgendaSet(surface_prop_agenda){
  Ignore(rte_los){}
  InterpAtmFieldToRteGps( surface_skin_t, t_field ){}
  surfaceBlackbody{
  }
}


AgendaSet(abs_scalar_gas_agenda){
  abs_scalar_gasExtractFromLookup{
}

}

AgendaSet( emission_agenda ){
   emissionPlanck{}
}

AgendaSet(opt_prop_gas_agenda){
  ext_matInit{
  }
  abs_vecInit{
  }
  ext_matAddGas{
  }
  abs_vecAddGas{
  }
}

AgendaSet(spt_calc_agenda){
  opt_prop_sptFromMonoData{
  }
}

AgendaSet(iy_space_agenda){
  Ignore(rte_pos){}
  Ignore(rte_los){}
  MatrixCBR(iy,f_grid){
  }
}

AgendaSet(ppath_step_agenda){
  ppath_stepGeometric{
    lmax=-1
  }
}

AgendaSet(rte_agenda){
  RteStd{
  }
}


#### Set some Monte Carlo parameters ###################################


ReadXML(mc_incoming){
filename="data/mc_incoming.xml"
}

IndexSet(stokes_dim){
value=4
}

MCSetSeedFromTime{}

Print(mc_seed){
level=1
}

AgendaSet( iy_cloudbox_agenda ){
  Ignore (rte_gp_p) {}
  Ignore (rte_gp_lat) {}
  Ignore (rte_gp_lon) {}
  ScatteringMonteCarlo{
    std_err=0.5
    max_time=20
    max_iter=-1
    incoming_lookup=1
    z_field_is_1D=1
  }

  #### Print calculated std. err., number of photons, #####

  VectorToTbByRJ(mc_error,mc_error) {}

  Print(mc_error){
    level=1
  }

  Print(mc_iteration_count){
    level=1
  }
}

#### Define viewing position and line of sight #########################

rte_losSet{
za=99.7841941981
aa=180
}

rte_posSet{
r_or_z=95000.1
lat=7.61968838781
lon=0
}

rte_posAddRgeoid{
}

Matrix1RowFromVector(sensor_pos,rte_pos){
}

Matrix1RowFromVector(sensor_los,rte_los){
}

#### Some other essentially pointless stuff ##########################

AntennaSet1D{
}

IndexSet(nelem){
value=1
}

VectorSet(mblock_za_grid){
value=0
}

sensorOff{
}

jacobianOff{
}

#### Perform RT calculation

RteCalc{
}

#### Print calculated Stokes vector #####

VectorToTbByRJ(y,y){
}

Print(y){
level=1
}

}