SUFI2_run.bat and SUFI2_post.bat run Error

[gebreme...@gmail.com]

Dear SWAT-CUP Users
I'm using SWATCUP 5.1.6.2 version for SWAT 2012 (ArcSWAT 2012.10_1.15) model and I am using in 32 bit processor. Arc SWAT was run successfully and didn't have any error and imported the files to database and saved the simulation so that it has displayed in the Scenarios "TablesOut" folder.
I am getting error in the SUFI2_run.bat and SUFI2_post.bat (forrt1: error <65> floating invalid) Even I left for almost 2 hour running the SUFI2_pre.batand it is still running.
I am trying to calibrate 16 parameters with 100 simulations in a daily basis. The parameters are:

r__CN2.mgt                                    35  ---   98.0
v__ALPHA_BF.gw                          0.0 --- 1.0
v__GW_DELAY.gw                        0.0---  500
v__REVAPMN.gw                          0.0---500
v__GW_REVAP.gw                        0.02----0.2
v__GWQMN.gw                            0.0---500
v__ESCO.hru                                  0.0---1.0
v__EPCO.hru                                 0.0---1.0
v__CANMX.hru                            0.0---10.0
r__SLSUBBSN.hru                       10.0---150.0
v__CH_N2.rte                             -0.01---0.3
v__CH_K2.rte                               0.0---150.0
r__SOL_AWC (1).sol                   0.0--1.0
r__SOL_K (1).sol                         0.0---2000
r__SOL_Z (1).sol                         0.0---3500
v__SURLAG.bsn                         1.0--

I want you also to make me clear which files are going to copy from the ArcSWAT project directory in to the SWAT-CUP directory so that I will have to achieve the link between the input/output of a calibration program and the model. (Because some doubts strike to me is which ArcSWAT files are to be linked in SWAT CUP project which is used for calibration). I herewith Attached you the Screen shots. I look forward your valuable suggestion

Best regards
Gebremedhin

[Abbaspour, Karim]

… please look at the swatcup manual at least once before running the calibration program. You should at least know what r__ or v__ mean or what the swat parameters mean.

Best, 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.

[zelale...@gmail.com]

Dear friend! I do have the same problem with SWAT-CUP. So, if your problem is resolved would you share your success stories. 

...

[Eeshan Kumar]

Hi,

You have specified wrong values for parameters. You may want to refer to the parameterization in swat cup section pg 41 in swat cup manual and that should solve the problem.

Best wishes,

Eeshan

[zelale...@gmail.com]

Dear Eshan Kumar,
Based on your comment I have checked SWAT-CUP manual pp 41. But, on this page there is only description of cods, there is no any value of parameter. So, would you please make more clear by mentioning the wrong value taken by him and indicating the correct value, so that we can identify or locate where our problems lays. Thank you very much in advance for your kind help. 
Regars,

[Eeshan Kumar]

For example: "r__CN2.mgt                                    35  ---   98.0" does not makes sense. "r" means that the existing parameter value is multiplied with by (1 + a given value) and if the CN2 is already say 70 then after multiplying it may reach to an unrealistic value. You may also refer to SWAT theoretical documentation to see what value a parameter may take.

Eeshan

[zelale...@gmail.com]

Dear Eeshan, Thank you very much for your valuable comment. Actually, the min and max value of all parameter is available with in the SWAT-CUP under par.inf.txt - Absolute-SWAT-Variables.txt. The min and max value of CN2 also taken from this data base, so where is the mistake lays, on top of this there is a warring in the same place which say that " If you choose "v" for spacial parameter such as soil parameter and CN2 then you will lose your spacial variable".  On top of this the calibration inputs- par.inf.txt needs min. and max. value of each parameter, which is already specified under the folder Absolute-SWAT- Variables.txt. For this reason I am confused with the comment given.Hoping to hear from you soon and also thank you in advance for your kind help. 

[gebreme...@gmail.com]

Dear Dr. Karim

I have read the user manual and changed the parameters range using the ranges from the Absolute_SWAT_Value.txt files and for the parameters specification I have used the r__for the spatial and v__for the non-spatial parameters. Besides, the SOL_K (1).sol, SOL_AWC (1).sol, SOL_Z (1).sol have changed to SOL_K ().sol, SOL_AWC ().sol, SOL_Z ().sol respectively for all the layers of all the HRU. The SUFI2_pre.bat and SUFI2_run.bat finish the process. However there is still getting error in the SUFI2_post.bat I herewith write the error displayed and attached the screen shot.

