Skip to content

Commit

Permalink
Prepare GTS snow depth observations for JEDI-based Land DA (#541)
Browse files Browse the repository at this point in the history
* Initiate to process the gts bufr data into the ioda format.

* Change the io paths for matching the global-workflow.

* Modify the ctest for the bufr2ioda yaml files.
  • Loading branch information
jiaruidong2017 authored Jul 24, 2023
1 parent 113ecc8 commit 7630832
Show file tree
Hide file tree
Showing 4 changed files with 97 additions and 3 deletions.
9 changes: 9 additions & 0 deletions parm/land/prep/prep_gts.yaml
Original file line number Diff line number Diff line change
@@ -0,0 +1,9 @@
gtsbufr:
mkdir:
- $(DATA)/obs
copy:
- [$(COM_OBS)/$(OPREFIX)adpsfc.tm00.bufr_d, $(DATA)/obs/]
gtsioda:
copy:
- [$(DATA)/$(OPREFIX)adpsfc_snow.nc4, $(COM_OBS)/$(OPREFIX)adpsfc_snow.nc4]

2 changes: 1 addition & 1 deletion test/land/CMakeLists.txt
Original file line number Diff line number Diff line change
Expand Up @@ -32,6 +32,6 @@ set_tests_properties(test_gdasapp_land_letkfoi_snowda

# Test convert BUFR to IODA
add_test(NAME test_gdasapp_convert_bufr_adpsfc_snow
COMMAND ${PROJECT_SOURCE_DIR}/test/soca/test_bufr2ioda.sh ${PROJECT_BINARY_DIR} bufr_adpsfc_snow
COMMAND ${PROJECT_SOURCE_DIR}/test/land/test_bufr2ioda.sh ${PROJECT_BINARY_DIR} bufr_adpsfc_snow
WORKING_DIRECTORY ${PROJECT_BINARY_DIR}/test/)

4 changes: 2 additions & 2 deletions test/testinput/bufr_adpsfc_snow.yaml
Original file line number Diff line number Diff line change
Expand Up @@ -7,7 +7,7 @@ observations:
- obs space:
name: bufr

obsdatain: "./testdata/gdas.t06z.adpsfc.tm00.bufr_d"
obsdatain: $(DATA)/obs/$(OPREFIX)adpsfc.tm00.bufr_d

exports:
variables:
Expand Down Expand Up @@ -41,7 +41,7 @@ observations:

ioda:
backend: netcdf
obsdataout: "./testoutput/gdas.t06z.adpsfc_snow.nc"
obsdataout: $(DATA)/$(OPREFIX)adpsfc_snow.nc4

variables:

Expand Down
85 changes: 85 additions & 0 deletions ush/land/bufr_adpsfc_snow.yaml
Original file line number Diff line number Diff line change
@@ -0,0 +1,85 @@
# (C) Copyright 2021-2022 NOAA/NWS/NCEP/EMC
#
# This software is licensed under the terms of the Apache Licence Version 2.0
# which can be obtained at http://www.apache.org/licenses/LICENSE-2.0.

observations:
- obs space:
name: bufr

obsdatain: "./testdata/gdas.t06z.adpsfc.tm00.bufr_d"

exports:
variables:
# MetaData
timestamp:
datetime:
year: "*/YEAR[1]"
month: "*/MNTH[1]"
day: "*/DAYS[1]"
hour: "*/HOUR[1]"
minute: "*/MINU[1]"
latitude:
query: "[*/CLAT, */CLATH]"
longitude:
query: "[*/CLON, */CLONH]"
stationIdentification:
query: "*/RPID"

height:
query: "[*/SELV, */HSMSL]"

# ObsValue
totalSnowDepth:
query: "[*/SNWSQ1/TOSD, */MTRMSC/TOSD, */STGDSNDM/TOSD]"
transforms:
- scale: 1000.0
filters:
- bounding:
variable: totalSnowDepth
upperBound: 10000000

ioda:
backend: netcdf
obsdataout: "./testoutput/gdas.t06z.adpsfc_snow.nc"

variables:

# MetaData
- name: "MetaData/dateTime"
coordinates: "longitude latitude"
source: variables/timestamp
longName: "Datetime"
units: "seconds since 1970-01-01T00:00:00Z"

- name: "MetaData/latitude"
coordinates: "longitude latitude"
source: variables/latitude
longName: "Latitude"
units: "degree_north"
range: [-90, 90]

- name: "MetaData/longitude"
coordinates: "longitude latitude"
source: variables/longitude
longName: "Longitude"
units: "degree_east"
range: [-180, 180]

- name: "MetaData/height"
coordinates: "longitude latitude"
source: variables/height
longName: "Height of Station"

- name: "MetaData/stationIdentification"
coordinates: "longitude latitude"
source: variables/stationIdentification
longName: "Identification of Observing Location"
units: "m"

# ObsValue
- name: "ObsValue/totalSnowDepth"
coordinates: "longitude latitude"
source: variables/totalSnowDepth
longName: "Total Snow Depth"
units: "mm"

0 comments on commit 7630832

Please sign in to comment.