#		NSL TECTUM MODEL - 8x8 COLUMN ARRAY

# network
nsl set system.model Tectum88Model

set t1 5.0

nsl set system.integration euler
nsl set system.simDelta 0.01
nsl set system.simEndTime $t1
nsl set system.diff.delta 0.01

nsl set Tectum88Model.tectum.gl_tm 2.34
nsl set Tectum88Model.tectum.lp_tm 0.3
nsl set Tectum88Model.tectum.sp_tm 0.9
nsl set Tectum88Model.tectum.sn_tm 1.625
nsl set Tectum88Model.tectum.py_tm 0.12
nsl set Tectum88Model.tectum.tp_tm 0.02

nsl set Tectum88Model.tectum.C 0.889

nsl set Tectum88Model.tectum.Wr2gl \
0.834 0.834 0.834 0.834 6.67 0.834 0.834 0.834 0.834
nsl set Tectum88Model.tectum.Wlpgl \
5.34 5.34 5.34 5.34 8.0 5.34 5.34 5.34 5.34
nsl set Tectum88Model.tectum.Wspgl \
0.67 0.67 0.67 0.67 0.67 0.67 0.67 0.67 0.67

nsl set Tectum88Model.tectum.Wr2lp \
0.125 0.125 0.125 0.125 1.0 0.125 0.125 0.125 0.125
nsl set Tectum88Model.tectum.Wgllp \
0 0 0 0 1 0 0 0 0
nsl set Tectum88Model.tectum.Wsplp \
0.6 0.6 0.6 0.6 0.8 0.6 0.6 0.6 0.6
nsl set Tectum88Model.tectum.Wsnlp \
8.2 8.2 8.2 8.2 8.0 8.2 8.2 8.2 8.2
nsl set Tectum88Model.tectum.Wtplp \
0 0 0 0 0.1 0 0 0 0

nsl set Tectum88Model.tectum.Wr2sp \
0.125 0.125 0.125 0.125 1.0 0.125 0.125 0.125 0.125
nsl set Tectum88Model.tectum.Wglsp \
0.5 0.5 0.5 0.5 1.0 0.5 0.5 0.5 0.5
nsl set Tectum88Model.tectum.Wsnsp \
0 0 0 0 20.0 0 0 0 0
nsl set Tectum88Model.tectum.Wtpsp \
0 0 0 0 0.1 0 0 0 0

nsl set Tectum88Model.tectum.Wlpsn \
5.25 5.25 5.25 5.25 5.25 5.25 5.25 5.25 5.25

nsl set Tectum88Model.tectum.Wr2py \
0.437 0.437 0.437 0.437 3.5 0.437 0.437 0.437 0.437
nsl set Tectum88Model.tectum.Wr3py \
0.0375 0.0375 0.0375 0.0375 0.3 0.0375 0.0375 0.0375 0.0375
nsl set Tectum88Model.tectum.Wr4py \
0.875 0.875 0.875 0.875 7.0 0.875 0.875 0.875 0.875
nsl set Tectum88Model.tectum.Wsppy \
0 0 0 0 2.0 0 0 0 0
nsl set Tectum88Model.tectum.Wlppy \
0.56 0.56 0.56 0.56 0.7 0.56 0.56 0.56 0.56
nsl set Tectum88Model.tectum.Wtppy \
0 0 0 0 0.9 0 0 0 0

nsl set Tectum88Model.tectum.Wr3tp \
0.0375 0.0375 0.0375 0.0375 0.3 0.0375 0.0375 0.0375 0.0375
nsl set Tectum88Model.tectum.Wr4tp \
0.625 0.625 0.625 0.625 5.0 0.625 0.625 0.625 0.625

# distance between the center of two squares (dx,dy) should be (0.5,0.5). 
# We set them to (1,1) for visual purposes. All stimuli values should multiplied by 4 accordingly
set indx 1
set indy 1

nsl set Tectum88Model.tectum.in.dx $indx     ;# 32x32 deg (32x32)
nsl set Tectum88Model.tectum.in.dy $indy     ;# 
nsl set Tectum88Model.tectum.in.xz 0		;#       ******* initial Position
nsl set Tectum88Model.tectum.in.yz 15		;#  

