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

LIS/Noah-MPv5.0 coupling #1666

Open
wants to merge 32 commits into
base: master
Choose a base branch
from
Open
Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
32 commits
Select commit Hold shift + click to select a range
f5b5a1d
coupling with NoahMP refactor 1
cenlinhe Jul 13, 2023
d4936fb
Merge branch 'master' of https://github.com/cenlinhe/LISF
cenlinhe Jul 13, 2023
f89d955
update lis parameter setup
cenlinhe Jul 14, 2023
6d9c3ea
update lis-noahmp driver 3
cenlinhe Jul 14, 2023
5692045
update integration 3
cenlinhe Jul 14, 2023
274706e
update integration for additional vars
cenlinhe Jul 14, 2023
ce277bd
update initialization
cenlinhe Jul 20, 2023
52012ea
update DA and sublsm
cenlinhe Jul 25, 2023
686e8d2
formating change
cenlinhe Jul 25, 2023
3f66017
bug fix for compilation error
cenlinhe Jul 26, 2023
9a10d75
bug fix
cenlinhe Jul 27, 2023
48697eb
add lis.config example for NoahMP v5
cenlinhe Jul 27, 2023
bab3bfb
bug fix
cenlinhe Jul 28, 2023
c79e5c9
sync with latest NoahMP v5 code
cenlinhe Nov 15, 2023
1f4801f
rename NoahMPnew to NoahMP50
cenlinhe Nov 17, 2023
3e56da4
Merge branch 'NASA-LIS:master' into LISF_NoahMPv5
cenlinhe Nov 17, 2023
e4950e6
change NoahMP.New in info message
cenlinhe Nov 18, 2023
8edb892
Fix Noah-MPv5 restart and model alarms for multiple nests following v…
cenlinhe Nov 30, 2023
ffd37b0
Merge branch 'NASA-LIS:master' into LISF_NoahMPv5
cenlinhe Nov 30, 2023
bb3694b
update DA TWS following NoahMPv4.0.1
cenlinhe Nov 30, 2023
7ddc550
remove redundancy in da_tws for snow depth variables
cenlinhe Nov 30, 2023
9e1639f
sync with latest noahmp github & bug fix
cenlinhe Nov 30, 2023
ab7d602
update lis.config for noahmp v5
cenlinhe Nov 30, 2023
7002100
Merge branch 'NASA-LIS:master' into LISF_NoahMPv5
cenlinhe Dec 6, 2023
2e6475d
Merge branch 'NASA-LIS:master' into LISF_NoahMPv5
cenlinhe Jan 5, 2024
1f82d43
bug fix for noahmp albedo & leafmass init & remove redundancy
cenlinhe Jan 5, 2024
ac29779
update lis.config for noahmp5
cenlinhe Jan 5, 2024
cb80970
update lis.config noahmp5 for AGRMET namelist updates
cenlinhe Jan 6, 2024
d9e9f95
bug fix for domain size dx dy input
cenlinhe Apr 23, 2024
750da6e
sync with latest noahmp bug fix
cenlinhe Apr 23, 2024
ce1405a
improve Noah-MPv5.0 naming convention and documentation
cenlinhe Dec 28, 2024
4dffcc7
update noahmp submodule link to official NCAR github
cenlinhe Dec 28, 2024
File filter

Filter by extension

Filter by extension

Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
4 changes: 4 additions & 0 deletions .gitmodules
Original file line number Diff line number Diff line change
@@ -0,0 +1,4 @@
[submodule "lis/surfacemodels/land/noahmp.5.0/phys/noahmp"]
path = lis/surfacemodels/land/noahmp.5.0/phys/noahmp
url = https://github.com/NCAR/noahmp
branch = release-v5.0-LIS
145 changes: 145 additions & 0 deletions lis/configs/Noah-MPv5.0/MODEL_OUTPUT_LIST_noahmp5.TBL
Original file line number Diff line number Diff line change
@@ -0,0 +1,145 @@
#short_name select? units signconv timeavg? min/max? std? vert.levels grib_id grib_scalefactor longname
#Energy balance components
Swnet: 1 W/m2 DN 0 0 0 1 111 10 # Net Shortwave Radiation (W/m2)
Lwnet: 1 W/m2 DN 0 0 0 1 112 10 # Net Longwave Radiation (W/m2)
Qle: 1 W/m2 UP 0 0 0 1 121 10 # Latent Heat Flux (W/m2)
Qh: 1 W/m2 UP 0 0 0 1 122 10 # Sensible Heat Flux (W/m2)
Qg: 1 W/m2 DN 0 0 0 1 155 10 # Ground Heat Flux (W/m2)
Qf: 0 W/m2 S2L 0 0 0 1 229 10 # Energy of fusion (W/m2)
Qv: 0 W/m2 S2V 0 0 0 1 134 10 # Energy of sublimation (W/m2)
Qa: 0 W/m2 DN 0 0 0 1 136 10 # Advective Energy (W/m2)
Qtau: 0 N/m2 DN 0 0 0 1 135 10 # Momentum flux (N/m2)
DelSurfHeat: 0 J/m2 INC 0 0 0 1 137 10 # Change in surface heat storage (J/m2)
DelColdCont: 0 J/m2 INC 0 0 0 1 138 10 # Change in snow cold content (J/m2)
BR: 0 - - 0 0 1 1 256 10 # Bowen ratio
EF: 0 - - 0 0 1 1 256 10 # Evaporative fraction

