Skip to content

Commit

Permalink
Merge pull request #1411 from cmclaug2/feature/testcaseSuiteUpdate
Browse files Browse the repository at this point in the history
Feature/testcase suite update
  • Loading branch information
jvgeiger authored Jul 31, 2024
2 parents 511888d + 29f6cf9 commit 26c3f9e
Show file tree
Hide file tree
Showing 57 changed files with 5,751 additions and 1 deletion.
22 changes: 22 additions & 0 deletions ldt/testcases/DAobs_MODISlai/README
Original file line number Diff line number Diff line change
@@ -0,0 +1,22 @@
MODIS MCD15A2H v006 LAI Data Assimilation Test Case

This testcase processes MODIS MCD15A2H v006 leaf area index (LAI) to generate
CDFs for 5km domain.

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.)

To run this test case:

* Generate the LDT executable.
* Run the LDT executable using the ldt.config file and the sample
input data.

Notes:

This feature was added in pull request #659
(https://github.com/NASA-LIS/LISF/pull/659).
230 changes: 230 additions & 0 deletions ldt/testcases/DAobs_MODISlai/ldt.config
Original file line number Diff line number Diff line change
@@ -0,0 +1,230 @@
#Overall driver options
LDT running mode: "DA preprocessing"
Processed LSM parameter filename: ./lis_input.d01.nc
LIS number of nests: 1
Number of surface model types: 1
Surface model types: LSM
Land surface model: "Noah-MP.4.0.1"
Routing model: none
Lake model: none
Water fraction cutoff value: 0.5
#Incorporate crop information: .false.
Number of met forcing sources: 0
Met forcing sources: none
Met spatial transform methods: bilinear
Topographic correction method (met forcing): "lapse-rate"
LDT diagnostic file: ldtlog
Mask-parameter fill diagnostic file: noahmp_MODIS.log
LDT output directory: ./
Undefined value: -9999.0
Number of ensembles per tile: 1
Number of processors along x: 1
Number of processors along y: 1

#LIS domain
Map projection of the LIS domain: latlon
Run domain lower left lat: 28.025
Run domain lower left lon: 34.025
Run domain upper right lat: 34.975
Run domain upper right lon: 39.975
Run domain resolution (dx): 0.05
Run domain resolution (dy): 0.05

#Landcover parameter inputs
Landcover data source: "MODIS_Native"
Landcover classification: "IGBPNCEP"
Landcover file: ./input/LS_PARAMETERS/noah_2dparms/igbp.bin
Landcover spatial transform: tile
Landcover fill option: none
Landcover map projection: latlon

#Landmask parameter inputs
Create or readin landmask: "create"
Landmask data source: "MODIS_Native"
Landmask file: none
Landmask spatial transform: none
Landmask map projection: latlon

#Soil parameter inputs
Soil fraction data source: none
Soils spatial transform: none
Soils fill option: none
Soils map projection: latlon
Porosity data source: none
Porosity map: none

#Soil texture inputs
Soil texture data source: ISRIC
Soil texture map: ./input/LS_PARAMETERS/soil_parms/ISRIC/v2017/TEXMHT_M_sl1_250m.tif
Soil texture spatial transform: mode
Soil texture fill option: neighbor
Soil texture fill radius: 5
Soil texture fill value: 6
Soil texture map projection: latlon


Soils spatial transform: none
Soils map projection: latlon


#Albedo inputs
Albedo data source: "NCEP_LIS"
Albedo map: ./input/LS_PARAMETERS/UMD/5KM/albedo_NCEP
Albedo climatology interval: monthly
Albedo spatial transform: none
Albedo fill option: neighbor
Albedo fill radius: 5
Albedo fill value: 0.14
Albedo map projection: latlon
Albedo lower left lat: -59.975
Albedo lower left lon: -179.975
Albedo upper right lat: 89.975
Albedo upper right lon: 179.975
Albedo resolution (dx): 0.05
Albedo resolution (dy): 0.05

#Maximum snow albedo inputs
Max snow albedo data source: "Barlage_Native"
Max snow albedo map: ./input/LS_PARAMETERS/noah_2dparms/maximum_snow_albedo.hdf
Max snow albedo spatial transform: budget-bilinear
Max snow albedo fill option: neighbor
Max snow albedo fill radius: 5
Max snow albedo fill value: 0.3
Max snow albedo map projection: latlon

#Greenness inputs
Greenness data source: "NCEP_LIS"
Greenness fraction map: ./input/LS_PARAMETERS/UMD/5KM/gvf_NCEP
Greenness climatology interval: monthly
Calculate min-max greenness fraction: .true.
Greenness spatial transform: none
Greenness fill option: neighbor
Greenness fill radius: 5
Greenness fill value: 0.20
Greenness maximum fill value: 0.80
Greenness minimum fill value: 0.05
Greenness map projection: latlon
Greenness lower left lat: -59.975
Greenness lower left lon: -179.975
Greenness upper right lat: 89.975
Greenness upper right lon: 179.975
Greenness resolution (dx): 0.05
Greenness resolution (dy): 0.05

#Slope type inputs
Slope type data source: NCEP_Native
Slope type map: ./input/LS_PARAMETERS/noah_2dparms/islope
Slope type spatial transform: neighbor
Slope type fill option: neighbor
Slope type fill radius: 5
Slope type fill value: 3.0
Slope type map projection: latlon

#Noah-MP LSM inputs
Noah-MP PBL Height Value: 900.

#Elevation map
GDAS forcing directory: ./input/MET_FORCING/GDAS
GDAS T126 elevation map: ./input/LS_PARAMETERS/metforcing_parms/GDAS/global_orography.t126.grb
GDAS T170 elevation map: ./input/LS_PARAMETERS/metforcing_parms/GDAS/global_orography.t170.grb
GDAS T254 elevation map: ./input/LS_PARAMETERS/metforcing_parms/GDAS/global_orography.t254.grb
GDAS T382 elevation map: ./input/LS_PARAMETERS/metforcing_parms/GDAS/global_orography.t382.grb
GDAS T574 elevation map: ./input/LS_PARAMETERS/metforcing_parms/GDAS/global_orography.t574.grb
GDAS T1534 elevation map: ./input/LS_PARAMETERS/metforcing_parms/GDAS/global_orography_uf.t1534.3072.1536.grb

# SRTM Elevation data entries:
Elevation data source: "SRTM_LIS"
Elevation metadata variable name: "HGT_M"
Elevation map: ./input/LS_PARAMETERS/topo_parms/SRTM/SRTM30/srtm_elev1km.1gd4r
Elevation number of bands: 1
Slope data source: "SRTM_LIS"
Slope map: ./input/LS_PARAMETERS/topo_parms/SRTM/SRTM30/srtm_slope1km.1gd4r
Slope number of bands: 1
Aspect data source: "SRTM_LIS"
Aspect map: ./input/LS_PARAMETERS/topo_parms/SRTM/SRTM30/srtm_aspect1km.1gd4r
Aspect number of bands: 1


Topography spatial transform: neighbor
Elevation fill option: neighbor
Elevation fill radius: 10
Elevation fill value: 10
Slope fill option: neighbor
Slope fill radius: 10
Slope fill value: 0.1
Aspect fill option: neighbor
Aspect fill radius: 10
Aspect fill value: 0

Topography map projection: latlon
Topography lower left lat: -59.995
Topography lower left lon: -179.995
Topography upper right lat: 89.995
Topography upper right lon: 179.995
Topography resolution (dx): 0.01
Topography resolution (dy): 0.01


#Bottom temperature inputs
Bottom temperature data source: "NCEP_LIS"
Bottom temperature map: ./input/LS_PARAMETERS/UMD/5KM/tbot_NCEP.1gd4r
Bottom temperature spatial transform: none
Bottom temperature fill option: neighbor
Bottom temperature fill radius: 3.
Bottom temperature fill value: 287.0
Bottom temperature topographic downscaling: "none"
Bottom temperature map projection: latlon
Bottom temperature lower left lat: -59.975
Bottom temperature lower left lon: -179.975
Bottom temperature upper right lat: 89.975
Bottom temperature upper right lon: 179.975
Bottom temperature resolution (dx): 0.05
Bottom temperature resolution (dy): 0.05

DA preprocessing method: "CDF generation"
DA observation source: "MCD15A2H LAI"
Name of the preprocessed DA file: "MCD15A2H_cdf_5km_2obs"

MCD15A2H LAI data directory: ./input/LS_PARAMETERS/MODIS/MCD15A2H.006
MCD15A2H LAI data version: "006"
MCD15A2H LAI apply climatological fill values: 1
MCD15A2H LAI apply QC flags: 1

Apply anomaly correction to obs: 0
Temporal resolution of CDFs: monthly
Number of bins to use in the CDF: 10
Observation count threshold: 2
Temporal averaging interval: "1da"
Apply external mask: 0
External mask directory: none

Search radius for openwater proximity detection: 1


Starting year: 2003
Starting month: 1
Starting day: 1
Starting hour: 0
Starting minute: 0
Starting second: 0
Ending year: 2005
Ending month: 1
Ending day: 1
Ending hour: 0
Ending minute: 0
Ending second: 0
LIS output timestep: "1da"

Maximum number of surface type tiles per grid: 1
Minimum cutoff percentage (surface type tiles): 0.05
Maximum number of soil texture tiles per grid: 1
Minimum cutoff percentage (soil texture tiles): 0.05
Maximum number of soil fraction tiles per grid: 1
Minimum cutoff percentage (soil fraction tiles): 0.05
Maximum number of elevation bands per grid: 1
Minimum cutoff percentage (elevation bands): 0.05
Maximum number of slope bands per grid: 1
Minimum cutoff percentage (slope bands): 0.05
Maximum number of aspect bands per grid: 1
Minimum cutoff percentage (aspect bands): 0.05

Original file line number Diff line number Diff line change
@@ -0,0 +1,34 @@
Focus: Downscaling 12 Member Ensemble Restart File to 8 Member Restart File

Description: This case involves the downscaling or generating an eight-member
restart file from a twelve-member restart file using the "random sampling"
sampling strategy.


This test case uses the following files:
(a) this README file
(b) the ldt.config input file used by LDT

To run such a case:

(a) Obtain an ensemble-based restart file from one of your LIS
model run directories, such as:
LIS_RST_NOAH33_201212312330.d01.nc

(b) Generate the LDT executable and copy to your working directory.

(c) Modify the ldt.config file provided for your case.

(d) Run the LDT executable with the provided ldt.config file by typing:
./LDT ldt.config

(e) View the NetCDF (*.nc) output using a tool, like
the ncview executable (if you have it installed).

Caveats:
(a) Please note that this is a more basic functional test for purpose
of learning major LDT features.

Notes:
This feature was added in pull request #853:
https://github.com/NASA-LIS/LISF/pull/853
Loading

0 comments on commit 26c3f9e

Please sign in to comment.