Skip to content

Commit

Permalink
Merge pull request #1678 from jvgeiger/feature/espc-d_testcase
Browse files Browse the repository at this point in the history
Add ESPC-D testcase files
  • Loading branch information
emkemp authored Jan 31, 2025
2 parents 8ea7b40 + 831c431 commit 8081e94
Show file tree
Hide file tree
Showing 5 changed files with 276 additions and 0 deletions.
19 changes: 19 additions & 0 deletions ldt/testcases/ESPC-D/README
Original file line number Diff line number Diff line change
@@ -0,0 +1,19 @@
ESPC-D Reader Test Case

This testcase processes Navy Earth System Prediction Capability (ESPC-D) data.

This directory contains:

* this README file.
* the ldt.config file used for this test case. (This file should be
edited to make sure that the locations of the parameter and
forcing files are specified correctly.)
* the input.ctl GrADS descriptor file for viewing the input parameter file.
* the output.ctl GrADS descriptor file for viewing the OUTPUT.
* the testcase.ctl GrADS descriptor file for viewing the TARGET_OUTPUT.

To run this test case:

* Generate the LDT executable.
* Run the LDT executable using the ldt.config file and the sample
input data.
107 changes: 107 additions & 0 deletions ldt/testcases/ESPC-D/input.ctl
Original file line number Diff line number Diff line change
@@ -0,0 +1,107 @@
dset ^lis_input.global.noah39.nc
dtype netcdf
options template
undef -9999
xdef 2560 linear -179.9297 0.140625
ydef 1920 linear -89.95312 0.09375
zdef 1 linear 1 1
* dummy tdef
tdef 1 linear 00z01dec2024 1hr
vars 96
DOMAINMASK=>DOMAINMASK 1 y,x description
LANDMASK=>LANDMASK 1 y,x description
SURFACETYPE=>SURFACETYPE1 0 0,y,x description
SURFACETYPE=>SURFACETYPE2 0 1,y,x description
SURFACETYPE=>SURFACETYPE3 0 2,y,x description
SURFACETYPE=>SURFACETYPE4 0 3,y,x description
SURFACETYPE=>SURFACETYPE5 0 4,y,x description
SURFACETYPE=>SURFACETYPE6 0 5,y,x description
SURFACETYPE=>SURFACETYPE7 0 6,y,x description
SURFACETYPE=>SURFACETYPE8 0 7,y,x description
SURFACETYPE=>SURFACETYPE9 0 8,y,x description
SURFACETYPE=>SURFACETYPE10 0 9,y,x description
SURFACETYPE=>SURFACETYPE11 0 10,y,x description
SURFACETYPE=>SURFACETYPE12 0 11,y,x description
SURFACETYPE=>SURFACETYPE13 0 12,y,x description
SURFACETYPE=>SURFACETYPE14 0 13,y,x description
SURFACETYPE=>SURFACETYPE15 0 14,y,x description
SURFACETYPE=>SURFACETYPE16 0 15,y,x description
SURFACETYPE=>SURFACETYPE17 0 16,y,x description
SURFACETYPE=>SURFACETYPE18 0 17,y,x description
SURFACETYPE=>SURFACETYPE19 0 18,y,x description
SURFACETYPE=>SURFACETYPE20 0 19,y,x description
SURFACETYPE=>SURFACETYPE21 0 20,y,x description
LANDCOVER=>LANDCOVER1 0 0,y,x description
LANDCOVER=>LANDCOVER2 0 1,y,x description
LANDCOVER=>LANDCOVER3 0 2,y,x description
LANDCOVER=>LANDCOVER4 0 3,y,x description
LANDCOVER=>LANDCOVER5 0 4,y,x description
LANDCOVER=>LANDCOVER6 0 5,y,x description
LANDCOVER=>LANDCOVER7 0 6,y,x description
LANDCOVER=>LANDCOVER8 0 7,y,x description
LANDCOVER=>LANDCOVER9 0 8,y,x description
LANDCOVER=>LANDCOVER10 0 9,y,x description
LANDCOVER=>LANDCOVER11 0 10,y,x description
LANDCOVER=>LANDCOVER12 0 11,y,x description
LANDCOVER=>LANDCOVER13 0 12,y,x description
LANDCOVER=>LANDCOVER14 0 13,y,x description
LANDCOVER=>LANDCOVER15 0 14,y,x description
LANDCOVER=>LANDCOVER16 0 15,y,x description
LANDCOVER=>LANDCOVER17 0 16,y,x description
LANDCOVER=>LANDCOVER18 0 17,y,x description
LANDCOVER=>LANDCOVER19 0 18,y,x description
LANDCOVER=>LANDCOVER20 0 19,y,x description
LANDCOVER=>LANDCOVER21 0 20,y,x description
TEXTURE=>TEXTURE1 0 0,y,x description
TEXTURE=>TEXTURE2 0 1,y,x description
TEXTURE=>TEXTURE3 0 2,y,x description
TEXTURE=>TEXTURE4 0 3,y,x description
TEXTURE=>TEXTURE5 0 4,y,x description
TEXTURE=>TEXTURE6 0 5,y,x description
TEXTURE=>TEXTURE7 0 6,y,x description
TEXTURE=>TEXTURE8 0 7,y,x description
TEXTURE=>TEXTURE9 0 8,y,x description
TEXTURE=>TEXTURE10 0 9,y,x description
TEXTURE=>TEXTURE11 0 10,y,x description
TEXTURE=>TEXTURE12 0 11,y,x description
TEXTURE=>TEXTURE13 0 12,y,x description
TEXTURE=>TEXTURE14 0 13,y,x description
TEXTURE=>TEXTURE15 0 14,y,x description
TEXTURE=>TEXTURE16 0 15,y,x description
ELEVFGRD=>ELEVFGRD 1 y,x description
ELEVATION=>ELEVATION 1 y,x description
SLOPEFGRD=>SLOPEFGRD 1 y,x description
SLOPE=>SLOPE 1 y,x description
ASPECTFGRD=>ASPECTFGRD 1 y,x description
ASPECT=>ASPECT 1 y,x description
GREENNESS=>GREENNESS1 0 0,y,x description
GREENNESS=>GREENNESS2 0 1,y,x description
GREENNESS=>GREENNESS3 0 2,y,x description
GREENNESS=>GREENNESS4 0 3,y,x description
GREENNESS=>GREENNESS5 0 4,y,x description
GREENNESS=>GREENNESS6 0 5,y,x description
GREENNESS=>GREENNESS7 0 6,y,x description
GREENNESS=>GREENNESS8 0 7,y,x description
GREENNESS=>GREENNESS9 0 8,y,x description
GREENNESS=>GREENNESS10 0 9,y,x description
GREENNESS=>GREENNESS11 0 10,y,x description
GREENNESS=>GREENNESS12 0 11,y,x description
SHDMIN=>SHDMIN 1 y,x description
SHDMAX=>SHDMAX 1 y,x description
ALBEDO=>ALBEDO1 0 0,y,x description
ALBEDO=>ALBEDO2 0 1,y,x description
ALBEDO=>ALBEDO3 0 2,y,x description
ALBEDO=>ALBEDO4 0 3,y,x description
ALBEDO=>ALBEDO5 0 4,y,x description
ALBEDO=>ALBEDO6 0 5,y,x description
ALBEDO=>ALBEDO7 0 6,y,x description
ALBEDO=>ALBEDO8 0 7,y,x description
ALBEDO=>ALBEDO9 0 8,y,x description
ALBEDO=>ALBEDO10 0 9,y,x description
ALBEDO=>ALBEDO11 0 10,y,x description
ALBEDO=>ALBEDO12 0 11,y,x description
MXSNALBEDO=>MXSNALBEDO 1 y,x description
TBOT=>TBOT 1 y,x description
lat=>lat 1 y,x description
lon=>lon 1 y,x description
endvars
118 changes: 118 additions & 0 deletions ldt/testcases/ESPC-D/ldt.config
Original file line number Diff line number Diff line change
@@ -0,0 +1,118 @@

