Here, the chosen internal pressure gradient method is selected and (in case that the compiler option SUBSTR_INI_PRESS is set), the initial pressure is calculated and subtracted from the updated internal pressure gradient.
If GETM is executed as slice model (compiler option SLICE_MODEL is set, the internal pressure gradient for is copied to . USES:
use getm_timers, only: tic, toc, TIM_INTPRESS IMPLICIT NONELOCAL VARIABLES:
integer :: i,j,k logical, save :: first=.true.