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

Examples: Running CBP with hspf and hsp2 #2

Open
rburghol opened this issue Sep 22, 2022 · 9 comments
Open

Examples: Running CBP with hspf and hsp2 #2

rburghol opened this issue Sep 22, 2022 · 9 comments

Comments

@rburghol
Copy link
Contributor

rburghol commented Sep 22, 2022

Running an hsp2 or hspf scenario

See also: how to set up a scenario HARPgroup/HARParchive#288 and HARPgroup/cbp_wsm#58

Overview

  • We run a simulation that uses hspf for the land simulation and uses hsp2 for the river
  • To get a list of land segments for a given river segment:
    • cbp get_landsegs PL3_5250_0001
  • If you see a message like ../config/seglists/PL3_5250_0001.land no such file or directory you should run the "basingen" command:
    • cbp basingen.csh vadeq_2021 5250

Check if Land is already Run

  • scenario=vadeq_2021
  • basin=PL3_5250_0001
  • cbp check_land $scenario $basin
  • Note: A river scenario may use a land scenario with a different scenario name - in fact, that is very common for our application where meteorology may change less often than scenarios. So, you must know the land scenario of your river scenario.\
    • cbp get_config [river scenario] river "LAND SCENARIO"
  • See also: documentation for check_met and check_river, check_met, check_prad and check_land: New Helper Functions cbp_wsm#84 give same guidance with regard to the scenario, i.e. MET config for land will tell you the met scenario to check for runoff, and the MET config for the river will tell you the met scenario used as direct precip/evap into river channels and impoundments.

Code 1: Using get_config to determine land scenario for river scenario.

cbp get_config $scenario river "LAND SCENARIO"
Loading /opt/model/p6/vadeq/hspf.config
Searching for LAND SCENARIO END LAND SCENARIO in /opt/model/p6/vadeq/config/control/river/vadeq_2021.con
all vadeq_2021

# scenario hsp2_2022
cbp get_config $scenario river "LAND SCENARIO"
# hsp2_2022

Check Meteorology/PRAD Status

  • Get the met dataset for the land (or river) simulation of interest
    • met_scen=`cbp get_config $scenario land MET`
  • Check the met WDM status:
    • cbp check_met $met_scen $basin
  • Get the prad dataset for the land (or river) simulation of interest
    • prad_scen=`cbp get_config $scenario land "PRECIP ATMOS DEPOSITION"`
  • Check the met WDM status:
    • cbp check_prad $prad_scen $basin

Code 2: Checking Land and river met scenario status.

