Dear Gorge
It is indeed possible to use AMR with two-phase-compressible.h from Dr. Daniel Fuster's sandbox. The trick is that you always stay at max level at the interface, you can do this by using a copy of scalar 'f' in the list which is an argument adapt wavelet function for example:
scalar fscalar[];
event adapt (i++) {
foreach(){
fscalar[] = f[];
}
boundary({fscalar});
adapt_wavelet((scalar *){fscalar},(double[]){0.0001,},maxlevel = MAXLEVEL);
}
Best
Mandeep