SUFI2 stopping at Variable 1...

252 views
Skip to first unread message

Bailey, David

unread,
Feb 10, 2015, 10:43:53 PM2/10/15
to swat...@googlegroups.com

Hello all:

I am simulating stream flow in a 1km2 watershed. I successfully simulated stream flow in ArcSWAT from 2005 to 2010 with a 2 year warm up so it’s actually from 2007 to 2010.  Now I am in SWAT CUP attempting calibration.

Every time I run calibration it stops at the following screen ( see attached image). Anyone have any idea what is wrong?

Here is more info:

Please help!

 

Master Watershed File: file.cio

Project Description:

General Input/Output section (file.cio):

2/10/2015 12:00:00 AM ARCGIS-SWAT interface AV

 

General Information/Watershed Configuration:

fig.fig

               6    | NBYR : Number of years simulated

            2005    | IYR : Beginning year of simulation

               1    | IDAF : Beginning julian day of simulation

             365    | IDAL : Ending julian day of simulation

Climate:

               0    | IGEN : Random number seed cycle code

               1    | PCPSIM : precipitation simulation code: 1=measured, 2=simulated

               0    | IDT : Rainfall data time step

               0    | IDIST : rainfall distribution code: 0 skewed, 1 exponential

           1.300    | REXP : Exponent for IDIST=1

               1    | NRGAGE: number of pcp files used in simulation

               1    | NRTOT: number of precip gage records used in simulation

               1    | NRGFIL: number of gage records in each pcp file

               1    | TMPSIM: temperature simulation code: 1=measured, 2=simulated

               1    | NTGAGE: number of tmp files used in simulation

               1    | NTTOT: number of temp gage records used in simulation

               1    | NTGFIL: number of gage records in each tmp file

               1    | SLRSIM : Solar radiation simulation Code: 1=measured, 2=simulated

               1    | NSTOT: number of solar radiation records in slr file

               1    | RHSIM : relative humidity simulation code: 1=measured, 2=simulated

               1    | NHTOT: number of relative humidity records in hmd file

               1    | WINDSIM : Windspeed simulation code: 1=measured, 2=simulated

               1    | NWTOT: number of wind speed records in wnd file

               0    | FCSTYR: beginning year of forecast period

               0    | FCSTDAY: beginning julian date of forecast period

               0    | FCSTCYCLES: number of time to simulate forecast period

Precipitation Files:

pcp1.pcp

 

 

Temperature Files:

tmp1.tmp

 

 

slr.slr             | SLRFILE: name of solar radiation file

hmd.hmd             | RHFILE: name of relative humidity file

wnd.wnd             | WNDFILE: name of wind speed file

cst.cst             | FCSTFILE: name of forecast data file

Watershed Modeling Options:

basins.bsn          | BSNFILE: name of basin input file

Database Files:

plant.dat            | PLANTDB: name of plant growth database file

till.dat            | TILLDB: name of tillage database file

pest.dat            | PESTDB: name of pesticide database file

fert.dat            | FERTDB: name of fertilizer database file

urban.dat           | URBANDB: name of urban database file

Special Projects:

               0    | ISPROJ: special project: 1=repeat simulation

               0    | ICLB: auto-calibration option: 0=no, 1=yes

                    | CALFILE: auto-calibration parameter file

Output Information:

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

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

               0    | ILOG: streamflow print code: 1=print log of streamflow

               0    | IPRP: print code for output.pst file: 1= print pesticide output

Starting Output Variable Section

Reach output variables:

   0   0   0   0   0   0   0   0   0   0   0   0   0   0   0   0   0   0   0   0

Subbasin output variables:

   0   0   0   0   0   0   0   0   0   0   0   0   0   0   0

HRU output variables:

   1   0   0   0   0   0   0   0   0   0   0   0   0   0   0   0   0   0   0   0

HRU data to be printed:

   1   0   0   0   0   0   0   0   0   0   0   0   0   0   0   0   0   0   0   0

ATMOSPERIC DEPOSITION

ATMO.ATM

               1    | IPHR: print code for hourly output 0=no 1=yes (hourq.out)

               1    | ISTO: print code for soil storage 0=no 1=yes (output.swr)

               0    | ISOL: Code for printing phosphorus/nitrogen in soil profile (output.snu)

               1    | I_SUBW: Code for routing headwaters

