# Atmospheric Radiation

To begin, we’ll load `foqat`.

    library(foqat)

## Calculate TUV in batch online

TUV models are available in both online and offline versions, but both require running on a daily basis (i.e., resetting parameters every day of the run).

`tuv_batch()` runs the online batch TUV by reading a time series of the proposed input parameters and their values, and aggregates the results to a new data frame.

`tuv_batch()` automates the online TUV calculation (https://www.acom.ucar.edu/Models/TUV/Interactive_TUV/).

Four output modes are supported:
Output mode 2: Molecular photolysis frequencies (109 photoreactions). Molecular photolysis frequencies (109 photoreactions).
Output mode 3: Weighted irradiance (27 weighting functions).
Output mode 4: Spectral actinic flux.
Output mode 5: Spectral irradiance.

    #Examples.
    #output molecular photolysis frequency (109 light reactions)
    #Time range: March 1, 2021
    #Resolution: hour
    #Ground elevation 1 km
    #Measured altitude 1.05 km
    #Longitude 109.747144 °
    #latitude 38.298267 °
    #Ozone column concentration 306 DU

    df=data.frame(
      time=seq(as.POSIXct("2021-03-01 00:00:00",tz="GMT"), as.POSIXct("2021-03-01 3:00:00",tz="GMT"), by="hour"),
      gAltitude=rep(1, 4),
      mAltitude=rep(1.05, 4),
      longitude=rep(109.747144, 4),
      latitude=rep(38.298267, 4),
      ozone=rep(306, 4)
    )
    tuv_df=tuv_batch(df, inputMode=0, outputMode=2, nStreams=-2)
    #> [1] "Running TUV on the inputs:"
    #>
      |
      |                                                                      |   0%
      |
      |==============                                                        |  20%
      |
      |============================                                          |  40%
      |
      |==========================================                            |  60%
      |
      |========================================================              |  80%
      |
      |======================================================================| 100%
    #Show columns 1 to 3
    head(tuv_df[,1:3])
    #>   time PHOTOLYSIS RATES    1 O2 -> O + O
    #> 1    0            1/sec                                      0.000E+00
    #> 2    1            1/sec                                      0.000E+00
    #> 3    2            1/sec                                      0.000E+00
    #> 4    3            1/sec                                      0.000E+00