#Water balance components
Snowf: 1 kg/m2s DN 0 0 0 1 161 10000 # Snowfall rate (kg/m2s)
Rainf: 1 kg/m2s DN 0 0 0 1 162 10000 # Rainfall rate (kg/m2s)
RainfConv: 0 kg/m2s DN 0 0 0 1 163 10000 # Convective Rainfall rate (kg/m2s)
TotalPrecip: 0 kg/m2s DN 0 0 0 1 164 10000 # Total Precipitation rate (kg/m2s)
Evap: 1 kg/m2s UP 0 0 0 1 57 10000 # Total Evapotranspiration (kg/m2s)
Qs: 1 kg/m2s OUT 0 0 0 1 235 10000 # Surface runoff (kg/m2s)
Qrec: 0 kg/m2s IN 0 0 0 1 143 10000 # Recharge (kg/m2s)
Qsb: 1 kg/m2s OUT 0 0 0 1 254 10000 # Subsurface runoff (kg/m2s)
Qsm: 0 kg/m2s S2L 0 0 0 1 99 10000 # Snowmelt (kg/m2s)
Qfz: 0 kg/m2s L2S 0 0 0 1 146 10000 # Refreezing of water in the snowpack (kg/m2s)
Qst: 0 kg/m2s - 0 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 0 0 0 1 150 1000 # Change in surface water storage (kg/m2)
DelIntercept: 0 kg/m2 INC 0 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 - 0 0 0 1 152 10 # Snow surface temperature (K)
VegT: 1 K - 0 0 0 1 153 10 # Vegetation canopy temperature (K)
BareSoilT: 1 K - 0 0 0 1 154 10 # Temperature of bare soil (K)
AvgSurfT: 0 K - 0 0 0 1 148 10 # Average surface temperature (K)
RadT: 1 K - 0 0 0 1 156 10 # Surface Radiative Temperature (K)
Albedo: 1 - - 0 0 0 1 84 100 # Surface Albedo (-)
SWE: 1 kg/m2 - 0 0 0 1 65 1000 # Snow Water Equivalent (kg/m2)
SWEVeg: 0 kg/m2 - 0 0 0 1 159 1000 # SWE intercepted by vegetation (kg/m2)
SurfStor: 0 kg/m2 - 0 0 0 1 160 1000 # Surface water storage (kg/m2)
TWS: 1 mm - 0 0 0 1 160 1000 # Terrestrial water storage (mm)
GWS: 0 mm - 0 0 0 1 176 100 # Ground water storage (mm)
WaterTableD: 0 m - 0 0 0 1 174 1 # Water table depth (m)
SWS: 0 mm - 0 0 0 1 333 10 # Surface water storage

#Subsurface state variables
SoilMoist: 1 m3/m3 - 0 0 0 4 86 1000 # Average layer soil moisture (kg/m2)
SoilTemp: 1 K - 0 0 0 4 85 1000 # Average layer soil temperature (K)
SmLiqFrac: 1 m3/m3 - 0 0 0 4 85 100 # Average layer fraction of liquid moisture (-)
SmFrozFrac: 0 - - 0 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 0 0 0 1 166 1 # Potential Evapotranspiration (kg/m2s)
ECanop: 1 kg/m2s UP 0 0 0 1 200 1 # Interception evaporation (kg/m2s)
TVeg: 1 kg/m2s UP 0 0 0 1 210 1 # Vegetation transpiration (kg/m2s)
ESoil: 1 kg/m2s UP 0 0 0 1 199 1 # Bare soil evaporation (kg/m2s)
EWater: 0 kg/m2s UP 0 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: 1 kg/m2 - 0 0 0 1 223 1000 # Total canopy water storage (kg/m2)
EvapSnow: 0 kg/m2s - 0 0 0 1 173 1000 # Snow evaporation (kg/m2s)
SubSnow: 1 kg/m2s - 0 0 0 1 198 1000 # Snow sublimation (kg/m2s)
SubSurf: 0 kg/m2s - 0 0 0 1 175 1000 # Sublimation of the snow free area (kg/m2s)

