(eng): We show that the stochastic 3D primitive equations with the Neumann boundary condition on the top, the lateral Dirichlet boundary condition and either the Dirichlet or the Neumann boundary condition on the bottom driven by multiplicative gradient-dependent white noise have unique maximal strong solutions both in the stochastic and PDE senses under certain assumptions on the growth of the noise. For the case of the Neumann boundary condition on the bottom, global existence is established by using the decomposition of the vertical velocity to the barotropic and baroclinic modes and an iterated stopping time argument. An explicit example of non-trivial infinite dimensional noise depending on the vertical average of the horizontal gradient of horizontal velocity is presented.