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

Xinanjiang #559

Merged
merged 6 commits into from
Jun 9, 2021
Merged

Xinanjiang #559

merged 6 commits into from
Jun 9, 2021

Conversation

prasanthvkrishna
Copy link
Contributor

@prasanthvkrishna prasanthvkrishna commented May 21, 2021

TYPE: New feature

KEYWORDS: A new option for Noah MP surface runoff estimation: The Xinanjiang Runoff Scheme (XAJ)

SOURCE: Prasanth Valayamkunnath, NCAR

DESCRIPTION OF CHANGES: Changes are made only on Noah MP codes.
M trunk/NDHMS/Land_models/NoahMP/IO_code/module_NoahMP_hrldas_driver.F
M trunk/NDHMS/Land_models/NoahMP/IO_code/module_hrldas_netcdf_io.F
M trunk/NDHMS/Land_models/NoahMP/phys/module_sf_noahmpdrv.F
M trunk/NDHMS/Land_models/NoahMP/phys/module_sf_noahmplsm.F
M trunk/NDHMS/Land_models/NoahMP/run/SOILPARM.TBL

ISSUE: The GitHub Issue that this PR addresses. For issue number 123, it would be:

Fixes #123

TESTS CONDUCTED: This scheme was developed and incorporated in the NWM runoff task. The model is tested for the selected 25 NWM calibration basins (V2.1) across the US, both with default and calibrated parameters. The streamflow results from this version of the model were discussed with the team members.

Notes: This scheme is part of the official Noah MP.

Checklist

Merging the PR depends on following checklist being completed. Add X between each of the square
brackets if they are completed in the PR itself. If a bullet is not relevant to you, please comment
on why below the bullet.

  • Closes issue #xxxx (An issue must exist or be created to be closed. The issue describes and documents the problem and general solution, the PR describes the technical details of the solution.)
  • Tests added (unit tests and/or regression/integration tests)
  • Backwards compatible (I think so, not sure!)
  • Requires new files? Yes, new soil_properties.nc file with additional variables, if this scheme is active.
  • Fully documented
  • Short description in the Development section of NEWS.md

@kafitzgerald
Copy link
Contributor

Looks like is passing CI at least (so w/out the new scheme active).

@prasanthvkrishna @aubreyd do either of you know if this code has been tested with the new scheme active or if there's a quick test we might run?

17, 11.55, 0.030, -10.472, 0.468, 0.454, 0.468, 9.74E-7, 1.12E-5, 0.030, 0.60, 'PLAYA'
18, 2.79, 0.006, -0.472, 0.200, 0.17, 0.069, 1.41E-4, 1.36E-4, 0.006, 0.52, 'LAVA'
19, 2.79, 0.01, -0.472, 0.339, 0.192, 0.069, 4.66E-5, 2.65E-5, 0.01, 0.92, 'WHITE SAND'
19,1 'BB DRYSMC F11 MAXSMC REFSMC SATPSI SATDK SATDW WLTSMC QTZ AXAJ BXAJ XXAJ '
Copy link
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

@aubreyd is this something that will break the soil properties script?

This shouldn't be going into the upcoming release anyway (so an issue for later), but just figured I'd check.

Copy link
Collaborator

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

I have a new version of create_soilproperties.R that handles these new params. Not sure where to stash though since this is not tied to a release (yet).

Copy link
Contributor

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Awesome!

Seems like these things might pile up a bit. I wonder if it make sense to just start a branch over on the repo so that we can start accruing changes in one spot (my memory is short wrt this kind of stuff)? Happy to do so or we can just wait if you'd rather.

Copy link
Contributor Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

@kafitzgerald @rcabell @aubreyd
I just made a working test case for the basin "06917000" to test the new scheme. I have also created the spatial parameter file. Please see this path "/glade/campaign/ral/hap/prasanth/Runoff_Task/Test_code/"
If you want to see the calibration directory, it is here "/glade/campaign/ral/hap/prasanth/Runoff_Task/CALIBRATION/XAJ". Please let me know if you have any questions.

@prasanthvkrishna
Copy link
Contributor Author

Looks like is passing CI at least (so w/out the new scheme active).

@prasanthvkrishna @aubreyd do either of you know if this code has been tested with the new scheme active or if there's a quick test we might run?

This code is developed in the Runoff Task. We ran, default, calibrated runs with this new scheme and compared the NWM performance with the original NWM version.

@rcabell
Copy link
Collaborator

rcabell commented May 24, 2021

Regression tests pass for CONUS on Cheyenne. Do we want to add new tests for this scheme before merging?

@aubreyd
Copy link
Collaborator

aubreyd commented May 24, 2021

Regression tests pass for CONUS on Cheyenne. Do we want to add new tests for this scheme before merging?

@prasanthvkrishna - Are all the new schemes in here or just Xinanjiang? May need documentation at some point... at least the new options and params.

@rcabell - Might be worth adding a runoff_option = 7 test since that is the proposed NWM option (should work out-of-the-box with the new TBL or I can provide new soil_properties file for either CONUS or Croton). We don't have automated testing for all of the NoahMP options now, so not sure how crazy we want to go.

@prasanthvkrishna
Copy link
Contributor Author

Regression tests pass for CONUS on Cheyenne. Do we want to add new tests for this scheme before merging?

@prasanthvkrishna - Are all the new schemes in here or just Xinanjiang? May need documentation at some point... at least the new options and params.

@rcabell - Might be worth adding a runoff_option = 7 test since that is the proposed NWM option (should work out-of-the-box with the new TBL or I can provide new soil_properties file for either CONUS or Croton). We don't have automated testing for all of the NoahMP options now, so not sure how crazy we want to go.

@aubreyd No, this PR has just the Xinanjiang. I can provide the documentation on these runoff schemes and updates in the parameter table.

@kafitzgerald
Copy link
Contributor

@aubreyd @rcabell

I think we'll want to update the soil properties files for the Croton and CONUS test cases and activate the new option in the relevant JSON file for the NWM tests. I think we've generally done this as a follow on PR in the past w/ the first one (similar to this) just ensuring we "do no harm."

@kafitzgerald
Copy link
Contributor

Thanks all!

Just went ahead and updated standard Croton test case in CI and activated new scheme activated in the NWM tests aside from Hawaii yesterday evening. Looks to be failing only where we'd expect differences with the new scheme (though another set of eyes wouldn't hurt).

And we have a CONUS parameter file update so hopefully we can get those running and merge here ASAP since folks are waiting on it.

That said, it'd be nice to get some documentation in the repo for this @prasanthvkrishna (this can wait for a follow on PR if needed, but should go in soon) in line with our poorly advertised contributing guidelines (see number 3 in the checklist for new features).

Something along the lines of this description for the channel loss functionality would be great, but ideally in the trunk/NDHMS/Doc directory as a single markdown file for now as described in the guidelines above. Happy to help if it moves things along.

@rcabell
Copy link
Collaborator

rcabell commented Jun 3, 2021

Croton CI testing shows only data drifts (no n-cores / restart failures). If someone from the science team can take a glance through and give a thumbs up/down on whether those drifts are acceptable, that would be great.

Copy link
Contributor

