A FORTRAN code for γγ → ZZ in SM and MSSM

Through the present paper, the code gamgamZZ is presented, which may be used to calculate all possible observables related to the process gamma gamma -> ZZ, in either the standard model, or the minimal super-symmetric standard model with real parameters.