HARPgroup / HARParchive

This repo houses HARP code development items, resources, and intermediate work products.
1 stars 0 forks source link

Update the hydrocode for all successfully created subwatersheds. #761

Open rburghol opened 1 year ago

rburghol commented 1 year ago

Info:

Tasks:

rburghol commented 1 year ago

Note. Some variation in FTABLE output, beyond what I expected, between the FTABLEs generated on January 4th, and ones I regenerated today. I think I am using the right script. When doing a diff, the script that I have in the model directory, and the script from HARP archive, are identical, with the exception that I added a few message statements for debugging and also re-added the model_version as an input to the script (don't think that broke it, but who knows?).

ls -lhrt 
-rw-rw-r-- 1 megpritch allmodelers 1021 Jan  4 18:09 /opt/model/p6/vadeq/input/param/river/vahydro_2022/ftables/YM2_6120_6430_ni_reservoir.ftable
-rw-rw-r-- 1 rob allmodelers 1021 Mar  8 20:30 /opt/model/p6/vadeq/input/param/river/vahydro_2022/ftables/YM2_6122_6120.ftable

The Tables are as follows: March 8th

FTABLE   6122
NOTE: FLOODPLAIN BASE = 5*BANKFULL WIDTH ***
      FLOODPLAIN SIDE-SLOPE = SAME AS CHANNEL'S ***
 ROWS COLS ***
   19    4
     DEPTH      AREA    VOLUME     DISCH  ***
      (FT)   (ACRES)   (AC-FT)     (CFS)  ***
     0.000     0.000      0.00      0.00
     0.598    18.779     10.68     13.89
     1.196    20.599     22.45     44.85
     1.793    22.420     35.31     89.83
     2.391    24.240     49.25    148.06
     2.989    26.060     64.29    219.36
     3.587    27.881     80.41    303.81
     4.184    29.701     97.62    401.63
     4.782    31.521    115.92    513.10
     5.380    33.341    135.30    638.58
     6.380   169.752    303.53    465.45
     8.272   175.515    630.23   1535.85
    10.165   181.278    967.83   3068.05
    12.057   187.041   1316.33   5009.80
    13.949   192.803   1675.75   7331.85
    15.842   198.566   2046.06  10015.68
    17.734   204.329   2427.29  13048.89
    19.627   210.092   2819.42  16422.96
    21.519   215.855   3222.45  20132.05

January 4th

FTABLE   6120
NOTE: FLOODPLAIN BASE = 5*BANKFULL WIDTH ***
      FLOODPLAIN SIDE-SLOPE = SAME AS CHANNEL'S ***
 ROWS COLS ***
   19    4
     DEPTH      AREA    VOLUME     DISCH  ***
      (FT)   (ACRES)   (AC-FT)     (CFS)  ***
     0.000     0.000      0.00      0.00
     0.471    70.254     31.51     31.96
     0.943    77.044     66.22    103.47
     1.414    83.834    104.13    207.73
     1.885    90.624    145.24    343.06
     2.356    97.414    189.55    509.11
     2.828   104.204    237.05    706.11
     3.299   110.994    287.76    934.60
     3.770   117.784    341.67   1195.30
     4.241   124.574    398.78   1489.04
     5.241   637.279   1028.86   1939.48
     6.707   658.395   1978.30   5639.79
     8.173   679.510   2958.69  10796.12
     9.638   700.625   3970.03  17260.05
    11.104   721.741   5012.31  24946.43
    12.569   742.856   6085.54  33800.83
    14.035   763.971   7189.71  43786.67
    15.500   785.087   8324.83  54878.92
    16.966   806.202   9490.90  67060.42
rburghol commented 1 year ago

Quick script to generate, and verify FTABLEs (uci should be todays date and 12k):


riverseg=PM7_4622_4620
cname="local_channel"

Rscript run/resegment/ftable_creation.R $riverseg "0. River Channel" vahydro-1.0 "/opt/model/p6/vadeq/input/param/river/vahydro_2022/ftables/"
# run the UCI maker
/opt/model/meta_model/run_model hsp2_cbp6 subsheds $riverseg tmp/scratch/rwb river prep 01_run_rug
ls -lhrt tmp/uci/river/subsheds/${riverseg}.uci