igen_theta0=0
& igen_theta0=iarr
if (igen_theta0.GT.0)
& call CTRL_MAP_GENARR3D(theta,igen_theta0,myThid)
integer igen_theta0, igen_salt0