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 to use in specfem3D, this file does not contain apparently the correct 2D elements at the free surface. I used the following boundary definition,

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)

Is there any suggestions? by the way, the volumes are also merged, so there must not been a problem with it.