@kafitzgerald kafitzgerald left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Diffs look reasonable to me (w/ the runoff scheme change we'd expect diffs in a large number of states and fluxes) though I don't have well calibrated expectations here.

Given the other testing that's been done, maybe we launch the CONUS test and pull some diff plots from that if there's interest in digging in? Would need to be a quick review though - this is needed for a few things.

@prasanthvkrishna
Copy link
Contributor Author

Thanks all!

Just went ahead and updated standard Croton test case in CI and activated new scheme activated in the NWM tests aside from Hawaii yesterday evening. Looks to be failing only where we'd expect differences with the new scheme (though another set of eyes wouldn't hurt).

And we have a CONUS parameter file update so hopefully we can get those running and merge here ASAP since folks are waiting on it.

That said, it'd be nice to get some documentation in the repo for this @prasanthvkrishna (this can wait for a follow on PR if needed, but should go in soon) in line with our poorly advertised contributing guidelines (see number 3 in the checklist for new features).

Something along the lines of this description for the channel loss functionality would be great, but ideally in the trunk/NDHMS/Doc directory as a single markdown file for now as described in the guidelines above. Happy to help if it moves things along.

@kafitzgerald Thanks for the review. I will provide the documentation and a summary text for the release note.

@rcabell
Copy link
Collaborator

rcabell commented Jun 8, 2021

CONUS testing with Long-Range physics produces only data differences, listed below:

channel_rt

     Variable    Count         Sum      Min      Max    Range   Mean  StdDev
0  streamflow  1771422  278007.274  -87.774  223.815  311.588  0.157   1.852
1       nudge     5027    -891.300 -538.997   13.790  552.787 -0.177   7.661
2   q_lateral  1056247   16134.686  -55.279   48.223  103.502  0.015   0.275
3    velocity  1729695   25568.012   -2.569    3.295    5.864  0.015   0.103
4        Head  1775503   39823.699 -739.615   34.292  773.907  0.022   0.588
     Variable    Count         Sum      Min       Max     Range   Mean  StdDev
0  streamflow  1992120  201595.561 -308.105  1579.162  1887.267  0.101   3.610
1       nudge     6003    -215.434  -41.990    92.167   134.156 -0.036   1.825
2   q_lateral  1397954    -862.603  -39.345    14.080    53.425 -0.001   0.120
3    velocity  1947369    5871.430   -3.048     3.294     6.342  0.003   0.074
4        Head  2000348   15385.224  -31.858   154.437   186.295  0.008   0.205

chanobs

     Variable  Count    Sum     Min     Max    Range   Mean  StdDev
0  streamflow    877  0.797 -63.669  41.337  105.006  0.001   3.265
     Variable  Count     Sum     Min     Max    Range   Mean  StdDev
0  streamflow   1155 -73.694 -60.428  58.991  119.419 -0.064    3.17

lakeout

         Variable  Count       Sum     Min      Max    Range   Mean  StdDev
0  water_sfc_elev   1875     5.740  -0.338    0.354    0.692  0.003   0.026
1          inflow   3711  2450.423 -80.473  397.558  478.030  0.660   8.372
2         outflow   1873   238.201 -16.966   27.736   44.703  0.127   1.476
         Variable  Count      Sum      Min     Max    Range   Mean  StdDev
0  water_sfc_elev   2587   12.762   -0.735   0.389    1.124  0.005   0.032
1          inflow   4288  992.072 -248.893  74.301  323.194  0.231   5.005
2         outflow   2585  496.421   -9.293  44.595   53.888  0.192   1.681

gwout

  Variable    Count       Sum    Min    Max  Range   Mean  StdDev
0   inflow  1012713   -92.766 -0.851  2.323  3.174 -0.000   0.006
1  outflow   857474   -43.421 -0.848  1.968  2.816 -0.000   0.005
2    depth   700957 -1358.921 -1.003  2.587  3.591 -0.002   0.048
  Variable    Count       Sum    Min    Max   Range   Mean  StdDev
0   inflow  1424876  -406.628 -0.905  4.416   5.322 -0.000   0.009
1  outflow  1286173  -218.043 -0.960  4.489   5.449 -0.000   0.008
2    depth  1080065 -8867.169 -4.617  5.521  10.138 -0.008   0.093

ldasout

     Variable    Count         Sum  ...      Range   Mean   StdDev
0        FVEG        2      -1.363  ...      0.080 -0.682    0.040
1       EMISS   372935      -0.127  ...      0.009 -0.000    0.000
2         FSA  2397794  -12079.591  ...     49.901 -0.005    0.065
3        FIRA  5587632  135661.758  ...     62.188  0.024    0.166
4      GRDFLX  5778605  188748.243  ...    575.617  0.033    0.725
..        ...      ...         ...  ...        ...    ...      ...
70    SOILSAT  5681883   -3342.028  ...      0.140 -0.001    0.002
71  SNOWT_AVG   512785 -431428.640  ...  20544.320 -0.841  111.111
72     ALBSND   967207      29.699  ...      1.870  0.000    0.010
73     ALBSNI   987993      29.109  ...      1.869  0.000    0.010
74      QRAIN   569952      -0.071  ...      0.000 -0.000    0.000

[75 rows x 8 columns]
     Variable    Count          Sum  ...      Range   Mean   StdDev
0        FVEG        5        2.167  ...      1.210  0.433    0.459
1       EMISS   582360       -0.949  ...      0.012 -0.000    0.000
2         FSA  3594055   -21707.886  ...     96.075 -0.006    0.094
3        FIRA  7201274   350709.868  ...     80.621  0.049    0.343
4      GRDFLX  7430236   165835.376  ...    518.198  0.022    1.181
..        ...      ...          ...  ...        ...    ...      ...
70    SOILSAT  7220626    -8316.211  ...      0.164 -0.001    0.004
71  SNOWT_AVG   591345 -1397007.949  ...  20544.320 -2.362  206.064
72     ALBSND  1706978      -37.675  ...      1.890 -0.000    0.012
73     ALBSNI  1740475      -36.396  ...      1.882 -0.000    0.012
74      QRAIN   578802       -0.097  ...      0.000 -0.000    0.000

[75 rows x 8 columns]

restart_hydro

None
      Variable    Count         Sum      Min      Max    Range   Mean  StdDev
0         stc1  4390930   24153.309   -4.814    9.290   14.104  0.006   0.051
1         smc1  6062675   -6493.181   -0.199    0.183    0.382 -0.001   0.005
2       sh2ox1  6062745   -6385.391   -0.199    0.183    0.382 -0.001   0.005
3         stc2  3553846    1310.857   -1.087    0.604    1.691  0.000   0.018
4         smc2  4966696   -3854.533   -0.065    0.075    0.140 -0.001   0.002
5       sh2ox2  4967454   -3822.898   -0.065    0.075    0.140 -0.001   0.002
6         stc3  1954730   -1365.851   -0.847    0.170    1.017 -0.001   0.007
7         smc3  4258046   -1149.084   -0.030    0.058    0.088 -0.000   0.001
8       sh2ox3  4262585   -1139.656   -0.030    0.058    0.088 -0.000   0.001
9         stc4   761669     -67.183   -0.172    0.168    0.340 -0.000   0.004
10        smc4  3322269    -307.950   -0.029    0.043    0.072 -0.000   0.001
11      sh2ox4  3330681    -307.887   -0.029    0.043    0.072 -0.000   0.001
12     infxsrt  2715959  268717.094  -17.290   19.323   36.613  0.099   0.353
13    soldrain  3299216    -547.560   -0.644    0.616    1.260 -0.000   0.006
14       hlink  1775503   39823.699 -739.615   34.292  773.907  0.022   0.588
15      qlink1  1771422  278007.274  -87.774  223.815  311.588  0.157   1.852
16      qlink2  1771422  278007.274  -87.774  223.815  311.588  0.157   1.852
17       resht     1875       5.740   -0.338    0.354    0.692  0.003   0.026
18      qlakeo     1873     238.201  -16.966   27.736   44.703  0.127   1.476
19      qlakei     3711    2450.423  -80.473  397.558  478.030  0.660   8.372
20  z_gwsubbas   702130      -1.359   -0.001    0.003    0.004 -0.000   0.000
      Variable    Count         Sum      Min       Max     Range   Mean  StdDev
