Boundary definition by GEOCUBIT for 2-layer model

Hello
I am trying to export the following mesh (with 2 irregular layers: volumes yellow and green) from cubit using geocubit library, but when i got the exported free_or_absorbing_surface_file_zmax file, this file contains the nodes in the interface as if they were in the free surface too. I used the following boundary condition, Is there any suggestions? by the way, the volumes are also merged, so there must not been a problem with it.

This is boundary_definition.py of GEOCUBIT

#… which extracts the bounding faces and defines them into blocks
boundary_definition.entities=[‘face’]
boundary_definition.define_bc(boundary_definition.entities,parallel=True)

thanks!