# intermittent satellite links running with/without HDLC in LLset val(chan)           Channel/Sat                 ;#Channel Typeset val(netif)          Phy/Sat            ;# network interface typeset val(mac)            Mac/Sat                 ;# MAC typeset val(ifq)            Queue/DropTail  ;# interface queue typeset val(ll)             LL/Sat/HDLC                    ;# link layer type#set val(ll)              LL/Sat# two types of eror models used#set val(err)            UrbanComplexMarkovErrorModelset val(err)            RuralComplexMarkovErrorModelset val(bw_up)          10Mbset val(bw_down)        10Mbset val(x)		250set val(y)		250set durlistA        "22.0 26.0 2.9 2.7"  ;# unblocked and blocked state duration avg value for urban modelset durlistB        "27.0 12.0 0.4 0.4"  ;# for rurral modelset BW 10; # in Mb/sset delay 254; # in ms# Set the queue length (in packets)#set qlen [expr round([expr ($BW / 8.0) * $delay * 2])]; # set buffer to pipe size# set sat link bandwidth and delayMac set bandwidth_ 10MbLL set delay_ 125msset val(bdp) [expr 10000/8.0 * 125/100]LL/Sat/HDLC set window_size_ $val(bdp) LL/Sat/HDLC set timeout_ 0.26LL/Sat/HDLC set max_timeouts_ 0# set selective repeat on; its GoBackN by defaultLL/Sat/HDLC set selRepeat_ 1set qlen $val(bdp)puts "queue len = $qlen"set val(ifqlen) $qlenproc start_time {} {return 0; }# proc start_time {} {return [expr ([eval ns-random] / 2147483647.0) * 0.1]}proc UniformErrorProc {} {	global val	set	errObj [new ErrorModel]	$errObj unit bit	#$errObj FECstrength $val(FECstrength) 	#$errObj datapktsize 512	#$errObj cntrlpktsize 80	return $errObj}proc MarkovErrorModel {} {	global durlistA		puts [lindex $durlistA 0]	puts [lindex $durlistA 1]	set errmodel [new ErrorModel/TwoStateMarkov $durlistA time]		#$errmodel drop-target [new Agent/Null]	return $errmodel}proc UrbanComplexMarkovErrorModel {} {	global durlistA	set errmodel [new ErrorModel/ComplexTwoStateMarkov $durlistA time]		$errmodel drop-target [new Agent/Null] 	return $errmodel}proc RuralComplexMarkovErrorModel {} {	global durlistB	set errmodel [new ErrorModel/ComplexTwoStateMarkov $durlistB time]		$errmodel drop-target [new Agent/Null] 	return $errmodel}proc stop {} {    global ns_ tracefd    $ns_ flush-trace    close $tracefd}proc setup {tcptype queuetype sources qlen duration file seed c_sources delay rate ecn run retries} {	global val traces ns_ tracefd	if {$retries == "NULL"} {		set val(ll) LL/Sat	} else {		set val(ll) LL/Sat/HDLC		LL/Sat/HDLC set max_timeouts_ $retries	}	# Initialize Global Variables	ns-random $seed	set ns_		[new Simulator]	set tracefd     [open out.tr w]	$ns_ trace-all $tracefd	$ns_ eventtrace-all	$ns_ rtproto Static		# Create sat n(0)		# configure node, please note the change below.	$ns_ node-config -satNodeType geo \	    -llType $val(ll) \	    -ifqType $val(ifq) \	    -ifqLen $val(ifqlen) \	    -macType $val(mac) \	    -phyType $val(netif) \	    -channelType $val(chan) \	    -downlinkBW $val(bw_down) \	    -wiredRouting ON	set n(0) [$ns_ node]	$n(0) set-position -95	$ns_ node-config -satNodeType terminal	set n(100) [$ns_ node]	$n(100) set-position 42.3 -71.1; # Boston	set n(101) [$ns_ node]	$n(101) set-position 37.9 -122.3; # Berkeley	$n(100) add-gsl geo $val(ll) $val(ifq) $val(ifqlen) $val(mac) $val(bw_up) \	    $val(netif) [$n(0) set downlink_] [$n(0) set uplink_]		$n(101) add-gsl geo $val(ll) $val(ifq) $val(ifqlen) $val(mac) $val(bw_up) \	    $val(netif) [$n(0) set downlink_] [$n(0) set uplink_]		$n(101) interface-errormodel [$val(err)] ;# intermittent downlink		$ns_ unset satNodeType_		# create regular routers	set n(99) [$ns_ node]    ;#....node 3	set n(102) [$ns_ node]   ;#....node 4		# determine the actual queutype and save the one asked for 	# assign queueing parameters -- simulate random drop with RED	if { [string match RandomDrop* $queuetype ] } then {		set queuetype [join \				   [concat "RED" \					[string range $queuetype [string length "RandomDrop"] end ]]\				   "" ]		Queue/$queuetype set thresh_ [expr 2* $qlen]		Queue/$queuetype set maxthresh_ [expr 2* $qlen]		# added for ECN comparison		Queue/$queuetype set setbit_ $ecn	} elseif { [string match RED* $queuetype ] } then {		#added by nirav		Queue/$queuetype set bytes_ false		Queue/$queuetype set queue_in_bytes_ false		# Queue/$queuetype set q_weight_ = -1		# Queue/$queuetype set max_p = 1		Queue/$queuetype set gentle_ true		#added by nirav		#Queue/$queuetype set thresh_ [expr $qlen * 0.5]		#Queue/$queuetype set maxthresh_ [expr $qlen * 0.75]		Queue/$queuetype set thresh_ 0		Queue/$queuetype set maxthresh_ 0		# added for ECN comparison		Queue/$queuetype set setbit_ $ecn	}		$ns_ duplex-link $n(99) $n(100) 10Mb 1ms $queuetype	$ns_ duplex-link $n(101) $n(102) 10Mb 1ms $queuetype# 	# TCP		Agent/$tcptype set window_ $val(bdp)# create nodes for sources and c_sources	puts "sources = $sources"	puts "c_sources = $c_sources"		for { set i 1 } {$i <= [expr 2* [expr $sources+$c_sources]]} {incr i} {	set n($i) [$ns_ node]	}		# connect sources and sinks to routers at 100 Mb/s	for {set i 1} {$i<=$sources} {incr i} {		$ns_ duplex-link $n($i) $n(99) 100Mb 1ms $queuetype		$ns_ duplex-link $n(102) $n([expr $i+$sources+$c_sources]) 100Mb 1ms $queuetype	}	# connect cross traffic sources and sinks	for {set i [expr $sources+1]} {$i<= [expr $c_sources+$sources]} {incr i} {		$ns_ duplex-link $n($i) $n(99) 10Mb 2ms $queuetype		$ns_ duplex-link $n(100) $n([expr $i+$sources+$c_sources]) 10Mb 2ms $queuetype	}		# Setup traffic flow between nodes	# Create source and sink Agents and connect them	for {set i 1} {$i<= [expr $sources+$c_sources]} {incr i} {	# Create the sender agent		set T($i) [new Agent/$tcptype]		# If Ecn is on set it in the sender		if { $ecn } { $T($i) set ecn_ 1 }		set hstcpflag 0		if {$hstcpflag==1} {			puts "creating high-speed connection"			set tcp($i) [$ns_ create-highspeed-connection TCP/Reno $n($i) TCPSink $n([expr $i+$c_sources+$sources]) $i]			if { $ecn } { $tcp($i) set ecn_ 1 }		} else {			puts "cerating normal $T($i) connection"			set tcp($i) [$ns_ create-connection $tcptype $n($i) TCPSink $n([expr $i+$c_sources+$sources]) $i]			if { $ecn } { $tcp($i) set ecn_ 1 }		}		if {$hstcpflag==1} {			[set tcp($i)] set windowOption_ 8			[set tcp($i)] set low_window_ 12.5			[set tcp($i)] set high_window_ 12500			[set tcp($i)] set high_p_ 0.0000096			[set tcp($i)] set high_decrease_ 0.1			puts "windowOption = [$tcp($i) set windowOption_]"			puts "low_window = [$tcp($i) set low_window_]"			puts "high_window = [$tcp($i) set high_window_]"			puts "high_p_ = [$tcp($i) set high_p_]"			puts "high_decrease = [$tcp($i) set high_decrease_]"		}		set F($i) [new Application/FTP]		$F($i) attach-agent $tcp($i)		$ns_ at [expr [start_time]+[expr ($i-1)*2]] "$F($i) start"	}	# tracing for sat links	$ns_ trace-all-satlinks $tracefd		# trace for tcp parameters	if { [lsearch $traces "TCP"] != -1 } {		for {set i 1} {$i<= [expr $sources+$c_sources]} {incr i} {			$tcp($i) trace cwnd_			$tcp($i) trace t_seqno_			$tcp($i) attach [open "$file.tcp_trace.$i" "w"]		}	}		# set queue size to pipe size	$ns_ queue-limit $n(99) $n(100) $val(ifqlen)	$ns_ queue-limit $n(101) $n(102) $val(ifqlen)		# trace queue	if { [lsearch $traces "queue"] != -1 } {		$ns_ trace-queue $n(99) $n(100) [open "$file.queue1" "w"]	}		# setup routing	set satrouteobject_ [new SatRouteObject]	$satrouteobject_ compute_routes	# Stop the sim at $duration	puts "duration = $duration\n"	$ns_ at $duration "stop"	$ns_ at $duration "puts \"NS EXITING...\" ; $ns_ halt"		# start up.	puts "Starting Simulation..."	$ns_ run}	# traces are the set of info to be put outset traces [ list "TCP" "queue" ]# this boring code sets the given parameters from teh command line in order or# assigns defaults if they're not given.if { [llength $argv]>0} then { set run [lindex $argv 0]}\else {set run 1 }if { [llength $argv]>1} then { set retries [lindex $argv 1]}\else {set retries 3 }if { [llength $argv]>2} then { set model [lindex $argv 2]}\else {set model "urban" }if { [llength $argv]>3} then { set sources [lindex $argv 3]}\else {set sources 1}if { [llength $argv]>4} then { set duration [lindex $argv 4]}\else {set duration 300}if { [llength $argv]>5} then { set file [lindex $argv 5]}\else {set file "case0"}if { [llength $argv]>6} then { set qlen [lindex $argv 6]}\else {set qlen 50 }if { [llength $argv]>7} then { set tcptype [lindex $argv 7]}\else {set tcptype "TCP/Newreno"}if { [llength $argv]>8} then { set queuetype [lindex $argv 8]}\else {set queuetype "RED"}if { [llength $argv]>9} then { set seed [lindex $argv 9]}\else {set seed 12345}if { [llength $argv]>10} then { set c_sources [lindex $argv 10]}\else {set  c_sources 0}if { [llength $argv]>11} then { set delay [lindex $argv 11]}\else {set delay "10ms"}if { [llength $argv]>12} then { set rate [lindex $argv 12]}\else {set rate "100k"}if { [llength $argv]>13} then { set ecn [lindex $argv 13]}\else {set ecn 1 }# Start the simputs "run=$run, retries=$retries, sources=$sources, duration=$duration, file=$file, qlen=$qlen, tcptype=$tcptype, queuetype=$queuetype, seed=$seed, c_sources=$c_sources"setup $tcptype $queuetype $sources $qlen $duration $file $seed $c_sources $delay $rate $ecn $run $retries