for(t2 = 0; t2 <= 4*an-4; t2++) {
  for(t4 = 0; t4 <= bm-1; t4++) {
    if (an >= t2+1) {
      s0(t2,t4,0);
      if (intMod(t2,4) == 0) {
        s1(t2/4,t4,0);
      }
    }
    else {
      if (intMod(t2,4) == 0) {
        s1(t2/4,t4,0);
      }
    }
    if (intMod(t2,4) == 0) {
      for(t6 = 1; t6 <= ambn-1; t6++) {
        s1(t2/4,t4,t6);
      }
    }
  }
}