set rdx [expr $indx*4]
set rdy [expr $indy*4]

nsl set Tectum88Model.tectum.r2.dx $rdx     ;# 
nsl set Tectum88Model.tectum.r2.dy $rdy     ;# 
nsl set Tectum88Model.tectum.r2.xz 0		;#       ******* initial Position
nsl set Tectum88Model.tectum.r2.yz 3		;#  

nsl set Tectum88Model.tectum.r3.dx $rdx     ;# 
nsl set Tectum88Model.tectum.r3.dy $rdy     ;# 
nsl set Tectum88Model.tectum.r3.xz 0		;#       ******* initial Position
nsl set Tectum88Model.tectum.r3.yz 3		;#  

nsl set Tectum88Model.tectum.r4.dx $rdx     ;# 
nsl set Tectum88Model.tectum.r4.dy $rdy     ;# 
nsl set Tectum88Model.tectum.r4.xz 0		;#       ******* initial Position
nsl set Tectum88Model.tectum.r4.yz 3		;#  

# All stimuli values should multiplied by 4 in accordance to distance between matrix elements
set scalex 4
set scaley 1
set scalev 1
set dx [expr 2*$scalex]    ;# Size, 2 4 8 16 32 
set dy [expr 2*$scaley]    ;# Size, 2 4 8 16 32
set vx [expr 7.6*$scalev]  ;# Speed, number of squares per second, 7.6 deg/sec (8 deg/sec).

nsl create block_stim stim -layer Tectum88Model.tectum.in \
	-spec_type NSL_CORNER -x0 [expr -$dx/2+$rdx/2] -y0 [expr -$dy/2+$rdy/2] -dx $dx -dy $dy -vx $vx
#	-spec_type NSL_CENTER -xc [expr -$dx/2] -yc 0 -dx $dx -dy $dy -vx $vx

#nsl create block_stim stim -layer Tectum88Model.tectum.in1 \
#	-spec_type NSL_CORNER -x0 [expr -$dx/2] -y0 [expr -$dy/2] -dx $dx -dy $dy -vx $vx
#	-spec_type NSL_CENTER -xc [expr -$dx/2] -yc 0 -dx $dx -dy $dy -vx $vx

#nsl create block_stim stim -layer Tectum88Model.tectum.in -val 1.5 \
#	-spec_type NSL_CORNER -x0 -2 -y0 -2 -dx 4 -dy 4 -vx 8

# -y0 15.0
# -dx 8
# -dy 2

# graphics
nsl create window_interface

nsl create display_frame .tect

set width 450
set height 450

#nsl create display_window .tect.in -graph area_level_graph -layer Tectum88Model.tectum.in \
#	-width $width -height $height -wymin 0 -wymax 1.5 -grid NSL_ON
nsl create display_window .tect.in -graph temporal_graph -layer Tectum88Model.tectum.in \
	-width $width -height $height -wymin 0 -wymax 1.5 -t1 $t1 -grid NSL_ON

nsl create display_frame .tect1

set width 600
set height 130

nsl create display_window .tect1.r2 -graph temporal_graph -layer Tectum88Model.tectum.r2 \
	-width $width -height $height -wymin 0 -wymax 1.5 -t1 $t1 -grid NSL_ON
nsl create display_window .tect1.r3 -graph temporal_graph -layer Tectum88Model.tectum.r3 \
	-width $width -height $height -wymin 0 -wymax 1.5 -t1 $t1 -grid NSL_ON
nsl create display_window .tect1.r4 -graph temporal_graph -layer Tectum88Model.tectum.r4 \
	-width $width -height $height -wymin 0 -wymax 1.5 -t1 $t1 -grid NSL_ON
nsl create display_window .tect1.py -graph temporal_graph -layer Tectum88Model.tectum.py \
	-width $width -height $height -wymin 0 -wymax 10 -t1 $t1 -grid NSL_ON
nsl create display_window .tect1.pyf -graph temporal_graph -layer Tectum88Model.tectum.PY \
	-width $width -height $height -wymin 0 -wymax 1.5 -t1 $t1 -grid NSL_ON -pulse NSL_ON

#nsl init