ERF
Energy Research and Forecasting: An Atmospheric Modeling Code
Classes | Functions | Variables
ERF_NoahMPVegetationTable.H File Reference

Noah-MP's canopy-resistance parameters, roughness length and monthly leaf area index by land-use category, for the two-stream surface energy balance (erf.radiation.seb_vegetation_type). More...

#include <AMReX_REAL.H>
#include <cstdio>
#include <string>
Include dependency graph for ERF_NoahMPVegetationTable.H:
This graph shows which files directly or indirectly include this file:

Go to the source code of this file.

Classes

struct  NoahMPVegetationParams
 

Functions

const NoahMPVegetationParams * noahmp_vegetation_params (int category) noexcept
 The parameters of MODIS land-use category category (1-20), as Noah-MP reads them. More...
 
amrex::Real noahmp_leaf_area_index (const NoahMPVegetationParams &params, amrex::Real day, amrex::Real days_in_year, bool southern) noexcept
 Leaf area index of params on day of year day, interpolated between the monthly values as Noah-MP does for its table-LAI options (PhenologyMainMod): month index 12 * day / days_in_year, linear between the two nearest months, wrapping December to January. In the southern hemisphere Noah-MP shifts the day by half a year first; pass southern = true for that. More...
 
bool noahmp_day_of_year (const std::string &datetime, amrex::Real &day, amrex::Real &days_in_year)
 The day of the year, as Noah-MP counts it, of a date "YYYY-MM-DD[ HH:MM[:SS]]" (ERF's start_datetime). More...
 

Variables

constexpr int noahmp_num_vegetation_categories = 20
 

Detailed Description

Noah-MP's canopy-resistance parameters, roughness length and monthly leaf area index by land-use category, for the two-stream surface energy balance (erf.radiation.seb_vegetation_type).

The 20 categories of the modified IGBP MODIS land-use dataset and their values in &noahmp_modis_parameters of Submodules/Noah-MP/parameters/NoahmpTable.TBL – the table Noah-MP itself reads – copied here so the balance can use them in builds without Noah-MP. The CTest NoahMPSoilTable_MatchesSubmodule (Tests/check_noahmp_soil_table.py) fails if the two drift apart.

Function Documentation

◆ noahmp_day_of_year()

bool noahmp_day_of_year ( const std::string &  datetime,
amrex::Real &  day,
amrex::Real &  days_in_year 
)
inline

The day of the year, as Noah-MP counts it, of a date "YYYY-MM-DD[ HH:MM[:SS]]" (ERF's start_datetime).

The day is Noah-MP's DayJulianInYear, which PhenologyMainMod takes as 0 <= day < days_in_year. It is 0 at 00:00 on 1 January, and the time of day is added as a fraction. days_in_year is 365 or 366.

Returns
false when the date cannot be read, the month is outside 1-12, the day is outside the month, or the time is outside 00:00:00-23:59:60.
132 {
133  int year = 0, month = 0, mday = 0, hour = 0, minute = 0;
134  double second = 0.0;
135  const int n = std::sscanf(datetime.c_str(), "%d-%d-%d %d:%d:%lf",
136  &year, &month, &mday, &hour, &minute, &second);
137  if (n < 3 || month < 1 || month > 12) { return false; }
138  if (n < 4) { hour = 0; }
139  if (n < 5) { minute = 0; }
140  if (n < 6) { second = 0.0; }
141  const bool leap = (year % 4 == 0 && year % 100 != 0) || year % 400 == 0;
142  const int month_days[12] = {31, leap ? 29 : 28, 31, 30, 31, 30, 31, 31, 30, 31, 30, 31};
143  if (mday < 1 || mday > month_days[month - 1] || hour < 0 || hour > 23 ||
144  minute < 0 || minute > 59 || second < 0.0 || second > 60.0) {
145  return false;
146  }
147  int whole = mday - 1;
148  for (int m = 0; m < month - 1; ++m) { whole += month_days[m]; }
149  day = static_cast<amrex::Real>(whole) +
150  static_cast<amrex::Real>((hour * 3600.0 + minute * 60.0 + second) / 86400.0);
151  days_in_year = leap ? amrex::Real(366) : amrex::Real(365);
152  return true;
153 }
Real
amrex::Real Real
Definition: ERF_ShocInterface.H:19
amrex
Definition: ERF_ConsoleIO.cpp:17

Referenced by RadChoice::init_params().

Here is the caller graph for this function:

◆ noahmp_leaf_area_index()

amrex::Real noahmp_leaf_area_index ( const NoahMPVegetationParams &  params,
amrex::Real  day,
amrex::Real  days_in_year,
bool  southern 
)
inlinenoexcept

Leaf area index of params on day of year day, interpolated between the monthly values as Noah-MP does for its table-LAI options (PhenologyMainMod): month index 12 * day / days_in_year, linear between the two nearest months, wrapping December to January. In the southern hemisphere Noah-MP shifts the day by half a year first; pass southern = true for that.

100 {
101  amrex::Real d = day;
102  if (southern) {
103  d = d + amrex::Real(0.5) * days_in_year;
104  while (d >= days_in_year) { d -= days_in_year; }
105  }
106  const amrex::Real month = amrex::Real(12.0) * d / days_in_year;
107  // Noah-MP: IntpMonth1 = MonthCurrent + 0.5 (an integer, so truncated),
108  // IntpMonth2 = IntpMonth1 + 1, weight of month 1 = (IntpMonth1 + 0.5) - MonthCurrent.
109  int m1 = static_cast<int>(month + amrex::Real(0.5));
110  int m2 = m1 + 1;
111  const amrex::Real w1 = (static_cast<amrex::Real>(m1) + amrex::Real(0.5)) - month;
112  const amrex::Real w2 = amrex::Real(1.0) - w1;
113  if (m1 < 1) { m1 = 12; }
114  if (m2 > 12) { m2 = 1; }
115  return w1 * params.lai[m1 - 1] + w2 * params.lai[m2 - 1];
116 }
NoahMPVegetationParams::lai
amrex::Real lai[12]
LAI_JAN ... LAI_DEC monthly leaf area index [m^2/m^2].
Definition: ERF_NoahMPVegetationTable.H:30

Referenced by RadChoice::init_params().

Here is the caller graph for this function:

◆ noahmp_vegetation_params()

const NoahMPVegetationParams* noahmp_vegetation_params ( int  category)
inlinenoexcept

The parameters of MODIS land-use category category (1-20), as Noah-MP reads them.

Returns
nullptr for a category outside 1-20
42 {
43  // {RS, RGL, HS, TOPT, RSMAX, Z0MVT, {LAI_JAN, ..., LAI_DEC}}
44  static constexpr NoahMPVegetationParams table[noahmp_num_vegetation_categories] = {
45  {amrex::Real(125.0), amrex::Real(30.0), amrex::Real(47.35), amrex::Real(298.0), amrex::Real(5000.0), amrex::Real(1.09), // 1 Evergreen Needleleaf Forest
46  {amrex::Real(4.0), amrex::Real(4.0), amrex::Real(4.0), amrex::Real(4.0), amrex::Real(4.0), amrex::Real(4.0), amrex::Real(4.0), amrex::Real(4.0), amrex::Real(4.0), amrex::Real(4.0), amrex::Real(4.0), amrex::Real(4.0)}},
47  {amrex::Real(150.0), amrex::Real(30.0), amrex::Real(41.69), amrex::Real(298.0), amrex::Real(5000.0), amrex::Real(1.10), // 2 Evergreen Broadleaf Forest
48  {amrex::Real(4.5), amrex::Real(4.5), amrex::Real(4.5), amrex::Real(4.5), amrex::Real(4.5), amrex::Real(4.5), amrex::Real(4.5), amrex::Real(4.5), amrex::Real(4.5), amrex::Real(4.5), amrex::Real(4.5), amrex::Real(4.5)}},
49  {amrex::Real(150.0), amrex::Real(30.0), amrex::Real(47.35), amrex::Real(298.0), amrex::Real(5000.0), amrex::Real(0.85), // 3 Deciduous Needleleaf Forest
50  {amrex::Real(0.0), amrex::Real(0.0), amrex::Real(0.0), amrex::Real(0.6), amrex::Real(1.2), amrex::Real(2.0), amrex::Real(2.6), amrex::Real(1.7), amrex::Real(1.0), amrex::Real(0.5), amrex::Real(0.2), amrex::Real(0.0)}},
51  {amrex::Real(100.0), amrex::Real(30.0), amrex::Real(54.53), amrex::Real(298.0), amrex::Real(5000.0), amrex::Real(0.80), // 4 Deciduous Broadleaf Forest
52  {amrex::Real(0.0), amrex::Real(0.0), amrex::Real(0.3), amrex::Real(1.2), amrex::Real(3.0), amrex::Real(4.7), amrex::Real(4.5), amrex::Real(3.4), amrex::Real(1.2), amrex::Real(0.3), amrex::Real(0.0), amrex::Real(0.0)}},
53  {amrex::Real(125.0), amrex::Real(30.0), amrex::Real(51.93), amrex::Real(298.0), amrex::Real(5000.0), amrex::Real(0.80), // 5 Mixed Forests
54  {amrex::Real(2.0), amrex::Real(2.0), amrex::Real(2.2), amrex::Real(2.6), amrex::Real(3.5), amrex::Real(4.3), amrex::Real(4.3), amrex::Real(3.7), amrex::Real(2.6), amrex::Real(2.2), amrex::Real(2.0), amrex::Real(2.0)}},
55  {amrex::Real(300.0), amrex::Real(100.0), amrex::Real(42.00), amrex::Real(298.0), amrex::Real(5000.0), amrex::Real(0.20), // 6 Closed Shrublands
56  {amrex::Real(0.0), amrex::Real(0.0), amrex::Real(0.3), amrex::Real(0.9), amrex::Real(2.2), amrex::Real(3.5), amrex::Real(3.5), amrex::Real(2.5), amrex::Real(0.9), amrex::Real(0.3), amrex::Real(0.0), amrex::Real(0.0)}},
57  {amrex::Real(170.0), amrex::Real(100.0), amrex::Real(39.18), amrex::Real(298.0), amrex::Real(5000.0), amrex::Real(0.06), // 7 Open Shrublands
58  {amrex::Real(0.0), amrex::Real(0.0), amrex::Real(0.2), amrex::Real(0.6), amrex::Real(1.5), amrex::Real(2.3), amrex::Real(2.3), amrex::Real(1.7), amrex::Real(0.6), amrex::Real(0.2), amrex::Real(0.0), amrex::Real(0.0)}},
59  {amrex::Real(300.0), amrex::Real(100.0), amrex::Real(42.00), amrex::Real(298.0), amrex::Real(5000.0), amrex::Real(0.60), // 8 Woody Savannas
60  {amrex::Real(0.2), amrex::Real(0.2), amrex::Real(0.4), amrex::Real(1.0), amrex::Real(2.4), amrex::Real(4.1), amrex::Real(4.1), amrex::Real(2.7), amrex::Real(1.0), amrex::Real(0.4), amrex::Real(0.2), amrex::Real(0.2)}},
61  {amrex::Real(70.0), amrex::Real(65.0), amrex::Real(54.53), amrex::Real(298.0), amrex::Real(5000.0), amrex::Real(0.50), // 9 Savannas
62  {amrex::Real(0.3), amrex::Real(0.3), amrex::Real(0.5), amrex::Real(0.8), amrex::Real(1.8), amrex::Real(3.6), amrex::Real(3.8), amrex::Real(2.1), amrex::Real(0.9), amrex::Real(0.5), amrex::Real(0.3), amrex::Real(0.3)}},
63  {amrex::Real(40.0), amrex::Real(100.0), amrex::Real(36.35), amrex::Real(298.0), amrex::Real(5000.0), amrex::Real(0.12), // 10 Grasslands
64  {amrex::Real(0.4), amrex::Real(0.5), amrex::Real(0.6), amrex::Real(0.7), amrex::Real(1.2), amrex::Real(3.0), amrex::Real(3.5), amrex::Real(1.5), amrex::Real(0.7), amrex::Real(0.6), amrex::Real(0.5), amrex::Real(0.4)}},
65  {amrex::Real(70.0), amrex::Real(65.0), amrex::Real(55.97), amrex::Real(298.0), amrex::Real(5000.0), amrex::Real(0.30), // 11 Permanent wetlands
66  {amrex::Real(0.2), amrex::Real(0.3), amrex::Real(0.3), amrex::Real(0.5), amrex::Real(1.5), amrex::Real(2.9), amrex::Real(3.5), amrex::Real(2.7), amrex::Real(1.2), amrex::Real(0.3), amrex::Real(0.3), amrex::Real(0.2)}},
67  {amrex::Real(40.0), amrex::Real(100.0), amrex::Real(36.25), amrex::Real(298.0), amrex::Real(5000.0), amrex::Real(0.15), // 12 Croplands
68  {amrex::Real(0.0), amrex::Real(0.0), amrex::Real(0.0), amrex::Real(0.0), amrex::Real(1.0), amrex::Real(2.0), amrex::Real(3.0), amrex::Real(3.0), amrex::Real(1.5), amrex::Real(0.0), amrex::Real(0.0), amrex::Real(0.0)}},
69  {amrex::Real(200.0), amrex::Real(999.0), amrex::Real(999.0), amrex::Real(298.0), amrex::Real(5000.0), amrex::Real(1.00), // 13 Urban and Built-Up
70  {amrex::Real(0.0), amrex::Real(0.0), amrex::Real(0.0), amrex::Real(0.0), amrex::Real(0.0), amrex::Real(0.0), amrex::Real(0.0), amrex::Real(0.0), amrex::Real(0.0), amrex::Real(0.0), amrex::Real(0.0), amrex::Real(0.0)}},
71  {amrex::Real(40.0), amrex::Real(100.0), amrex::Real(36.25), amrex::Real(298.0), amrex::Real(5000.0), amrex::Real(0.14), // 14 cropland/natural vegetation mosaic
72  {amrex::Real(0.2), amrex::Real(0.3), amrex::Real(0.3), amrex::Real(0.4), amrex::Real(1.1), amrex::Real(2.5), amrex::Real(3.2), amrex::Real(2.2), amrex::Real(1.1), amrex::Real(0.3), amrex::Real(0.3), amrex::Real(0.2)}},
73  {amrex::Real(999.0), amrex::Real(999.0), amrex::Real(999.0), amrex::Real(298.0), amrex::Real(5000.0), amrex::Real(0.00), // 15 Snow and Ice
74  {amrex::Real(0.0), amrex::Real(0.0), amrex::Real(0.0), amrex::Real(0.0), amrex::Real(0.0), amrex::Real(0.0), amrex::Real(0.0), amrex::Real(0.0), amrex::Real(0.0), amrex::Real(0.0), amrex::Real(0.0), amrex::Real(0.0)}},
75  {amrex::Real(999.0), amrex::Real(999.0), amrex::Real(999.0), amrex::Real(298.0), amrex::Real(5000.0), amrex::Real(0.00), // 16 Barren or Sparsely Vegetated
76  {amrex::Real(0.0), amrex::Real(0.0), amrex::Real(0.0), amrex::Real(0.0), amrex::Real(0.0), amrex::Real(0.0), amrex::Real(0.0), amrex::Real(0.0), amrex::Real(0.0), amrex::Real(0.0), amrex::Real(0.0), amrex::Real(0.0)}},
77  {amrex::Real(100.0), amrex::Real(30.0), amrex::Real(51.75), amrex::Real(298.0), amrex::Real(5000.0), amrex::Real(0.00), // 17 Water
78  {amrex::Real(0.0), amrex::Real(0.0), amrex::Real(0.0), amrex::Real(0.0), amrex::Real(0.0), amrex::Real(0.0), amrex::Real(0.0), amrex::Real(0.0), amrex::Real(0.0), amrex::Real(0.0), amrex::Real(0.0), amrex::Real(0.0)}},
79  {amrex::Real(150.0), amrex::Real(100.0), amrex::Real(42.00), amrex::Real(298.0), amrex::Real(5000.0), amrex::Real(0.30), // 18 Wooded Tundra
80  {amrex::Real(1.0), amrex::Real(1.0), amrex::Real(1.1), amrex::Real(1.3), amrex::Real(1.7), amrex::Real(2.1), amrex::Real(2.1), amrex::Real(1.8), amrex::Real(1.3), amrex::Real(1.1), amrex::Real(1.0), amrex::Real(1.0)}},
81  {amrex::Real(150.0), amrex::Real(100.0), amrex::Real(42.00), amrex::Real(298.0), amrex::Real(5000.0), amrex::Real(0.20), // 19 Mixed Tundra
82  {amrex::Real(0.6), amrex::Real(0.6), amrex::Real(0.7), amrex::Real(0.8), amrex::Real(1.2), amrex::Real(1.8), amrex::Real(1.8), amrex::Real(1.3), amrex::Real(0.8), amrex::Real(0.7), amrex::Real(0.6), amrex::Real(0.6)}},
83  {amrex::Real(200.0), amrex::Real(100.0), amrex::Real(42.00), amrex::Real(298.0), amrex::Real(5000.0), amrex::Real(0.03), // 20 Barren Tundra
84  {amrex::Real(0.0), amrex::Real(0.0), amrex::Real(0.0), amrex::Real(0.0), amrex::Real(0.0), amrex::Real(0.0), amrex::Real(0.0), amrex::Real(0.0), amrex::Real(0.0), amrex::Real(0.0), amrex::Real(0.0), amrex::Real(0.0)}},
85  };
86  if (category < 1 || category > noahmp_num_vegetation_categories) { return nullptr; }
87  return &table[category - 1];
88 }
noahmp_num_vegetation_categories
constexpr int noahmp_num_vegetation_categories
Definition: ERF_NoahMPVegetationTable.H:34
NoahMPVegetationParams
Definition: ERF_NoahMPVegetationTable.H:23

Referenced by RadChoice::init_params(), TwoStreamRadiation::seb_surface_moisture(), and ERF::set_surface_layer_skin().

Here is the caller graph for this function:

Variable Documentation

◆ noahmp_num_vegetation_categories

constexpr int noahmp_num_vegetation_categories = 20
inlineconstexpr

Number of MODIS land-use categories; categories are numbered from 1, as in Noah-MP.

Referenced by noahmp_vegetation_params().