Skip to content

Commit

Permalink
Merge pull request #1464 from jvgeiger/feature/distbinaryrst_testcase
Browse files Browse the repository at this point in the history 
Add testcase files for distributed binary restart testcase
  • Loading branch information
@jvgeiger
jvgeiger committed Dec 6, 2023
2 parents ef9eb8e + 879d742 commit d648d72
Show file tree
Hide file tree
Showing 8 changed files with 11,271 additions and 0 deletions.
153 changes: 153 additions & 0 deletions lis/testcases/misc/dist_binary_rst/NOAHMP_OUTPUT_LIST.TBL
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,153 @@
#short_name select? units signconv timeavg? min/max? std? vert.levels grib_id grib_scalefactor longname

#Energy balance components
Swnet: 1 W/m2 DN 1 0 0 1 111 10 # Net Shortwave Radiation (W/m2)
Lwnet: 1 W/m2 DN 1 0 0 1 112 10 # Net Longwave Radiation (W/m2)
Qle: 1 W/m2 UP 1 0 0 1 121 10 # Latent Heat Flux (W/m2)
Qh: 0 W/m2 UP 1 0 0 1 122 10 # Sensible Heat Flux (W/m2)
Qg: 0 W/m2 DN 1 0 0 1 155 10 # Ground Heat Flux (W/m2)
Qf: 0 W/m2 S2L 1 0 0 1 229 10 # Energy of fusion (W/m2)
Qv: 0 W/m2 S2V 1 0 0 1 134 10 # Energy of sublimation (W/m2)
Qa: 0 W/m2 DN 1 0 0 1 136 10 # Advective Energy (W/m2)
Qtau: 0 N/m2 DN 1 0 0 1 135 10 # Momentum flux (N/m2)
DelSurfHeat: 0 J/m2 INC 1 0 0 1 137 10 # Change in surface heat storage (J/m2)
DelColdCont: 0 J/m2 INC 1 0 0 1 138 10 # Change in snow cold content (J/m2)
BR: 0 - - 1 0 1 1 256 10 # Bowen ratio
EF: 0 - - 1 0 1 1 256 10 # Evaporative fraction

#Water balance components
Snowf: 1 kg/m2s DN 1 0 0 1 161 10000 # Snowfall rate (kg/m2s)
Rainf: 0 kg/m2s DN 1 0 0 1 162 10000 # Rainfall rate (kg/m2s)
RainfConv: 0 kg/m2s DN 1 0 0 1 163 10000 # Convective Rainfall rate (kg/m2s)
TotalPrecip: 1 kg/m2s DN 1 0 0 1 164 10000 # Total Precipitation rate (kg/m2s)
Evap: 0 kg/m2s UP 1 0 0 1 57 10000 # Total Evapotranspiration (kg/m2s)
Qs: 1 kg/m2s OUT 1 0 0 1 235 10000 # Surface runoff (kg/m2s)
Qrec: 0 kg/m2s IN 1 0 0 1 143 10000 # Recharge (kg/m2s)
Qsb: 1 kg/m2s OUT 1 0 0 1 254 10000 # Subsurface runoff (kg/m2s)
Qsm: 1 kg/m2s S2L 1 0 0 1 99 10000 # Snowmelt (kg/m2s)
Qfz: 0 kg/m2s L2S 1 0 0 1 146 10000 # Refreezing of water in the snowpack (kg/m2s)
Qst: 0 kg/m2s - 1 0 0 1 147 10000 # Snow throughfall (kg/m2s)
DelSoilMoist: 0 kg/m2 INC 0 0 0 1 148 10000 # Change in soil moisture (kg/m2)
DelSWE: 0 kg/m2 INC 0 0 0 1 149 1000 # Change in snow water equivalent (kg/m2)
DelSurfStor: 0 kg/m2 INC 1 0 0 1 150 1000 # Change in surface water storage (kg/m2)
DelIntercept: 0 kg/m2 INC 1 0 0 1 151 1000 # Change in interception storage (kg/m2)
RHMin: 0 - - 0 0 0 1 51 10 # Minimum 2 meter relative humidity (-)

#Surface state variables
SnowT: 0 K - 1 0 0 1 152 10 # Snow surface temperature (K)
VegT: 1 K - 1 0 0 1 153 10 # Vegetation canopy temperature (K)
BareSoilT: 0 K - 1 0 0 1 154 10 # Temperature of bare soil (K)
AvgSurfT: 1 K - 1 0 0 1 148 10 # Average surface temperature (K)
RadT: 1 K - 1 0 0 1 156 10 # Surface Radiative Temperature (K)
Albedo: 0 - - 1 0 0 1 84 100 # Surface Albedo (-)
DTB_10H_36H: 0 K - 0 0 0 1 65 1000 # TB difference between 10H and 36H (K)
DTB_10V_36V: 0 K - 0 0 0 1 65 1000 # TB difference between 10V and 36V (K)
DTB_18H_36H: 0 K - 0 0 0 1 65 1000 # TB difference between 18H and 36H (K)
DTB_18V_36V: 0 K - 0 0 0 1 65 1000 # TB difference between 18V and 36V (K)
TB_10H: 0 K - 0 0 0 1 65 1000 # TB at 10H (K)
TB_10V: 0 K - 0 0 0 1 65 1000 # TB at 10V (K)
TB_18H: 0 K - 0 0 0 1 65 1000 # TB at 18H (K)
TB_18V: 0 K - 0 0 0 1 65 1000 # TB at 18V (K)
TB_36H: 0 K - 0 0 0 1 65 1000 # TB at 36H (K)
TB_36V: 0 K - 0 0 0 1 65 1000 # TB at 36V (K)
SWE: 1 kg/m2 - 1 0 0 1 65 1000 # Snow Water Equivalent (kg/m2)
SWEVeg: 0 kg/m2 - 1 0 0 1 159 1000 # SWE intercepted by vegetation (kg/m2)
SurfStor: 0 kg/m2 - 1 0 0 1 160 1000 # Surface water storage (kg/m2)
TWS: 0 mm - 1 0 0 1 160 1000 #

