volatile uint8_t cnt=0,s0=1,s1=1;　//　グローバル変数を指定します。

　/*    タイマー０の割り込みルーチン  */
SIGNAL(SIG_OVERFLOW0) {     /*   SIG_OVERFLOW0はデバイスにより違います。シンボルを探してください。*/
  TCNT0=216;                /* 4MHz/(1024*40) = 100Hzを作ります。タイマーのカウンタセット　TCNT0はデバイスにより違います。　*/
  s0=s1;                    // swを調べる　新データを旧データ変数に入れて
  s1=(PIND & 1);            // 新データを読み込む　bit0にswがあるから off=1 on=0　bit3だったら　PIND & 0b00001000 とする。
  if(s0>s1){                //　旧＝1、新＝0　のとき（押された）
    cnt++;                  //　swのカウンタをインクリメント
    if(cnt==4){cnt=0;}      // 3までにする　変更でそれ以上にもできる
  }
}

/* メインルーチン内の設定 */
  DDRD=0b11111110;          /* PORTD スイッチポートは入力モード */
  PORTD=0b00000001;         /* PORTD　スイッチポートは pull up    */
  TCCR0=5;                  /* タイマ0の分周比＝1/1024    */
  TIMSK=(1<<TOIE0);         /* タイマ0　割り込み許可   */
  TCNT0=216;                /* タイマーのカウンタセット */
  sei();                    /* 割り込み総元許可 */

/**************************************************************************************************************
 pwm4.c  ATmega48 高速PWMの実験  050815 050814  im
 
 ２つの８ビットタイマ/カウンタ0,2を使って高速のPWM信号を４ポート( OC0A,OC0B,OC2A,OC2B )に取り出す。
 
 クロック８MHzで分周比１なら31.25kHzのPWM周波数が得られる。１サイクル約30μsである。 8分周でも、3.9kHzとなる。
 
 この例では、１つのLEDを 暗−明−暗 と等比級数的に変化させ（視覚にあわせるため）、かつ、４つのLEDがずれて点滅
 するように設定してある。
 
 メインルーチンでセンサによるコントロールを行えばバリエーションが可能になる。
 
***************************************************************************************************************/

#include <avr/io.h>

void wait(void){                            /* 変化を遅くするための時間稼ぎ */
  uint16_t j,k;
  for(j=1;j<60000;j++){k=j;}
}

int main (void) {                           /* メインルーチン */
  uint8_t a1,a2,a3,a4;
  uint8_t d[57]={0,1,2,3,3,4,4,5,6,7,9,11,13,15,18,22,27,32,38,46,55,66,
                 79,95,114,137, 164,197,236,197,164,137,114,95,79,
                 66,55,46,38,32,27,22,18,15,13,11,9,7,6,5,4,4,3,3,2,1,0}; /* 明るさ変化のため等比数列 */
  
  DDRB=0xff;                                /* 全ポートを出力としたが、必要なのは PD6,PD5,PB3,PD3 である */
  DDRC=0xff;                                /*                〃                                         */
  DDRD=0xff;                                /*                〃                                         */
  
  TCCR0B=0x01;                              /* 1=分周無し 2=8分周 3=64分周 4=256分周 5=1024分周 */
  TCCR0A=0b10100011;                        /* 1010=ABとも非反転  00=予約  11=高速PWM(TCCR0BのWGM02=0と) */
  TCCR2B=0x01;                              /* 1=分周無し 2=8分周 3=32分周 4=64分周 5=128分周 */
  TCCR2A=0b10100011;
  
  a1=1;                                      /* LEDの点滅をずらすためのオフセット */
  a2=5;
  a3=10;
  a4=15;
  
  for(;;){                                   /* メインの無限ループ PWMに関係なく自由に使える */
    wait();wait();                           /* 時間稼ぎ */
    a1++;if(a1==57){a1=0;}                   /* 明るさ番号をインクリメント */
    a2++;if(a2==57){a2=0;}
    a3++;if(a3==57){a3=0;}
    a4++;if(a4==57){a4=0;}
    
    OCR0A=d[a1];                                /* 明るさ番号に応じたデューティ比較値を設定 */
    OCR0B=d[a2];
    OCR2A=d[a3];
    OCR2B=d[a4];
  }
}
      /* if(PORTD&_BV(4)){PORTD&=~_BV(4);}else{PORTD|=_BV(4);} トグル出力・ポートチェック用 */