#Cold season processes
Snowcover: 1 - - 0 0 0 1 66 100 # Snow Cover (-)
SnowDepth: 1 m - 0 0 0 1 66 1000 # Snow Depth (m)
SLiqFrac: 0 - - 0 0 0 1 65 1000 # Fraction of SWE in the liquid phase
SnowTProf: 1 K - 0 0 0 3 239 1000 # Temperature of the snow pack (K)

#Forcings
Wind_f: 1 m/s - 0 0 0 1 177 10 # Near Surface Wind (m/s)
Rainf_f: 1 kg/m2s DN 0 0 0 1 162 1000 # Average rainfall rate
Snowf_f: 0 kg/m2s DN 0 0 0 1 161 1000 # Average snowfall rate
Tair_f: 1 K - 0 0 0 1 11 10 # Near surface air temperature
Qair_f: 1 kg/kg - 0 0 0 1 51 1000 # Near surface specific humidity
Psurf_f: 1 Pa - 0 0 0 1 1 10 # Surface pressure
SWdown_f: 1 W/m2 DN 0 0 0 1 204 10 # Surface incident shortwave radiation
LWdown_f: 1 W/m2 DN 0 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: 1 - - 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: 1 - - 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 - - 0 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 - - 0 0 0 1 190 100 # LAI
SAI: 1 - - 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: 1 - - 0 0 0 1 87 100 # Greenness
Tempbot: 0 - - 0 0 0 1 194 10 # Bottom soil temperature

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

#Carbon variables
GPP: 0 kg/m2s2 DN 1 0 0 1 256 1 # Gross Primary Production
NPP: 0 kg/m2s2 DN 1 0 0 1 256 1 # Net Primary Production
NEE: 0 kg/m2s2 UP 1 0 0 1 256 1 # Net Ecosystem Exchange

#Irrigation
Irrigated water: 0 kg/m2s - 0 0 0 1 333 10 #Irrigation amount

#Temperature
VegGrndT: 0 K - 0 0 0 1 152 10 #vegetated_ground_surface_temperature
VegCanopT: 1 K - 0 0 0 1 152 10 #canopy_air_temperature
AvgGrndT: 1 K - 0 0 0 1 152 10 #average_ground_surface_temperature
VegT2m: 1 K - 0 0 0 1 152 10 #2-m_air_temperature_over_vegetated_part
BARE2MT: 1 K - 0 0 0 1 152 10 #2-m_air_temperature_over_bareground_part
22 changes: 22 additions & 0 deletions lis/configs/Noah-MPv5.0/forcing_variables.txt
Original file line number Diff line number Diff line change
@@ -0,0 +1,22 @@
#ALMA Name select vlevels units
Tair: 1 1 K # Near Surface Air Temperature
Qair: 1 1 kg/kg # Near Surface Specific Humidity
SWdown: 1 1 W/m2 # Incident Shortwave Radiation
SWdirect: 0 1 W/m2 # Incident Shortwave Radiation
SWdiffuse: 0 1 W/m2 # Incident Shortwave Radiation
LWdown: 1 1 W/m2 # Incident Longwave Radiation
Wind_E: 1 1 W/m2 # Eastward Wind
Wind_N: 1 1 m/s # Northward Wind
Psurf: 1 1 Pa # Surface Pressure
Rainf: 1 1 kg/m2s # Rainfall Rate
Snowf: 0 1 kg/m2s # Snowfall Rate
CRainf: 1 1 kg/m2s # Convective Rainfall Rate
Forc_Hgt: 0 1 m # Height of Forcing Variables
Ch: 0 1 - # Surface Exchange Coefficient for Heat
Cm: 0 1 - # Surface Exchange Coefficient for Momentum
Q2sat: 0 1 - # Saturated Mixing Ratio
Emiss: 0 1 - # Surface Emissivity
Cosz: 0 1 - # Cosine of Zenith Angle
Albedo: 0 1 - # Surface Albedo
LPressure: 0 1 Pa # Level pressure
O3: 0 1 - # Ozone concentration
Loading