############################################################
##  40m model
##  31 Dec, 2010 by K. Arai
#
############################################################

# modulation frequency
const fmod 11.065399M
const fmod3 33.196197M
const fmod5 55.326995M
const loss 0.000050 # 50ppm per reflection

############################################################
###############
# 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    0         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.756                    nPRC5 nPRC6

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

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

# X Arm
s    sMIX  2.013                      nMIX nPOX
m1   ITMX  0.013846 $loss    0        nPOX nARMX1

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

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

# Y Arm
s    sMIY  1.987                     nMIY nPOY
m1   ITMY  0.013846 $loss    0       nPOY nARMY1

### Arm length
### 2011/Dec Value:                  37.54 meter
### complete anti reso for fmod:     33.866
### phi(rcav_fmod)=phi(rcav_fmod5):  37.25258
s    sARMY 37.54                     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.4125                     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

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


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

### TEST1 PRC phase scan
run1: xaxis PRM phi lin -5 5 1000
run1: ad adPRCRF0  0       nPRC1
run1: ad adPRCRF1  $fmod   nPRC1
run1: ad adPRCRF5  $fmod5  nPRC1
run1: gnuterm pdf model_40m_PRC_length_test_run1.pdf
run1: yaxis db

### TEST2 PRC macro length scan
run2: xaxis sPR3 L lin 2.6 3 1000
run2: ad adPRCRF0  0       nPRC1
run2: ad adPRCRF1  $fmod   nPRC1
run2: ad adPRCRF5  $fmod5  nPRC1
run2: gnuterm pdf model_40m_PRC_length_test_run2.pdf
run2: yaxis db

### TEST3 SRC macro length scan
run3: xaxis sSR3 L lin 1 2 1000
run3: ad adPRCRF5  $fmod5  nPRC1
run3: ad addarkRF5  $fmod5  nDARK
run3: ad adPOXRF5  $fmod5   nPOX
run3: ad adPOYRF5  $fmod5   nPOY
run3: gnuterm pdf model_40m_PRC_length_test_run3.pdf

### TEST4 SRC & Schnupp asymmetry scan
run4: xaxis sSR3 L lin 1 2 1000 # 1.4125
run4: x2axis sTRXdummy L lin 0.01 0.1 9
run4: func sMIXcal = 2 - $x2/2
run4: func sMIYcal = 2 + $x2/2
run4: put sMIX L $sMIXcal
run4: put sMIY L $sMIYcal
run4: ad addarkRF5  $fmod5  nDARK
run4: ad adPOXRF5  $fmod5   nPOX
run4: ad adPOYRF5  $fmod5   nPOY

### TEST5 PRC & Schnupp asymmetry scan
run5: xaxis sPR3 L lin 2 3 1000
run5: x2axis sTRXdummy L lin 0.01 0.1 9
run5: func sMIXcal = 2 - $x2/2
run5: func sMIYcal = 2 + $x2/2
run5: put sMIX L $sMIXcal
run5: put sMIY L $sMIYcal
run5: ad adPRCRF1  $fmod   nPRC1
run5: ad adPRCRF5  $fmod5  nPRC1

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