# Vizualization on adjacent meshes with freefem

**URL:** https://community.freefem.org/t/vizualization-on-adjacent-meshes-with-freefem/1363
**Category:** General Discussion
**Created:** [December 2, 2021, 12:49pm UTC](https://community.freefem.org/t/vizualization-on-adjacent-meshes-with-freefem/1363 "2021-12-02T12:49:04Z")
**Posts on this page:** 1
**Page:** 1

<div class="post-metadata">

### Author: ![Loic](https://avatars.discourse-cdn.com/v4/letter/l/e9c0ed/32.png) [@Loic](https://community.freefem.org/u/Loic)
#### Post date: [December 2, 2021, 12:49pm UTC](https://community.freefem.org/t/vizualization-on-adjacent-meshes-with-freefem/1363/1 "2021-12-02T12:49:04Z")

</div>

Hello all,

I have different velocity fields **vk[]** define on adjacent (non-overlapping) meshes **ThK** and on FE space **VhK** (each mesh shares one adjacent boundary edge with its neighbour)

How can I gather all these velocity fields, to plot them on the same figure with FreeFem ? What is the best strategy ?

I have tried to add all the meshes **ThK** to re-create the mesh **TH** and then interpolate the results on the global mesh TH :

> //TH=sum(ThK)  
> fespace VH(TH, P1)  
> VH UG, uloc  
> // loop on index k //  
> matrix interpol = interpolate(VH, VhK, inside=1);  
> uloc= interpol\*vk;  
> UG=UG+uloc;

But by doing this, the visualization is very bad and in some cases I can hardly see something.

Thank you,

Best regards,

Loïc,