// ############################## START FUNC ############################################
function recognize_features(patterns,MW)
{

     var multipath=patterns.split("%");
     var collected_string="";
     var group="";
     var stringparts=new Array(0);
     var d=0;
     var xt=0;
     var yt=0;
     var zt=0;
     var count1;
     var count2;
     var count3;
     var count4;
     var NH=0;
     var NHstr="";

       for(n=0;n<multipath.length;n++){
       group="";
       stringparts=multipath[n].split("£");
       stringparts[0]=stringparts[0].substring(1,stringparts[0].length);
        if(multipath[n]!="#" && stringparts[5].indexOf("@@@")==-1 && stringparts[4]!="complex"){      // =========== if starts ==============
        stringparts[5]="#"+stringparts[5];
switch(stringparts[5]){

  case "#*7**6**6**7* r *7**6**6**7* r *7**6**6**7* r *7**6**6**7* r #":
    group="tertiary_amine";
    break;
  case "#*7**6**6**7* r *7**6**6**7* r #":
// This is inaccurate. See later in the list.
    group="azo_group";
    break;
  case "#*16**6**6**16* r *16**6**6**16* r #":
    group="sulfide";
    break;
  case "#*16**6**6**16* r *16**1**1**16* r #":
    group="mercaptan";
    break;
  case "#*6**8* d *8**6* r  d *6**7**7**6**6**7* r *7**1**1**7* r *7**6* r #":
    group="primary_amide";
     if(stringparts[4]=="i") group="primary_lactam";
    break;
  case "#*6**8* d *8**6* r  d *6**7**7**6**6**7* r *7**6**6**7* r *7**6* r #":
    group="secondary_amide";
     if(stringparts[4]=="i") group="secondary_lactam";
    break;
  case "#*6**8* d *8**6* r d *6**7* d *7**6* d *6**7* r d *7**6* r d #":
    group="secondary_amide";
     if(stringparts[4]=="i") group="secondary_lactam";
    break;
  case "#*6**8**8**6**6**8* r *8**6* r *6**8* d *8**6* r  d #":
    group="esther";
     if(stringparts[4]=="i") group="lactone";
    break;
  case "#*7**6**6**7* r *7**1**1**7* r *7**1**1**7* r #":
    group="primary_amine";
    break;
  case "#*7**6**6**7* r #":
    group="primary_amine";
    break;
  case "#*6**8**8**1**1**8* r *8**6* r *6**8* d *8**6* r  d #":
    group="carbonyl_acid";
    break;
  case "#*6**8* d *8**6* r  d *6**8* d *8**6* r  d #":   // If deprotonated acid
    group="carbonyl_acid";
    break;
  case "#*7**6**6**7* r *7**6* d *6**7* r  d #":
    group="imine";
    break;
  case "#*7**8**8**7* r *7**8**8**7* r *7**6**6**7* r #":
    group="nitro_group";
    break;
  case "#*7**6* d *6**7* r  d *7**6* d *6**7* r  d #": // imine with neigbour C=C carbon
    group="imine";
    break;
// Because the N below has no double bond, all neighbour double bond are ignored
  case "#*7**6**6**7* r *7**6**6**7* r *7**6**6**7* r #": 
    group="tertiary_amine";
    break;
  case "#*8**6**6**8* r *8**6**6**8* r #":
    group="ether";
    break;
  case "#*6**8* d *8**6* r *6**1**1**6* r #":
    group="aldehyde";
    break;
  case "#*6**8* d *8**6* r  d *6**1**1**6* r #":
    group="aldehyde";
    break;
  case "#*8**6**6**8* r *8**1**1**8* r #":
    group="alcohole";
     if(stringparts[4]=="y") group="aromatic_alcohole";
    break;
  case "#*8**6* d *6**8* r  d #":
    group="ketone";
    break;
  case "#*7**7* d *7**6**6**7* r *7**7* r  d *7**6**6**7* r #":
    group="azo_group";                         // Will be doubled first and fixed in the function "complete".
    break;
  default:
    if(stringparts[5].length>69) group="electronegative_complex";
}
        }                                            // =========== if else ==============
        else{
         if(multipath[n]!="#" && stringparts[4]!="complex"){
         d=parseInt(stringparts[3]);
         group="big_fat_group";
         if(d<18) group="medium_fat_group";
         if(d<12) group="small_fat_group";
           if(stringparts[4]=="r" && d<18) group="";  
           // Ignore medium/small fat groups on fat rings above
           if(MW/d<3) group="";
           // Ignore fat complex which involves big parts of the molecule.
         }
         if(stringparts[4]=="complex"){
         NH=0;
         if(stringparts[5].indexOf("*1*")!=-1) NH=2;                // 2 if NH present in ring, else 0
         count1=(stringparts[5].split("*").length-1-NH*4)/8;
         count2=(stringparts[5].split("*7*").length-1-NH)/4;
         count3=(stringparts[5].split("*8*").length-1)/4;
         count4=(stringparts[5].split("*16*").length-1)/4;
         NHstr="";
         if(NH==2) NHstr="_NH";

         // So count1=size of ring, count2 numbers of N in ring, count3 numbers of O, count4 S
         // Remember - this is a aromatic ring
         
         group="Compl_"+count1+"_N"+count2+"_O"+count3+"_S"+count4+NHstr;
         }
        }                                            // =========== if ends ==============
     d=stringparts[3];
     xt=parseInt(1000*stringparts[0]/d)/1000;
     yt=parseInt(1000*stringparts[1]/d)/1000;
     zt=parseInt(1000*stringparts[2]/d)/1000;
      if(group!="") collected_string=collected_string+" "+group+"  "+ xt+ " "+yt+" "+zt+"#";
       }

     collected_string=collected_string.substring(0,collected_string.length-1);

return collected_string;

}
// ############################### end func ############################################