← Back to branch summary

~fluidity-core/fluidity/sediment

« back to all changes in this revision

Viewing changes to preprocessor/Boundary_Conditions_From_Options.F90

  • Committer: Sam Parkinson
  • Date: 2011-12-01 15:18:13 UTC
  • mfrom: (3413.1.30 sediment_restructure)
  • Revision ID: s.parkinson11@imperial.ac.uk-20111201151813-9vygndld22tdzepx
merge with sediment_restructure branch

Show diffs side-by-side

added added

removed removed

Lines of Context:
preprocessor/Boundary_Conditions_From_Options.F90
50
50
use integer_set_module !for iceshelf
51
51
use fields_base !for iceshelf
52
52
use fields ! for iceshelf
53
 
use sediment, only: set_sediment_reentrainment, set_sediment_bc_id
 
53
use sediment, only: set_sediment_reentrainment
54
54
use halos_numbering
55
55
use halos_base
56
56
use fefields
223
223
 
224
224
       ! Same thing for sediments. It's of type sediment_reentrainment
225
225
       if (trim(bc_type) .eq. "sediment_reentrainment") then
 
226
          ewrite(2,*) "Changing sediment_reentrainment BC type to neumann"
226
227
          bc_type = "neumann"
227
 
          ! set which ID it is as sediment classes are created onthe fly so
228
 
          ! aren't in the options tree as normal fields
229
 
          call set_sediment_bc_id(field%name, i+1)
230
228
       end if
231
229
 
232
230
       ! Same thing for k_epsilon turbulence model.
705
703
       end if
706
704
 
707
705
       if (trim(bc_type) .eq. "sediment_reentrainment") then
708
 
            ! skip sediment boundareis - done seperately
709
 
            ! see assemble/Sediment.F90
710
 
            cycle boundary_conditions
 
706
          ! skip sediment boundaries - done seperately
 
707
          ! see assemble/Sediment.F90
 
708
          cycle boundary_conditions
711
709
       end if
712
710
 
713
711
       if (trim(bc_type) .eq. "k_epsilon") then
725
723
            surface_element_list=surface_element_list)
726
724
 
727
725
       if((surface_mesh%shape%degree==0).and.(bc_type=="dirichlet")) then
 
726
 
728
727
         ! if the boundary condition is on a 0th degree mesh and is of type strong dirichlet
729
728
         ! then the positions used to calculate the bc should be body element centred not
730
729
         ! surface element centred