        mask=(img2==22000 | (img2==12000 & xxm<60));
        lset=sqrt((xxm-43).^2*1.5+(yym-62).^2/2)-8;
        nb=(lset<19 & lset>-2); delt=0.5; niter=50;
        if fuserf==0,lset=lset.*(1-nb)+lset1;end
        evolve;

        mask=(img2==22000 | (img2==12000 & xxm<60));
        lset=sqrt((xxm-43).^2+(yym-62).^2*2)-26;
        nb=(lset<10 & lset>-24); delt=0.5; niter=50;
        if fuserf==0,lset=lset.*(1-nb)+lset1;end
        evolve;

        mask=(img2==22000 | (img2==12000 & xxm>60));
        lset=sqrt((xxm-98).^2*1.5+(yym-62).^2/2)-8;
        nb=(lset<19 & lset>-2); delt=0.5; niter=50;
        if fuserf==0,lset=lset.*(1-nb)+lset2;end
        evolve;

        mask=(img2==22000 | (img2==12000 & xxm>60));
        lset=sqrt((xxm-93).^2+(yym-62).^2*1.5)-23;
        nb=(lset<9 & lset>-22); delt=0.5; niter=50;
        if fuserf==0,lset=lset.*(1-nb)+lset2;end
        evolve;
