SAWT CUP SUFI 2 simulation stuck at "Execution successully completed variable #1"

[Vaishali Sharda]

Hello,
I am running SWAT_CUP Sufi 2 to calibrate my swat model for flow using 3 variables for 2006 to 2011 with 2004 and 2005 as warm up period. The run gets stuck at the figure attached (SWATCUP_error1.jpeg). I have also attached other files. Can someone please point out what could be going on?
Thanks
VAS

[Karim Abbaspour]

Make sure there is a space between the fields in observed.txt

Karim

--
You received this message because you are subscribed to the Google Groups "SWAT-CUP" group.
To unsubscribe from this group and stop receiving emails from it, send an email to swat-cup+u...@googlegroups.com.
For more options, visit https://groups.google.com/d/optout.

<file.cio>

<observed.txt>

<echo_extract_rch.txt>

<echo_goal_fn.txt>

<SWATCUP_error1.jpg>

[Vaishali Sharda]

Thanks Dr. Abbaspour.
I am sure there is space between the fields. Am I missing something that you can point to?
Best Regards
VAS

[Karim Abbaspour]

That error is usually caused by something wrong in the observed file, may be the timing is wrong. Pls look at the files in Echo, there may be a clue.

Karim

[Vaishali Sharda]

I tried but cannot find anything. Can you please look at the echo files attached in the first post. I will really appreciate it.
Best Regards
vas

[Karim Abbaspour]

in the echo_observed_rch you have

variable no. = 1 
     1     2     3     4     5     6     7     0     0     0     0     0     0     0     0     0     0     0

Don't you see anything wrong here with time array? There must be something wrong with your obsereved_rch.txt file.

Also, what does v__REVAPMN of -85 mean?

Karim 

To unsubscribe from this group and stop receiving emails from it, send an email to swat-cup+unsubscribe@googlegroups.com.

[Yuan Lifeng]

Hello, Vas,

     Have you solved this issue?  I meet the same problem when I run SWAT-CUP, I use 13 parameters and 11 observation to calibrate model. SWAT-CUP stuck at 2003 year. calibration period is 1995-2008, warm-up is 5-yr, beginning year is 1990, running period is 19 year. 

在 2018年2月15日星期四 UTC-6上午11:13:57，Vaishali Sharda写道：

[Vaishali Sharda]

Hi Yuan,
Not yet. Please see Dr. Abbaspour replies above. I am still trying to figure out what is wron with my observed_rch.txt.

[Vincent Wu]

My guess is your number of year (NBYR) in file.cio is wrong. Try to set it as 24 (19 yrs plus 5 yrs warm-up). Hope this helps.

On Friday, February 16, 2018 at 8:53:40 AM UTC-6, Yuan Lifeng wrote:

[Vaishali Sharda]

Change your IYR from 1990 to 1995

On Friday, February 16, 2018 at 8:53:40 AM UTC-6, Yuan Lifeng wrote:

[Karim Abbaspour]

11              : number of reaches (subbasins) to get for the first variable
1  7  16  22  50   51  61  82  17  24  71  : reach (subbasin) numbers for the first variable

These subbasins should be in increasing order everywhere, also in the observation files.

r__CN2.mgt    35   98

look in the manual to learn what r__ means.

Karim

--

[Yuan Lifeng]

Hello, Vincent,

     Thank you for your suggestion, if I set it as 24, it will stuck at 2005 year, so I change it back to 19 yrs, right now it stuck at 2003 year.  My calibration period is 1995-2008. Anyway, I will try and see what will happen.

    Thanks again

Lifeng

在 2018年2月16日星期五 UTC-6上午10:29:00，Vincent Wu写道：

[Yuan Lifeng]

Thank you, Vaishali. we do the same things right now, I am checking my Observed_rch.txt, it's hard to find something wrong there. I just export these data from excel file and keep all necessary space. but still stuck .

在 2018年2月16日星期五 UTC-6上午9:47:15，Vaishali Sharda写道：

[Yuan Lifeng]