0         stc1  5959483   27857.593   -5.508     8.707    14.215  0.005   0.110
1         smc1  7619366  -11926.100   -0.229     0.141     0.370 -0.002   0.008
2       sh2ox1  7619434  -11317.587   -0.229     0.141     0.370 -0.001   0.007
3         stc2  5357715  -55540.317   -1.879     1.317     3.196 -0.010   0.089
4         smc2  6399762   -9984.222   -0.169     0.086     0.255 -0.002   0.005
5       sh2ox2  6400469   -9433.546   -0.169     0.086     0.255 -0.001   0.005
6         stc3  3934915  -11588.045   -0.798     0.188     0.986 -0.003   0.021
7         smc3  5705017   -3317.550   -0.070     0.063     0.133 -0.001   0.002
8       sh2ox3  5710520   -3241.136   -0.070     0.063     0.133 -0.001   0.002
9         stc4  1925804    -348.172   -0.260     0.161     0.421 -0.000   0.007
10        smc4  4779425    -984.753   -0.022     0.049     0.071 -0.000   0.001
11      sh2ox4  4787179    -984.763   -0.022     0.049     0.071 -0.000   0.001
12     infxsrt  2803224  430719.561  -17.413    42.691    60.104  0.154   0.493
13    soldrain  4773306   -2685.476   -0.671     0.722     1.394 -0.001   0.009
14       hlink  2000348   15385.224  -31.858   154.437   186.295  0.008   0.205
15      qlink1  1992120  201595.561 -308.105  1579.162  1887.267  0.101   3.610
16      qlink2  1992120  201595.561 -308.105  1579.162  1887.267  0.101   3.610
17       resht     2587      12.762   -0.735     0.389     1.124  0.005   0.032
18      qlakeo     2585     496.421   -9.293    44.595    53.888  0.192   1.681
19      qlakei     4288     992.072 -248.893    74.301   323.194  0.231   5.005
20  z_gwsubbas  1081243      -8.867   -0.005     0.006     0.010 -0.000   0.000

restart_lsm

None
     Variable     Count          Sum      Min      Max     Range   Mean  StdDev
0      SOIL_T  10661175    24031.131   -4.814    9.290    14.104  0.002   0.035
1      SNOW_T     15923   -21300.182 -273.160  273.160   546.320 -1.338  23.873
2         SMC  18609686   -11804.748   -0.199    0.183     0.382 -0.001   0.003
3        SH2O  18623465   -11655.833   -0.199    0.183     0.382 -0.001   0.003
4       ZSNSO   8103612       16.623   -0.095    0.109     0.204  0.000   0.000
5       SNICE   3329326     -568.063  -43.434   48.506    91.940 -0.000   0.134
6       SNLIQ   3503201     -195.276   -7.983    7.851    15.834 -0.000   0.023
7       QSNOW   1070851       -0.000   -0.000    0.000     0.000 -0.000   0.000
8        FWET   2164674    -4514.862   -1.000    1.000     2.000 -0.002   0.014
9      SNEQVO   1305155      -80.678  -11.617   13.975    25.592 -0.000   0.049
10        EAH   5363762 -3963288.347 -864.984  651.231  1516.215 -0.739   4.698
11        TAH   4691285    21851.118  -20.416    3.308    23.724  0.005   0.029
12         CM   5143790       41.979   -0.012    0.015     0.028  0.000   0.000
13         CH   5167777      170.691   -0.046    0.048     0.094  0.000   0.000
14      ISNOW       122       78.000   -1.000    1.000     2.000  0.639   0.769
15     CANLIQ   2106793     -446.478   -0.377    0.431     0.808 -0.000   0.001
16     CANICE   1024527       -0.666   -0.003    0.004     0.008 -0.000   0.000
17      SNEQV   1345979      -61.737   -8.026    7.881    15.906 -0.000   0.039
18      SNOWH   1388304       -1.306   -0.009    0.012     0.021 -0.000   0.000
19         TV   4791195    18861.101  -20.770    3.520    24.290  0.004   0.029
20         TG   4335779    42280.137   -6.028    5.542    11.571  0.010   0.042
21     ACMELT   1514800      596.111   -7.887    8.045    15.932  0.000   0.037
22      TAUSS    823870       86.089   -2.097    3.406     5.503  0.000   0.016
23       QSFC   4916689       17.435   -0.005    0.003     0.008  0.000   0.000
24  SFCRUNOFF   6983022  2836288.413  -82.647   42.219   124.866  0.406   1.796
25   UDRUNOFF   2738557    -4068.586   -5.678    8.488    14.167 -0.001   0.069
26    ACCECAN   3136568     2286.333   -1.173    0.446     1.620  0.001   0.003
27    ACCEDIR   5454535   -31171.760   -2.490    1.193     3.684 -0.006   0.027
28   ACCETRAN   4967658      249.365   -0.188    0.172     0.360  0.000   0.001
     Variable     Count           Sum  ...     Range   Mean  StdDev
0      SOIL_T  17177917 -3.961894e+04  ...    14.215 -0.002   0.082
1      SNOW_T     25222 -4.641882e+04  ...   546.320 -1.840  40.148
2         SMC  24503570 -2.621262e+04  ...     0.370 -0.001   0.005
3        SH2O  24517602 -2.497703e+04  ...     0.370 -0.001   0.005
4       ZSNSO   9944561 -3.444700e+01  ...     0.220 -0.000   0.001
5       SNICE   4093359  7.153740e+03  ...    85.986  0.002   0.311
6       SNLIQ   4025255 -4.593365e+03  ...    16.150 -0.001   0.045
7       QSNOW    813085 -0.000000e+00  ...     0.000 -0.000   0.000
8        FWET   1915509 -5.845052e+03  ...     2.000 -0.003   0.036
9      SNEQVO   1600132  2.883680e+03  ...    26.930  0.002   0.092
10        EAH   6928566 -1.014604e+07  ...  2216.532 -1.464  14.860
11        TAH   6220330  5.743055e+04  ...    33.214  0.009   0.056
12         CM   6871650  1.132810e+02  ...     0.050  0.000   0.000
13         CH   6901979  3.917870e+02  ...     0.184  0.000   0.000
14      ISNOW       546  1.700000e+02  ...     2.000  0.311   0.950
15     CANLIQ   1925580 -5.963490e+02  ...     0.828 -0.000   0.003
16     CANICE    813136 -2.675000e+00  ...     0.033 -0.000   0.000
17      SNEQV   1689341  3.678763e+03  ...    18.513  0.002   0.079
18      SNOWH   1707943  7.304000e+00  ...     0.031  0.000   0.000
19         TV   6324650  5.211214e+04  ...    34.561  0.008   0.055
20         TG   5965586  9.248912e+04  ...    14.178  0.016   0.083
21     ACMELT   1924534 -1.640632e+03  ...    18.535 -0.001   0.077
22      TAUSS   1327895  3.640700e+01  ...     4.836  0.000   0.015
23       QSFC   6577981  3.910900e+01  ...     0.012  0.000   0.000
24  SFCRUNOFF   8614704  7.314143e+06  ...   157.706  0.849   3.160
25   UDRUNOFF   4316377 -5.069207e+04  ...    44.493 -0.012   0.246
26    ACCECAN   4551926  6.835422e+03  ...     3.771  0.002   0.007
27    ACCEDIR   7279834 -1.098942e+05  ...    10.365 -0.015   0.071
28   ACCETRAN   6747108 -3.259390e+02  ...     1.259 -0.000   0.004

