← Back to branch summary

~fluidity-core/fluidity/sea-ice-branch

« back to all changes in this revision

Viewing changes to assemble/Advection_Diffusion_DG.F90

  • Committer: Simon Mouradian
  • Date: 2012-03-30 11:11:47 UTC
  • mfrom: (3520.32.239 fluidity)
  • Revision ID: simon.mouradian06@imperial.ac.uk-20120330111147-y4n4y4th52x2dqcf
merge from lp:fluidity

Show diffs side-by-side

added added

removed removed

Lines of Context:
assemble/Advection_Diffusion_DG.F90
528
528
    call allocate(rhs_diff, T%mesh, trim(field_name)//" Diffusion RHS")
529
529
   
530
530
    call construct_advection_diffusion_dg(matrix, rhs, field_name, state, &
531
 
         mass, matrix_diff, rhs_diff, semidiscrete=.true., &
 
531
         mass=mass, diffusion_m=matrix_diff, diffusion_rhs=rhs_diff, semidiscrete=.true., &
532
532
         velocity_name=velocity_name)
533
533
 
534
 
    call get_dg_inverse_mass_matrix(inv_mass, mass)
 
534
    ! mass has only been assembled only for owned elements, so we can only compute
 
535
    ! its inverse for owned elements
 
536
    call get_dg_inverse_mass_matrix(inv_mass, mass, only_owned_elements=.true.)
535
537
    
536
538
    ! Note that since theta and dt are module global, these lines have to
537
539
    ! come after construct_advection_diffusion_dg.