code: def point_stokes(self, context):
          
          stokes = np.empty(context.shape, context.dtype)
          stokes[:,:,0] = 1
          stokes[:,:,1:4] = 0
          return stokes
docs: Return a stokes parameter array to montblanc