[29 rows x 8 columns]

restart_nudging

          Variable   Count       Sum      Min      Max     Range   Mean  StdDev
0  model_discharge  354311  9604.451 -564.075  644.436  1208.511  0.027   2.826
1            nudge    5027  -891.300 -538.997   13.790   552.787 -0.177   7.661
          Variable   Count        Sum  ...     Range   Mean  StdDev
0  model_discharge  896851  40279.691  ...  1208.511  0.045   2.621
1            nudge    6003   -215.434  ...   134.156 -0.036   1.825

[2 rows x 8 columns]

@rcabell
Copy link
Collaborator

rcabell commented Jun 8, 2021

CONUS testing with full routing physics produces only data differences, listed below:

channel_rt

     Variable    Count         Sum     Min     Max    Range   Mean  StdDev
0  streamflow  1656905  131531.655 -23.307  98.586  121.893  0.079   1.112
1       nudge     5838    -130.184 -26.226   4.903   31.130 -0.022   0.539
2   q_lateral  1526981    6489.239 -11.007  57.974   68.981  0.004   0.122
3    velocity  1622369    8971.582  -1.306   2.948    4.254  0.006   0.052
4        Head  1667818   18162.786  -9.237  13.775   23.012  0.011   0.130
     Variable    Count        Sum      Min      Max    Range   Mean  StdDev
0  streamflow  1842065  93468.630 -170.194  180.546  350.740  0.051   1.237
1       nudge     6122    176.569  -71.091  278.789  349.879  0.029   3.756
2   q_lateral  1738249   -741.764  -29.198    7.867   37.065 -0.000   0.069
3    velocity  1806278   2533.764   -2.621    2.651    5.272  0.001   0.043
4        Head  1854713   6778.192   -6.392  136.372  142.764  0.004   0.143

chanobs

     Variable  Count     Sum     Min    Max  Range   Mean  StdDev
0  streamflow   1013 -83.981 -10.896  8.924  19.82 -0.083   0.873
     Variable  Count     Sum    Min     Max   Range   Mean  StdDev
0  streamflow   1100 -52.068 -9.548  12.473  22.021 -0.047   1.001

lakeout

         Variable  Count      Sum     Min     Max    Range   Mean  StdDev
0  water_sfc_elev   1976    3.819  -0.153   0.308    0.461  0.002   0.017
1          inflow   4102  592.821 -35.750  96.399  132.150  0.145   2.669
2         outflow   1974   96.229  -5.391  18.059   23.451  0.049   0.715
         Variable  Count      Sum     Min     Max    Range   Mean  StdDev
0  water_sfc_elev   2454    7.761  -0.187   0.352    0.540  0.003   0.023
1          inflow   4378  126.551 -91.830  25.285  117.115  0.029   2.042
2         outflow   2451  125.190 -10.628  11.860   22.487  0.051   0.734

gwout

  Variable    Count       Sum     Min     Max   Range   Mean  StdDev
0   inflow  1621188   -77.624  -2.356   0.899   3.255 -0.000   0.006
1  outflow  1522144   -47.303  -0.869   0.826   1.695 -0.000   0.004
2    depth  1438037  2332.733 -39.692  37.363  77.055  0.002   0.245
  Variable    Count       Sum     Min     Max    Range   Mean  StdDev
0   inflow  1798552   -74.190  -1.314   0.850    2.165 -0.000   0.006
1  outflow  1736941   -48.611  -1.214   0.720    1.934 -0.000   0.005
2    depth  1649763  7522.898 -62.541  72.505  135.046  0.005   0.405

rtout

       Variable      Count           Sum  ...      Range   Mean  StdDev
0    zwattablrt    1278878  2.086635e+05  ...      4.000  0.163   0.691
1  sfcheadsubrt   52420025  1.589939e+07  ...  10162.759  0.303   3.724
2    QSTRMVOLRT    1735801  1.117775e+07  ...  14392.184  6.440  42.912
3       QBDRYRT       1741 -1.030273e+03  ...    103.361 -0.592   3.828
4        SOIL_M  315381587 -7.505986e+04  ...      0.842 -0.000   0.004

[5 rows x 8 columns]
       Variable      Count           Sum  ...      Range   Mean  StdDev
0    zwattablrt    1801619  1.536477e+05  ...      4.000  0.085   0.704
1  sfcheadsubrt   53599189  2.069484e+07  ...  13941.218  0.386   4.426
2    QSTRMVOLRT    2439924  2.155969e+07  ...  22221.239  8.836  60.114
3       QBDRYRT       2529 -6.827520e+03  ...    261.449 -2.700  11.295
4        SOIL_M  395619026 -1.129229e+05  ...      0.887 -0.000   0.004

[5 rows x 8 columns]

ldasout

     Variable    Count         Sum  ...      Range   Mean  StdDev
0       EMISS   269254      -0.014  ...      0.008 -0.000   0.000
1         FSA  1520951   -3752.242  ...     31.705 -0.002   0.069
2        FIRA  5603488   50444.655  ...     86.280  0.009   0.205
3      GRDFLX  5936729   97426.942  ...    434.096  0.016   1.027
4         HFX  5901139  113666.683  ...    465.800  0.019   1.278
..        ...      ...         ...  ...        ...    ...     ...
68    SOILSAT  5639699   -1903.321  ...      0.333 -0.000   0.003
69  SNOWT_AVG   463026 -184897.705  ...  20544.320 -0.399  97.832
70     ALBSND   725227      62.796  ...      1.872  0.000   0.010
71     ALBSNI   741857      61.663  ...      1.866  0.000   0.010
72      QRAIN   562330      -0.108  ...      0.000 -0.000   0.000

[73 rows x 8 columns]
     Variable    Count         Sum  ...      Range   Mean   StdDev
0        FVEG        1       0.692  ...      0.000  0.692    0.000
1       EMISS   416207      -0.505  ...      0.009 -0.000    0.000
2         FSA  2411159   -3288.130  ...     59.891 -0.001    0.087
3        FIRA  6917686   27901.449  ...     93.884  0.004    0.235
4      GRDFLX  7361525   40654.541  ...    393.687  0.006    1.050
..        ...      ...         ...  ...        ...    ...      ...
69    SOILSAT  6894244   -2933.963  ...      0.372 -0.000    0.003
70  SNOWT_AVG   564281 -873129.034  ...  20544.320 -1.547  150.412
71     ALBSND  1353339      34.251  ...      1.883  0.000    0.011
72     ALBSNI  1382878      33.442  ...      1.883  0.000    0.011
73      QRAIN   559973      -0.070  ...      0.000 -0.000    0.000

[74 rows x 8 columns]

restart_hydro

None
        Variable     Count           Sum  ...      Range   Mean  StdDev
