& momPressureForcing, metricTerms, useNHMTerms,
LOGICAL metricTerms
IF ( usingSphericalPolarGrid .AND. metricTerms ) THEN
IF ( usingCylindricalGrid .AND. metricTerms ) THEN
CALL WRITE_0D_L( metricTerms, INDEX_NONE, 'metricTerms =',
metricTerms = .FALSE.
& metricTerms, vectorInvariantMomentum,
metricTerms = .TRUE.
IF ( metricTerms ) THEN