#Subsurface state variables
SoilMoist: 1 m3/m3 - 1 0 0 4 86 1000 # Average layer soil moisture (kg/m2)
SoilTemp: 1 K - 1 0 0 4 85 1000 # Average layer soil temperature (K)
SmLiqFrac: 0 - - 1 0 0 4 85 100 # Average layer fraction of liquid moisture (-)
SmFrozFrac: 0 - - 1 0 0 4 85 100 # Average layer fraction of frozen moisture (-)
SoilWet: 0 - - 0 0 0 1 85 100 # Total soil wetness (-)
RelSMC: 0 m3/m3 - 0 0 0 1 86 1000 # Relative soil moisture
RootTemp: 0 K - 0 0 0 1 85 1000 # Rootzone temperature (K)

#Evaporation components
PotEvap: 0 kg/m2s UP 1 0 0 1 166 1 # Potential Evapotranspiration (kg/m2s)
ECanop: 0 kg/m2s UP 1 0 0 1 200 1 # Interception evaporation (kg/m2s)
TVeg: 0 kg/m2s UP 1 0 0 1 210 1 # Vegetation transpiration (kg/m2s)
ESoil: 0 kg/m2s UP 1 0 0 1 199 1 # Bare soil evaporation (kg/m2s)
EWater: 0 kg/m2s UP 1 0 0 1 170 1 # Open water evaporation (kg/m2s)
RootMoist: 0 kg/m2 - 0 0 0 1 171 1 # Root zone soil moisture (kg/m2)
CanopInt: 0 kg/m2 - 1 0 0 1 223 1000 # Total canopy water storage (kg/m2)
EvapSnow: 0 kg/m2s - 1 0 0 1 173 1000 # Snow evaporation (kg/m2s)
SubSnow: 0 kg/m2s - 1 0 0 1 198 1000 # Snow sublimation (kg/m2s)
SubSurf: 0 kg/m2s - 1 0 0 1 175 1000 # Sublimation of the snow free area (kg/m2s)
ACond: 0 m/s - 1 0 0 1 179 100000 # Aerodynamic conductance
CCond: 0 m/s - 1 0 0 1 179 1000000 # Canopy conductance
VegCanopT: 0 K - 1 0 0 1 179 1000000 # Canopy air temperature

#Cold season processes
Snowcover: 0 - - 1 0 0 1 66 100 # Snow Cover (-)
SnowDepth: 1 m - 1 0 0 1 66 1000 # Snow Depth (m)
SLiqFrac: 0 - - 0 0 0 1 65 1000 # Fraction of SWE in the liquid phase
SnowIce: 0 mm - 0 0 0 3 65 1000 # Snow-layer ice (mm)
SnowAge: 0 - - 1 0 0 1 65 1000 # Snow again term (-)
SAlbedo: 0 - - 1 0 0 1 65 1000 # Snow albedo at last time step (-)
SnowTProf: 0 K - 0 0 0 3 65 1000 # Snow layer temperature (K)
SnowLiq: 1 mm - 0 0 0 3 65 1000 # Snow layer liquid water (mm)
ActSnowNL: 0 - - 1 0 0 1 65 1000 # actual number of snow layers (-)
z_snow: 0 m - 0 0 0 3 65 1000 # Snow layer-bottom depth from snow surface (m)
z_soil: 0 m - 0 0 0 4 65 1000 # Soil layer-bottom depth from snow surface (m)

#Forcings
Wind_f: 0 m/s - 1 0 0 1 177 10 # Near Surface Wind (m/s)
Rainf_f: 1 kg/m2s DN 1 0 0 1 162 1000 # Average rainfall rate
Snowf_f: 0 kg/m2s DN 1 0 0 1 161 1000 # Average snowfall rate
Tair_f: 1 K - 1 0 0 1 11 10 # Near surface air temperature
Qair_f: 0 kg/kg - 1 0 0 1 51 1000 # Near surface specific humidity
Psurf_f: 0 Pa - 1 0 0 1 1 10 # Surface pressure
SWdown_f: 0 W/m2 DN 1 0 0 1 204 10 # Surface incident shortwave radiation
LWdown_f: 0 W/m2 DN 1 0 0 1 205 10 # Surface incident longwave radiation

