utoronto-comp-aero - Sept 2013 Merge Notes

Click here to return to top level Jetstream Merge Notes page.

Important change for PTC RANS cases: A bug was found in the time stepping. This is for PTC RANS only.
The following code:
if ( (diablo% nk_time == 17) .and. &
- (jac_app .and. (niter == diablo%strtup%krylv_max) ) .or. &
- ((.not. jac_app) .and. (niter == diablo%newton%krylv_max)) ) then
- deltat_ref = 0.5d0*deltat_ref
- hold = .true.
else
- hold = .false.
end if
has been corrected to the following:
if ( (diablo% nk_time == 17) .and. ( &
- (jac_app .and. (niter == diablo%strtup%krylv_max) ) .or. &
- ((.not. jac_app) .and. (niter == diablo%newton%krylv_max)) ) ) then
- deltat_ref = 0.5d0*deltat_ref
- hold = .true.
else
- hold = .false.
end if
Note the different brackets. If (and only if) you have a cse that runs RANS PTC, you will need to rerun your main unmerged branch with the above changes to update the benchmark for your test case to get your results to match.
Wing-Body-Tail optimization capability via Region Variables implemented by T.Reist
New "least-squares" approach for enforcing a prescribed wing lift distribution as a constraint in an aerodynamic shape optimization case
If you want to use PETSC for the parallel meshmove, you no longer have to uncomment the "#define PETSC" in the movegrid/FEM_cpp.f90 file (and this line is in fact removed). Instead, in the main jetstream Makefile set "petsc = on" and this will pass the required preprocessor directives at compile-time. Setting "petsc = off" will compile the code without PETSC capability.

Sept 2013 Merge Notes