mahmood19227
/
SimulationCore2
1
0
Fork 0
Code Issues 0 Pull Requests 0 Releases 0 Wiki Activity
Simulation Core
You can not select more than 25 topics Topics must start with a letter or number, can include dashes ('-') and can be up to 35 characters long.
19 Commits
2 Branches
134 MiB
 
 
 
 
 
 
Tree: 8d59cc0a67
Branches Tags
${ item.name }
Create branch ${ searchTerm }
from '8d59cc0a67'
${ noResults }
SimulationCore2/Equipments/MudSystem/Plot_Final_Mud_Elements.f90

413 lines
20 KiB
Raw Blame History 

  1. subroutine PlotFinalMudElements    !  is called in subroutine CirculationCodeSelect

  2.     

  3.         Use GeoElements_FluidModule

  4.         USE CMudPropertiesVariables

  5.         USE MudSystemVARIABLES

  6.         USE Pumps_VARIABLES

  7.         Use TD_StringConnectionData

  8.         !USE CHOKEVARIABLES

  9.         !USE CDataDisplayConsoleVariables ,  StandPipePressureDataDisplay=>StandPipePressure

  10.         !use CManifolds

  11.         use CDrillWatchVariables

  12.         !use CHOKEVARIABLES

  13.         !use CChokeManifoldVariables

  14.         !use CTanksVariables, TripTankVolume2 => DrillingWatch%TripTankVolume, TripTankDensity2 => TripTankDensity

  15.         USE sROP_Other_Variables

  16.         USE sROP_Variables

  17.         Use KickVariables

  18.         USE CKellyConnectionEnumVariables

  19.         USE UTUBEVARS

  20.         use CLog1

  21.         Use CError

  22.         Use , intrinsic :: IEEE_Arithmetic

  23.    

  24.    implicit none

  25.    

  26.    integer jelement, jmud, jsection,ielement,i

  27.    integer jopelement,jopmud,jopsection

  28.    character(len=120) :: temp1, temp2

  29. 

  30.    

  31.    

  32.    

  33.                     if (ChokeControlPanel%ChokePanelStrokeResetSwitch == 1) then

  34.   write(*,*) 'well cap=' , sum(PipeSection_VolumeCapacity(F_StringIntervalCounts+1:NoPipeSections)) + sum(OpSection_VolumeCapacity(1:F_BottomHoleIntervalCounts)) 

  35.                         DeltaWellCap= sum(PipeSection_VolumeCapacity(F_StringIntervalCounts+1:NoPipeSections)) + sum(OpSection_VolumeCapacity(1:F_BottomHoleIntervalCounts)) - WellCapOld

  36.                         WellCapOld= sum(PipeSection_VolumeCapacity(F_StringIntervalCounts+1:NoPipeSections)) + sum(OpSection_VolumeCapacity(1:F_BottomHoleIntervalCounts))

  37.                         write(*,*) 'cap_reset,DeltaWellCap=' , DeltaWellCap

  38.                     endif 

  39.    

  40.   

  41.    

  42. 

  43. !========================ANNULUS END================= 

  44. if ((Ann_Mud_Forehead_X%Last() - BopStackSpecification%AboveAnnularHeight) > 0.8 .or. Ann_Density%Last()==0.0) then ! for Line (BellNippleToWell-NonFullWell) 

  45.     WellisNOTFull= .true.

  46. else

  47.     WellisNOTFull= .false.

  48. endif

  49. 

  50. !WRITE(*,*) 'Ann_Mud_Forehead_X%Last() , KillHeight', Ann_Mud_Forehead_X%Last() , KillHeight

  51. if ((Ann_Mud_Forehead_X%Last() - BopStackSpecification%KillHeight)>0.8 .or. Ann_Density%Last()==0.0) then ! for Line j4 , WellToChokeManifold(Through 26)

  52.     ChokeLineNOTFull= .true.

  53. else

  54.     ChokeLineNOTFull= .false.

  55. endif

  56. 

  57. 

  58. 

  59. !=========================================================                

  60.         

  61.         jmud= 1

  62.         jsection= 1

  63.         jelement= 0     ! number of final mud elements

  64. 

  65.         

  66.         

  67.         call Xend_MudElement%Empty()

  68.         call TVDend_MudElement%Empty()

  69.         call Density_MudElement%Empty()

  70.         call MudGeoType%Empty()

  71.         call PipeID_MudElement%Empty()

  72.         call PipeOD_MudElement%Empty()

  73.         !call Angle_MudElement%Empty()

  74.         call MudType_MudElement%Empty()

  75.         

  76. 

  77.         

  78.         DO WHILE(jmud <= Hz_Mud_Forehead_X%Length() .and. jsection<=1)

  79.             

  80.                 jelement= jelement+1

  81.                 TrueMinValue= min(Hz_Mud_Forehead_X%Array(jmud), Xend_PipeSection(jsection))

  82. 

  83.                 call Xend_MudElement%Add(TrueMinValue)

  84.                 call TVD_Calculator(TrueMinValue,MudCircVerticalDepth)

  85.                 call TVDend_MudElement%Add(MudCircVerticalDepth)

  86.                 call Density_MudElement%Add(Hz_Density%Array(jmud))

  87.                 call PipeID_MudElement%Add(ID_PipeSectionInch(jsection))

  88.                 call PipeOD_MudElement%Add(OD_PipeSectionInch(jsection))

  89.                 !call Angle_MudElement%Add(Angle_PipeSection(jsection))

  90.                 call MudType_MudElement%Add(Hz_MudOrKick%Array(jmud))

  91.                 

  92.                 

  93.                 if (Xend_MudElement%Array(jelement)== Hz_Mud_Forehead_X%Array(jmud)) then

  94.                     jmud= jmud+1

  95.                 else

  96.                     jsection= jsection+1

  97.                 endif

  98.         

  99.         ENDDO 

  100.         

  101.         NoHorizontalMudElements= jelement   

  102.         

  103. 

  104.         

  105.         

  106.         jmud= 1        

  107.         jsection= 2        

  108.         

  109.         DO WHILE(jmud <= St_Mud_Forehead_X%Length() .and. jsection<=F_StringIntervalCounts)

  110.             

  111.                 jelement= jelement+1

  112.                 TrueMinValue= min(St_Mud_Forehead_X%Array(jmud), Xend_PipeSection(jsection))

  113. 

  114.                 call Xend_MudElement%Add(TrueMinValue)

  115.                 call TVD_Calculator(TrueMinValue,MudCircVerticalDepth)

  116.                 call TVDend_MudElement%Add(MudCircVerticalDepth)                

  117.                 call Density_MudElement%Add(St_Density%Array(jmud))

  118.                 call PipeID_MudElement%Add(ID_PipeSectionInch(jsection))

  119.                 call PipeOD_MudElement%Add(OD_PipeSectionInch(jsection))

  120.                 !call Angle_MudElement%Add(Angle_PipeSection(jsection))

  121.                 call MudType_MudElement%Add(St_MudOrKick%Array(jmud))

  122.                 

  123.                 

  124.                 if (Xend_MudElement%Array(jelement)== St_Mud_Forehead_X%Array(jmud)) then

  125.                     jmud= jmud+1

  126.                 else

  127.                     jsection= jsection+1

  128.                 endif

  129.         

  130.         ENDDO

  131. 

  132.         NoStringMudElements= jelement- NoHorizontalMudElements

  133.         

  134.         

  135. 

  136.         

  137.         

  138.         jmud= 1

  139.         jsection= F_StringIntervalCounts+1

  140.         DO WHILE(jmud<=  Ann_Mud_Forehead_X%Length() .and. jsection<=NoPipeSections)

  141.             

  142.                 jelement= jelement+1

  143.                 TrueMinValue= max(Ann_Mud_Forehead_X%Array(jmud), Xend_PipeSection(jsection))

  144. 

  145.                 call Xend_MudElement%Add(TrueMinValue)

  146.                 call TVD_Calculator(TrueMinValue,MudCircVerticalDepth)

  147.                 call TVDend_MudElement%Add(MudCircVerticalDepth)                

  148.                 call Density_MudElement%Add(Ann_Density%Array(jmud))

  149.                 call PipeID_MudElement%Add(ID_PipeSectionInch(jsection))

  150.                 call PipeOD_MudElement%Add(OD_PipeSectionInch(jsection))

  151.                 !call Angle_MudElement%Add(Angle_PipeSection(jsection))

  152.                 call MudType_MudElement%Add(Ann_MudOrKick%Array(jmud))

  153.                 

  154.                 

  155.                 if (Xend_MudElement%Array(jelement)== Ann_Mud_Forehead_X%Array(jmud)) then

  156.                     jmud= jmud+1

  157.                 else

  158.                     jsection= jsection+1

  159.                 endif

  160.         

  161.         ENDDO        

  162.         

  163.                 do i= 2, Xend_MudElement%Length()

  164.                     if ( i== NoHorizontalMudElements+NoStringMudElements+1) then

  165.                         call Xstart_MudElement%Add (Ann_Mud_Backhead_X%Array(1))        ! start of annulus

  166.                         call TVD_Calculator(Ann_Mud_Backhead_X%Array(1),MudCircVerticalDepth)

  167.                         call TVDstart_MudElement%Add(MudCircVerticalDepth) 

  168.                     elseif ( i== NoHorizontalMudElements+1 ) then

  169.                         call Xstart_MudElement%Add (St_Mud_Backhead_X%Array(1))        ! start of stirng

  170.                         call TVD_Calculator(St_Mud_Backhead_X%Array(1),MudCircVerticalDepth)

  171.                         call TVDstart_MudElement%Add(MudCircVerticalDepth)                         

  172.                     else

  173.                         call Xstart_MudElement%Add(Xend_MudElement%Array(i-1))          ! normal calculation

  174.                         call TVDstart_MudElement%Add(TVDend_MudElement%Array(i-1))          ! normal calculation

  175.                     endif

  176.                     

  177.                 enddo

  178. 

  179.             NoCasingMudElements = jelement- NoStringMudElements- NoHorizontalMudElements

  180.             

  181.             

  182. !=========================For Torque and Drag========================            

  183.         if (allocated(TDXstart_MudElementArray)) deallocate(TDXstart_MudElementArray)

  184.         allocate(TDXstart_MudElementArray(NoHorizontalMudElements+NoStringMudElements+NoCasingMudElements))

  185.         if (allocated(TDXend_MudElementArray)) deallocate(TDXend_MudElementArray)

  186.         allocate(TDXend_MudElementArray(NoHorizontalMudElements+NoStringMudElements+NoCasingMudElements))

  187.         if (allocated(TDDensity_MudElementArray)) deallocate(TDDensity_MudElementArray)

  188.         allocate(TDDensity_MudElementArray(NoHorizontalMudElements+NoStringMudElements+NoCasingMudElements))

  189.         

  190.         TDNoHorizontalMudElements= NoHorizontalMudElements

  191.         TDNoStringMudElements= NoStringMudElements

  192.         TDNoCasingMudElements= NoCasingMudElements 

  193.         

  194. 

  195.            TDXstart_MudElementArray(:) = Xstart_MudElement%Array(:)

  196.            TDXend_MudElementArray(:) = Xend_MudElement%Array(:)

  197.            TDDensity_MudElementArray(:) = Density_MudElement%Array(:)

  198. !=====================================================================            

  199.           

  200.            

  201.         !do i=NoHorizontalMudElements+1, NoHorizontalMudElements+NoStringMudElements            ! 2-string elements

  202.         !        write(*,333) 'STRING:', i,'Xstart\=', Xstart_MudElement%Array(i), 'Xend=' , Xend_MudElement%Array(i), 'Density=' , Density_MudElement%Array(i), 'MudType=' , MudType_MudElement%Array(i)

  203.         !enddo           

  204.            

  205.            

  206.            

  207.  !================================================================

  208.             

  209.  !                   Open Hole Mud Elements    

  210.         jopmud= 1

  211.         jopsection= 1

  212.         jopelement= 0     ! number of final mud elements

  213. 

  214.         

  215.         call Xend_OpMudElement%Empty()

  216.         call TVDend_OpMudElement%Empty()

  217.         call Density_OpMudElement%Empty()

  218.         call PipeID_OpMudElement%Empty()

  219.         call PipeOD_OpMudElement%Empty()

  220.         !call Angle_OpMudElement%Empty()

  221.         call MudTypeOp_MudElement%Empty()

  222.         

  223.             

  224. 

  225.         DO WHILE(jopmud<=  Op_Mud_Forehead_X%Length() .and. jopsection<=F_BottomHoleIntervalCounts)

  226.             

  227.                 jopelement= jopelement+1

  228.                 TrueMinValue= max(Op_Mud_Forehead_X%Array(jopmud), Xend_OpSection(jopsection))

  229.                 call Xend_OpMudElement%Add(TrueMinValue)

  230.                 call TVD_Calculator(TrueMinValue,MudCircVerticalDepth)

  231.                 call TVDend_OpMudElement%Add(MudCircVerticalDepth)                 

  232.                 call Density_OpMudElement%Add(Op_Density%Array(jopmud))

  233.                 call PipeID_OpMudElement%Add(ID_OpSectionInch(jopsection))

  234.                 call PipeOD_OpMudElement%Add(OD_OpSectionInch(jopsection))

  235.                 !call Angle_MudElement%Add(Angle_PipeSection(jopsection))

  236.                 call MudTypeOp_MudElement%Add(Op_MudOrKick%Array(jopmud))

  237.                 

  238.                 

  239.                 if (Xend_OpMudElement%Array(jopelement)== Op_Mud_Forehead_X%Array(jopmud)) then

  240.                     jopmud= jopmud+1

  241.                 else

  242.                     jopsection= jopsection+1

  243.                 endif

  244.         

  245.         ENDDO        

  246.         

  247.                 do i= 2, Xend_OpMudElement%Length()

  248.                     call Xstart_OpMudElement%Add(Xend_OpMudElement%Array(i-1))

  249.                     call TVDstart_OpMudElement%Add(TVDend_OpMudElement%Array(i-1))

  250.                 enddo

  251. 

  252.             NoBottomHoleMudElements = jopelement

  253.             

  254.             

  255.   !================================================================

  256.              

  257.  

  258.         

  259.         if(allocated(StringMudElement)) deallocate(StringMudElement)

  260.         allocate(StringMudElement(NoStringMudElements))

  261.         

  262.         if(allocated(CasingMudElement)) deallocate(CasingMudElement)

  263.         allocate(CasingMudElement(NoCasingMudElements+NoBottomHoleMudElements))        

  264.         

  265.         istring=0    

  266.         icasing=0

  267.         

  268.         BitMudDensity= Density_MudElement%Array(NoHorizontalMudElements+NoStringMudElements)  ! (for ROP module)

  269.   !================================================================   

  270.         

  271.   !============================ UTUBE =============================

  272.         

  273.              !IF (UtubePossibility== .true. .and. Get_KellyConnection() /= KELLY_CONNECTION_STRING .and. WellHeadIsOpen) THEN

  274.         IF (UtubePossibility== .true. .and. TD_FluidStringConnectionMode==0 .and. WellHeadIsOpen .AND. NoGasPocket == 0) THEN

  275.                 CALL WellPressureDataTransfer

  276.                 !WRITE (*,*) ' U-Tube Done 1'

  277.                 CALL Utube

  278.                 !WRITE (*,*) ' U-Tube Done 2'

  279.                 if (QUtubeInput> 0.0) call Utube1_and_TripIn

  280.                 if (QUtubeOutput> 0.0) call Utube2_and_TripIn

  281.         END IF

  282.        

  283.   !========================== UTUBE- end ========================= 

  284.        

  285.         ! do imud=1, st_MudDischarged_Volume%Length()

  286.         !    write(*,*) 'st-plot:', imud, St_MudDischarged_Volume%Array(imud), St_Mud_Backhead_X%Array(imud) ,St_Mud_Forehead_X%Array(imud)

  287.         !enddo          

  288. !==================== Display ========================

  289.         !do i=1, St_MudOrKick%Length()

  290.         !    write(*,555) i,'St_Volume(i), type=' ,St_MudDischarged_Volume%Array(i),St_MudOrKick%Array(i)

  291.         !

  292.         !    IF (IEEE_Is_NaN(St_MudDischarged_Volume%Array(i)))      call ErrorStop('NaN in St Volume-Plot')

  293.         !    IF (St_MudDischarged_Volume%Array(i)<0.)      call ErrorStop('St Volume <0' , St_MudDischarged_Volume%Array(i))

  294.         !enddo              

  295.              

  296.              

  297.         IF (ANY(IEEE_Is_NaN(Op_MudDischarged_Volume%Array(:))) .OR. ANY(Op_MudDischarged_Volume%Array(:) <= 0.0)) THEN

  298.             do i = 1 , Op_MudOrKick%Length()

  299.                 write(*,555) i,'Op_Volume(i), type=' ,Op_MudDischarged_Volume%Array(i) , Op_MudOrKick%Array(i) , Op_Density%Array(i)

  300.             end do        

  301.             call ErrorStop('NaN in Op Volume-Plot or Op Volume <=0')

  302.         END IF

  303.         

  304.             

  305.         IF (ANY(IEEE_Is_NaN(Ann_MudDischarged_Volume%Array(:))) .OR. ANY(Ann_MudDischarged_Volume%Array(:) <= 0.0)) THEN

  306.             do i = 1 , Ann_MudOrKick%Length()

  307.                 write(*,555) i,'Ann_Volume(i), type=' ,Ann_MudDischarged_Volume%Array(i) , Ann_MudOrKick%Array(i) , Ann_Density%Array(i)

  308.             end do        

  309.             call ErrorStop('NaN in Ann Volume-Plot or Ann Volume <=0')

  310.         END IF

  311.         

  312.         !do i=1, Ann_MudOrKick%Length()

  313.         !    !write(*,555) i,'Ann_Volume(i), type=' ,Ann_MudDischarged_Volume%Array(i),Ann_MudOrKick%Array(i),Ann_Density%Array(i)

  314.         !    

  315.         !    IF (IEEE_Is_NaN(Ann_MudDischarged_Volume%Array(i)))      call ErrorStop('NaN in Ann Volume-Plot')

  316.         !    IF (Ann_MudDischarged_Volume%Array(i)<=0.)      call ErrorStop('Ann Volume <=0' , Ann_MudDischarged_Volume%Array(i))

  317.         !enddo

  318.         

  319. 555     FORMAT(I3,5X,A42,(f12.5),5X,I3,5X,(f12.5))         

  320.         

  321.         NoStringMudElementsForPlot= NoStringMudElements

  322.         

  323.         !   1-Horizontal Mud Elements are not shown

  324.                 !write(*,333) 'Horiz:', 1,'Xstart\=', Xstart_MudElement%Array(1), 'Xend=' , Xend_MudElement%Array(1), 'Density=' , Density_MudElement%Array(1), 'MudType=' , MudType_MudElement%Array(1)

  325. 

  326.         do i=NoHorizontalMudElements+1, NoHorizontalMudElements+NoStringMudElements            ! 2-string elements

  327.                 if (Xend_MudElement%Array(i) <= 0.0) then

  328.                     NoStringMudElementsForPlot= NoStringMudElementsForPlot-1

  329.                     cycle

  330.                 endif

  331.                 istring= istring+1

  332.                 StringMudElement(istring)%StartMd = Xstart_MudElement%Array(i)

  333.                 StringMudElement(istring)%EndMd = Xend_MudElement%Array(i)

  334.                 !StringMudElement(istring)%Id = PipeID_MudElement%Array(i)

  335.                 !StringMudElement(istring)%Od = PipeOD_MudElement%Array(i)

  336.                 StringMudElement(istring)%Density = Density_MudElement%Array(i)

  337. 

  338.                 if (MudType_MudElement%Array(i) == 104) then

  339.                     MudType_MudElement%Array(i)= 4  ! air

  340.                 elseif (MudType_MudElement%Array(i) > 0 .and. MudType_MudElement%Array(i) < 100) then       ! all kicks

  341.                     MudType_MudElement%Array(i)= 1  ! gas kick

  342.                 endif

  343. 

  344.                 StringMudElement(istring)%MudType = MudType_MudElement%Array(i)

  345.                 !write(*,333) 'STRING:', i,'Xstart\=', Xstart_MudElement%Array(i), 'Xend=' , Xend_MudElement%Array(i), 'Density=' , Density_MudElement%Array(i), 'MudType=' , MudType_MudElement%Array(i)

  346.         enddo

  347.         

  348.        

  349. 

  350.         do i=Xend_MudElement%Length(), NoHorizontalMudElements+NoStringMudElements+1  , -1               ! 3-casing elements  

  351.                 icasing= icasing+1

  352.                 CasingMudElement(icasing)%StartMd = Xend_MudElement%Array(i)

  353.                 CasingMudElement(icasing)%EndMd = Xstart_MudElement%Array(i)

  354.                 !CasingMudElement(icasing)%Id = PipeID_MudElement%Array(i)

  355.                 !CasingMudElement(icasing)%Od = PipeOD_MudElement%Array(i)

  356.                 !write(*,333) 'CASING:', i,'Xstart\=', Xstart_MudElement%Array(i), 'Xend=' , Xend_MudElement%Array(i), 'Density=' , Density_MudElement%Array(i), 'MudType=' , MudType_MudElement%Array(i)

  357.                 !call Log_1(temp1)

  358.                 !write(*,444) 'CASING:', i,'Xstart\=', Xstart_MudElement%Array(i), 'Xend=' , Xend_MudElement%Array(i), 'PipeID_MudElement%Array(i)=' , PipeID_MudElement%Array(i), 'PipeOD_MudElement%Array(i)=' , PipeOD_MudElement%Array(i)

  359.                 CasingMudElement(icasing)%Density = Density_MudElement%Array(i)

  360.                 

  361.                 if (MudType_MudElement%Array(i) == 104) then

  362.                     MudType_MudElement%Array(i)= 4  ! air

  363.                 elseif (MudType_MudElement%Array(i) > 0 .and. MudType_MudElement%Array(i) < 100) then

  364.                     MudType_MudElement%Array(i)= 1  ! gas kick

  365.                 endif

  366.                 

  367.                 CasingMudElement(icasing)%MudType = MudType_MudElement%Array(i)

  368.                 

  369.         enddo

  370.         

  371.         do i= NoBottomHoleMudElements, 1  , -1               ! 4-open hole elements 

  372.                 icasing= icasing+1

  373.                 CasingMudElement(icasing)%StartMd = Xend_OpMudElement%Array(i)

  374.                 CasingMudElement(icasing)%EndMd = Xstart_OpMudElement%Array(i)

  375.                 !CasingMudElement(icasing)%Id = PipeID_OpMudElement%Array(i)

  376.                 !CasingMudElement(icasing)%Od = PipeOD_OpMudElement%Array(i)

  377.                 !write(*,333) 'OpenHole:',i,'Xstart\=', Xstart_OpMudElement%Array(i), 'Xend=' , Xend_OpMudElement%Array(i), 'Density=' , Density_OpMudElement%Array(i), 'MudType=' , MudTypeOp_MudElement%Array(i)

  378.                 !call Log_1(temp2)

  379.                 !write(*,444) 'OpenHole:',i,'Xstart\=', Xstart_OpMudElement%Array(i), 'Xend=' , Xend_OpMudElement%Array(i), 'PipeID_MudElement%Array(i)=' , PipeID_MudElement%Array(i), 'PipeOD_MudElement%Array(i)=' , PipeOD_MudElement%Array(i)

  380.                 CasingMudElement(icasing)%Density = Density_OpMudElement%Array(i)

  381.                 

  382.                 if (MudTypeOp_MudElement%Array(i) == 104) then

  383.                     MudTypeOp_MudElement%Array(i)= 4  ! air

  384.                 elseif (MudTypeOp_MudElement%Array(i) > 0 .and. MudTypeOp_MudElement%Array(i) < 100) then

  385.                     MudTypeOp_MudElement%Array(i)= 1  ! gas kick

  386.                 endif                

  387.                 

  388.                 CasingMudElement(icasing)%MudType = MudTypeOp_MudElement%Array(i)

  389.         enddo        

  390.         

  391. 

  392.         

  393. 333     FORMAT(A10,I3,5X,A8,(f12.5),5X,A8,(f12.5),5X,A8,(f12.5),5X,A8,I3)  

  394. 444     FORMAT(A10,I2,5X,A8,(f12.3),5X,A8,(f12.3),5X,A8,(f12.3),5X,A8,(f12.3))        

  395.         

  396.         

  397.         

  398.         ! shomare gozari be tartib HZ mud, ST mud, Casing

  399.         ! shomare gzari OpenHole jodagane ast az 1

  400.         

  401.         call SetStringFluids(NoStringMudElementsForPlot, StringMudElement)      !for data display in string

  402.         call SetAnnalusFluids(NoCasingMudElements+NoBottomHoleMudElements, CasingMudElement)       !for data display in casing

  403. 

  404. 

  405. !===========================================================================================================================

  406. !===========================================================================================================================                 

  407. 

  408. 

  409.     end subroutine PlotFinalMudElements

  410.     

  411.     

  412.