cbp check_met $met_scen $basin
Looking for land segments: N51059 N51061 N51107 N51153 N51600 N51683 N51685`
Missing 0 land-segment met: 

Set Up Model Environment

MODEL_ROOT=/opt/model/p6/vadeq
META_MODEL_ROOT=/opt/model/meta_model
IGNORE_PROBLEMS=0 # set this to 1 if you want to plow forward regardless of errors, not advisable but sometimes...
export MODEL_ROOT META_MODEL_ROOT IGNORE_PROBLEMS

Prepare Meteorology

Note: This is only done in the early spring each year, and during drought re-analysis, or if an error was found in old downloaded met data. Otherwise, this condition can be skipped.

segs=`cbp get_landsegs $basin`
for i in $segs; do
  #only need to create a WDM for each land seg
  # this script reads the file /etc/hspf.config to get directories.
  wdm_pm_one $i 1984010100 2020123123 nldas2 harp2021 $met_scen $prad_scen
done

Run Land Simulation

Run Land with HSPF (fast)

# now use new sbatch technique to run land simultaneously
segs=`cbp get_landsegs $basin`
location="p6_vadeq" # crucial to avoid collisions if running multiple installs
read -r scope land_scenario < <(cbp get_config $scenario river "LAND SCENARIO")
for i in $segs; do
  job_name="${location}_${scenario}_${i}"
  echo "sbatch --job-name=\"$job_name\" /opt/model/meta_model/run_model hspf_cbp6 $land_scenario $i auto land "
  sbatch --job-name="$job_name" /opt/model/meta_model/run_model hspf_cbp6 $land_scenario $i auto land 
done

Run Land with HSP2 (slow, but detailed output)

Phase 6
  • Note: the first 2 steps to set up the scenario only need to be done once
# (do once) Created duplicate met scenario for hsp2 testing: 
cbp copy_scenario hsp2_2022 hsp2_met 
# (do once) Set up land directories: 
cbp make_land_directories.csh hsp2_met 

river="JL1_7170_6800"
scenario="hsp2_met"
segs=`cbp get_landsegs $river`
for i in $segs; do
  echo "sbatch /opt/model/meta_model/run_model $model $scenario $i auto land "
  sbatch /opt/model/meta_model/run_model $model $scenario $i auto land
done

Phase 5
i=A51081
model=hsp2_cbp5
scenario=p532sova_2021
/opt/model/meta_model/run_model $model $scenario $i auto land init

Run a Full Basin River Simulation

Manual Method

segs=`cbp get_riversegs $basin`
for i in $segs; do
  /opt/model/meta_model/run_model hsp2_cbp6 $scenario $i auto river 
done

Run with mm_run_basin

  • Wraps the steps in a backgroundable nohup-able process.
nohup ./run/mm/mm_run_basin $scenario $basin &

Run Single river with auto temp dir

  • land from another scenario
  • see hsp2 results: om-model-info/6893975
  • see hspf results: om-model-info/6894028
MODEL_ROOT=/opt/model/p6/vadeq
META_MODEL_ROOT=/opt/model/meta_model
export MODEL_ROOT META_MODEL_ROOT
# run in hsp2
/opt/model/meta_model/run_model hsp2_cbp6 hsp2_2022 PS2_5550_5560 auto river 
# run in hspf
/opt/model/meta_model/run_model hspf_cbp6 vadeq_2021 PS2_5550_5560 auto river


Run a headwater with specific temp dir

  • note: Uses runoff data from another, previously run scenario, so only need to run the river module
/opt/model/meta_model/run_model hsp2_cbp6 hsp2_2022 JL1_7080_7190 /opt/model/p6/vadeq/tmp/rob-scratch/h2_test river run

Run the next downstream segment:

/opt/model/meta_model/run_model hsp2_cbp6 hsp2_2022 JL1_7190_7250 /opt/model/p6/vadeq/tmp/rob-scratch/h2_test river`

Run a small basin

  • note: Uses runoff data from another, previously run scenario, so only need to run the river module
segs=`cbp get_riversegs JL2_7250_7090`
for i in $segs; do
  /opt/model/meta_model/run_model hsp2_cbp6 hsp2_2022 $i /opt/model/p6/vadeq/tmp/rob-scratch/h2_test river
done

Run a basin with slurm dependency

@rburghol
Copy link
Contributor Author

rburghol commented Sep 23, 2022

Run hspf version of p6

Prep met

Run a test in the non-tidal James, Rockfish

# generate phase 6 meteorology WDMs
cd /opt/model/p6/vadeq/
# do that subwatershed
segs=`cbp get_landsegs JL2_7250_7090`
for i in $segs; do
  #only need to create a WDM for each land seg
  # this script reads the file /etc/hspf.config to get directories.
  wdm_pm_one $i 1984010100 2020123123 nldas2 harp2021 nldas1221 p20211221
done

Run land segments with sbatch

# now use new sbatch technique to run land simultaneously
for i in $segs; do
  echo "sbatch /opt/model/meta_model/run_model hspf_cbp6 vadeq_2021 $i auto land "
  sbatch /opt/model/meta_model/run_model hspf_cbp6 vadeq_2021 $i auto land 
done

Run a single segment

/opt/model/meta_model/run_model  hspf_cbp6 vadeq_2021 JL1_7080_7190 auto river run 