/**********************************************************************************
  pwmh1.c     tiny2313  高速PWM   0508015 im

tiny2313の16ビットカウンタを 8ビット高速PWMとして使用する。

クロックは内蔵RC８MHz。PWMクロックは1/8分周。したがってPWM周波数は3.9kHz。

TCCR1A: COM1A1 COM1A0 COM1B1 COM1B0 ------ ------ WGM11 WGM10     1010 00 01     1010=一致でlow(AB)
TCCR1B: ICNC1  ICES1  -----  WGM13  WGM12  CS12   CS11  CS10      000 00 010     010=1/8分周
TCCR1C: FOC1A  FOC1B  -      -      -      -      -     -         0000 0000
 （WGM13  WGM12  WGM11 WGM10 = 0001 ・・8ビットPWM動作）

TCCR0A: COM0A1 COM0A0 COM0B1 COM0B0 -      -      WGM01 WGM00     1010 00 11     1010=一致でlow(AB)
TCCR0B: FOC0A  FOC0B  -      -      WGM02  CS02   CS01  CS00      00 00 0 010    010=1/8分周
  (WGM=011 高速PWM)

*********************************************************************************** */

#include <avr/io.h>

void wait(void){                            /* 変化を遅くするための時間稼ぎ */
  uint16_t j,k;
  for(j=1;j<60000;j++){k=j;}
}

/*   メインルーチン   */
int main (void) {
  uint8_t a1,a2,a3,a4;
  uint8_t d[57]={0,1,2,3,3,4,4,5,6,7,9,11,13,15,18,22,27,32,38,46,55,66,
                 79,95,114,137, 164,197,236,197,164,137,114,95,79,
                 66,55,46,38,32,27,22,18,15,13,11,9,7,6,5,4,4,3,3,2,1,0}; /* 明るさ変化のため等比数列 */
  
  DDRB=0xff;                                 /* 必要なのは PB2 PB3 PB4 */  
  DDRD=0xff;                                 /* 必要なのは PD5 */  

  TCCR1A=0b10100001;
  TCCR1B=0b00000010;
  TCCR1C=0b00000000;
  TCCR0A=0b10100011;
  TCCR0B=0b00000010;
  
  a1=1;                                      /* LEDの点滅をずらすためのオフセット */
  a2=5;
  a3=10;
  a4=15;

  for(;;) {                                  /* 無限ループ  この中でデューティの値k1〜k4を自由に決められる */
    wait();wait();                           /* 時間稼ぎ */
    a1++;if(a1==57){a1=0;}                   /* 明るさ番号をインクリメント */
    a2++;if(a2==57){a2=0;}
    a3++;if(a3==57){a3=0;}
    a4++;if(a4==57){a4=0;}
    
    OCR1A=d[a1];                              /* 明るさ番号に応じたデューティ比較値を設定 */
    OCR1B=d[a2];
    OCR0A=d[a3];
    OCR0B=d[a4];
    
  }/* forの終わり */
}
/*   メインルーチン終わり   */

MAXIM　RTC　DS3234を基準時間とした周波数カウンタのバラックが完成しましたので報告します。

このRTCはDS3231のSPI版で秋月でも取り扱っています。
DS3231とほぼ同じ構成で、１HzをSQW端子から取り出してゲート時間を作りました。

16MHzの水晶発振器の４分周出力を測定したところ、
秋月の12.8MHz基準のカウンタで　3.999954〜55Hz
DS3231基準のカウンタで　　　　 3.999942〜43Hz
DS3234基準のカウンタで　　　　 3.999945〜46Hz
となりました。

秋月12.8MHzを基準にしますと、
DS3231　は　-12Hz　（3ppmクロックが速い）
DS3234　は　- 9Hz　（2ppm強クロックが速い）
となります。

真の値は不明ですが、これらの中間にあるとすれば±2ppm程度の正確さで使えるのではないでしょうか。
素人の基準を持たないケースではこの程度の誤差は問題になることはないようですから、いずれを基準に
しても実用になろうかと思います。
一般の水晶は50ppm程度と聞きますから一桁あまり正確であると言えるでしょうか。

なお、このDS3231は長時間の平均誤差が1.5ppm（クロックが速い）と報告した時計の物とは別の個体です。