0           stc1   4155504  1.257013e+04  ...     12.661  0.003   0.055
1           smc1   6397832 -1.646223e+03  ...      0.356 -0.000   0.005
2         sh2ox1   6400260 -1.578064e+03  ...      0.356 -0.000   0.005
3           stc2   3409136 -1.621658e+03  ...      2.454 -0.000   0.020
4           smc2   5600881 -1.792645e+03  ...      0.256 -0.000   0.002
5         sh2ox2   5605120 -1.764950e+03  ...      0.256 -0.000   0.002
6           stc3   1848038 -7.737650e+02  ...      0.716 -0.000   0.006
7           smc3   4875043 -9.362590e+02  ...      0.162 -0.000   0.002
8         sh2ox3   4880039 -9.338050e+02  ...      0.162 -0.000   0.002
9           stc4    752234 -4.093400e+01  ...      0.642 -0.000   0.005
10          smc4   4066710 -4.143590e+02  ...      0.190 -0.000   0.001
11        sh2ox4   4067561 -4.144210e+02  ...      0.190 -0.000   0.001
12       infxsrt   4449981  1.107588e+06  ...    680.838  0.249   2.023
13      soldrain   4044056 -6.254760e+02  ...     10.262 -0.000   0.009
14     sfcheadrt   4640890  9.937117e+05  ...    630.702  0.214   1.747
15       QBDRYRT      1741 -1.030273e+03  ...    103.361 -0.592   3.828
16      infxswgt  49348306 -1.000000e-03  ...      2.000 -0.000   0.078
17  sfcheadsubrt  52420025  1.589939e+07  ...  10162.759  0.303   3.724
18      sh2owgt1  90281947 -5.000000e-03  ...      8.171 -0.000   0.010
19      sh2owgt2  79600390 -6.000000e-03  ...      6.020 -0.000   0.007
20      sh2owgt3  70132807  2.000000e-03  ...      4.229  0.000   0.006
21      sh2owgt4  60943401  3.000000e-03  ...      1.742  0.000   0.004
22    qstrmvolrt   1735801  1.117775e+07  ...  14392.184  6.440  42.912
23         hlink   1667818  1.816279e+04  ...     23.012  0.011   0.130
24        qlink1   1656905  1.315317e+05  ...    121.893  0.079   1.112
25        qlink2   1656905  1.315317e+05  ...    121.893  0.079   1.112
26         resht      1976  3.819000e+00  ...      0.461  0.002   0.017
27        qlakeo      1974  9.622900e+01  ...     23.451  0.049   0.715
28        qlakei      4102  5.928210e+02  ...    132.150  0.145   2.669
29    z_gwsubbas   1439211  2.333000e+00  ...      0.077  0.000   0.000

[30 rows x 8 columns]
        Variable      Count           Sum  ...      Range   Mean  StdDev
0           stc1    5273633 -1.789976e+04  ...     14.700 -0.003   0.083
1           smc1    7730887 -2.519125e+03  ...      0.423 -0.000   0.005
2         sh2ox1    7732982 -2.197212e+03  ...      0.423 -0.000   0.005
3           stc2    4727691 -3.878239e+04  ...      4.690 -0.008   0.073
4           smc2    6817537 -3.365201e+03  ...      0.276 -0.000   0.003
5         sh2ox2    6821352 -3.062172e+03  ...      0.276 -0.000   0.003
6           stc3    3299094 -5.026133e+03  ...      1.107 -0.002   0.016
7           smc3    6117464 -1.183650e+03  ...      0.165 -0.000   0.002
8         sh2ox3    6120757 -1.154704e+03  ...      0.165 -0.000   0.002
9           stc4    1633206 -1.267330e+02  ...      0.623 -0.000   0.007
10          smc4    5317469 -6.435400e+02  ...      0.185 -0.000   0.001
11        sh2ox4    5315771 -6.435920e+02  ...      0.185 -0.000   0.001
12       infxsrt    4538121  1.437196e+06  ...    841.718  0.317   2.188
13      soldrain    5307142 -1.171276e+03  ...      9.183 -0.000   0.009
14     sfcheadrt    4772567  1.293428e+06  ...    796.637  0.271   1.887
15       QBDRYRT       2529 -6.827520e+03  ...    261.449 -2.700  11.295
16      infxswgt   51394779 -1.000000e-03  ...      2.000 -0.000   0.087
17  sfcheadsubrt   53599189  2.069484e+07  ...  13941.218  0.386   4.426
18      sh2owgt1  111376732  1.000000e-03  ...      8.596  0.000   0.011
19      sh2owgt2   99388714  1.000000e-02  ...      7.548  0.000   0.009
20      sh2owgt3   89777373  5.000000e-03  ...      6.605  0.000   0.008
21      sh2owgt4   80933163  9.000000e-03  ...      1.774  0.000   0.006
22    qstrmvolrt    2439924  2.155969e+07  ...  22221.239  8.836  60.114
23         hlink    1854713  6.778192e+03  ...    142.764  0.004   0.143
24        qlink1    1842065  9.346863e+04  ...    350.740  0.051   1.237
25        qlink2    1842065  9.346863e+04  ...    350.740  0.051   1.237
26         resht       2454  7.761000e+00  ...      0.540  0.003   0.023
27        qlakeo       2451  1.251900e+02  ...     22.487  0.051   0.734
28        qlakei       4378  1.265510e+02  ...    117.115  0.029   2.042
29    z_gwsubbas    1651025  7.523000e+00  ...      0.135  0.000   0.000

[30 rows x 8 columns]

restart_lsm

None
     Variable     Count           Sum  ...      Range   Mean  StdDev
0      SOIL_T  10164912  1.013377e+04  ...     12.661  0.001   0.037
1      SNOW_T     17925 -8.464254e+03  ...    546.320 -0.472  20.499
2         SMC  20940466 -4.789486e+03  ...      0.356 -0.000   0.003
3        SH2O  20952980 -4.691241e+03  ...      0.356 -0.000   0.003
4       ZSNSO   6893592 -1.868000e+00  ...      0.180 -0.000   0.000
5       SNICE   2891041  1.537817e+03  ...     70.371  0.001   0.141
6       SNLIQ   3133684 -4.495490e+02  ...     14.307 -0.000   0.020
7       QSNOW   1057778 -0.000000e+00  ...      0.000 -0.000   0.000
8        FWET   2152468 -6.700677e+03  ...      2.000 -0.003   0.029
9      SNEQVO   1139430  1.445860e+03  ...     26.292  0.001   0.043
10        EAH   5386943 -1.856760e+06  ...   1763.032 -0.345   7.304
11        TAH   4426819  7.844991e+03  ...      6.554  0.002   0.038
12         CM   5031846  1.413500e+01  ...      0.037  0.000   0.000
13         CH   5064121  6.157000e+01  ...      0.126  0.000   0.000
14      ISNOW       101  3.100000e+01  ...      2.000  0.307   0.952
15     CANLIQ   2114952 -5.767530e+02  ...      0.936 -0.000   0.002
16     CANICE   1012229 -9.700000e-02  ...      0.004 -0.000   0.000
17      SNEQV   1187143  1.526421e+03  ...     14.353  0.001   0.035
18      SNOWH   1229689  1.134000e+00  ...      0.036  0.000   0.000
19         TV   4575224  3.851161e+03  ...      7.974  0.001   0.039
20         TG   4011593  1.721964e+04  ...     17.878  0.004   0.056
21     ACMELT   1436524 -1.205435e+03  ...     14.301 -0.001   0.032
22      TAUSS    650038  4.337000e+01  ...      5.283  0.000   0.012
23       QSFC   4564595  8.694000e+00  ...      0.015  0.000   0.000
24  SFCRUNOFF   7053628  1.890102e+07  ...  11459.825  2.680  32.308
25   UDRUNOFF   3555243 -8.689312e+03  ...    218.862 -0.002   0.166
26    ACCECAN   3225073  4.910385e+03  ...      2.462  0.002   0.009
27    ACCEDIR   5630095 -2.137223e+04  ...      7.177 -0.004   0.048
28   ACCETRAN   5223602  2.247180e+02  ...      0.528  0.000   0.002