E:\ArcSWAT_1973\test1973.Sufi2.SwatCup>SUFI2_goal_fn.exe

Could Not Find E:\ArcSWAT_1973\test1973.Sufi2.SwatCup\SUFI2.OUT\beh_*.*

R²=0.022312

! ERROR Simulation number 2 in variable file 1 is missing or has a problem

 

E:\ArcSWAT_1973\test1973.Sufi2.SwatCup>SUFI2_new_pars.exe

Reading goal file

Simulation no 3 is missing from the goal.txt file

E:\ArcSWAT_1973\test1973.Sufi2.SwatCup>SUFI2_95ppu.exe

Finished reading observation file

Reading FLOW_OUT_1

Error! Simulation 1 has problems

 

E:\ArcSWAT_1973\test1973.Sufi2.SwatCup>SUFI2_95ppu_beh.exe

Cannot open input file /SUFI2.OUT/no_beh_sims.txt

THIS IS NOT AN ERROR!!

There were no solutions with the objective function value >=the behavioral threshold

 

E:\ArcSWAT_1973\test1973.Sufi2.SwatCup>rem 95ppu_NO_Obs.exe

 

Here are the Parameters range taken from Absolute_SWAT_Value.txt

r__CN2.mgt        35  98

v__ALPHA_BF.gw    0    1

v__GW_DELAY.gw    0   500

v__REVAPMN.gw     0   500

v__GW_REVAP.gw    0.02 0.2

v__GWQMN.gw       0   5000

v__ESCO.hru       0   1

v__EPCO.hru       0   1

v__CANMX.hru      0  100

r__SLSUBBSN.hru   10  150

v__CH_N2.rte     -0.01  0.3

v__CH_K2.rte     -0.01  500

r__SOL_AWC().sol  0    1

r__SOL_K().sol    0    2000

r__SOL_Z().sol    0    3500

v__SURLAG.bsn     0.05  24

Thanking in advance for your valuable corrections and kind help to fix my problem.

Best, Gebremedhin

...

[Abbaspour, Karim]

When you use r__ it means a relative change of a few percentages. Therefore you should not use a range of 35 to 98. It should be something around +/- 2 to 3 percents, meaning -0.2  to  +0.2

But this may not be your problem either. You have to look at the timing carefully, the time swat simulates as stated in file.cio and the time you specify in the extract file. They should match.

Best, Karim

--

[gebreme...@gmail.com]

Dear Dr. Karim,
Thank you very much for your fast response and kind help! I understand the range now. When I have changed the ranges  I am still getting the error Massage.

E:\ArcSWAT_1973\test1973.Sufi2.SwatCup>SUFI2_goal_fn.exe

Could Not Find E:\ArcSWAT_1973\test1973.Sufi2.SwatCup\SUFI2.OUT\beh_*.*

R²=0.022312

! ERROR Simulation number 2 in variable file 1 is missing or has a problem

Is timing mean the time series of the observed values? If so I have checked and it is okay!
Herewith I have attached the file.cio and the SUFI2_extract.rch file to have a look.

Thanking you,

Best, Gebremedhin

On Wednesday, January 28, 2015 at 4:41:56 PM UTC+3, gebreme...@gmail.com wrote:

...

[gebreme...@gmail.com]

Dear Dr. Karim,
I thank you very much for your endless support! I have solved the problems running the SUFI2_post.bat with your great help but when I run 16 parameters with 100 simulations I got R2=0.13 and NSE= -0.03. I need your valuable suggestion

Thanking you
Gebremedhin

...

[gebreme...@gmail.com]

Dear Dr. Karim,
After second iteration with 500 simulations having the parameters of second iteration below I have got R2 and NSE<0.5. i need your help.
r__CN2.mgt                0.245786     0.399240
v__ALPHA_BF.gw      -0.005631     0.195599
v__GW_DELAY.gw       349.815155   524.267395
v__REVAPMN.gw        200.659470   254.622787
v__GW_REVAP.gw       0.035327     0.072417
v__GWQMN.gw          1966.364624  2694.270264
v__ESCO.hru              0.648635     0.739781
r__SLSUBBSN.hru      0.642318     0.773152
v__CH_N2.rte             0.118820     0.183924
v__CH_K2.rte             109.668709   133.195999
r__SOL_AWC().sol    -0.247897     0.032307
r__SOL_K().sol          0.371455     0.774718
v__SURLAG.bsn        7.801311    14.019550

