How to set "No Observation" file - Problem with climate change scenario

799 views
Skip to first unread message

Carolina N

unread,
Jan 23, 2017, 11:50:50 AM1/23/17
to SWAT-CUP
Hello all,

I've calibrated my model monthly for a  baseperiod (1977-1987) and now I need to simulate climate change for 2020-2099.

I've run ArcSWAt again for my climate change scenario (2020-2099) and started the simulation in SWAT-CUP putting the parameters range calibrated in Par_inf.

As I read in the swat cup manual and in posts in this group. I've set "No observation" file as follow (I have already change thousand of times these sets but I still get an error, the 95_ppu_No_Obs plot is not created).


Extract_rch_No_Obs.def

SUFI2          : SWAT-CUP program: SUFI2, GLUE, ParaSol, PSO, MCMC
output.rch     : swat output file name
1              : number of variables to get (such as: discharge, sediment, ET etc.)
7               : variable column number(s) in the swat output file (as many as the number above)
R-FLOW       :Names of variables

27            : total number of subbasins in the project

1             : number of subbasins to get for the first variable
1            : subbasin numbers for the first variable. Write "All" if equal to the total number of HRUs



2020           : beginning year of simulation not including warm up period
2099           : end year of simulation

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


Extract_hru_No_Obs.def

SUFI2          : SWAT-CUP program: SUFI2, GLUE, ParaSol, PSO, MCMC
output.hru     : swat output file name
1              : number of variables to get (such as: discharge, sediment, ET etc.)
7         : variable column number(s) in the swat output file (as many as the number above)
R-FLOW     :variable names


27              : total number of subbasins in the project

308              : total number of HRUs in the project

1              : number of HRUs to get for the first variable
1            : HRU numbers for the first variable. Write "All" if equal to the total number of HRUs



2020           : beginning year of simulation not including warm up period
2099           : end year of simulation

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



Extract_sub_No_Obs.def

SUFI2          : SWAT-CUP program: SUFI2, GLUE, ParaSol, PSO, MCMC
output.sub     : swat output file name
1              : number of variables to get (such as: discharge, sediment, ET etc.)
7        : variable column number(s) in the swat output file (as many as the number above)
R-FLOW  : names of variables


27: total number of subbasins in the project

1             : number of subbasins to get for the first variable
1            : HRU numbers for the first variable. Write "All" if equal to the total number of HRUs



2020          : beginning year of simulation not including warm up period
2099           : end year of simulation

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

95ppu_No_Obs.def

SUFI2    :name of optimization program: SUFI2, Para_Sol, GLUE, PSO, MCMC
1       : number of variables to calculate 95PPU for without any observations

 :variable names, which should be the same as simulation file names:
R-FLOW.txt


960            : number of data points or simulation time steps

Observed Function (Do I need to set this file? And if I leave in blank the space for observed values? If I delet this part, the program doesn't run)

Observed

1     : number of observed variables
5     : Objective function type, 1=mult,2=sum,3=r2,4=chi2,5=NS,6=br2,7=ssqr,8=PBIAS,9=KGE,10=RSR,11=MNS
0.5   : min value of objective function threshold for the behavioral solutions
1     : if objective function is 11=MNS (modified NS),indicate the power, p.

FLOW_OUT_1   : this is the name of the variable and the subbasin number to be included in the objective function
1     : weight of the variable in the objective function
-1    : Dynamic flow separation. Not considered if -1. If 1, then values should be added in the forth column below after observations
-1    : constant flow separation, threshold value. (not considered if -1)
1     : if separation of signal is considered, this is weight of the smaller values in the objective function
1     : if separation of signal is considered, this is weight of the larger values in the objective function
10    : percentage of measurement error
1068   : number of data points for this variable as it follows below. First column is a sequential number from beginning
      : of the simulation, second column is variable name and date (format arbitrary), third column is variable value.



Thanks in advance

Carolina 

Abbaspour, Karim

unread,
Jan 23, 2017, 1:31:58 PM1/23/17
to swat...@googlegroups.com

No need to change 1000 times! Just do it once correctly. The format is very clear. If you are just interested in flow, then you just need the rch extraction file.

Look in the batch files  in the Executable Files, which let you decide what to run and what not to run.

 

In the SUFI2_extract.bat you determine from which file you want to extract.

In the SUFI2_Post.bat, if you have no observation, then only check 95ppu_NO_Obs.exe (use text view to rem others). No need to run the programs for when you have observation.

 

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.

Carolina N

unread,
Jan 23, 2017, 3:23:53 PM1/23/17
to SWAT-CUP
Hi Dr. Abbaspour,

Thank you for answering. I have forgotten of checking No obs in the executable files.

Now, I am with an error as follow in the images about 'std::bad_alloc' (Contact the application's support team for more information).

Could you help me?
error1.png
error2.png

Abbaspour, Karim

unread,
Jan 24, 2017, 10:13:14 AM1/24/17
to swat...@googlegroups.com

In SUFI2_post.bat remark (rem) everything except 95ppu_NO_Obs.exe.

Reply all
Reply to author
Forward
Message has been deleted
0 new messages