Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Observer End-to-End Test for OMPS NP #76

Open
4 tasks done
emilyhcliu opened this issue Jul 3, 2023 · 2 comments
Open
4 tasks done

Observer End-to-End Test for OMPS NP #76

emilyhcliu opened this issue Jul 3, 2023 · 2 comments
Assignees

Comments

@emilyhcliu
Copy link

emilyhcliu commented Jul 3, 2023

End-to-End Tests for OMPS NP Processing

  • Create BUFR converter - python script
  • Test BUFR converted IODA obs input to H(x)
  • Test BUFR converted IODA obs input to H(x) + initial error assignment + read routine QC
  • Test BUFR converted IODA obs input to H(x) + initial error assignment + read routine QC + setup routine QC

Issues found in the existing data filtering yaml (omi_aura.yaml) are documented here:

  • Need to exclude total ozone data from the profile data (add a filter to remove total ozone data)
  • Gross check needs to be updated - obs and gross errors changed for top 5 layers
  • Top five layers are assimilated
  • Quality control variable types do not comply with IODA convention --- they need to be integers
  • Quality control variable names do not comply with IODA convention --- should be camel case

Updates made and tested

  • Added filter to remove total ozone data
  • Updated observation error assignment
  • Updated gross check
  • Removed the filter, which excludes the top 5 layers
  • The quality control flags (toq and poq) are integers. Modified the filters accordingly in the YAML
  • The quality control flags names do not comply IODA convention (should be camel case).
    SBUVTOQ (toq) - total ozone quality ---> totalOzoneQuality
    SBUVPOQ(poa) - profile ozone quality ---> profileOzoneQuality

QC Procedures: (Data thinning not applied)

read routine QC

  • Remove total ozone data from profile data
  • Observation value sanity check
  • Total Ozone quality check
  • Profile Ozone quality check

setup routine QC

  • Background departure gross check |O-F| > min(10*obserr, 300)
@emilyhcliu emilyhcliu self-assigned this Jul 9, 2023
@emilyhcliu
Copy link
Author

emilyhcliu commented Jul 9, 2023

HofX (with total column removed) Comparison: UFO vs. GSI for each pressure layer

Pressure (hPa) from High to Low (1013, 640, 400, 250, 160, 101, 64, 40, 25, 16, 10, 6, 5, 2.5, 1.6, 1, 0.5, 0,4, 0.25, 0.16, 0.1]

UFO O-F GSI O-F
ufo_ompsnp_npp_omf_ozoneLayer_totozchk_101325 gsi_ompsnp_npp_omf_ozoneLayer_noqc_101325
ufo_ompsnp_npp_omf_ozoneLayer_totozchk_63936 gsi_ompsnp_npp_omf_ozoneLayer_noqc_63936
ufo_ompsnp_npp_omf_ozoneLayer_totozchk_40327 gsi_ompsnp_npp_omf_ozoneLayer_noqc_40327
ufo_ompsnp_npp_omf_ozoneLayer_totozchk_25432 gsi_ompsnp_npp_omf_ozoneLayer_noqc_25432
ufo_ompsnp_npp_omf_ozoneLayer_totozchk_16009 gsi_ompsnp_npp_omf_ozoneLayer_noqc_16009
ufo_ompsnp_npp_omf_ozoneLayer_totozchk_10132 gsi_ompsnp_npp_omf_ozoneLayer_noqc_10132
ufo_ompsnp_npp_omf_ozoneLayer_totozchk_6393 gsi_ompsnp_npp_omf_ozoneLayer_noqc_6393
ufo_ompsnp_npp_omf_ozoneLayer_totozchk_4032 gsi_ompsnp_npp_omf_ozoneLayer_noqc_4032
ufo_ompsnp_npp_omf_ozoneLayer_totozchk_2543 gsi_ompsnp_npp_omf_ozoneLayer_noqc_2543
ufo_ompsnp_npp_omf_ozoneLayer_totozchk_1600 gsi_ompsnp_npp_omf_ozoneLayer_noqc_1600
ufo_ompsnp_npp_omf_ozoneLayer_totozchk_1013 gsi_ompsnp_npp_omf_ozoneLayer_noqc_1013
ufo_ompsnp_npp_omf_ozoneLayer_totozchk_639 gsi_ompsnp_npp_omf_ozoneLayer_noqc_639
ufo_ompsnp_npp_omf_ozoneLayer_totozchk_403 gsi_ompsnp_npp_omf_ozoneLayer_noqc_403
ufo_ompsnp_npp_omf_ozoneLayer_totozchk_254 gsi_ompsnp_npp_omf_ozoneLayer_noqc_254
ufo_ompsnp_npp_omf_ozoneLayer_totozchk_160 gsi_ompsnp_npp_omf_ozoneLayer_noqc_160
ufo_ompsnp_npp_omf_ozoneLayer_totozchk_101 gsi_ompsnp_npp_omf_ozoneLayer_noqc_101
ufo_ompsnp_npp_omf_ozoneLayer_totozchk_63 gsi_ompsnp_npp_omf_ozoneLayer_noqc_63
ufo_ompsnp_npp_omf_ozoneLayer_totozchk_40 gsi_ompsnp_npp_omf_ozoneLayer_noqc_40
ufo_ompsnp_npp_omf_ozoneLayer_totozchk_25 gsi_ompsnp_npp_omf_ozoneLayer_noqc_25
ufo_ompsnp_npp_omf_ozoneLayer_totozchk_16 gsi_ompsnp_npp_omf_ozoneLayer_noqc_16
ufo_ompsnp_npp_omf_ozoneLayer_totozchk_10 gsi_ompsnp_npp_omf_ozoneLayer_noqc_10

@emilyhcliu
Copy link
Author

emilyhcliu commented Jul 9, 2023

HofX (after QC) Comparison: UFO vs. GSI for each pressure layer

Pressure (hPa) from High to Low (1013, 640, 600, 250, 160, 101, 64, 40, 25, 16, 10, 6, 5, 2.5, 1.6, 1, 0.5, 0,4, 0.25, 0.16, 0.1]

UFO O-F vs. GSI O-F @ 100 hPa
ufo_ompsnp_npp_omf_ozoneLayer_qc_10132
gsi_ompsnp_npp_omf_ozoneLayer_qc_10132

Plots for other layers are here

@emilyhcliu emilyhcliu changed the title UFO Observer End-to-End Test for OMPS NP Observer End-to-End Test for OMPS NP Aug 2, 2023
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant