# 1 "/home/admin/SimulationCore2/TorqueDrag/TD_DrillingSubs/TD_DrillStemConfiguration.f90"
subroutine TD_DrillStemConfiguration
Use CStringConfigurationVariables
use SoftwareInputsVariables
! Use CNearFloorConnection
use CDataDisplayConsole
use SimulationVariables
use SimulationVariables !@
use UnitySignalVariables
use UnitySignalsModule
! use UnitySignalsModuleVariables
use OperationScenariosModule
Use TD_DrillStemComponents
use SimulationVariables !@
use SimulationVariables !@
Use TD_StringConnectionData
Use sROP_Variables
use DownHoleModule
implicit none
Integer :: i , j , k , kk
!====================================================
! Drill Stem Components Data Modification
!====================================================
!=====> Drill Stem Total Length&Weight Calculation
data%State%TD_String%DrillStemTotalLength = 0.0d0
data%State%TD_String%DrillStemTotalLengthIni = 0.0d0
data%State%TD_String%DrillStemTotalWeight = 0.0d0
Do i= 1,data%State%TD_String%DrillStemComponentsNumbs
data%State%TD_String%DrillStemTotalLength = data%State%TD_String%DrillStemTotalLength+data%State%TD_DrillStems(i)%Length
data%State%TD_String%DrillStemTotalLengthIni = data%State%TD_String%DrillStemTotalLengthIni+data%State%TD_DrillStems(i)%LengthIni
data%State%TD_String%DrillStemTotalWeight = data%State%TD_String%DrillStemTotalWeight+data%State%TD_DrillStems(i)%Weight
End Do
!=====> Top&Down Depth Calculation Of Initial Drill Stem Components (Graphic)
data%State%TD_DrillStems(1)%DownDepth = data%State%TD_String%DrillStemTotalLength-data%State%TD_StConn%ConnectionHeight
data%State%TD_String%DrillStemBottom = data%State%TD_DrillStems(1)%DownDepth
if ( data%State%TD_DrillStems(1)%DownDepth>=data%State%TD_WellGeneral%WellTotalLength ) then
!if ( data%State%TD_StConn%HookHeight>=data%State%TD_StConn%HookHeightOld .and. data%State%ROP_Bit%RateOfPenetration==0. ) then
! data%State%TD_StConn%GRigidConnectionHeight = data%State%TD_StConn%GRigidConnectionHeight
!else
data%State%TD_StConn%GRigidConnectionHeight = data%State%TD_StConn%GRigidConnectionHeight-(data%State%TD_WellGeneral%WellTotalLength-data%State%TD_DrillStems(1)%DownDepthIniG)
data%State%TD_DrillStems(1)%DownDepthIniG = data%State%TD_WellGeneral%WellTotalLength !???????????
! !data%State%TD_StConn%GRigidConnectionHeight = data%State%TD_StConn%GRigidConnectionHeight
!end if
else
data%State%TD_StConn%GRigidConnectionHeight = data%State%TD_StConn%ConnectionHeight
data%State%TD_DrillStems(1)%DownDepthIniG = data%State%TD_String%DrillStemTotalLengthIni-data%State%TD_StConn%GRigidConnectionHeight
end if
!data%State%TD_DrillStems(1)%DownDepthIniG = data%State%TD_String%DrillStemTotalLengthIni-data%State%TD_StConn%GRigidConnectionHeight
!if ( data%State%TD_DrillStems(1)%DownDepthIniG>data%State%TD_WellGeneral%WellTotalLength ) then
! data%State%TD_StConn%GRigidConnectionHeight = data%State%TD_StConn%GRigidConnectionHeight+(data%State%TD_DrillStems(1)%DownDepthIniG-data%State%TD_WellGeneral%WellTotalLength)
! data%State%TD_DrillStems(1)%DownDepthIniG = data%State%TD_WellGeneral%WellTotalLength !???????????
!end if
data%State%TD_DrillStems(1)%TopDepthIniG = data%State%TD_DrillStems(1)%DownDepthIniG-data%State%TD_DrillStems(1)%LengthIni
Do i = 2,data%State%TD_String%DrillStemComponentsNumbs
data%State%TD_DrillStems(i)%TopDepthIniG = data%State%TD_DrillStems(i-1)%TopDepthIniG-data%State%TD_DrillStems(i)%LengthIni
data%State%TD_DrillStems(i)%DownDepthIniG = data%State%TD_DrillStems(i-1)%DownDepthIniG-data%State%TD_DrillStems(i-1)%LengthIni
End Do
!=====> Top&Down Depth Calculation Of Initial Drill Stem Components (for fluid module)
if ( data%State%TD_DrillStems(1)%DownDepthIniG>=(data%State%TD_WellGeneral%WellTotalLength-.1d0) ) then
data%State%TD_StConn%RigidConnectionHeight = data%State%TD_StConn%GRigidConnectionHeight+(.1d0-(data%State%TD_WellGeneral%WellTotalLength-data%State%TD_DrillStems(1)%DownDepthIniG))
else
data%State%TD_StConn%RigidConnectionHeight = data%State%TD_StConn%GRigidConnectionHeight
end if
data%State%TD_DrillStems(1)%DownDepthIni = data%State%TD_String%DrillStemTotalLengthIni-data%State%TD_StConn%RigidConnectionHeight
data%State%TD_DrillStems(1)%TopDepthIni = data%State%TD_DrillStems(1)%DownDepthIni-data%State%TD_DrillStems(1)%LengthIni
Do i = 2,data%State%TD_String%DrillStemComponentsNumbs
data%State%TD_DrillStems(i)%TopDepthIni = data%State%TD_DrillStems(i-1)%TopDepthIni-data%State%TD_DrillStems(i)%LengthIni
data%State%TD_DrillStems(i)%DownDepthIni = data%State%TD_DrillStems(i-1)%DownDepthIni-data%State%TD_DrillStems(i-1)%LengthIni
End Do
!=====> Top&Down Depth Calculation Of Drill Stem Components
data%State%TD_DrillStems(1)%TopDepth = data%State%TD_String%DrillStemTotalLength-data%State%TD_DrillStems(1)%Length-data%State%TD_StConn%GRigidConnectionHeight
data%State%TD_DrillStems(1)%DownDepth = data%State%TD_String%DrillStemTotalLength-data%State%TD_StConn%GRigidConnectionHeight
Do i = 2,data%State%TD_String%DrillStemComponentsNumbs
data%State%TD_DrillStems(i)%TopDepth = data%State%TD_DrillStems(i-1)%TopDepth-data%State%TD_DrillStems(i)%Length
data%State%TD_DrillStems(i)%DownDepth = data%State%TD_DrillStems(i-1)%DownDepth-data%State%TD_DrillStems(i-1)%Length
End Do
!data%State%TD_String%DrillStemBottom = data%State%TD_DrillStems(1)%DownDepth
!=====> Hole Type & Inclination Determination Of Drill Stem Components
Do i = 1,data%State%TD_String%DrillStemComponentsNumbs
if (data%State%TD_DrillStems(i)%TopDepth .ge. 0.d0) then
Do j = 1,data%State%TD_WellGeneral%WellIntervalsCount
if (data%State%TD_DrillStems(i)%TopDepth>data%State%TD_WellGeo(j)%TopDepth .and. data%State%TD_DrillStems(i)%TopDepth<data%State%TD_WellGeo(j)%DownDepth) then
data%State%TD_DrillStems(i)%HoleType = data%State%TD_WellGeo(j)%HoleType
if (data%State%TD_DrillStems(i)%HoleType == 1) then
data%State%TD_DrillStems(i)%StartAngle = data%State%TD_WellGeo(j)%StartAngle+(((data%State%TD_DrillStems(i)%TopDepth-data%State%TD_WellGeo(j)%TopDepth)/data%State%TD_WellGeo(j)%RCurvature)) ![rad]
data%State%TD_DrillStems(i)%EndAngle = data%State%TD_WellGeo(j)%StartAngle+(((data%State%TD_DrillStems(i)%DownDepth-data%State%TD_WellGeo(j)%TopDepth)/data%State%TD_WellGeo(j)%RCurvature))
data%State%TD_DrillStems(i)%RCurvature = data%State%TD_WellGeo(j)%RCurvature
else if (data%State%TD_DrillStems(i)%HoleType == 2) then
data%State%TD_DrillStems(i)%StartAngle = data%State%TD_WellGeo(j)%StartAngle-(((data%State%TD_DrillStems(i)%TopDepth-data%State%TD_WellGeo(j)%TopDepth)/data%State%TD_WellGeo(j)%RCurvature))
data%State%TD_DrillStems(i)%EndAngle = data%State%TD_WellGeo(j)%StartAngle-(((data%State%TD_DrillStems(i)%DownDepth-data%State%TD_WellGeo(j)%TopDepth)/data%State%TD_WellGeo(j)%RCurvature))
data%State%TD_DrillStems(i)%RCurvature = data%State%TD_WellGeo(j)%RCurvature
else
data%State%TD_DrillStems(i)%StartAngle = data%State%TD_WellGeo(j)%EndAngle !Straight
data%State%TD_DrillStems(i)%EndAngle = data%State%TD_DrillStems(i)%StartAngle
end if
end if
End Do
else
data%State%TD_DrillStems(i)%HoleType = 0 !Straight
data%State%TD_DrillStems(i)%StartAngle = 0.0d0
data%State%TD_DrillStems(i)%EndAngle = 0.0d0
end if
End Do
!=====> Hole(Well) Diameter Determination Of Drill Stem Components
Do i = 1,data%State%TD_String%DrillStemComponentsNumbs
If (data%State%TD_DrillStems(i)%TopDepth>data%State%TD_Casing(1)%TopDepth .and. data%State%TD_DrillStems(i)%TopDepth<data%State%TD_Casing(1)%DownDepth) then
data%State%TD_DrillStems(i)%HoleDiameter = data%State%TD_Casing(1)%Id
Else if (data%State%TD_Liner(1)%Length/=0.d0 .and. data%State%TD_DrillStems(i)%TopDepth>data%State%TD_Liner(1)%TopDepth .and. data%State%TD_DrillStems(i)%TopDepth<data%State%TD_Liner(1)%DownDepth) then
data%State%TD_DrillStems(i)%HoleDiameter = data%State%TD_Liner(1)%Id
Else if (data%State%TD_OpenHole(1)%Length/=0.d0 .and. data%State%TD_DrillStems(i)%TopDepth>data%State%TD_OpenHole(1)%TopDepth .and. data%State%TD_DrillStems(i)%TopDepth<data%State%TD_OpenHole(1)%DownDepth) then
data%State%TD_DrillStems(i)%HoleDiameter = data%State%TD_OpenHole(1)%Id
Else if (data%State%TD_ROPHole(1)%Length/=0.d0 .and. data%State%TD_DrillStems(i)%TopDepth>data%State%TD_ROPHole(1)%TopDepth .and. data%State%TD_DrillStems(i)%TopDepth<data%State%TD_ROPHole(1)%DownDepth) then
data%State%TD_DrillStems(i)%HoleDiameter = data%State%TD_ROPHole(1)%Id !????????? check
End If
!=====> Viscosity Correction Coefficient Calculation
if (data%State%TD_DrillStems(i)%HoleType == 1 .or. data%State%TD_DrillStems(i)%HoleType == 2) then
data%State%TD_DrillStems(i)%MudVisCorrectCoef = (2.0d0/pi)*atan(sqrt(data%State%TD_DrillStems(i)%HoleDiameter**2-(data%State%TD_DrillStems(i)%Od**2)) &
/ (data%State%TD_DrillStems(i)%Od))*(4.0d0/pi-1.0d0) + 1.0d0
else if (data%State%TD_DrillStems(i)%HoleType == 0) then
data%State%TD_DrillStems(i)%MudVisCorrectCoef = 1.d0
end if
End Do
!=====> Out_of_Well Length Of DrillStem
if (data%State%TD_DrillStems(data%State%TD_String%DrillStemComponentsNumbs)%TopDepth<0.) then
data%State%TD_String%OutOfWellLength = abs(data%State%TD_DrillStems(data%State%TD_String%DrillStemComponentsNumbs)%TopDepth)
else
data%State%TD_String%OutOfWellLength = 0.d0
end if
!=====> Near Floor Connection
Do i = data%State%TD_String%DrillStemComponentsNumbs, 1, -1
if ( data%State%TD_DrillStems(i)%TopDepth<=0.d0 ) then
data%State%TD_String%NearFloorConnectionNo = i
data%State%TD_String%NearFloorConnectionHeight = abs(data%State%TD_DrillStems(i)%TopDepth)
end if
if (data%State%TD_DrillStems(i)%TopDepth>0.d0) then
exit
end if
End Do
Call Set_NearFloorConnection(real(data%State%TD_String%NearFloorConnectionHeight))
!=====> Graphic Output
kk = 0
k = 0
Do i= 1,data%State%TD_String%StringConfigurationCount
k = k+1
if ( data%State%TD_DrillStems(k)%DownDepthIniG>0.d0 ) then
kk = kk+1
if (data%State%TD_DrillStems(k+data%State%TD_DrillStem(i)%Numbs-1)%TopDepthIniG<=0.d0) then
exit
end if
k = k+data%State%TD_DrillStem(i)%Numbs-1
end if
End Do
if (allocated(data%State%G_StringElement)) deallocate(data%State%G_StringElement)
allocate(data%State%G_StringElement(kk))
k = 0
Do i= 1,kk
data%State%G_StringElement(kk-i+1)%Id = data%State%TD_DrillStem(i)%Id
data%State%G_StringElement(kk-i+1)%Od = data%State%TD_DrillStem(i)%Od
if ( data%State%TD_DrillStem(i)%ComponentType>4 ) then
data%State%G_StringElement(kk-i+1)%ComponentType = 3
else
data%State%G_StringElement(kk-i+1)%ComponentType = data%State%TD_DrillStem(i)%ComponentType
end if
data%State%G_StringElement(kk-i+1)%TopDepth = data%State%TD_DrillStems(k+data%State%TD_DrillStem(i)%Numbs)%TopDepthIniG
data%State%G_StringElement(kk-i+1)%DownDepth = data%State%TD_DrillStems(k+1)%DownDepthIniG
k = k+data%State%TD_DrillStem(i)%Numbs
End Do
if (data%State%G_StringElement(1)%TopDepth<=0.d0) then
data%State%G_StringElement(1)%TopDepth = 0.d0
end if
Call SetString(kk, data%State%G_StringElement)
if(data%State%TD_DrillStems(1)%DownDepth<0.d0) then
data%Equipments%DrillingWatch%BitPosition = data%State%TD_DrillStems(1)%DownDepth
else
data%Equipments%DrillingWatch%BitPosition = data%State%G_StringElement(kk)%DownDepth !???????????
!if ( data%State%G_StringElement(kk)%DownDepth>data%State%TD_WellGeneral%WellTotalLength ) then
! data%Equipments%DrillingWatch%BitPosition = data%State%TD_WellGeneral%WellTotalLength !???????????
!end if
end if
!if(print_log) print*, 'kk=' , kk
!Do i=1,kk
! if(print_log) print*, 'Number/' , 'downDepth/' , 'topDepth/'
! if(print_log) print*, 'ID/' , 'OD/' , 'Type/'
! if(print_log) print*, i , data%State%G_StringElement(i)%DownDepth , data%State%G_StringElement(i)%TopDepth
! if(print_log) print*, data%State%G_StringElement(i)%Id , data%State%G_StringElement(i)%Od , data%State%G_StringElement(i)%ComponentType
!End Do
!if(print_log) print*, 'Top Joint Height before=' , data%State%TD_String%TopJointHeight
!=====> Top Joint Height
if ( Get_OperationCondition()==OPERATION_DRILL .and. Get_KellyConnection() == KELLY_CONNECTION_STRING ) then
data%State%TD_String%TopJointHeight = abs(data%State%TD_DrillStems(data%State%TD_String%DrillStemComponentsNumbs-1)%TopDepth)
!if(print_log) print*, 'TopDepth1=' , data%State%TD_DrillStems(data%State%TD_String%DrillStemComponentsNumbs-1)%TopDepth
!if(print_log) print*, 'numb1=' , data%State%TD_String%DrillStemComponentsNumbs-1
!if(print_log) print*, 'data%State%TD_String%TopJointHeight=' , data%State%TD_String%TopJointHeight
else
data%State%TD_String%TopJointHeight = abs(data%State%TD_DrillStems(data%State%TD_String%DrillStemComponentsNumbs)%TopDepth)
!if(print_log) print*, 'TopDepth2=' , data%State%TD_DrillStems(data%State%TD_String%DrillStemComponentsNumbs)%TopDepth
!if(print_log) print*, 'numb2=' , data%State%TD_String%DrillStemComponentsNumbs
!if(print_log) print*, 'data%State%TD_String%TopJointHeight2=' , data%State%TD_String%TopJointHeight
end if
!if(print_log) print*, 'data%State%TD_String%DrillStemComponentsNumbs=' , data%State%TD_String%DrillStemComponentsNumbs
!if(print_log) print*, 'data%State%TD_String%TopJointHeight_=' , data%State%TD_String%TopJointHeight
!if(print_log) print*, 'data%State%TD_StConn%ConnectionHeight_=' , data%State%TD_StConn%ConnectionHeight
!if(print_log) print*, 'data%State%TD_StConn%StringConnectionMode_=' , data%State%TD_StConn%StringConnectionMode
!
!if(print_log) print*, 'Top Joint Height after=' , data%State%TD_String%TopJointHeight
end subroutine