module profiles save c double precision, allocatable:: psi(:) double precision, allocatable:: theta(:) double precision, allocatable:: psibig(:) double precision, allocatable:: pprof(:) double precision, allocatable:: qprof(:) double precision, allocatable:: gprof(:) c double precision:: rma,zma,psifactor end module profiles