Thanks a lot, Dr. Abbaspour, i will change CN2 parameter and order the number of subbasin. I appreciate your help much. 

Best,

Lifeng

在 2018年2月16日星期五 UTC-6上午10:40:52，Abbaspour写道：

To unsubscribe from this group and stop receiving emails from it, send an email to swat-cup+u...@googlegroups.com.

[Vaishali Sharda]

Dr. Abbaspour,
Thanks for taking time every time to respond. I have gone through the observed_rch.txt again and I do not understand why it is getting stuck at 7. I have created it in excel and copied to SWAT CUP. I did the same for monthly calibration and it worked fine.
I have changed the REVAPMN upper and lower limit and still getting the same error.
Please let me know if you have any suggestions.
Best Regards
VAS

[Karim Abbaspour]

Lets see the observed_rch file. What do you mean stuck at 7?

Karim

[Yuan Lifeng]

Dr. Abbaspour, 

     Thank you very much. CUP run successfully right now.  

Goal_type= Nash_Sutcliff    No_sims= 1     Best_sim_no= 1    Best_goal = 8.965662e-002

Variable           p-factor  r-factor  R2    NS     bR2      MSE       SSQR       PBIAS  KGE  RSR   MNS   VOL_FR  ---  Mean_sim(Mean_obs)   StdDev_sim(StdDev_obs)

FLOW_OUT_1:        0.13      0.00      0.31  0.19   0.1477   1.6e-001  1.8e-002   -11.6 0.52  0.90  0.13  0.90            0.46(0.42)           0.38(0.44)

FLOW_OUT_7:        0.12      0.00      0.25  0.03   0.1130   1.8e+000  1.7e-001   -22.5 0.44  0.99  0.06  0.82            1.51(1.23)           1.26(1.37)

FLOW_OUT_16:       0.07      0.00      0.33  0.19   0.1681   8.8e+000  1.0e+000   -21.2 0.51  0.90  0.13  0.83            3.63(3.00)           2.95(3.29)

FLOW_OUT_17:       0.11      0.00      0.26  0.21   0.0997   4.5e-001  8.0e-002   -1.8  0.45  0.89  0.18  0.98            0.68(0.67)           0.56(0.76)

FLOW_OUT_22:       0.13      0.00      0.45  0.28   0.2901   2.4e+002  5.1e+001   -30.7 0.55  0.85  0.17  0.77           23.72(18.16)          17.50(18.31)

FLOW_OUT_24:       0.14      0.00      0.41  0.24   0.2483   2.1e+002  4.2e+001   -28.4 0.54  0.87  0.12  0.78           21.18(16.49)          15.81(16.78)

FLOW_OUT_50:       0.11      0.00      0.41  0.27   0.2326   5.9e+002  1.2e+002   -30.2 0.52  0.85  0.17  0.77           35.35(27.15)          25.16(28.53)

FLOW_OUT_51:       0.10      0.00      0.36  0.13   0.1935   1.1e+002  2.9e+001   -46.3 0.38  0.93  0.08  0.68           13.95(9.53)           9.99(11.31)

FLOW_OUT_61:       0.08      0.00      0.35  -0.08  0.1939   3.5e+000  1.2e+000   -72.3 0.17  1.04  -0.08 0.58            2.32(1.34)           1.72(1.79)

FLOW_OUT_71:       0.08      0.00      0.24  -0.29  0.1166   1.2e+001  2.9e+000   -58.7 0.22  1.13  -0.18 0.63            3.91(2.46)           3.07(3.00)

FLOW_OUT_82:       0.05      0.00      0.30  -0.18  0.1736   2.2e+003  5.6e+002   -49.1 0.33  1.09  -0.10 0.67           64.85(43.50)          44.32(42.70)

Goal_type= Nash_Sutcliff    No_sims= 1     Best_sim_no= 1    Best_goal = 8.965662e-002

 Parameter_Name           Fitted_Value     Min_value      Max_value

1:R__CN2.mgt            0.000000        -0.200000      0.200000

2:V__ALPHA_BF.gw        0.500000        0.000000       1.000000