# Overall driver options
LDT running mode: "USAFSI analysis"
Processed LSM parameter filename: lis_input.global.noah39.nc # When using LIS-Noah
LIS number of nests: 1
Number of surface model types: 0
Number of met forcing sources: 0
LDT diagnostic file: ldtlog.2024120100
LDT output directory: LDTOUT
Undefined value: -9999.0

# Rotated GALWEM 10-km domain
Map projection of the LIS domain: latlon
Run domain lower left lat: -89.9531250
Run domain lower left lon: -179.9296875
Run domain upper right lat: 89.9531250
Run domain upper right lon: 179.9296875
Run domain resolution (dx): 0.140625
Run domain resolution (dy): 0.093750

# ------------------------
# Section for USAFSI

# *** Former Environment Variables ***
USAFSI netcdf filename prefix: USAFSI
USAFSI valid date (YYYYMMDDHH): 2024120100
USAFSI fractional snow data directory: ./USAFSIIN/snofrac/
USAFSI modified data directory: ./USAFSIIN/modified/ # not used
USAFSI surface obs data format: 1
USAFSI surface obs data directory: ./USAFSIIN/sfcobs/
# note: options for 2: XCAL GMI and 3: AMSR2 are still in development or test state
USAFSI brightness temperature data option: 1 # 1: SSMIS, 2: XCAL GMI, 3: AMSR2
USAFSI XCAL GMI data directory: ./GMI_retrievals/
USAFSI AMSR2 data directory: ./AMSR2_retrievals/
USAFSI SSMIS data directory: ./SSMIS_retrievals/
USAFSI FNMOC SST GRIB1 data directory: ./USAFSIIN/FNMOC_SST/
USAFSI surface temperature data directory: ./USAFSIIN/sfctmps/
USAFSI static data directory: ./USAFSIIN/static/
USAFSI unmodified data directory: ./USAFSIIN/unmodified/ # not used
USAFSI VIIRS data directory: ./USAFSIIN/viirs

