FLIPS (Fortran Linear Inverse Problem Solver) is a Fortran95 module for solving large scale statistical linear inverse problems of formm = Ax + ewhere m is called the measurement, A is the theory matrix, x is the unknown and e is the measurement error. Usually, the measurement and the theory matrix are given and some statistical information about the error is known. The task is to extract as much information about the unknown as possible. FLIPS is able to calculate the maximum a posteriori estimate (MAP estimate) and the posteriori covariance matrix of the unknown. If the error and (possible) a priori information are assumed to be Gaussian, so is the posteriori distribution of the unknown, and the MAP estimate together with the posteriori covariance matrix are enough to determine it.FLIPS solves the inverse problem by first transforming it into a equivalent (overdetermined) least squares problem. It is then further transformed into a simple upper triangular system using Givens rotations which is then easy to solve using back substitution. The Givens rotations are made row-by-row which makes it possible to feed the problem data (measurements, theory matrix rows & error variance/error covariance matrix) into FLIPS in small pieces, thus decreasing the memory footprint in the computer.