#Additional forcings
DirectSW_f: 0 W/m2 - 0 0 0 1 256 10 # Surface direct incident shortwave radiation
DiffuseSW_f: 0 W/m2 - 0 0 0 1 256 10 # Surface diffuse incident shortwave radiation
NWind_f: 0 m/s N 0 0 0 1 256 10 # Northward wind
EWind_f: 0 m/s E 0 0 0 1 256 10 # Eastward wind
FHeight_f: 0 m - 0 0 0 1 256 10 # Height of forcing variables
CH_f: 0 - - 0 0 0 1 256 10 # Surface exchange coefficient for heat
CM_f: 0 - - 0 0 0 1 256 10 # Surface Exchange Coefficient for momentum
Emiss_f: 0 - - 0 0 0 1 256 10 # Surface emissivity
MixRatio_f: 0 kg/kg - 0 0 0 1 256 10 # Surface mixing ration
CosZenith_f: 0 - - 0 0 0 1 256 10 # Cosine of zenith angle
Albedo_f: 0 - - 0 0 0 1 256 10 # Surface albedo

#Parameters
Landmask: 0 - - 0 0 0 1 81 1 # Land Mask (0 - Water, 1- Land)
Landcover: 0 - - 0 0 0 1 186 1 # Land cover
Soiltype: 0 - - 0 0 0 1 187 1 # Soil type
SandFrac: 0 - - 0 0 0 1 999 1 # Sand fraction
ClayFrac: 0 - - 0 0 0 1 999 1 # Clay fraction
SiltFrac: 0 - - 0 0 0 1 999 1 # Silt fraction
Porosity: 0 - - 3 0 0 1 999 1 # Porosity
Soilcolor: 0 - - 0 0 0 1 188 1 # Soil color
Elevation: 0 m - 0 0 0 1 189 10 # Elevation
Slope: 0 - - 0 0 0 1 999 10 # Slope
LAI: 1 - - 1 0 0 1 190 100 # LAI
SAI: 0 - - 0 0 0 1 191 100 # SAI
Snfralbedo: 0 - - 0 0 0 1 192 100 # Snow fraction albedo
Mxsnalbedo: 0 - - 0 0 0 1 192 100 # Maximum snow albedo
Greenness: 0 - - 0 0 0 1 87 100 # Greenness
Tempbot: 0 - - 0 0 0 1 194 10 # Bottom soil temperature

#Routing
Streamflow: 0 m3/s - 1 0 0 1 333 10 #Streamflow
RiverStor: 0 m3 - 1 0 0 1 333 10 #RiverStorage
RiverDepth: 0 m - 1 0 0 1 333 10 #RiverDepth
RiverVelocity: 0 m/s - 1 0 0 1 333 10 #RiverVelocity
FloodQ: 0 m3/s - 1 0 0 1 333 10 #FloodDischarge
FloodEvap: 0 m3 - 1 0 0 1 333 10 #FloodEvap
FloodStor: 0 m3 - 1 0 0 1 333 10 #FloodStorage
FloodDepth: 0 m - 1 0 0 1 333 10 #FloodDepth
FloodVelocity: 0 m/s - 1 0 0 1 333 10 #FloodVelocity
FloodedFrac: 0 - - 1 0 0 1 333 10 #FloodedFrac
FloodedArea: 0 m2 - 1 0 0 1 333 10 #FloodedArea
SurfElev: 0 m - 1 0 0 1 333 10 #SurfElev
RunoffStor: 0 m3 - 1 0 0 1 333 10 #RunoffStor
BaseflowStor: 0 m3 - 1 0 0 1 333 10 #BaseflowStor

#Irrigation
Irrigated water: 0 kg/m2s - 1 0 0 1 333 10 #Irrigation amount
34 changes: 34 additions & 0 deletions lis/testcases/misc/dist_binary_rst/README
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,34 @@
Distributed Binary Restart Testcase

This test demonstrates writing distributed binary for both LIS_RST and LIS_DAPERT files with a assimilation example over the Tuolumne watershed.

This is a testcase that uses:
(a) The MERRA2 forcing
(b) Bilinear Interpolation
(c) The NoahMP 4.0.1 land surface model
(d) Tuolumne watershed
(e) A time period from 0z 29 Oct 2015 to 0z 3 Nov 2015

This directory contains:
(a) This README file
(b) The "lis.config" file used for this testcase
(c) The "MODEL_OUTPUT_LIST.TBL" file used by the lis.config file
to select the output variables
(d) An "output.ctl" file used to display the output in GrADS from
the LIS run

Note that the "lis.config" file should be edited
to make sure that the locations of the parameter and forcing files
are specified correctly.

To run this testcase:
(a) Generate the LIS executable
(b) Run the LIS executable using the "lis.config" file with the
testcase input data
(c) View the binary output using the testcase GrADS files

Caveats:
(a) Please note that this is a simple functional test and the output
from the testcase is not expected to be used for any scientific
evaluation.

Loading

0 comments on commit d648d72

Please sign in to comment.