Run all segments

segs=`cbp get_riversegs JL2_7250_7090`
for i in $segs; do
  /opt/model/meta_model/run_model hspf_cbp6 vadeq_2021 $i auto river
done

@rburghol
Copy link
Contributor Author

rburghol commented Sep 23, 2022

Run James headwaters

# generate phase 6 meteorology WDMs
cd /opt/model/p6/vadeq/
# do that subwatershed
segs=`cbp get_landsegs JL2_7250_7090`
for i in $segs; do
  #only need to create a WDM for each land seg
  # this script reads the file /etc/hspf.config to get directories.
  wdm_pm_one $i 1984010100 2020123123 nldas2 harp2021 nldas1221 p20211221
done

@rburghol
Copy link
Contributor Author

rburghol commented Oct 4, 2022

Running Occoquan

Headwater in Occoquan, PL2_4970_5250:

  • Check if land has already run successfully
    • cbp basingen.csh vadeq_2021 4970
    • cbp check_land vadeq_2021 PL2_4970_5250
      • Looking for land segments: N51059 N51061 N51107 N51153 N51600 N51683 N51685
      • Found 43 land-uses, and 7 land-segments
      • Missing land segment-landuse combos:
      • NONE
  • Since NONE is result of check_land, we can try to run this single segment (note: we do NOT need to run basingen.csh to run a single segment now)
    • /opt/model/meta_model/run_model hspf_cbp6 vadeq_2021 PL2_4970_5250 auto river

Run the whole basin in hspf

segs=`cbp get_riversegs PL3_5250_0001`
for i in $segs; do   
  /opt/model/meta_model/run_model hspf_cbp6 vadeq_2021 $i auto river
done

Test headwater in hsp2: /opt/model/meta_model/run_model hsp2_cbp6 hsp2_2022 PL2_4970_5250 auto river

Run the whole basin in hsp2:

segs=`cbp get_riversegs PL3_5250_0001`
for i in $segs; do   
  /opt/model/meta_model/run_model hsp2_cbp6 hsp2_2022 $i auto river
done

@rburghol
Copy link
Contributor Author

rburghol commented Jan 10, 2023

step by step to debug new subsheds:

MODEL_ROOT=/opt/model/p6/vadeq
META_MODEL_ROOT=/opt/model/meta_model
export MODEL_ROOT META_MODEL_ROOT

temp_dir=" /opt/model/p6/vadeq/tmp/scratch/11677"
/opt/model/meta_model/run_model hsp2_cbp6 subsheds PS2_5568_5560 $temp_dir river prep 00_init
/opt/model/meta_model/run_model hsp2_cbp6 subsheds PS2_5568_5560 $temp_dir river prep 01_run_rug
/opt/model/meta_model/run_model hsp2_cbp6 subsheds PS2_5568_5560 $temp_dir river prep 03_special_action
/opt/model/meta_model/run_model hsp2_cbp6 subsheds PS2_5568_5560 $temp_dir river prep 04_convert_hsp2
/opt/model/meta_model/run_model hsp2_cbp6 subsheds PS2_5568_5560 $temp_dir river prep 05_ann_load
/opt/model/meta_model/run_model hsp2_cbp6 subsheds PS2_5568_5560 $temp_dir river prep 06_ave_load
/opt/model/meta_model/run_model hsp2_cbp6 subsheds PS2_5568_5560 $temp_dir river prep 08_etm

