marine_gdas_plots.yaml

# template YAML file to create EVA YAML files
# based on obs spaces listed in JEDI YAML files
datasets:
  - name: experiment
    type: IodaObsSpace
    filenames:
      - @FILENAME@
    @CHANNELSKEY@
    groups:
      - name: ObsValue
        variables: &variables @VARIABLES@
      - name: ObsError
      - name: ombg
      - name: oman
      - name: hofx0
      - name: EffectiveQC0
      - name: MetaData
      - name: PreQC
transforms:

  # bkg
  - transform: arithmetic
    new name: experiment::bkg::${variable}
    equals: experiment::ObsValue::${variable}-experiment::ombg::${variable}
    for:
      variable: *variables

  # Generate omb that passed QC for JEDI
  - transform: accept where
    new name: experiment::OmBQC::${variable}
    starting field: experiment::ombg::${variable}
    where:
      - experiment::EffectiveQC0::${variable} == 0
    for:
      variable: *variables

  # Generate oma that passed QC for JEDI
  - transform: accept where
    new name: experiment::OmAQC::${variable}
    starting field: experiment::oman::${variable}
    where:
      - experiment::EffectiveQC0::${variable} == 0
    for:
      variable: *variables

  # Generate oma that passed QC for JEDI
  - transform: accept where
    new name: experiment::hofxQC::${variable}
    starting field: experiment::hofx0::${variable}
    where:
      - experiment::EffectiveQC0::${variable} == 0
    for:
      variable: *variables





graphics:

  plotting_backend: Emcpy
  figure_list:

  # ---------- Map Plots ----------
  # Map plot of OmBQC
  # --------

  - batch figure:
      variables: *variables
      @CHANNELSKEY@
    dynamic options:
      - type: vminvmaxcmap
        data variable: experiment::OmBQC::${variable}
    figure:
      layout: [1,1]
      figure size: [20,10]
      title: 'OmB post QC | @NAME@ @CYCLE@ | ${variable_title}'
      output name: map_plots/@NAME@/${variable}/@CHANNELVAR@/@NAME@_${variable}@CHANNELVAR@OmBQC.png
      tight_layout: true
    plots:
      - mapping:
          projection: plcarr
          domain: global
        add_map_features: ['coastline']
        add_grid:
        add_colorbar:
           label: '${variable}'
        layers:
        - type: MapScatter
          longitude:
            variable: experiment::MetaData::longitude
          latitude:
            variable: experiment::MetaData::latitude
          data:
            variable: experiment::OmBQC::${variable}
            @CHANNELKEY@
          markersize: 0.01
          label: '$(variable)'
          colorbar: true
          # below may need to be edited/removed
          cmap: 'seismic'
          vmin: ${dynamic_vmin}
          vmax: ${dynamic_vmax}


  # Histogram plots
  # ---------------

  # OmA pre and post QC
  - batch figure:
      variables: *variables
    figure:
      layout: [1,1]
      title: 'OmA pre- and post QC | @NAME@ | ${variable_title}'
      output name: histograms/@NAME@/${variable}/oma_pre_post_qc_${variable}.png
    plots:
      - add_xlabel: 'Observation minus anl pre- and post-QC'
        add_ylabel: 'Count'
        add_legend:
          loc: 'upper left'
        statistics:
          fields:
            - field_name: experiment::OmAQC::${variable}
              xloc: 0.5
              yloc: -0.10
              kwargs:
                fontsize: 6
          statistics_variables:
          - n
          - min
          - mean
          - max
          - std
        layers:
        - type: Histogram
          data:
            variable: experiment::OmAQC::${variable}
          color: 'blue'
          label: 'OmA (post QC)'
          bins: 100
          alpha: 0.5
        - type: Histogram
          data:
            variable: experiment::oman::${variable}
          color: 'red'
          label: 'OmA (pre QC)'
          bins: 100
          alpha: 0.5

  # diff between OmA and OmB
  - batch figure:
      variables: *variables
    figure:
      layout: [1,1]
      title: 'OmA and OmB, post-QC | @NAME@ | ${variable_title}'
      output name: histograms/@NAME@/${variable}/oma_omb_histogram_${variable}.png
    plots:
      - add_xlabel: 'OmA and OmB post-QC'
        add_ylabel: 'Count'
        add_legend:
          loc: 'upper left'
        statistics:
          fields:
            - field_name: experiment::OmAQC::${variable}
              xloc: 0.5
              yloc: -0.10
              kwargs:
                fontsize: 6
          statistics_variables:
          - n
          - min
          - mean
          - max
          - std
        layers:
        - type: Histogram
          data:
            variable: experiment::OmAQC::${variable}
          color: 'blue'
          label: 'OmA (post QC)'
          bins: 100
          alpha: 0.5
        - type: Histogram
          data:
            variable: experiment::OmBQC::${variable}
          color: 'red'
          label: 'OmB (post QC)'
          bins: 100
          alpha: 0.5

  - batch figure:
      variables: *variables
      @CHANNELSKEY@
    figure:
      layout: [1,1]
      title: 'Observations vs. JEDI h(x) | @NAME@ @CYCLE@ | ${variable_title}'
      output name: observation_scatter_plots/jedi_hofx_vs_obs_@CYCLE@_@NAME@_${variable}@CHANNELVAR@.png
    plots:
      - add_xlabel: 'Observation Value'
        add_ylabel: 'JEDI h(x)'
        add_grid:
        add_legend:
          loc: 'upper left'
        layers:
        - type: Scatter
          x:
            variable: experiment::ObsValue::${variable}
          y:
            variable: experiment::hofx0::${variable}
          @CHANNELKEY@
          markersize: 1
          color: 'black'
          label: 'JEDI h(x) versus obs (all obs)'
        - type: Scatter
          x:
            variable: experiment::ObsValue::${variable}
          y:
            variable: experiment::hofxQC::${variable}
          @CHANNELKEY@
          markersize: 1
          color: 'red'
          label: 'JEDI h(x) versus obs (passed QC in JEDI)'