subroutine TD_DrillStemStartUp
    
    Use CCasingLinerChokeVariables
    Use CStringConfigurationVariables
    Use CSimulationVariables
    Use TD_DrillStemComponents
    Use TD_WellElements
    Use TD_WellGeometry
    Use TD_GeneralData
    Use TD_StringConnectionData
   

   implicit none
   
   
  
   
   TD_DrillStemTotalLength    = 500.d0     !???????????????????
   TD_DrillStemTotalLengthIni = 500.d0     !???????????????????
   TD_OutOfWellLength         = 3.48d0
   TD_DrillStemAxialVelocity  = 0.d0
   TD_DrillStemRotVelocity    = 0.d0
   TD_WeightOnBit             = 0.d0
   TD_BitTorque               = 0.d0
   TD_TotalTorque             = 0.d0
   TD_StaticHookLoad          = 0.d0
   TD_DlMax                   = 0.d0
   TD_DlTotal                 = 0.d0
   TD_DlTouch                 = 1.d0           !(=/0.) because of TD_StringConnectionModes subroutine  (initial value for first iteration in if condition)
   TD_HookLoad                = 0.d0
   TD_StringTorque            = 0.d0
   TD_ToolJointRange          = 0.4005d0*3.28  ![m]*3.28=[ft]  ,  ToolJointLength=400.5mm
   TD_RemoveVolume            = 0.0d0
   
   TD_DrillStemBottom           = TD_DrillStemTotalLength-TD_OutOfWellLength   !????????????
   !TD_DrillStems%TopDepth      = 
   !TD_DrillStems%DownDepth     = 
   !TD_DrillStems%StartAngle    = 
   !TD_DrillStems%EndAngle      = 
   TD_DrillStems%Area           = 0.d0
   TD_DrillStems%RCurvature     = 0.d0
   TD_DrillStems%RtoolJoint     = 0.d0
   TD_DrillStems%ToolJointRange = 0.d0
   TD_DrillStems%HoleDiameter   = 0.d0
   
   TD_DrillStems%Force1       = 0.d0
   TD_DrillStems%Force2       = 0.d0
   TD_DrillStems%Torque       = 0.d0
   TD_DrillStems%Drag         = 0.d0
   TD_DrillStems%CombVelRatio = 0.d0
   
   TD_DrillStems%MudDensityIn      = 10.d0*7.48051948d0
   TD_DrillStems%MudDensityOut     = 10.d0*7.48051948d0
   TD_DrillStems%MudViscosity      = 0.16d0   !????????????????
   TD_DrillStems%MudVisCorrectCoef = 1.d0     !????????????????
   TD_DrillStems%MudWeight         = 10.d0
   !TD_DrillStems%MudPlasticVis    =   !????????????????
   !TD_DrillStems%MudYieldPoint    =   !????????????????
   
   !TD_DrillStems%FricFactor = 
   !TD_DrillStems%           =
   !TD_DrillStems%           =
   
   !TD_DrillStems%StaticHookLoad = 
   !TD_DrillStems%TotalSHookLoad = 
   !TD_DrillStems%Dl             = 
   !TD_DrillStems%DlTotal        = 
   
   !if (allocated(TD_FluidMudDensity))  deallocate(TD_FluidMudDensity)
   !if (allocated(TD_FluidMudEndX))     deallocate(TD_FluidMudEndX)
   !if (allocated(TD_FluidMudStartX))   deallocate(TD_FluidMudStartX)
   !Allocate (TD_FluidMudDensity(1) , TD_FluidMudEndX(1) , TD_FluidMudStartX(1))
   !TD_FluidMudDensity(:) = 10.d0
   !TD_FluidMudStartX(:)  = -268.d0     !??????????
   !TD_FluidMudEndX(:)    = 9990.d0     !??????????
   
   
   
   
  
end subroutine