[29 rows x 8 columns]
     Variable     Count           Sum  ...      Range   Mean  StdDev
0      SOIL_T  14933624 -6.183501e+04  ...     14.700 -0.004   0.065
1      SNOW_T     28106 -2.375703e+04  ...    546.320 -0.845  29.542
2         SMC  25983357 -7.711516e+03  ...      0.423 -0.000   0.003
3        SH2O  25990862 -7.057680e+03  ...      0.423 -0.000   0.003
4       ZSNSO   9121507 -2.704100e+01  ...      0.238 -0.000   0.000
5       SNICE   3814263  4.220058e+03  ...     88.866  0.001   0.262
6       SNLIQ   3800638 -2.780178e+03  ...     16.093 -0.001   0.035
7       QSNOW    794037 -0.000000e+00  ...      0.000 -0.000   0.000
8        FWET   1879438 -3.404921e+03  ...      2.000 -0.002   0.029
9      SNEQVO   1487216  1.886331e+03  ...     18.431  0.001   0.073
10        EAH   6638913 -1.384643e+06  ...   2447.127 -0.209  11.783
11        TAH   5451981  6.505352e+03  ...     21.279  0.001   0.046
12         CM   6343083  2.035200e+01  ...      0.039  0.000   0.000
13         CH   6388617  6.542600e+01  ...      0.142  0.000   0.000
14      ISNOW       329  8.700000e+01  ...      2.000  0.264   0.964
15     CANLIQ   1903710 -3.287160e+02  ...      0.866 -0.000   0.003
16     CANICE    794934 -1.548000e+00  ...      0.040 -0.000   0.000
17      SNEQV   1574861  2.277380e+03  ...     16.154  0.001   0.063
18      SNOWH   1606776  5.584000e+00  ...      0.041  0.000   0.000
19         TV   5632356  4.185322e+03  ...     24.145  0.001   0.046
20         TG   5194369  1.126166e+04  ...     16.410  0.002   0.062
21     ACMELT   1879018 -1.042148e+03  ...     16.207 -0.001   0.060
22      TAUSS   1160710  6.292700e+01  ...      4.368  0.000   0.014
23       QSFC   5826894  5.597000e+00  ...      0.024  0.000   0.000
24  SFCRUNOFF   8484484  3.575350e+07  ...  26910.438  4.214  59.588
25   UDRUNOFF   4990634 -3.233155e+04  ...    407.296 -0.006   0.306
26    ACCECAN   4615756  8.588254e+03  ...      6.112  0.002   0.012
27    ACCEDIR   7254991 -3.758447e+04  ...     16.946 -0.005   0.071
28   ACCETRAN   6782408  1.265810e+02  ...      1.212  0.000   0.003

[29 rows x 8 columns]

restart_nudging

          Variable   Count       Sum     Min     Max   Range   Mean  StdDev
0  model_discharge  485653  3007.212 -22.760  26.562  49.321  0.006   0.401
1            nudge    5838  -130.184 -26.226   4.903  31.130 -0.022   0.539
          Variable    Count        Sum  ...     Range   Mean  StdDev
0  model_discharge  1050361  13568.083  ...  1036.691  0.013   1.715
1            nudge     6122    176.569  ...   349.879  0.029   3.756

[2 rows x 8 columns]

@aubreyd aubreyd requested a review from gochis June 8, 2021 15:56
@prasanthvkrishna
Copy link
Contributor Author

prasanthvkrishna commented Jun 8, 2021

CONUS testing with full routing physics produces only data differences, listed below:

channel_rt

     Variable    Count         Sum     Min     Max    Range   Mean  StdDev
0  streamflow  1656905  131531.655 -23.307  98.586  121.893  0.079   1.112
1       nudge     5838    -130.184 -26.226   4.903   31.130 -0.022   0.539
2   q_lateral  1526981    6489.239 -11.007  57.974   68.981  0.004   0.122
3    velocity  1622369    8971.582  -1.306   2.948    4.254  0.006   0.052
4        Head  1667818   18162.786  -9.237  13.775   23.012  0.011   0.130
     Variable    Count        Sum      Min      Max    Range   Mean  StdDev
0  streamflow  1842065  93468.630 -170.194  180.546  350.740  0.051   1.237
1       nudge     6122    176.569  -71.091  278.789  349.879  0.029   3.756
2   q_lateral  1738249   -741.764  -29.198    7.867   37.065 -0.000   0.069
3    velocity  1806278   2533.764   -2.621    2.651    5.272  0.001   0.043
4        Head  1854713   6778.192   -6.392  136.372  142.764  0.004   0.143

chanobs

     Variable  Count     Sum     Min    Max  Range   Mean  StdDev
0  streamflow   1013 -83.981 -10.896  8.924  19.82 -0.083   0.873
     Variable  Count     Sum    Min     Max   Range   Mean  StdDev
0  streamflow   1100 -52.068 -9.548  12.473  22.021 -0.047   1.001

lakeout

         Variable  Count      Sum     Min     Max    Range   Mean  StdDev
0  water_sfc_elev   1976    3.819  -0.153   0.308    0.461  0.002   0.017
1          inflow   4102  592.821 -35.750  96.399  132.150  0.145   2.669
2         outflow   1974   96.229  -5.391  18.059   23.451  0.049   0.715
         Variable  Count      Sum     Min     Max    Range   Mean  StdDev
0  water_sfc_elev   2454    7.761  -0.187   0.352    0.540  0.003   0.023
1          inflow   4378  126.551 -91.830  25.285  117.115  0.029   2.042
2         outflow   2451  125.190 -10.628  11.860   22.487  0.051   0.734

gwout

  Variable    Count       Sum     Min     Max   Range   Mean  StdDev
0   inflow  1621188   -77.624  -2.356   0.899   3.255 -0.000   0.006
1  outflow  1522144   -47.303  -0.869   0.826   1.695 -0.000   0.004
2    depth  1438037  2332.733 -39.692  37.363  77.055  0.002   0.245
  Variable    Count       Sum     Min     Max    Range   Mean  StdDev
0   inflow  1798552   -74.190  -1.314   0.850    2.165 -0.000   0.006
1  outflow  1736941   -48.611  -1.214   0.720    1.934 -0.000   0.005
2    depth  1649763  7522.898 -62.541  72.505  135.046  0.005   0.405

rtout

       Variable      Count           Sum  ...      Range   Mean  StdDev
0    zwattablrt    1278878  2.086635e+05  ...      4.000  0.163   0.691
1  sfcheadsubrt   52420025  1.589939e+07  ...  10162.759  0.303   3.724
2    QSTRMVOLRT    1735801  1.117775e+07  ...  14392.184  6.440  42.912
3       QBDRYRT       1741 -1.030273e+03  ...    103.361 -0.592   3.828
4        SOIL_M  315381587 -7.505986e+04  ...      0.842 -0.000   0.004