3:V__GW_DELAY.gw        250.000000      0.000000       500.000000

4:V__GWQMN.gw           2500.000000     0.000000       5000.000000

5:V__GW_REVAP.gw        1.010000        0.020000       2.000000

6:R__ESCO.hru           0.505000        0.010000       1.000000

7:R__SOL_K(..).sol      1200.000000     900.000000     1500.000000

8:R__SOL_AWC(..).sol    0.500000        0.000000       1.000000

9:R__SOL_BD(..).sol     1.700000        0.900000       2.500000

10:V__CH_N2.rte         0.145000        -0.010000      0.300000

11:V__CH_K2.rte         249.995010      -0.010000      500.000000

12:R__CANMX.hru         50.000000       0.000000       100.000000

13:V__SURLAG.bsn        12.025001       0.050000       24.000000

    0.000000     0.500000   250.000000  2500.000000     1.010000     0.505000  1200.000000     0.500000     1.700000     0.145000   249.995010    50.000000    12.025001

r__CN2.mgt                    0.000000

v__ALPHA_BF.gw                0.500000

v__GW_DELAY.gw                250.000000

v__GWQMN.gw                   2500.000000

v__GW_REVAP.gw                1.010000

r__ESCO.hru                   0.505000

r__SOL_K().sol                1200.000000

r__SOL_AWC().sol              0.500000

r__SOL_BD().sol               1.700000

v__CH_N2.rte                  0.145000

v__CH_K2.rte                  249.995010

r__CANMX.hru                  50.000000

v__SURLAG.bsn                 12.025001

在 2018年2月16日星期五 UTC-6上午10:40:52，Abbaspour写道：

11              : number of reaches (subbasins) to get for the first variable

To unsubscribe from this group and stop receiving emails from it, send an email to swat-cup+u...@googlegroups.com.

[Vaishali Sharda]

Thanks. The files are attached.
Best Regards
VAS

[Yuan Lifeng]

Hello, Vas, 

      I think you should use 2008 as first year output in observed_rch txt, not 2006.  2006 is beginning year of simulation, warm-up is 2 years(2006,2007), so beginning year is 2008 in observed_rch.txt, end year is 2011. 

      Best,

      Lifeng

在 2018年2月16日星期五 UTC-6上午11:31:41，Vaishali Sharda写道：

[Vaishali Sharda]

Hi Lifeng,
2004-2005 is 2 year warm-up and simulation (SWAT-CUP) starts at 2006. This is what I have in mind.
Thanks
Vaishali.

[Karim Abbaspour]

               8    | NBYR : Number of years simulated
            2006    | IYR : Beginning year of simulation

You skip two years. SWAT does not give output for 2006 and 2007. Therefore you cannot start the observed files from 2006.

Karim

To unsubscribe from this group and stop receiving emails from it, send an email to swat-cup+unsubscribe@googlegroups.com.

[Yuan Lifeng]

Hello, Vaishali,

      Why do you write 2006 for IYR in your file.cio as below：

General Information/Watershed Configuration:

fig.fig

               8    | NBYR : Number of years simulated

            2006    | IYR : Beginning year of simulation

      I think it should be : 2004 | IYR : Beginning year of simulation.   Good luck. 

 Lifeng

在 2018年2月16日星期五 UTC-6上午11:54:12，Vaishali Sharda写道：

[Vaishali Sharda]

Thanks Lifeng and Dr. Abbaspour.
I changed the IYR to 2004. But still the same problem.
VAS

[Yuan Lifeng]

"total reaches in the project= 15" is there really 15 subbasins in your project? 

you'd better to enter the number of subbasin. How many subbasin in your project? you can get it from watershed shp file. Good luck.

Best

Lifeng

在 2018年2月16日星期五 UTC-6下午12:04:32，Vaishali Sharda写道：

[Vaishali Sharda]

Yes, That is right. 15 reaches.
I have already completed a project with calibrating this same watershed at a monthly timescale. Now working on a daily time step and that is where I am getting this error.
Thanks
VAS