# *** PMW snow depth retrieval algorithm
USAFSI XCAL GMI raw data directory: ./GMI_L1CR_web/
USAFSI AMSR2 raw data directory: ./AMSR2_L1R/
USAFSI SSMIS raw data directory: ./USAFSIIN/OBSI_SSMISUPP/
# note: 4: Kelly, 2009 algorithm only works for XCAL GMI and AMSR2 datasets; under development
USAFSI PMW snow depth retrieval algorithm option: 3 # 1: Hollinger,1991, 2: Chang et al.,1987
# 3:Foster et al.,1997, 4: Kelly, 2009
USAFSI forest fraction file: ./USAFSIIN/static/ForestFraction_0p25deg.nc # for algorithm 3 and 4
USAFSI forest density file: ./MOD44B/PTC_MOD44B_0p25.nc # for algorithm 4

# *** Snow climatology
USAFSI Snow Climatology: 1 # 1: legacy 2: updated 10-km climo

# *** Former Namelist Variables ***
USAFSI decimal fraction adjustment of snow depth towards climo: 0.1
USAFSI default snow depth (m) when actual depth unknown: 0.1
USAFSI minimum snow depth (m) for which to print a diagnostic: 0.025
USAFSI maximum number of surface observations allowed: 200000
USAFSI SSMIS shallow snow depth threshold (m): 0.15
# There should be three latitudes listed below, defining three latitude bands
# for sanity checking snow reports in different seasons.
# Strongly recommend 4000 3000 2000
USAFSI latitudes (deg * 100) for summer climo check: 4000 3000 2000
# There should be four elevations listed below, defining four elevation
# levels for sanity checking snow reports in different seasons.
# Strongly recommend 1000 1500 1000 1000
USAFSI elevations (m) for summer climo check: 1000 1500 1000 1000
USAFSI temperature (deg K * 10) above which no snow is allowed: 2820
USAFSI max reported temperature (deg K * 10) allowed around poles: 3030
USAFSI minimum ice concentration (%) needed to set ice flag: 15
# There should 24 latitudes defined below, each row definining a hemisphere
# (NH first row, SH second row), and each column defining a month
# (J F M A M J J A S O N D). These are used to define "high latitudes"
# in the SSMIS-based sea ice analysis.
USAFSI high latitude thresholds (deg) for sea ice::
81.0 81.0 81.0 81.0 81.0 81.0 81.0 82.0 82.0 82.0 81.0 81.0
90.0 90.0 78.0 78.0 68.0 68.0 67.0 67.0 66.0 66.0 68.0 70.0
::
# There should 24 latitudes defined below, each row definining a hemisphere
# (NH first row, SH second row), and each column defining a month
# (J F M A M J J A S O N D). These are used to define "low latitudes"
# in the SSMIS-based sea ice analysis.
USAFSI low latitude thresholds (deg) for sea ice::
45.0 44.0 44.0 45.0 51.0 52.0 55.0 65.0 65.0 62.0 53.0 50.0
60.0 62.0 62.0 60.0 58.0 56.0 55.0 55.0 53.0 54.0 55.0 57.0
::
USAFSI max age of VIIRS pixels to use: 3
USAFSI min VIIRS fraction to mark point as bare ground: 0.6
USAFSI min VIIRS/CDFS-II fraction to mark point as snow: 0.4
USAFSI use CDFS-II fractional snow data: .true.
USAFSI use VIIRS snow mask: .true.
# Bratseth settings
USAFSI observation error variance (m^2): 0.0006
USAFSI background error variance (m^2): 0.0010
USAFSI background error horizontal correlation length (m): 55555.
USAFSI background error vertical correlation length (m): 800.
USAFSI elevQC difference threshold (m): 400.
USAFSI skewed backQC snow depth threshold (m): 0.4