/opt/model/meta_model/run_model hsp2_cbp6 subsheds PS2_5568_5560 $temp_dir river run 00_init
/opt/model/meta_model/run_model hsp2_cbp6 subsheds PS2_5568_5560 $temp_dir river run 02_confluence
/opt/model/meta_model/run_model hsp2_cbp6 subsheds PS2_5568_5560 $temp_dir river run 03_upstream
/opt/model/meta_model/run_model hsp2_cbp6 subsheds PS2_5568_5560 $temp_dir river run 04_ps_sep_div
/opt/model/meta_model/run_model hsp2_cbp6 subsheds PS2_5568_5560 $temp_dir river run 05_compare
/opt/model/meta_model/run_model hsp2_cbp6 subsheds PS2_5568_5560 $temp_dir river run 07_get_uci
/opt/model/meta_model/run_model hsp2_cbp6 subsheds PS2_5568_5560 $temp_dir river run 08_run_river
/opt/model/meta_model/run_model hsp2_cbp6 subsheds PS2_5568_5560 $temp_dir river run 09_hydr_for_sed
/opt/model/meta_model/run_model hsp2_cbp6 subsheds PS2_5568_5560 $temp_dir river run 10_hydr_to_wdm
/opt/model/meta_model/run_model hsp2_cbp6 subsheds PS2_5568_5560 $temp_dir river run 12_res_typeII
/opt/model/meta_model/run_model hsp2_cbp6 subsheds PS2_5568_5560 $temp_dir river run 13_scour_org
/opt/model/meta_model/run_model hsp2_cbp6 subsheds PS2_5568_5560 $temp_dir river run 99_tidy_up

@rburghol rburghol mentioned this issue Jan 26, 2023
@rburghol
Copy link
Contributor Author

MODEL_ROOT=/opt/model/p6/vadeq
META_MODEL_ROOT=/opt/model/meta_model
IGNORE_PROBLEMS=0 # set this to 1 if you want to plow forward regardless of errors, not advisable but sometimes...
export MODEL_ROOT META_MODEL_ROOT IGNORE_PROBLEMS

# check the land sim
cbp check_land hsp2_2022 JL2_6850_6890
# run land
segs=`cbp get_landsegs JL2_6850_6890`
land_scenario="hsp2_2022"
for i in $segs; do
  echo "sbatch /opt/model/meta_model/run_model hspf_cbp6 $land_scenario $i auto land "
  sbatch /opt/model/meta_model/run_model hspf_cbp6 $land_scenario $i auto land 
done

segs=`cbp get_riversegs JL2_6850_6890`
for i in $segs; do
  /opt/model/meta_model/run_model hsp2_cbp6 hsp2_2022 $i auto river 
done

@rburghol
Copy link
Contributor Author

rburghol commented Jan 26, 2023

  • James River land sim running as of 5:00pm 1/26/2023
cbp check_land hsp2_2022 JL7_7070_0001
segs=`cbp get_landsegs JL7_7070_0001`
# use read since LAND SCENARIO is in the form "all hsp2_2022"
read -r scope land_scenario < <(cbp get_config hsp2_2022 river "LAND SCENARIO")
for i in $segs; do
  echo "sbatch /opt/model/meta_model/run_model hspf_cbp6 $land_scenario $i auto land "
  sbatch /opt/model/meta_model/run_model hspf_cbp6 $land_scenario $i auto land 
done

segs=`cbp get_riversegs JL7_7070_0001`
for i in $segs; do
  /opt/model/meta_model/run_model hsp2_cbp6 hsp2_2022 $i auto river 
done


Debug a single watershed

. hspf_config # loads the MODEL_ROOT and META_MODEL_ROOT config now
basin="JU2_6602_6600" # JU1_6290_6590
segs=`cbp get_landsegs $basin`
cbp check_met nldas1221 $basin
read -r scope land_scenario < <(cbp get_config hsp2_2022 river "LAND SCENARIO")
for i in $segs; do
   echo "sbatch /opt/model/meta_model/run_model hspf_cbp6 $land_scenario $i auto land ";
   sbatch /opt/model/meta_model/run_model hspf_cbp6 $land_scenario $i auto land ; 
done

/opt/model/meta_model/run_model hsp2_cbp6 hsp2_2022 $basin auto river 

Do a group