[5 rows x 8 columns]
       Variable      Count           Sum  ...      Range   Mean  StdDev
0    zwattablrt    1801619  1.536477e+05  ...      4.000  0.085   0.704
1  sfcheadsubrt   53599189  2.069484e+07  ...  13941.218  0.386   4.426
2    QSTRMVOLRT    2439924  2.155969e+07  ...  22221.239  8.836  60.114
3       QBDRYRT       2529 -6.827520e+03  ...    261.449 -2.700  11.295
4        SOIL_M  395619026 -1.129229e+05  ...      0.887 -0.000   0.004

[5 rows x 8 columns]

ldasout

     Variable    Count         Sum  ...      Range   Mean  StdDev
0       EMISS   269254      -0.014  ...      0.008 -0.000   0.000
1         FSA  1520951   -3752.242  ...     31.705 -0.002   0.069
2        FIRA  5603488   50444.655  ...     86.280  0.009   0.205
3      GRDFLX  5936729   97426.942  ...    434.096  0.016   1.027
4         HFX  5901139  113666.683  ...    465.800  0.019   1.278
..        ...      ...         ...  ...        ...    ...     ...
68    SOILSAT  5639699   -1903.321  ...      0.333 -0.000   0.003
69  SNOWT_AVG   463026 -184897.705  ...  20544.320 -0.399  97.832
70     ALBSND   725227      62.796  ...      1.872  0.000   0.010
71     ALBSNI   741857      61.663  ...      1.866  0.000   0.010
72      QRAIN   562330      -0.108  ...      0.000 -0.000   0.000

[73 rows x 8 columns]
     Variable    Count         Sum  ...      Range   Mean   StdDev
0        FVEG        1       0.692  ...      0.000  0.692    0.000
1       EMISS   416207      -0.505  ...      0.009 -0.000    0.000
2         FSA  2411159   -3288.130  ...     59.891 -0.001    0.087
3        FIRA  6917686   27901.449  ...     93.884  0.004    0.235
4      GRDFLX  7361525   40654.541  ...    393.687  0.006    1.050
..        ...      ...         ...  ...        ...    ...      ...
69    SOILSAT  6894244   -2933.963  ...      0.372 -0.000    0.003
70  SNOWT_AVG   564281 -873129.034  ...  20544.320 -1.547  150.412
71     ALBSND  1353339      34.251  ...      1.883  0.000    0.011
72     ALBSNI  1382878      33.442  ...      1.883  0.000    0.011
73      QRAIN   559973      -0.070  ...      0.000 -0.000    0.000

[74 rows x 8 columns]

restart_hydro

None
        Variable     Count           Sum  ...      Range   Mean  StdDev
0           stc1   4155504  1.257013e+04  ...     12.661  0.003   0.055
1           smc1   6397832 -1.646223e+03  ...      0.356 -0.000   0.005
2         sh2ox1   6400260 -1.578064e+03  ...      0.356 -0.000   0.005
3           stc2   3409136 -1.621658e+03  ...      2.454 -0.000   0.020
4           smc2   5600881 -1.792645e+03  ...      0.256 -0.000   0.002
5         sh2ox2   5605120 -1.764950e+03  ...      0.256 -0.000   0.002
6           stc3   1848038 -7.737650e+02  ...      0.716 -0.000   0.006
7           smc3   4875043 -9.362590e+02  ...      0.162 -0.000   0.002
8         sh2ox3   4880039 -9.338050e+02  ...      0.162 -0.000   0.002
9           stc4    752234 -4.093400e+01  ...      0.642 -0.000   0.005
10          smc4   4066710 -4.143590e+02  ...      0.190 -0.000   0.001
11        sh2ox4   4067561 -4.144210e+02  ...      0.190 -0.000   0.001
12       infxsrt   4449981  1.107588e+06  ...    680.838  0.249   2.023
13      soldrain   4044056 -6.254760e+02  ...     10.262 -0.000   0.009
14     sfcheadrt   4640890  9.937117e+05  ...    630.702  0.214   1.747
15       QBDRYRT      1741 -1.030273e+03  ...    103.361 -0.592   3.828
16      infxswgt  49348306 -1.000000e-03  ...      2.000 -0.000   0.078
17  sfcheadsubrt  52420025  1.589939e+07  ...  10162.759  0.303   3.724
18      sh2owgt1  90281947 -5.000000e-03  ...      8.171 -0.000   0.010
19      sh2owgt2  79600390 -6.000000e-03  ...      6.020 -0.000   0.007
20      sh2owgt3  70132807  2.000000e-03  ...      4.229  0.000   0.006
21      sh2owgt4  60943401  3.000000e-03  ...      1.742  0.000   0.004
22    qstrmvolrt   1735801  1.117775e+07  ...  14392.184  6.440  42.912
23         hlink   1667818  1.816279e+04  ...     23.012  0.011   0.130
24        qlink1   1656905  1.315317e+05  ...    121.893  0.079   1.112
25        qlink2   1656905  1.315317e+05  ...    121.893  0.079   1.112
26         resht      1976  3.819000e+00  ...      0.461  0.002   0.017
27        qlakeo      1974  9.622900e+01  ...     23.451  0.049   0.715
28        qlakei      4102  5.928210e+02  ...    132.150  0.145   2.669
29    z_gwsubbas   1439211  2.333000e+00  ...      0.077  0.000   0.000

[30 rows x 8 columns]
        Variable      Count           Sum  ...      Range   Mean  StdDev
0           stc1    5273633 -1.789976e+04  ...     14.700 -0.003   0.083
1           smc1    7730887 -2.519125e+03  ...      0.423 -0.000   0.005
2         sh2ox1    7732982 -2.197212e+03  ...      0.423 -0.000   0.005
3           stc2    4727691 -3.878239e+04  ...      4.690 -0.008   0.073
4           smc2    6817537 -3.365201e+03  ...      0.276 -0.000   0.003
5         sh2ox2    6821352 -3.062172e+03  ...      0.276 -0.000   0.003
6           stc3    3299094 -5.026133e+03  ...      1.107 -0.002   0.016
7           smc3    6117464 -1.183650e+03  ...      0.165 -0.000   0.002
8         sh2ox3    6120757 -1.154704e+03  ...      0.165 -0.000   0.002
9           stc4    1633206 -1.267330e+02  ...      0.623 -0.000   0.007
10          smc4    5317469 -6.435400e+02  ...      0.185 -0.000   0.001
11        sh2ox4    5315771 -6.435920e+02  ...      0.185 -0.000   0.001
12       infxsrt    4538121  1.437196e+06  ...    841.718  0.317   2.188
13      soldrain    5307142 -1.171276e+03  ...      9.183 -0.000   0.009
14     sfcheadrt    4772567  1.293428e+06  ...    796.637  0.271   1.887
15       QBDRYRT       2529 -6.827520e+03  ...    261.449 -2.700  11.295
16      infxswgt   51394779 -1.000000e-03  ...      2.000 -0.000   0.087
17  sfcheadsubrt   53599189  2.069484e+07  ...  13941.218  0.386   4.426
18      sh2owgt1  111376732  1.000000e-03  ...      8.596  0.000   0.011
19      sh2owgt2   99388714  1.000000e-02  ...      7.548  0.000   0.009
20      sh2owgt3   89777373  5.000000e-03  ...      6.605  0.000   0.008
21      sh2owgt4   80933163  9.000000e-03  ...      1.774  0.000   0.006
22    qstrmvolrt    2439924  2.155969e+07  ...  22221.239  8.836  60.114
23         hlink    1854713  6.778192e+03  ...    142.764  0.004   0.143
24        qlink1    1842065  9.346863e+04  ...    350.740  0.051   1.237
25        qlink2    1842065  9.346863e+04  ...    350.740  0.051   1.237
26         resht       2454  7.761000e+00  ...      0.540  0.003   0.023
27        qlakeo       2451  1.251900e+02  ...     22.487  0.051   0.734
28        qlakei       4378  1.265510e+02  ...    117.115  0.029   2.042
29    z_gwsubbas    1651025  7.523000e+00  ...      0.135  0.000   0.000

