%test_ssi
for m = [1 7]
   for d = [10 20]
      for t = [ 9 12]
         for L = [ 0 40]
            sdn = datenum(1998, m, d, t, 0, 0);
            Z = solarzenith(sdn, L, 0);
            E = ssi(sdn, L, 0);
            E2 = frouin(sdn, L, 0);
            fprintf('%4i %02i %02i %6.3f %6.3f %6.3f %4.1f %6.3f %6.3f %6.3f %6.3f %6.3f %6.3f\n', ...
               1998, m , d, t, L, 0, Z, E(1), E(2), E(3), E(4), E(5), E(6));
            fprintf('%4i %02i %02i %6.3f %6.3f %6.3f %4.1f %6.3f %6.3f %6.3f %6.3f %6.3f %6.3f\n', ...
               1998, m , d, t, L, 0, Z, E2(1), E2(2), E2(3), E2(4), E2(5), E2(6));
         end
      end
   end
end

   