basin="JU4_7330_7000"
cbp check_land nldas1221 $basin
segs=`cbp get_landsegs $basin`
for i in $segs; do
   echo "sbatch /opt/model/meta_model/run_model hspf_cbp6 $land_scenario $i auto land ";
   sbatch /opt/model/meta_model/run_model hspf_cbp6 $land_scenario $i auto land ; 
done


segs=`cbp get_riversegs $basin`
for i in $segs; do
  /opt/model/meta_model/run_model hsp2_cbp6 hsp2_2022 $i auto river 
done

@jdkleiner
Copy link
Member

jdkleiner commented Feb 1, 2023

Full CIA Run - Potomac

PM7_4820_0001

jdk@deq2:/opt/model/p6/vadeq$ cbp check_land hsp2_2022 PM7_4820_0001
   Loading /opt/model/p6/vadeq/hspf.config
   Loading /opt/model/p6/vadeq/hspf.config
   Loading /opt/model/p6/vadeq/hspf.config
   /opt/model/p6/vadeq/config/seglists/PM7_4820_0001.land: No such file or directory.
   Looking for land segments:
   Found 43 land-uses, and 0 land-segments
   Missing land segment-landuse combos:

   Land segments with at least one LU missing:
jdk@deq2:/opt/model/p6/vadeq$ cbp basingen.csh hsp2_2022 4820
   Loading /opt/model/p6/vadeq/hspf.config
jdk@deq2:/opt/model/p6/vadeq$ cbp check_land hsp2_2022 PM7_4820_0001
   Loading /opt/model/p6/vadeq/hspf.config
   Loading /opt/model/p6/vadeq/hspf.config
   Loading /opt/model/p6/vadeq/hspf.config
   Looking for land segments: N24031 N51061 N51107 N24013 N24021 H24021 N51069 N51171 N42001 N24043 N42055 H51015 N51015 N51820 H51165 N51165 N54031 N51187 N51043 N51840 N51660 N42041 N54003 N54037 N51139 L51015 N51790 N42099 N54065 N42057 H54031 N54027 N24001 N42009 H51139 L54071 N51091 N54071 H54023 L54023 N54023 N54057 H54071 L54031 N42111 N24023 H24023 N54077 N54093 H54057 N51059 N51600 N11001
   Found 43 land-uses, and 53 land-segments
   Missing land segment-landuse combos:
   aopN24031.wdm aopN24013.wdm aopN24021.wdm aopH24021.wdm aopN51069.wdm aopN51171.wdm aopN42001.wdm aopN24043.wdm aopN42055.wdm aopN51820.wdm aopH51165.wdm aopN51165.wdm...
   Land segments with at least one LU missing:
   N24031 N24013 N24021 H24021 N51069 N51171 N42001 N24043 N42055 N51820 H51165 N51165 N54031 N51187 N51043 N51840 N51660 N42041 N54003 N54037 N51139 L51015 N51790 N42099 N54065 N42057 H54031 N54027 N24001 N42009 H51139 L54071 N54071 H54023 L54023 N54023 N54057 H54071 L54031 N42111 N24023 H24023 N54077 N54093 H54057 N11001

^ This is saying the land simulation needs to be run

jdk@deq2:/opt/model/p6/vadeq$ cbp check_met nldas1221 PM7_4820_0001
   Loading /opt/model/p6/vadeq/hspf.config
   Loading /opt/model/p6/vadeq/hspf.config
   Loading /opt/model/p6/vadeq/hspf.config
   Looking for land segments: N24031 N51061 N51107 N24013 N24021 H24021 N51069 N51171 N42001 N24043 N42055 H51015 N51015 N51820 H51165 N51165 N54031 N51187 N51043 N51840 N51660 N42041 N54003 N54037 N51139 L51015 N51790 N42099 N54065 N42057 H54031 N54027 N24001 N42009 H51139 L54071 N51091 N54071 H54023 L54023 N54023 N54057 H54071 L54031 N42111 N24023 H24023 N54077 N54093 H54057 N51059 N51600 N11001
   Found 53 land-segments
   Found 53 land-segments with met WDMs
   Missing 0 land-segment met:
   Missing land segment-landuse combos:

^ This is saying the meteorology is there, just need to run the land simulation

MODEL_ROOT=/opt/model/p6/vadeq
META_MODEL_ROOT=/opt/model/meta_model
IGNORE_PROBLEMS=0 # set this to 1 if you want to plow forward regardless of errors, not advisable but sometimes...
export MODEL_ROOT META_MODEL_ROOT IGNORE_PROBLEMS

# now use new sbatch technique to run land simultaneously
segs=`cbp get_landsegs PM7_4820_0001`
# land_scenario="vadeq_2021"
read -r scope land_scenario < <(cbp get_config vadeq_2021 river "LAND SCENARIO")
for i in $segs; do
  echo "sbatch /opt/model/meta_model/run_model hspf_cbp6 $land_scenario $i auto land "
  sbatch /opt/model/meta_model/run_model hspf_cbp6 $land_scenario $i auto land 
done

# to monitor as its going
squeue

^ note land_scenario is just the scenario name, hspf_cbp6 is the model

Next step: run the river simulation:

MODEL_ROOT=/opt/model/p6/vadeq
META_MODEL_ROOT=/opt/model/meta_model
IGNORE_PROBLEMS=0 # set this to 1 if you want to plow forward regardless of errors, not advisable but sometimes...
export MODEL_ROOT META_MODEL_ROOT IGNORE_PROBLEMS

segs=`cbp get_riversegs PM7_4820_0001`
for i in $segs; do
  /opt/model/meta_model/run_model hsp2_cbp6 hsp2_2022 $i auto river 
done

Results

  • Example: om-model-info/6963628
  • Monitor progress during model run:
SELECT rseg.hydrocode AS "rseg.hydrocode",
	   model.pid AS "model.pid", 
	   model.propname AS "model.propname", 
	   model.propcode AS "model.propcode", 
	   scenario.pid AS "scenario.pid", 
	   scenario.propname AS "scenario.propname", 
	   scenario.propcode AS "scenario.propcode",
	   Qout.propvalue AS "Qout"
FROM dh_properties as model
LEFT OUTER JOIN dh_properties as scenario 
  ON (model.pid = scenario.featureid)
LEFT OUTER JOIN dh_properties as Qout 
  ON (scenario.pid = Qout.featureid)
LEFT OUTER JOIN dh_feature as rseg 
  ON (model.featureid = rseg.hydroid)   
WHERE model.propcode = 'cbp-6.0' AND scenario.propname = 'hsp2_2022' AND Qout.propname = 'Qout'
ORDER BY model.modified DESC
LIMIT 10;

Count scenarios run

SELECT COUNT(*)
FROM dh_properties as model
LEFT OUTER JOIN dh_properties as scenario 
  ON (model.pid = scenario.featureid)
LEFT OUTER JOIN dh_properties as Qout 
  ON (scenario.pid = Qout.featureid)
LEFT OUTER JOIN dh_feature as rseg 
  ON (model.featureid = rseg.hydroid)   
WHERE model.propcode = 'cbp-6.0' AND scenario.propname = 'hsp2_2022' AND Qout.propname = 'Qout' AND rseg.hydrocode LIKE '%P%';
 count
-------
    39
(1 row)

count number of segs:

jdk@deq2:/opt/model/p6/vadeq$ printf '%s\n' $segs | wc -w
126

@jdkleiner
Copy link
Member

jdkleiner commented Feb 15, 2023

Potomac Seg QA 2/15/23:

  • Headwater segment: Seneca Creek
    • vahydrosw_wshed_PU1_5520_5210
  • om-model-info/67800/dh_feature
jdk@deq2:/opt/model/p6/vadeq$

MODEL_ROOT=/opt/model/p6/vadeq
META_MODEL_ROOT=/opt/model/meta_model
IGNORE_PROBLEMS=0 # set this to 1 if you want to plow forward regardless of errors, not advisable but sometimes...
export MODEL_ROOT META_MODEL_ROOT IGNORE_PROBLEMS