[30 rows x 8 columns]

restart_lsm

None
     Variable     Count           Sum  ...      Range   Mean  StdDev
0      SOIL_T  10164912  1.013377e+04  ...     12.661  0.001   0.037
1      SNOW_T     17925 -8.464254e+03  ...    546.320 -0.472  20.499
2         SMC  20940466 -4.789486e+03  ...      0.356 -0.000   0.003
3        SH2O  20952980 -4.691241e+03  ...      0.356 -0.000   0.003
4       ZSNSO   6893592 -1.868000e+00  ...      0.180 -0.000   0.000
5       SNICE   2891041  1.537817e+03  ...     70.371  0.001   0.141
6       SNLIQ   3133684 -4.495490e+02  ...     14.307 -0.000   0.020
7       QSNOW   1057778 -0.000000e+00  ...      0.000 -0.000   0.000
8        FWET   2152468 -6.700677e+03  ...      2.000 -0.003   0.029
9      SNEQVO   1139430  1.445860e+03  ...     26.292  0.001   0.043
10        EAH   5386943 -1.856760e+06  ...   1763.032 -0.345   7.304
11        TAH   4426819  7.844991e+03  ...      6.554  0.002   0.038
12         CM   5031846  1.413500e+01  ...      0.037  0.000   0.000
13         CH   5064121  6.157000e+01  ...      0.126  0.000   0.000
14      ISNOW       101  3.100000e+01  ...      2.000  0.307   0.952
15     CANLIQ   2114952 -5.767530e+02  ...      0.936 -0.000   0.002
16     CANICE   1012229 -9.700000e-02  ...      0.004 -0.000   0.000
17      SNEQV   1187143  1.526421e+03  ...     14.353  0.001   0.035
18      SNOWH   1229689  1.134000e+00  ...      0.036  0.000   0.000
19         TV   4575224  3.851161e+03  ...      7.974  0.001   0.039
20         TG   4011593  1.721964e+04  ...     17.878  0.004   0.056
21     ACMELT   1436524 -1.205435e+03  ...     14.301 -0.001   0.032
22      TAUSS    650038  4.337000e+01  ...      5.283  0.000   0.012
23       QSFC   4564595  8.694000e+00  ...      0.015  0.000   0.000
24  SFCRUNOFF   7053628  1.890102e+07  ...  11459.825  2.680  32.308
25   UDRUNOFF   3555243 -8.689312e+03  ...    218.862 -0.002   0.166
26    ACCECAN   3225073  4.910385e+03  ...      2.462  0.002   0.009
27    ACCEDIR   5630095 -2.137223e+04  ...      7.177 -0.004   0.048
28   ACCETRAN   5223602  2.247180e+02  ...      0.528  0.000   0.002

[29 rows x 8 columns]
     Variable     Count           Sum  ...      Range   Mean  StdDev
0      SOIL_T  14933624 -6.183501e+04  ...     14.700 -0.004   0.065
1      SNOW_T     28106 -2.375703e+04  ...    546.320 -0.845  29.542
2         SMC  25983357 -7.711516e+03  ...      0.423 -0.000   0.003
3        SH2O  25990862 -7.057680e+03  ...      0.423 -0.000   0.003
4       ZSNSO   9121507 -2.704100e+01  ...      0.238 -0.000   0.000
5       SNICE   3814263  4.220058e+03  ...     88.866  0.001   0.262
6       SNLIQ   3800638 -2.780178e+03  ...     16.093 -0.001   0.035
7       QSNOW    794037 -0.000000e+00  ...      0.000 -0.000   0.000
8        FWET   1879438 -3.404921e+03  ...      2.000 -0.002   0.029
9      SNEQVO   1487216  1.886331e+03  ...     18.431  0.001   0.073
10        EAH   6638913 -1.384643e+06  ...   2447.127 -0.209  11.783
11        TAH   5451981  6.505352e+03  ...     21.279  0.001   0.046
12         CM   6343083  2.035200e+01  ...      0.039  0.000   0.000
13         CH   6388617  6.542600e+01  ...      0.142  0.000   0.000
14      ISNOW       329  8.700000e+01  ...      2.000  0.264   0.964
15     CANLIQ   1903710 -3.287160e+02  ...      0.866 -0.000   0.003
16     CANICE    794934 -1.548000e+00  ...      0.040 -0.000   0.000
17      SNEQV   1574861  2.277380e+03  ...     16.154  0.001   0.063
18      SNOWH   1606776  5.584000e+00  ...      0.041  0.000   0.000
19         TV   5632356  4.185322e+03  ...     24.145  0.001   0.046
20         TG   5194369  1.126166e+04  ...     16.410  0.002   0.062
21     ACMELT   1879018 -1.042148e+03  ...     16.207 -0.001   0.060
22      TAUSS   1160710  6.292700e+01  ...      4.368  0.000   0.014
23       QSFC   5826894  5.597000e+00  ...      0.024  0.000   0.000
24  SFCRUNOFF   8484484  3.575350e+07  ...  26910.438  4.214  59.588
25   UDRUNOFF   4990634 -3.233155e+04  ...    407.296 -0.006   0.306
26    ACCECAN   4615756  8.588254e+03  ...      6.112  0.002   0.012
27    ACCEDIR   7254991 -3.758447e+04  ...     16.946 -0.005   0.071
28   ACCETRAN   6782408  1.265810e+02  ...      1.212  0.000   0.003

[29 rows x 8 columns]

restart_nudging

          Variable   Count       Sum     Min     Max   Range   Mean  StdDev
0  model_discharge  485653  3007.212 -22.760  26.562  49.321  0.006   0.401
1            nudge    5838  -130.184 -26.226   4.903  31.130 -0.022   0.539
          Variable    Count        Sum  ...     Range   Mean  StdDev
0  model_discharge  1050361  13568.083  ...  1036.691  0.013   1.715
1            nudge     6122    176.569  ...   349.879  0.029   3.756

[2 rows x 8 columns]

@aubreyd @rcabell Results looks reasonable to me.

@rcabell rcabell merged commit 8a61c8f into NCAR:master Jun 9, 2021
jdmattern pushed a commit to jdmattern/wrf_hydro_nwm_public that referenced this pull request Jun 28, 2021
* Xinanjiang runoff scheme added
* spatial soil XAJ params added
* change NWM tests to use new runoff scheme

Co-authored-by: Katelyn FitzGerald <katelynw@ucar.edu>
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

Successfully merging this pull request may close these issues.

None yet

4 participants