#DEFINITIONS:  -*-sh-*-

Arts2 {

StringCreate( datapath )
ArrayOfStringCreate( filenames )
ArrayOfStringCreate( fieldnames )


# ECMWF data (chevallier and similar)
ArrayOfStringSet( filenames, [
  "chevallierl91_all_ccol.xml.gz",
  "chevallierl91_all_q.xml.gz",
  "chevallierl91_all_t.xml.gz",
  "chevallierl91_all_oz.xml.gz",
  "chevallierl91_all_rcol.xml.gz",
] )

ArrayOfStringSet( fieldnames, [
                 "T", "z",
                 "scat_species-LWC-mass_density", "scat_species-IWC-mass_density",
                 "scat_species-RR-mass_flux", "scat_species-SR-mass_flux",
                 "abs_species-H2O", "abs_species-O3"
] )




###########################################

StringCreate( fileloc )
StringCreate( filename )
Copy( filename, datapath )
IndexCreate( ncases )
nelemGet( ncases, filenames )

# Some prep settings
#-------------------
output_file_formatSetZippedAscii
#output_file_formatSetAscii
IndexSet (atmosphere_dim, 1)
ArrayOfMatrixCreate( arrayofmatrix_1 )



AgendaSet( forloop_agenda ){
  Extract( fileloc, filenames, forloop_index )
  Append( filename, fileloc )

  # Read ArrayOfMatrix batch data
  #------------------------------
  ReadXML( arrayofmatrix_1, filename )


  # Convert to batch_atm_fields_compact
  # -----------------------------------
  batch_atm_fields_compactFromArrayOfMatrix(
    atmospheres_fields=arrayofmatrix_1,
    field_names=fieldnames )

  WriteXML( in=batch_atm_fields_compact, filename=filename )
}

IndexStepDown( ncases, ncases )
ForLoop( forloop_agenda, 0, ncases, 1 )

} # End of Main