basin="PU1_5520_5210"

cbp check_met nldas1221 $basin
cbp check_land hsp2_2022 $basin
# cbp basingen.csh hsp2_2022 5520

# run the land simulation 
segs=`cbp get_landsegs $basin`
read -r scope land_scenario < <(cbp get_config hsp2_2022 river "LAND SCENARIO")
for i in $segs; do
  echo "sbatch /opt/model/meta_model/run_model hspf_cbp6 $land_scenario $i auto land "
  sbatch /opt/model/meta_model/run_model hspf_cbp6 $land_scenario $i auto land 
done

# run the river simulation
segs=`cbp get_riversegs $basin`
for i in $segs; do
  /opt/model/meta_model/run_model hsp2_cbp6 hsp2_2022 $i auto river 
done

result:
Model_Map_vahydro-1 0runid_11_cbp-6 0hsp2_2022_potomac

@jdkleiner
Copy link
Member

jdkleiner commented Feb 28, 2023

Potomac Seg QA

2/28/23:

  • Headwater segment: South Branch Potomac
    • vahydrosw_wshed_PU0_6080_5620
cd /opt/model/p6/vadeq

# set up model environment
MODEL_ROOT=/opt/model/p6/vadeq
META_MODEL_ROOT=/opt/model/meta_model
IGNORE_PROBLEMS=0
export MODEL_ROOT META_MODEL_ROOT IGNORE_PROBLEMS

# input arguments
rseg=PU0_6080_5620
rmodel=hsp2_cbp6
rscenario=hsp2_2022
lmodel=hspf_cbp6
lscenario=vadeq_2021

# check meteorology status
metscenario=`cbp get_config $rscenario land MET`
cbp check_met $metscenario $rseg 
cbp check_land $rscenario $rseg

# run the land simulation (sbatch uses slurm to run land simultaneously)
lsegs=`cbp get_landsegs $rseg`
for i in $lsegs; do
  echo "sbatch $META_MODEL_ROOT/run_model $lmodel $lscenario $i auto land"
  sbatch $META_MODEL_ROOT/run_model $lmodel $lscenario $i auto land 
done

# run the river simulation
rsegs=`cbp get_riversegs $rseg`
for i in $rsegs; do
  echo "$META_MODEL_ROOT/run_model $rmodel $rscenario $i auto river"
  $META_MODEL_ROOT/run_model $rmodel $rscenario $i auto river 
done

results:

  • land sim completed in ~8.5mins
  • river sim ran into errors:
 ... ... ... SEDM         2019   0.00000000
 ... ... ... SLON         2020   0.00000000
 ... ... ... SEDM         2020   0.00000000
DONE sho for segment PU0_6080_5620 land scenario hsp2_2022

 PROBLEM FILE WRITTEN

Error: opening wdm=   -2
../../../tmp/wdm/land/sho/hsp2_2022/sho_L54071_PU0_6080_5620.wdm


Error: Problem found in step hsp2_cbp6 => river / prep / 05_ann_load
       See /opt/model/p6/vadeq/tmp/scratch/12796/problem

3/1/23:

  • revised setting of lscenario
    • read -r scope lscenario < <(cbp get_config $rscenario river "LAND SCENARIO")
  • had to rerun land simulation using the correct lscenario
  • land ran in ~8mins
  • river run completed no errors
  • d.dh/om-model-info/605598/dh_feature

3/6/23:

  • Ran another test seg to test HARP analysis script updates
  • Headwater segment: Seneca Creek
    • PU1_5520_5210
    • land ran in ~7:50mins
    • river run completed with errors
    • model container created: /d.dh/om-model-info/605610/dh_feature
  • glenn seg test: PL3_5250_0001
    • land ran in ~10mins
  • single test seg: PL2_4970_5250

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

No branches or pull requests

2 participants