############################################################
##  40m model
##  1 Aug, 2010 by K. Arai
#
############################################################

# modulation frequency
const fmod 11.065399M
const fmod3 33.196197M
const fmod5 55.326995M
const fmod15 165.980985M
const loss 0.000050 # 50ppm

############################################################
###############
# IOO section
###############

#    name Pow     f_ofs              node  
l    L1   1       0                  nL         # Laser

#    name f_mod   depth order pm/am  nodes
mod  EOM1 $fmod   .1     3     pm     nL  n1	# EOM
mod  EOM2 $fmod5  .1     3     pm     n1  n2	# EOM

s    s1   1.0                        n2  nREFL

###############
# PRC section
###############
# PRM
#    name T     Loss     phi         node1 node2
m1   PRM  0.05637 $loss    90         nREFL nPRC1
### NO PRM ###
#m1   PRM  1 0 0         nREFL nPRC1

s    sPR1 1.0                        nPRC1 nPRC2

# bs name T        loss    phi     alpha     node1 node2 node3 node4
bs1  PR2  0.000025 $loss   0       5         nPRC2 nPRC3 nPR2F nPR2B

s    sPR2 1.0                        nPRC3 nPRC4

# bs name T        loss    phi     alpha     node1 node2 node3 node4
bs1  PR3  0.000025 $loss   0       5         nPRC4 nPRC5 nPR3F nPR3B

### add 4 meters to this to get actual PRC length ###
s    sPR3 2.7538                    nPRC5 nPRC6

###############
# Michelson section
###############

# Main BS
# bs name T   loss    phi     alpha     node1 node2 node3 node4
bs1  BS   0.5 $loss   63.6396105 45        nPRC6 nMIY  nMIX  nSRC6

# X Arm
s    sMIX  2.0171                     nMIX nPOX1

bs2  POX   0.00041692 $loss   0 1     nPOX1 nPOXF  nPOX2  nPOXB

s    sPOX  0                          nPOX2 nPOX3

m1   ITMX  0.013846 $loss    0        nPOX3 nARMX1

### Arm length
### 2011/Dec Value:                  37.54 meter
### complete anti reso for fmod:     33.866
### phi(rcav_fmod)=phi(rcav_fmod5):  37.25258
### phi(rcav_fmod)=phi(rcav_fmod5)/5:  37.7974
s    sARMX 37.7974                     nARMX1 nARMX2

m1   ETMX  0.0000137 $loss   0       nARMX2 nTRX
s    sTRXdummy 1	    	     nTRX nTRXdummy

# Y Arm
s    sMIY  1.9829                    nMIY nPOY1

bs2  POY   0.00041692 $loss   0 1    nPOY1 nPOYF nPOY2 nPOYB

s    sPOY  0                         nPOY2 nPOY3

m1   ITMY  0.013846 $loss    0       nPOY3 nARMY1

### Arm length
### 2011/Dec Value:                  37.54 meter
### complete anti reso for fmod:     33.866
### phi(rcav_fmod)=phi(rcav_fmod5):  37.25258
### phi(rcav_fmod)=phi(rcav_fmod5)/5:  37.7974
s    sARMY 37.7974                     nARMY1 nARMY2

m1   ETMY  0.0000137 $loss   0       nARMY2 nTRY
s    sTRYdummy 1	    	     nTRY nTRYdummy

###############
# SRC section
###############

### add 4 meters to this to get actual SRC length ###
s    sSR3 1.39915                     nSRC6 nSRC5

# bs name T        loss    phi     alpha     node1 node2 node3 node4
bs1  SR3  0.000025 $loss   0       5         nSRC5 nSRC4 nSR3F nSR3B

s    sSR2 1.0                        nSRC4 nSRC3

# bs name T        loss    phi     alpha     node1 node2 node3 node4
bs1  SR2  0.000025 $loss   0       5         nSRC3 nSRC2 nSR2F nSR2B

s    sSR1 1.0                        nSRC2 nSRC1

# SRM
#m1   SRM  0.09903      $loss   90       nSRC1 nDARK # with SRM
m1   SRM  0.09903      $loss   0       nSRC1 nDARK # with SRM

### NO SRM ###
#m1   SRM  1 0  95       nSRC1 nDARK            # no SRM


#####################################
###  This is the end of the model ###
#####################################
############################################################

##   Test 1 Carrier Power Recycling
run1: xaxis L1 f log 0.01 1M 10000

run1: ad adREFL  0   nREFL
run1: put adREFL f $x1

run1: ad adPRC1  0   nPRC1
run1: put adPRC1 f $x1

run1: ad adDARK  0   nDARK
run1: put adDARK f $x1

run1: ad adMIX  0   nMIX* #power returning from the arm
run1: put adMIX f $x1

run1: ad adARMX  0   nARMX1* #power inside the arm
run1: put adARMX f $x1

run1: ad adTRX  0 nTRX
run1: put adTRX f $x1

run1: gnuterm pdf model_40m_power_budget_test_run1.pdf
run1: yaxis log abs

##   Test 2 Refl RF structure
run2: xaxis L1 f lin -60M 60M 10000

run2: ad adREFL  0   nREFL
run2: put adREFL f $x1

run2: ad adPRC1  0   nPRC1
run2: put adPRC1 f $x1

run2: ad adDARK  0   nDARK
run2: put adDARK f $x1

run2: gnuterm pdf model_40m_power_budget_test_run2.pdf
run2: yaxis log abs

##   Test 3 Power (Transfer Function) to each port

run3: xaxis L1 f lin 0 $fmod15 15

run3: ad adREFL  0   nREFL
run3: put adREFL f $x1

run3: ad adPRC1f  0   nPRC1
run3: put adPRC1f f $x1

run3: ad adPRC1b  0   nPRC1*
run3: put adPRC1b f $x1

run3: ad adPOP1  0   nPR2B
run3: put adPOP1 f $x1

run3: ad adDARK  0   nDARK
run3: put adDARK f $x1

run3: ad adMIXf  0   nMIX # power towards the arm
run3: put adMIXf f $x1

run3: ad adMIXb  0   nMIX* # power returns from the arm
run3: put adMIXb f $x1

run3: ad adPOX  0 nPOXB
run3: put adPOX f $x1

run3: ad adARMX  0   nARMX1* # power inside the arm
run3: put adARMX f $x1

run3: ad adTRX  0   nTRX # arm transmittion
run3: put adTRX f $x1

run3: yaxis log abs

run3: gnuterm pdf model_40m_power_budget_test_run3.pdf

#yaxis abs
#yaxis log abs:deg

#gnuterm pdf model_40m_power_budget_test.pdf
#gnuterm windows
#gnuterm x11

pause -1
GNUPLOT
set grid
#set dgrid3d
#set nosurface
#set contour
#set cntrparam levels 10
END