##################################################                                       
#### Script to configure the follwing Variances (NOT Standard deviations)
#### based on user inputed switches
#### SORO, RFSP, Tuf, Sop, CCP, Vpf, Vdd, mdd

### (Nov 29, 2018) M.Secor - Added call to scaling factors file
source('Inputs//ScalingFactors.R')

#### Reading user inputs
source('Inputs//SetVariances-Standard.R')

##################################################
#### SORO
if(Ssoro == 1){
    set.SORO.vc = Vsoro.variance.common           		#Exchangeable random
    set.SORO.vs = Vsoro.variance.specific      		    #Exchangeable common
}else if(Ssoro == 3){
    set.SORO.v = Vsoro.variance.common+
                 Vsoro.variance.specific  	#Exchangeable diagonal
    set.SORO.l = 10.30 	       		        #Exchangeable length-scale
}

##################################################
#### RFSP
#### Standard Deviations 1.78%/1.54% (All Epistemic ?)
if(Srfsp == 2){
    set.RFSP.vs2 = Vrfsp.variance.cop.specific           	#Partitioned Exchangeable
    set.RFSP.vc2 = Vrfsp.variance.cop.common
    set.RFSP.vo2 = Vrfsp.variance.copfop.common
    set.RFSP.ws2 = Vrfsp.variance.fop.specific
    set.RFSP.wc2 = Vrfsp.variance.fop.common
}else if(Srfsp == 4){
    set.RFSP.vc4 = Vrfsp.variance.cop.specific+Vrfsp.variance.cop.common   	#Partitioned Spatial
    set.RFSP.vo4 = Vrfsp.variance.copfop.common
    set.RFSP.wc4 = Vrfsp.variance.fop.specific+Vrfsp.variance.fop.common
    set.RFSP.l4  = 10.30
}


##################################################
#### Tuf
if(Stuf == 1){#Exchangeable
	set.Tuf.vs1.n = Vtuf.variance.specific.north
	set.Tuf.vs1.s = Vtuf.variance.specific.south
	set.Tuf.vc1.n = Vtuf.variance.common.north
	set.Tuf.vc1.s = Vtuf.variance.common.south
    set.Vtuf.grand = Vtuf.variance.cross
}else if(Stuf == 3){#Spatial
	set.Tuf.v3.n = Vtuf.variance.specific.north+Vtuf.variance.common.north
	set.Tuf.v3.s = Vtuf.variance.specific.south+Vtuf.variance.common.south
	set.Tuf.l3 = 10.30
}


##################################################
#### Sop
#### This model is rather odd, and we don't have the data for it.
#### Zsop: A-variance of SOP^J
if(Ssop == 2){
    set.Zsop.vs2 = Zsop.variance.cop.specific
	set.Zsop.vc2 = Zsop.variance.cop.common
	set.Zsop.vo2 = Zsop.variance.copfop.common
	set.Zsop.ws2 = Zsop.variance.fop.specific
	set.Zsop.wc2 = Zsop.variance.fop.common
}else if(Ssop == 4){
    set.Zsop.vc4 = Zsop.variance.cop.specific+Zsop.variance.cop.common
	set.Zsop.vo4 = Zsop.variance.copfop.common
	set.Zsop.wc4 = Zsop.variance.fop.specific+Zsop.variance.fop.common
    set.Zsop.l4 =  10.30
}

##################################################
#### CCP
#### Zccp  
#### This is the error variance for CCP
#### In Dumitru's documents, this is, CCP Error 'Operating Conditions'
## North can be read from a file
North<-c(6:10,18:24,33:40,50:58,69:78,90:100,112:122,134:144,157:168,
         181:192,205:216,229:240,253:264,277:288,301:312,325:336,
         348:358,370:380,392:402,413:422,432:440,449:456,464:470,476:480)
South<-setdiff(1:Ni,North)

if(Sccp == 1){ #Exchangeable
	set.Zccp.vs1.n = Zccp.variance.specific.north
	set.Zccp.vs1.s = Zccp.variance.specific.south
	set.Zccp.vc1.n = Zccp.variance.common.north
	set.Zccp.vc1.s = Zccp.variance.common.south
	set.Zccp.grand = Zccp.variance.cross

}else if(Sccp == 3){ #Spatial
	set.Zccp.v3.n = Zccp.variance.specific.north+Zccp.variance.common.north
	set.Zccp.v3.s = Zccp.variance.specific.south+Zccp.variance.common.south
	set.Zccp.l3 = 10.30
}


##################################################
#### Vpf
#### See SUMMARY OF UNCERTAINTY STRUCTURE document
Vpf =Vpf.variance

##################################################
#### Vdd, mdd
Vddt = Vnop.variance

if(Msec == 1){
    mdd<-SDS1.DD.mean
    Vdd = SDS1.DD.variance
    Vdd11 = SDS1.buffer.variance
    Vdd12 = SDS1.dsc.variance
    Vdde = Vddt + Vdd11 + Vdd12
}else if(Msec == 2){
    mdd<-SDS2.DD.mean
    Vdd = SDS1.DD.variance
    Vdd21 = SDS2.buffer.variance2
    Vdd22 =SDS2.dsc.variance
    Vdde = Vddt + Vdd21 + Vdd22
}

## EOF