Goal_type= R2    No_sims= 30     Best_sim_no= 16    Best_goal = 2.206211e-001

Variable        p-factor  r-factor    R2      NS       bR2       MSE      SSQR    PBIAS     KGE      RSR    VOL_FR  ---  Mean_sim(Mean_obs)      StdDev_sim(StdDev_obs)
FLOW_OUT_1         0.78     0.46     0.22     0.19     0.0604   1.8e+003  4.5e+002 34.5     0.24     0.90   1.53        13.24(20.21)               27.12(46.50)

Best, Gebremedhin
On Wednesday, January 28, 2015 at 4:41:56 PM UTC+3, gebreme...@gmail.com wrote:

...

[Carlos Montellano]

Hi,

Would you please share how to solve this error that comes when run Post_bat. 

All ranges in my files and timings for my 490 observation points (discontinuous observations) are Ok. Just to track where the error is, I tried to run the project with other observation dataset (actually it's the same to the current one but with less observed points, just 455) and works good with this file.

 

I know SWAT CUP accepts a wide range of observations points, therefore I'm not sure what I did wrong...Below my file.cio and SUFI_extract_rch.def:

- File.cio:

General Information/Watershed Configuration:

fig.fig

              13    | NBYR : Number of years simulated

            1992    | IYR : Beginning year of simulation

               1    | IDAF : Beginning julian day of simulation

             365    | IDAL : Ending julian day of simulation

............

............

Output Information:

               1    | IPRINT: print code (month, day, year)

               3    | NYSKIP: number of years to skip output printing/summarization

-SUFI2_extract_rch.def

output.rch     : swat output file name

1              : number of variables to get

6         : variable column number(s) in the swat output file (as many as the above number)

38              : total number of reaches (subbasins) in the project

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

36        : reach (subbasin) numbers for the first variable

1995           : beginning year of simulation not including the warm up period

2001           : end year of simulation

1              : time step (1=daily, 2=monthly, 3=yearly)

Best,

Carlos

[Carlos Montellano]

Hi,

I found out the error from my previous post was because of a duplicate observation in my dataset. 

Best,

Carlos

[yeshiwas mekonnen]

dawetm...@gmail.com

Dear SWAT-CUP Users
I'm using SWATCUP 5.1.6.2 version for SWAT 2012 (ArcSWAT 2012.10_1.15) model and I am using in 64 bit processor. Arc SWAT was run successfully and didn't have any error and imported the files to database and saved the simulation so that it has displayed in the Scenarios "TablesOut" folder.
I am getting error in the SUFI2_run. GIS   MON     AREAkm2  FLOW_INcms FLOW_OUTcms     EVAPcms    TLOSScms  SED_INtons SED_OUTtonsSEDCONCmg/kg   ORGN_INkg  ORGN_OUTkg   ORGP_INkg  ORGP_OUTkg    NO3_INkg   NO3_OUTkg    NH4_INkg   NH4_OUTkg    NO2_INkg   NO2_OUTkg   MINP_INkg  MINP_OUTkg   CHLA_INkg  CHLA_OUTkg   CBOD_INkg  CBOD_OUTkg  DISOX_INkg DISOX_OUTkg SOLPST_INmgSOLPST_OUTmg SORPST_INmgSORPST_OUTmg  REACTPSTmg    VOLPSTmg  SETTLPSTmgRESUSP_PSTmgDIFFUSEPSTmgREACBEDPSTmg   BURYPSTmg   BED_PSTmg BACTP_OUTctBACTLP_OUTct  CMETAL#1kg  CMETAL#2kg  CMETAL#3kg     TOT Nkg     TOT Pkg NO3ConcMg/l

I want you also to make me clear which files are going to copy from the ArcSWAT project directory in to the SWAT-CUP directory so that I will have to achieve the link . I herewith Attached you the Screen shots. I look forward your valuable suggestion

Best regards
yeshiwas mekonnen