septwq.dat

               1    | IA_B: Code for binary output of files (.rch, .sub, .hru files only)

               0    | IHUMUS: Print watqual.out file 0=no 1=yes (output.wql)

               1    | ITEMP: 0=print no file(s) 1=print output.vel/output.dep

               1    | ISNOW: 0=do not print snowband.out; 1=print output.snw

               1    | IMGT: 0=do not print output.mtg; 1=print output.mgt

               1    | IWTR: Code for printing output.pot and output.wtr files

0         | ICALEN: Code for printing out calendar or julian dates to .rch, .sub and .hru files

 

18  : Number of Parameters (the program only reads the first 4 parameters or any number indicated here)

200  : number of simulations

r__CN2.mgt                 -0.2       0.2

v__ALPHA_BF.gw             0.0       1.0

v__GW_DELAY.gw                     30.0      450.0

v__GWQMN.gw                0.0       2.0

v__GW_REVAP.gw             0.0       0.2

v__ESCO.hru                0.8       1.0

v__CH_N2.rte              0.0       0.3

v__CH_K2.rte               5.0       130.0

v__ALPHA_BNK.rte           0.0       1.0

r__SOL_AWC(1).sol         -0.2       0.4

r__SOL_K(1).sol                         -0.8       0.8

r__SOL_BD(1).sol          -0.5       0.6

v__SFTMP.bsn              -5.0       5.0

r__CANMX.hru    0   100

r__BLAI{42}.plant.dat    0.5   10

r__TLAPS.sub    -10   10

r__TIMP.bsn    0   1

r__SURLAG.bsn    0.05   24

 

1     : number of observed variables

FLOW_OUT_17   : this is the name of the variable and the subbasin number to be included in the objective function

1461   : 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.

 

1              FLOW_OUT_001_2007   0.0023

2              FLOW_OUT_002_2007   0.0022

3              FLOW_OUT_003_2007   0.0023

4              FLOW_OUT_004_2007   0.0054

5              FLOW_OUT_005_2007   0.0037

6              FLOW_OUT_006_2007   0.0029

7              FLOW_OUT_007_2007   0.0027

8              FLOW_OUT_008_2007   0.0026

9              FLOW_OUT_009_2007   0.0025

10           FLOW_OUT_010_2007   0.0025

11           FLOW_OUT_011_2007   0.0024

12           FLOW_OUT_012_2007   0.0022

13           FLOW_OUT_013_2007   0.0022

14           FLOW_OUT_014_2007   0.0022

15           FLOW_OUT_015_2007   0.0022

16           FLOW_OUT_016_2007   0.0022

17           FLOW_OUT_017_2007   0.0021

……..

1443       FLOW_OUT_347_2010   0.0100

1444       FLOW_OUT_348_2010   0.0100

1445       FLOW_OUT_349_2010   0.0132

1446       FLOW_OUT_350_2010   0.0120

1447       FLOW_OUT_351_2010   0.0351

1448       FLOW_OUT_352_2010   0.1243

1449       FLOW_OUT_353_2010   0.3705

1450       FLOW_OUT_354_2010   0.1251

1451       FLOW_OUT_355_2010   0.0674

1452       FLOW_OUT_356_2010   0.0550

1453       FLOW_OUT_357_2010   0.0486

1454       FLOW_OUT_358_2010   0.0445

1455       FLOW_OUT_359_2010   0.0402

1456       FLOW_OUT_360_2010   0.0372

1457       FLOW_OUT_361_2010   0.0351

1458       FLOW_OUT_362_2010   0.0369

1459       FLOW_OUT_363_2010   0.0391

1460       FLOW_OUT_364_2010   0.0330

1461       FLOW_OUT_365_2010   0.0302


Thanks!

--
- David J. Bailey

SWATCUP_Stopping.JPG

Abbaspour, Karim

unread,
Feb 11, 2015, 12:43:53 AM2/11/15
to <swat-cup@googlegroups.com>
You tell file.cio you want to do monthly simulation, then you give daily data! 


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

--
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.
<SWATCUP_Stopping.JPG>

David Bailey

unread,
Feb 11, 2015, 2:11:56 AM2/11/15
to swat...@googlegroups.com
Dr. Abbaspour:

Thanks for catching that! It appears to be working now.

Thanks!

- David 

Sent from my iPhone
Reply all
Reply to author
Forward
0 new messages