← Back to branch summary

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

« back to all changes in this revision

Viewing changes to assemble/Geostrophic_Pressure.F90

  • Committer: Simon Mouradian
  • Date: 2012-05-30 09:08:23 UTC
  • mfrom: (3520.32.274 fluidity)
  • Revision ID: simon.mouradian06@imperial.ac.uk-20120530090823-pq0fgemjpqn72j8d
merge from lp:fluidity

Show diffs side-by-side

added added

removed removed

Lines of Context:
assemble/Geostrophic_Pressure.F90
110
110
    !! Mass option path
111
111
    character(len = OPTION_PATH_LEN) :: mass_option_path
112
112
    !! Divergence matrix
113
 
    type(block_csr_matrix) :: ct_m
 
113
    type(block_csr_matrix), pointer :: ct_m
114
114
    !! RHS terms from integrating the divergence operator by parts
115
115
    type(scalar_field) :: ct_rhs
116
116
    !! Mass matrix. Only used when not lumping mass for continuous u_mesh.
785
785
    ewrite(2, *) "Boundary conditions field: ", trim(lbcfield%name)
786
786
    
787
787
    ct_sparsity => get_csr_sparsity_firstorder(state, matrices%p_mesh, matrices%u_mesh)
 
788
    allocate(matrices%ct_m)
788
789
    call allocate(matrices%ct_m, ct_sparsity, blocks = (/1, dim/), name = "CT")
789
790
    call allocate(matrices%ct_rhs, matrices%p_mesh, "CTRHS")
790
791
    
1197
1198
      end select
1198
1199
    end if
1199
1200
    call deallocate(matrices%ct_m)
 
1201
    deallocate(matrices%ct_m)
1200
1202
    call deallocate(matrices%ct_rhs)
1201
1203
    
1202
1204
    if(matrices%have_cmc_m) then
1434
1436
    !!< Compute the divergence of a field
1435
1437
  
1436
1438
    type(vector_field), intent(in) :: field
1437
 
    type(block_csr_matrix), intent(in) :: ct_m
 
1439
    type(block_csr_matrix), pointer :: ct_m
1438
1440
    type(csr_matrix), intent(in) :: mass
1439
1441
    type(scalar_field), intent(inout) :: div
1440
1442