def init_interface(field):
interface_shape = np.array(field.shape); interface_shape[-<NUM_LIT:1>] += <NUM_LIT:1><EOL>interfaces = np.tile(False,len(interface_shape)); interfaces[-<NUM_LIT:1>] = True<EOL>interface_zero = Field(np.zeros(interface_shape), domain=field.domain, interfaces=interfaces)<EOL>return interface_zero<EOL>
Return a Field object defined at the vertical interfaces of the input Field object.