~fluidity-core/fluidity/adjoint

« back to all changes in this revision

Viewing changes to adjoint/Adjoint_Python_Fortran.F90

Merge the fluidity_revolve branch.

Show diffs side-by-side

added added

removed removed

Lines of Context:
129
129
          call adj_chkierr(ierr)
130
130
          cycle
131
131
        end if
 
132
 
 
133
        if (index(trim(field_name), "MassMatrix") /= 0) then
 
134
          assert(.not. have_option("/mesh_adaptivity")) ! your coordinate is no longer auxiliary, because you compute it
 
135
          ierr = adj_variable_set_auxiliary(result(j), .true.)
 
136
          call adj_chkierr(ierr)
 
137
          cycle
 
138
        end if
132
139
 
133
140
        ! now we try scalar/vector/tensor
134
141
        if (have_option("/material_phase::" // trim(material_phase_name) // "/scalar_field::" // trim(field_name))) then