# Other new settings
USAFSI bogus climatology snow depth value (m): 0.2
USAFSI source of ocean data: "ESPC-D"
USAFSI GOFS SST data directory: ./USAFSIIN/GOFS
USAFSI GOFS CICE data directory: ./USAFSIIN/GOFS
USAFSI ESPC-D SST data directory: ./USAFSIIN/ESPC-D
USAFSI ESPC-D CICE data directory: ./USAFSIIN/ESPC-D
USAFSI LIS GRIB2 data directory: ./USAFSIIN/LIS_T2
USAFSI LIS GRIB2 security class: U
USAFSI LIS GRIB2 data category: C
USAFSI LIS GRIB2 data resolution: C0P09DEG
USAFSI LIS GRIB2 area of data: GLOBAL

USAFSI GALWEM root directory: ./USAFSIIN/GALWEM
USAFSI GALWEM subdirectory: GALWEM
USAFSI GALWEM use timestamp directories: 1
USAFSI GALWEM nominal resolution (km): 17

16 changes: 16 additions & 0 deletions ldt/testcases/ESPC-D/output.ctl
Original file line number Diff line number Diff line change
@@ -0,0 +1,16 @@
dset ^USAFSI_2024120100.nc
dtype netcdf
options template
undef -1
xdef 2560 linear -179.9297 0.140625
ydef 1920 linear -89.95312 0.09375
zdef 1 linear 1 1
* dummy tdef
tdef 1 linear 00z01dec2024 1hr
vars 5
snoanl=>snoanl1 0 0,y,x description
snoage=>snoage1 0 0,y,x description
icecon=>icecon1 0 0,y,x description
icemask=>icemask1 0 0,y,x description
iceage=>iceage1 0 0,y,x description
endvars
16 changes: 16 additions & 0 deletions ldt/testcases/ESPC-D/testcase.ctl
Original file line number Diff line number Diff line change
@@ -0,0 +1,16 @@
dset ^TARGET_OUTPUT/USAFSI_2024120100.nc
dtype netcdf
options template
undef -1
xdef 2560 linear -179.9297 0.140625
ydef 1920 linear -89.95312 0.09375
zdef 1 linear 1 1
* dummy tdef
tdef 1 linear 00z01dec2024 1hr
vars 5
snoanl=>snoanl1 0 0,y,x description
snoage=>snoage1 0 0,y,x description
icecon=>icecon1 0 0,y,x description
icemask=>icemask1 0 0,y,x description
iceage=>iceage1 0 0,y,x description
endvars

0 comments on commit 8081e94

Please sign in to comment.