prep_years_data.js

const SunCalc = require("./views/suncalc")
d2r = Math.PI / 180;
r2d = 180 / Math.PI;
//var zn = ["Mesha", "Vrushabha", "Mithuna", "Karkataka", "Simha", "Kanya", "Tula", "Vrushchika", "Dhanu", "Makara", "Kumbha", "Meena"];
var zn = ["మేషా", "వృషభ", "మిథున", "కర్కాటక", "సింహా", "కన్యా", "తుల", "వృశ్చిక", "ధను", "మకర", "కుంభ", "మీన"];
var wd = ["ఆదివారం", "సోమవారం", "మంగళవారం", "బుధవారం", "గురువారం", "శుక్రవారం", "శనివారం"];
var range = [1, 31, 0, 0, -3000, 4000, 0, 23, 0, 59, -12, 12, 0, 59];
var naks = ["అశ్విని", "భరణి", "కృతిక", "రోహిణి", "మృగశిర", "ఆరుద్ర", "పునర్వసు", "పుష్యమి", "ఆశ్లేష", "మఖా", "పుభ", "ఉత్తర", "హస్త", "చిత్తా", "స్వాతి", "విశాఖ", "అనురాధ", "జ్యేస్టా", "మూల", "పూర్వాషాడ", "ఉత్తరాషాడ", "శ్రవణం", "ధనిష్ట", "శతభిషం", "పూర్వాభద్ర", "ఉత్తరాభద్ర", "రేవతి"];
var tith = ["పాడ్యమి", "విధియా", "తదియ", "చవితి", "పంచమి", "షష్టి", "సప్తమి", "అష్టమి", "నవమి", "దశమి", "ఏకాదశి", "ద్వాదశి", "త్రయోదశి", "చతుర్దశి", "పౌర్ణమి", "పాడ్యమి", "విధియా", "తదియ", "చవితి", "పంచమి", "షష్టి", "సప్తమి", "అష్టమి", "నవమి", "దశమి", "ఏకాదశి", "ద్వాదశి", "త్రయోదశి", "చతుర్దశి", "అమావాస్య"];
var kar = ["బావా", "బాలవ", "కౌలవ", "తైతులా", "గరజ", "వనిజ", "విష్టి", "శకునా", "చతుష్పాద", "నాగవ", "కింస్తుఘన"];
var yog = ["విష్కంభ", "పృథి", "ఆయుష్మాన్", "సౌభాగ్య", "శోభన", "అతిగండ", "సుకర్మాన్", "ధృతి", "సూల","గండ", "వృద్ధి", "ధృవ", "వ్యాఘాత", "హర్షణ", "వజ్ర", "సిద్ధి", "వ్యతిపాత", "వరియన్", "పరిఘ", "శివ","సిద్ధ", "సాధ్య", "శుభ", "శుక్ల", "బ్రహ్మ", "ఇంద్ర", "వైధృతి"];
var tipnaks = [2, 5, 6, 0, 1, 4, 3, 2, 4, 5, 5, 0, 2, 1, 3, 6, 1, 4, 4, 5, 0, 3, 3, 3, 5, 0, 1];
var Lmoon, Lsun, skor, LmoonYoga, LsunYoga, dt;
var ayanamsa = 0;
var masa = ["చైత్ర", "వైశాఖ", "జ్యయిష్ఠ", "ఆషాఢ", "శ్రావణ", "భాద్రపద", "ఆశ్వయుజ", "కార్తీక", "మార్గశీర్ష", "పుష్య", "మాఘ", "ఫాల్గుణ"];
var masa_en = ["Chaitra", "Vaishaakha", "Jyaishtha", "Aashaadha", "Shraavana", "Bhaadrapada", "Aashvayuja", "Kaartika", "Maargashiirsha", "Pushya", "Maagha", "Phaalguna"];
var samvatsara = ["ప్రభవ", "విభవ", "శుక్ల", "ప్రమోద", "ప్రజోత్పత్తి", "ఆంగిరస", "శ్రీముఖ", "భావ", "యువన్", "ధాతు", "ఈశ్వర", "బహుధాన్య", "ప్రమతిన్", "విక్రమ", "వృష", "చిత్రభాను", "స్వభాను", "తారణ", "పార్థివ", "వ్యాయ", "సర్వజిత్", "సర్వధారిన్", "విరోధిన్", "వికృత", "ఖర", "నందన", "విజయ", "జయ", "మన్మథ", "దుర్ముఖ", "హేమలంబి", "విళంబి", "వికారిన్", "శార్వరి", "ప్లావ", "శుభకృత్", "శోభకృతు", "క్రోధిన్","విశ్వావసు", "పరాభవ", "ప్లవంగ", "కీలక", "సౌమ్య", "సాధారణ", "విరోధీకృత", "పరితాపిన్", "ప్రమాదిన్", "ఆనంద", "రాక్షస", "నల", "పింగళ", "కాళయుక్తి", "సిద్ధార్థిన్", "రౌద్ర", "దుర్మతి", "దుందుభి", "రుధిరోద్గారిన్", "రక్తాక్ష", "క్రోధన", "అక్షయ"];
function isNewMoon(_xdt, logFlag = false) {
   var dt = new Date(_xdt);
   dt.setHours(0); dt.setMinutes(0); dt.setSeconds(0);
   let moonI = 0.0140;
   for (var i = 0; i < 60 * 24 - 1; i++) {
     dt.setMinutes(dt.getMinutes() + 1);
     let x = SunCalc.getMoonIllumination(dt).fraction;
     if (logFlag) {
       console.log("moon illumination", dt, x);
     }
     if (x <= moonI) {
       return true;
     }
   }

   return false;
}

 function isFullMoon(_xdt, logFlag = false) {
   var dt = new Date(_xdt);
   dt.setHours(0); dt.setMinutes(0); dt.setSeconds(0);
   let moonI = 0.99;
   for (var i = 0; i < 60 * 24 - 1; i++) {
     dt.setMinutes(dt.getMinutes() + 1);
     let x = SunCalc.getMoonIllumination(dt).fraction;
     
     if (moonI > x) {
       moonI = x;
     }
   }

   if (logFlag) {
       console.log("moon illumination", dt, moonI);
   }
   
   return false;
 }
 function getPrevNewMoon(_curDate) {
   var newMoonFound = false;
   var d = new Date(_curDate);
   while (!newMoonFound) {
     if (isNewMoon(d)) {
       newMoonFound = true;
     } else {
       d.setDate(d.getDate() - 1);
     }
   }
   return d;
 }

 function getNextNewMoon(_curDate, logFlag=false) {
   var newMoonFound = false;
   var d = new Date(_curDate);
   while (!newMoonFound) {
     if (isNewMoon(d, logFlag)) {
       newMoonFound = true;
     } else {
       d.setDate(d.getDate() + 1);
     }
   }
//    var prevMonth = d.getMonth() - 1;
//    if (prevMonth == 0 || prevMonth == 2 || prevMonth == 4 || prevMonth == 6 || prevMonth == 7 || prevMonth == 9 || prevMonth == 11) {
//     d.setDate(d.getDate() + 1);
//    }

   if (d.getHours() > 12) {
    d.setDate(d.getDate() + 1);
   }
   return d;
 }
//---------------------------------------------------------------------------
// Data on the Moon outrage in length.
//---------------------------------------------------------------------------
function corr(mlcor, mscor, fcor, dcor, lcor) {
   this.mlcor = mlcor;
   this.mscor = mscor;
   this.fcor = fcor;
   this.dcor = dcor;
   this.lcor = lcor;
}

function corr2(l, ml, ms, f, d) {
   this.l = l;
   this.ml = ml;
   this.ms = ms;
   this.f = f;
   this.d = d;
}

var corrMoon = new Array(); // main
i = 0;
// ml, ms, f, d, l
corrMoon[i++] = new corr(0, 0, 0, 4, 13.902);
corrMoon[i++] = new corr(0, 0, 0, 2, 2369.912);
corrMoon[i++] = new corr(1, 0, 0, 4, 1.979);
corrMoon[i++] = new corr(1, 0, 0, 2, 191.953);
corrMoon[i++] = new corr(1, 0, 0, 0, 22639.500);
corrMoon[i++] = new corr(1, 0, 0, -2, -4586.465);
corrMoon[i++] = new corr(1, 0, 0, -4, -38.428);
corrMoon[i++] = new corr(1, 0, 0, -6, -0.393);
corrMoon[i++] = new corr(0, 1, 0, 4, -0.289);
corrMoon[i++] = new corr(0, 1, 0, 2, -24.420);
corrMoon[i++] = new corr(0, 1, 0, 0, -668.146);
corrMoon[i++] = new corr(0, 1, 0, -2, -165.145);
corrMoon[i++] = new corr(0, 1, 0, -4, -1.877);
corrMoon[i++] = new corr(0, 0, 0, 3, 0.403);
corrMoon[i++] = new corr(0, 0, 0, 1, -125.154);
corrMoon[i++] = new corr(2, 0, 0, 4, 0.213);
corrMoon[i++] = new corr(2, 0, 0, 2, 14.387);
corrMoon[i++] = new corr(2, 0, 0, 0, 769.016);
corrMoon[i++] = new corr(2, 0, 0, -2, -211.656);
corrMoon[i++] = new corr(2, 0, 0, -4, -30.773);
corrMoon[i++] = new corr(2, 0, 0, -6, -0.570);
corrMoon[i++] = new corr(1, 1, 0, 2, -2.921);
corrMoon[i++] = new corr(1, 1, 0, 0, -109.673);
corrMoon[i++] = new corr(1, 1, 0, -2, -205.962);
corrMoon[i++] = new corr(1, 1, 0, -4, -4.391);
corrMoon[i++] = new corr(1, -1, 0, 4, 0.283);
corrMoon[i++] = new corr(1, -1, 0, 2, 14.577);
corrMoon[i++] = new corr(1, -1, 0, 0, 147.687);
corrMoon[i++] = new corr(1, -1, 0, -2, 28.475);
corrMoon[i++] = new corr(1, -1, 0, -4, 0.636);
corrMoon[i++] = new corr(0, 2, 0, 2, -0.189);
corrMoon[i++] = new corr(0, 2, 0, 0, -7.486);
corrMoon[i++] = new corr(0, 2, 0, -2, -8.096);
corrMoon[i++] = new corr(0, 0, 2, 2, -5.741);
corrMoon[i++] = new corr(0, 0, 2, 0, -411.608);
corrMoon[i++] = new corr(0, 0, 2, -2, -55.173);
corrMoon[i++] = new corr(0, 0, 2, -4, 0.025);
corrMoon[i++] = new corr(1, 0, 0, 1, -8.466);
corrMoon[i++] = new corr(1, 0, 0, -1, 18.609);
corrMoon[i++] = new corr(1, 0, 0, -3, 3.215);
corrMoon[i++] = new corr(0, 1, 0, 1, 18.023);
corrMoon[i++] = new corr(0, 1, 0, -1, 0.560);
corrMoon[i++] = new corr(3, 0, 0, 2, 1.060);
corrMoon[i++] = new corr(3, 0, 0, 0, 36.124);
corrMoon[i++] = new corr(3, 0, 0, -2, -13.193);
corrMoon[i++] = new corr(3, 0, 0, -4, -1.187);
corrMoon[i++] = new corr(3, 0, 0, -6, -0.293);
corrMoon[i++] = new corr(2, 1, 0, 2, -0.290);
corrMoon[i++] = new corr(2, 1, 0, 0, -7.649);
corrMoon[i++] = new corr(2, 1, 0, -2, -8.627);
corrMoon[i++] = new corr(2, 1, 0, -4, -2.740);
corrMoon[i++] = new corr(2, -1, 0, 2, 1.181);
corrMoon[i++] = new corr(2, -1, 0, 0, 9.703);
corrMoon[i++] = new corr(2, -1, 0, -2, -2.494);
corrMoon[i++] = new corr(2, -1, 0, -4, 0.360);
corrMoon[i++] = new corr(1, 2, 0, 0, -1.167);
corrMoon[i++] = new corr(1, 2, 0, -2, -7.412);
corrMoon[i++] = new corr(1, 2, 0, -4, -0.311);
corrMoon[i++] = new corr(1, -2, 0, 2, 0.757);
corrMoon[i++] = new corr(1, -2, 0, 0, 2.580);

corrMoon[i++] = new corr(1, -2, 0, -2, 2.533);
corrMoon[i++] = new corr(0, 3, 0, -2, -0.344);
corrMoon[i++] = new corr(1, 0, 2, 2, -0.992);
corrMoon[i++] = new corr(1, 0, 2, 0, -45.099);
corrMoon[i++] = new corr(1, 0, 2, -2, -0.179);
corrMoon[i++] = new corr(1, 0, -2, 2, -6.382);
corrMoon[i++] = new corr(1, 0, -2, 0, 39.528);
corrMoon[i++] = new corr(1, 0, -2, -2, 9.366);
corrMoon[i++] = new corr(0, 1, 2, 0, 0.415);
corrMoon[i++] = new corr(0, 1, 2, -2, -2.152);
corrMoon[i++] = new corr(0, 1, -2, 2, -1.440);
corrMoon[i++] = new corr(0, 1, -2, -2, 0.384);
corrMoon[i++] = new corr(2, 0, 0, 1, -0.586);
corrMoon[i++] = new corr(2, 0, 0, -1, 1.750);
corrMoon[i++] = new corr(2, 0, 0, -3, 1.225);
corrMoon[i++] = new corr(1, 1, 0, 1, 1.267);
corrMoon[i++] = new corr(1, -1, 0, -1, -1.089);
corrMoon[i++] = new corr(0, 0, 2, -1, 0.584);
corrMoon[i++] = new corr(4, 0, 0, 0, 1.938);
corrMoon[i++] = new corr(4, 0, 0, -2, -0.952);
corrMoon[i++] = new corr(3, 1, 0, 0, -0.551);
corrMoon[i++] = new corr(3, 1, 0, -2, -0.482);
corrMoon[i++] = new corr(3, -1, 0, 0, 0.681);
corrMoon[i++] = new corr(2, 0, 2, 0, -3.996);
corrMoon[i++] = new corr(2, 0, 2, -2, 0.557);
corrMoon[i++] = new corr(2, 0, -2, 2, -0.459);
corrMoon[i++] = new corr(2, 0, -2, 0, -1.298);
corrMoon[i++] = new corr(2, 0, -2, -2, 0.538);
corrMoon[i++] = new corr(1, 1, -2, -2, 0.426);
corrMoon[i++] = new corr(1, -1, 2, 0, -0.304);
corrMoon[i++] = new corr(1, -1, -2, 2, -0.372);
corrMoon[i++] = new corr(0, 0, 4, 0, 0.418);
corrMoon[i++] = new corr(2, -1, 0, -1, -0.352);


var corrMoon2 = new Array(); // additional
i = 0;
// l, ml, ms, f, d
corrMoon2[i++] = new corr2(0.127, 0, 0, 0, 6);
corrMoon2[i++] = new corr2(-0.151, 0, 2, 0, -4);
corrMoon2[i++] = new corr2(-0.085, 0, 0, 2, 4);
corrMoon2[i++] = new corr2(0.150, 0, 1, 0, 3);
corrMoon2[i++] = new corr2(-0.091, 2, 1, 0, -6);
corrMoon2[i++] = new corr2(-0.103, 0, 3, 0, 0);
corrMoon2[i++] = new corr2(-0.301, 1, 0, 2, -4);
corrMoon2[i++] = new corr2(0.202, 1, 0, -2, -4);
corrMoon2[i++] = new corr2(0.137, 1, 1, 0, -1);
corrMoon2[i++] = new corr2(0.233, 1, 1, 0, -3);
corrMoon2[i++] = new corr2(-0.122, 1, -1, 0, 1);
corrMoon2[i++] = new corr2(-0.276, 1, -1, 0, -3);
corrMoon2[i++] = new corr2(0.255, 0, 0, 2, 1);
corrMoon2[i++] = new corr2(0.254, 0, 0, 2, -3);
corrMoon2[i++] = new corr2(-0.100, 3, 1, 0, -4);
corrMoon2[i++] = new corr2(-0.183, 3, -1, 0, -2);
corrMoon2[i++] = new corr2(-0.297, 2, 2, 0, -2);
corrMoon2[i++] = new corr2(-0.161, 2, 2, 0, -4);
corrMoon2[i++] = new corr2(0.197, 2, -2, 0, 0);
corrMoon2[i++] = new corr2(0.254, 2, -2, 0, -2);
corrMoon2[i++] = new corr2(-0.250, 1, 3, 0, -2);
corrMoon2[i++] = new corr2(-0.123, 2, 0, 2, 2);
corrMoon2[i++] = new corr2(0.173, 2, 0, -2, -4);
corrMoon2[i++] = new corr2(0.263, 1, 1, 2, 0);
corrMoon2[i++] = new corr2(0.130, 3, 0, 0, -1);
corrMoon2[i++] = new corr2(0.113, 5, 0, 0, 0);
corrMoon2[i++] = new corr2(0.092, 3, 0, 2, -2);



function daysInMonth(m, y) {
   var g_days = new Array(31, 28, 31, 30, 31, 30, 31, 31, 30, 31, 30, 31);
   if (((y % 4 == 0) && (y % 100 != 0)) || (y % 400 == 0)) g_days[1] = 29;
   return g_days[m];
}


//-----------------------------------------------------------------------------------
// Calculating geotsent p avoid longitude Moon and angular sector p News.
// (2 sec accuracy. longitude)
//-----------------------------------------------------------------------------------
function moon(jd) {
   // days from 1900
   tdays = jd - 2415020;
   t = tdays / 36525;
   t2 = t * t;
   t3 = t * t * t;

   // slope travels to the equator
   ob = 23.452294 - 0.0130125 * t - 0.00000164 * t2 + 0.000000503 * t3;
   // the average length moon
   l = 270.4337361 + 13.176396544528099 * tdays - 5.86 * t2 / 3600 + 0.0068 * t3 / 3600;
   // the difference medium length Moon and the Sun (the averageElongation Moon):
   d = 350.7374861110581 + 445267.1142166667 * t - t2 * 1.436111132303874e-3 + 0.0000018888889 * t3;
   // Perigee moon
   pe = 334.329556 + 14648522.52 * t / 3600 - 37.17 * t2 / 3600 - 0.045 * t3 / 3600;
   // the average anomoly sun
   ms = 358.4758333333334 + 35999.04974999958 * t - t2 * 1.500000059604645e-4 - t3 * 3.3333333623078e-6;
   // The average anomoloy moon
   //ml = 296.1046083333757 + 477198.8491083336*t + 0.0091916667090522*t2 + 0.0000143888893*t3;
   ml = fix360(l - pe);
   // Rising length node orbit the moon:
   om = 259.183275 - 6962911.23 * t / 3600 + 7.48 * t2 / 3600 + 0.008 * t3 / 3600;
   // the average length Moon, measured from the bottom up hub orbit:

   f = fix360(l - om);

   with (Math) {
      // periodic revisions
      r2rad = 360.0 * d2r;
      tb = tdays * 1e-12; // *10^12
      t2c = tdays * tdays * 1e-16; // *10^16
      a1 = sin(r2rad * (0.53733431 - 10104982 * tb + 191 * t2c));
      a2 = sin(r2rad * (0.71995354 - 147094228 * tb + 43 * t2c));
      c2 = cos(r2rad * (0.71995354 - 147094228 * tb + 43 * t2c));
      a3 = sin(r2rad * (0.14222222 + 1536238 * tb));
      a4 = sin(r2rad * (0.48398132 - 147269147 * tb + 43 * t2c));
      c4 = cos(r2rad * (0.48398132 - 147269147 * tb + 43 * t2c));
      a5 = sin(r2rad * (0.52453688 - 147162675 * tb + 43 * t2c));
      a6 = sin(r2rad * (0.84536324 - 11459387 * tb));
      a7 = sin(r2rad * (0.23363774 + 1232723 * tb + 191 * t2c));
      a8 = sin(r2rad * (0.58750000 + 9050118 * tb));
      a9 = sin(r2rad * (0.61043085 - 67718733 * tb));

      dlm = 0.84 * a3 + 0.31 * a7 + 14.27 * a1 + 7.261 * a2 + 0.282 * a4 + 0.237 * a6;
      dpm = -2.1 * a3 - 2.076 * a2 - 0.840 * a4 - 0.593 * a6;
      dkm = 0.63 * a3 + 95.96 * a2 + 15.58 * a4 + 1.86 * a5;
      dls = -6.4 * a3 - 0.27 * a8 - 1.89 * a6 + 0.20 * a9;
      dgc = (-4.318 * c2 - 0.698 * c4) / 3600.0 / 360.0;
      dgc = (1.000002708 + 139.978 * dgc);

      ml = d2r * (ml + (dlm - dpm) / 3600.0); //Average anomoly moon
      ms = d2r * (ms + dls / 3600.0); //Average anomoly sun
      f = d2r * (f + (dlm - dkm) / 3600.0);
      d = d2r * (d + (dlm - dls) / 3600.0); //avg elongation moon

      lk = 0; lk1 = 0; sk = 0; sinp = 0; nib = 0; g1c = 0;
      i1corr = 1.0 - 6.8320e-8 * tdays;
      i2corr = dgc * dgc;

      for (i = 0; i < 93; i++) { // outrage at length
         arg = corrMoon[i].mlcor * ml + corrMoon[i].mscor * ms + corrMoon[i].fcor * f + corrMoon[i].dcor * d;
         sinarg = sin(arg);
         if (corrMoon[i].mscor != 0) {
            sinarg *= i1corr;
            if (corrMoon[i].mscor == 2 || corrMoon[i].mscor == -2) sinarg *= i1corr;
         }
         if (corrMoon[i].fcor != 0) sinarg *= i2corr;
         lk += corrMoon[i].lcor * sinarg;
      }
      for (i = 0; i < 27; i++) { // outrage at length additional
         arg = corrMoon2[i].ml * ml + corrMoon2[i].ms * ms + corrMoon2[i].f * f + corrMoon2[i].d * d;
         sinarg = sin(arg);
         lk1 += corrMoon2[i].l * sinarg;
      }

      // resentments of the planets
      dlid = 0.822 * sin(r2rad * (0.32480 - 0.0017125594 * tdays));
      dlid += 0.307 * sin(r2rad * (0.14905 - 0.0034251187 * tdays));
      dlid += 0.348 * sin(r2rad * (0.68266 - 0.0006873156 * tdays));
      dlid += 0.662 * sin(r2rad * (0.65162 + 0.0365724168 * tdays));
      dlid += 0.643 * sin(r2rad * (0.88098 - 0.0025069941 * tdays));
      dlid += 1.137 * sin(r2rad * (0.85823 + 0.0364487270 * tdays));
      dlid += 0.436 * sin(r2rad * (0.71892 + 0.0362179180 * tdays));
      dlid += 0.327 * sin(r2rad * (0.97639 + 0.0001734910 * tdays));

      l = l + nutation(jd) + (dlm + lk + lk1 + dlid) / 3600.0;
      LmoonYoga = l;
      //alert("Lmoon="+l);
      l = fix360(l);

      // angular velocity of the moon on ecliptic (deg/day):
      vl = 13.176397;
      vl = vl + 1.434006 * cos(ml);
      vl = vl + .280135 * cos(2 * d);
      vl = vl + .251632 * cos(2 * d - ml);
      vl = vl + .09742 * cos(2 * ml);
      vl = vl - .052799 * cos(2 * f);
      vl = vl + .034848 * cos(2 * d + ml);
      vl = vl + .018732 * cos(2 * d - ms);
      vl = vl + .010316 * cos(2 * d - ms - ml);
      vl = vl + .008649 * cos(ms - ml);
      vl = vl - .008642 * cos(2 * f + ml);
      vl = vl - .007471 * cos(ms + ml);
      vl = vl - .007387 * cos(d);
      vl = vl + .006864 * cos(3 * ml);
      vl = vl + .00665 * cos(4 * d - ml);
      vl = vl + .003523 * cos(2 * d + 2 * ml);
      vl = vl + .003377 * cos(4 * d - 2 * ml);
      vl = vl + .003287 * cos(4 * d);
      vl = vl - .003193 * cos(ms);
      vl = vl - .003003 * cos(2 * d + ms);
      vl = vl + .002577 * cos(ml - ms + 2 * d);
      vl = vl - .002567 * cos(2 * f - ml);
      vl = vl - .001794 * cos(2 * d - 2 * ml);
      vl = vl - .001716 * cos(ml - 2 * f - 2 * d);
      vl = vl - .001698 * cos(2 * d + ms - ml);
      vl = vl - .001415 * cos(2 * d + 2 * f);
      vl = vl + .001183 * cos(2 * ml - ms);
      vl = vl + .00115 * cos(d + ms);
      vl = vl - .001035 * cos(d + ml);
      vl = vl - .001019 * cos(2 * f + 2 * ml);
      vl = vl - .001006 * cos(ms + 2 * ml);
   }
   skor = vl;
   //l += ay;
   //if(l < 0.0)l += 360.0;
   return l;
}


//----------------------------------------------------------------------
// Calculating geotsent p avoid longitude Sun.
// (the acuracy of 1 sec . longitude)
//----------------------------------------------------------------------
function sun(jd) {
   // days frm 1900:
   tdays = jd - 2415020;

   t = tdays / 36525;
   t2 = t * t;
   t3 = t * t * t;

   // the avg len sun

   ls = 279.696678 + 0.9856473354 * tdays + 1.089 * t2 / 3600;
   // perigee sun
   pes = 101.220833 + 6189.03 * t / 3600 + 1.63 * t2 / 3600 + 0.012 * t3 / 3600;
   // avg anomoly sun

   ms = fix360(ls - pes + 180);
   g = ms + (0.266 * Math.sin((31.8 + 119.0 * t) * d2r) + 6.40 * Math.sin((231.19 + 20.2 * t) * d2r) + (1.882 - 0.016 * t) * Math.sin((57.24 + 150.27 * t) * d2r)) / 3600.0;
   // Rising sun node len
   oms = 259.18 - 1934.142 * t;
   // eccentricity orbit sun
   ex = 0.01675104 - 0.0000418 * t - 0.000000126 * t2;
   // avg length moon
   l = 270.4337361 + 13.176396544528099 * tdays - 5.86 * t2 / 3600 + 0.0068 * t3 / 3600;
   // avg anomaly moon
   ml = 296.1046083333757 + 477198.8491083336 * t + 0.0091916667090522 * t2 + 0.0000143888893 * t3;
   // avg len earth
   le = 99.696678 + 0.9856473354 * tdays + 1.089 * t2 / 3600;

   om = 259.183275 - 6962911.23 * t / 3600 + 7.48 * t2 / 3600 + 0.008 * t3 / 3600;

   // eccentric anomoloy calculation iterative method
   u = kepler(g, ex, 0.0000003);

   with (Math) {
      // cal true anomaly sun
      b = sqrt((1 + ex) / (1 - ex));
      if (abs(Math.PI - u) < 1.0e-10) truanom = u;
      else truanom = 2.0 * atan(b * tan(u / 2));
      truanom = fix360(truanom * r2d);

      //corrections for cal of longitude and radius vector
      u1 = (153.23 + 22518.7541 * t) * d2r;
      u2 = (216.57 + 45037.5082 * t) * d2r;
      u3 = (312.69 + 32964.3577 * t) * d2r;
      u4 = (350.74 + 445267.1142 * t - 0.00144 * t2) * d2r;
      u6 = (353.4 + 65928.71550000001 * t) * d2r;
      u5 = (315.6 + 893.3 * t) * d2r;

      dl = 0.00134 * cos(u1);
      dl += 0.00154 * cos(u2);
      dl += 0.002 * cos(u3);
      dl += 0.00179 * sin(u4);
      dl += 0.202 * sin(u5) / 3600;

      dr = 0.00000543 * sin(u1);
      dr += 0.00001575 * sin(u2);
      dr += 0.00001627 * sin(u3);
      dr += 0.00003076 * cos(u4);
      dr += 9.26999999e-06 * sin(u6);

      // true len of sun (deg)
      il = ls + dl + truanom - ms;

      // corrections to abberations links
      r1 = 1.0000002 * (1 - ex * ex) / (1 + ex * cos(truanom * d2r));
      rs = r1 + dr; // radius vector
      ab = (20.496 * (1 - ex * ex) / rs) / 3600;
      ls = il + nutation(jd) - ab; // app len sun
      LsunYoga = ls;

      ls = fix360(ls);
   }
   return ls;
}

//----------------------------------------------------------------------------
// cal start and end of tithi (len = 12)and karana (len = 6)
//----------------------------------------------------------------------------
function tithi(jd, n1, tzone, len) {
   var s_t = {};
   var flag;
   jdt = jd;
   knv = Math.floor(((jd - 2415020) / 365.25) * 12.3685);

   for (itit = n1; itit < (n1 + 2); ++itit) {
      aspect = len * itit; // sun n moon in the early tithi
      flag = 0;
      if (aspect == 0) { jdt = novolun(jd, knv); flag = 1; }
      if (aspect == 360) { jdt = novolun(jd, (knv + 1)); flag = 1; }
      while (flag < 1) {
         Lsun0 = sun(jdt);
         Lmoon0 = moon(jdt);
         a = fix360(Lsun0 + aspect); // pt should be where luna
         asp1 = a - Lmoon0; // assymptots of the moon to ur point
         if (asp1 > 180) asp1 -= 360;
         if (asp1 < -180) asp1 += 360;
         flag = 1;

         if (Math.abs(asp1) > 0.001) { jdt += (asp1 / (skor - 1)); flag = 0; }
      }
      if (itit == n1) s_t.start = calData(jdt + (tzone - dt) / 24);
      if (itit == (n1 + 1)) s_t.end = calData(jdt + (tzone - dt) / 24);
   }
   return s_t;
}


//----------------------------------------------------------------------------
// cal entry and exit moon in nakshatra
//----------------------------------------------------------------------------
function nakshatra(jd, n_naksh, tzone) {
   var s_t = {};
   var flag;
   jdt = jd;

   for (inak = n_naksh; inak < (n_naksh + 2); ++inak) {
      n1 = fix360(inak * 80 / 6); // co-ordinate start of nakshatra
      flag = 0;
      while (flag < 1) {
         Lmoon0 = fix360(moon(jdt) + ayanamsa);
         asp1 = n1 - Lmoon0; // distance frm moon before nakshatra(degree)
         if (asp1 > 180) asp1 -= 360;
         if (asp1 < -180) asp1 += 360;
         flag = 1;
         if (Math.abs(asp1) > 0.001) { jdt += (asp1 / skor); flag = 0; }
      }
      if (inak == n_naksh) s_t.start = calData(jdt + (tzone - dt) / 24);
      if (inak == (n_naksh + 1)) s_t.end = calData(jdt + (tzone - dt) / 24);
   }
   return s_t;
}


//----------------------------------------------------------------------------
// cal begin and end of yoga
//----------------------------------------------------------------------------
function yoga(jd, zyoga, tzone) {
   var s_t = {};
   var flag;
   jdt = jd;
   z = zyoga;
   var nn_yoga = new Array(2);
   nn_yoga[0] = Math.floor(z * 6 / 80) * 80 / 6;
   nn_yoga[1] = (Math.floor(z * 6 / 80) + 1) * 80 / 6;
   for (iyog = 0; iyog < 2; ++iyog) {
      flag = 0;
      while (flag < 1) {
         Lsun0 = sun(jdt);
         Lmoon0 = moon(jdt);
         dmoonYoga = (LmoonYoga + ayanamsa - 491143.07698973856);
         dsunYoga = (LsunYoga + ayanamsa - 36976.91240579201);
         //alert(LmoonYoga+"\r"+LsunYoga+"\r"+ayanamsa);
         z = dmoonYoga + dsunYoga;
         asp1 = nn_yoga[iyog] - z;
         flag = 1;
         if (Math.abs(asp1) > 0.001) { jdt += (asp1 / (skor + 1.0145616633)); flag = 0; }
         //if (Math.abs(asp1) > 0.001) {jdt += (asp1 / skor) + (58.13 * Math.sin(asp1*d2r)); flag = 0;}
      }
      if (iyog == 0) s_t.start = calData(jdt + (tzone - dt) / 24);
      if (iyog == 1) s_t.end = calData(jdt + (tzone - dt) / 24);
   }
   return s_t;
}


//-----------------------------------------------------------------------------
//cal time in the near past novoluna (err less then 2 min)
//-----------------------------------------------------------------------------
function novolun(jd, knv) {

   t = (jd - 2415020) / 36525;
   t2 = t * t;
   t3 = t * t * t;

   with (Math) {
      jdnv = 2415020.75933 + 29.53058868 * knv + 0.0001178 * t2 - 0.000000155 * t3;
      jdnv += 0.00033 * sin((166.56 + 132.87 * t - 0.009173 * t2) * d2r);
      m = 359.2242 + 29.10535608 * knv - 0.0000333 * t2 - 0.00000347 * t3;
      ml = 306.0253 + 385.81691806 * knv + 0.0107306 * t2 + 0.00001236 * t3;
      f = 21.2964 + 390.67050646 * knv - 0.0016528 * t2 - 0.00000239 * t3;
      m *= d2r;
      ml *= d2r;
      f *= d2r;

      djd = (0.1734 - 0.000393 * t) * sin(m);
      djd += 0.0021 * sin(2 * m);
      djd -= 0.4068 * sin(ml);
      djd += 0.0161 * sin(2 * ml);
      djd -= 0.0004 * sin(3 * ml);
      djd += 0.0104 * sin(2 * f);
      djd -= 0.0051 * sin(m + ml);
      djd -= 0.0074 * sin(m - ml);
      djd += 0.0004 * sin(2 * f + m);
      djd -= 0.0004 * sin(2 * f - m);
      djd -= 0.0006 * sin(2 * f + ml);
      djd += 0.001 * sin(2 * f - ml);
      djd += 0.0005 * sin(m + 2 * ml);

      jdnv += djd;
   }
   return jdnv;
}

//-----------------------------------------------------
// decision equation kepler (in rad)
//-----------------------------------------------------
function kepler(m, ex, err) {
   //val u0, delta;

   m *= d2r;
   u0 = m;
   err *= d2r;
   delta = 1;
   while (Math.abs(delta) >= err) {
      delta = (m + ex * Math.sin(u0) - u0) / (1 - ex * Math.cos(u0));
      u0 += delta;
   }
   return u0;
}

//-----------------------------------------------------
// cal nutation in len
//-----------------------------------------------------
function nutation(jd) {

   t = (jd - 2415020) / 36525;
   t2 = t * t;

   // avg len sun
   ls = 279.6967 + 36000.7689 * t + 0.000303 * t2;
   // avg len moon
   l = 270.4341639 + 481267.8831417 * t - 0.0011333333 * t2;
   // avg anomaly sun
   ms = 358.4758333333334 + 35999.04974999958 * t - t2 * 1.500000059604645e-4;
   // avg anomaly moon
   ml = 296.1046083333757 + 477198.8491083336 * t + 0.0091916667090522 * t2;
   // the diff medium len of moon and sun (avg elongation moon)
   d = 350.7374861110581 + 445267.1142166667 * t - t2 * 1.436111132303874e-3;

   om = 259.1832750002543 - 1934.142008333206 * t + .0020777778 * t2;
   ls *= d2r; l *= d2r; ms *= d2r; ml *= d2r; d *= d2r; om *= d2r;
   d2 = d * d; l2 = l * l; ls2 = ls * ls;

   with (Math) {
      nut = (-17.2327 - 0.01737 * t) * sin(om);
      nut += 0.2088 * sin(2.0 * om);
      nut += 0.0675 * sin(ml);
      nut -= 0.0149 * sin(ml - d2);
      nut -= 0.0342 * sin(l2 - om);
      nut += 0.0114 * sin(l2 - ml);
      nut -= 0.2037 * sin(l2);
      nut -= 0.0261 * sin(l2 + ml);
      nut += 0.0124 * sin(ls2 - om);
      nut += 0.0214 * sin(ls2 - ms);
      nut -= 1.2729 * sin(ls2);
      nut -= 0.0497 * sin(ls2 + ms);
      nut += 0.1261 * sin(ms);
      nut = nut / 3600.0;
   }
   return nut;
}

//-----------------------------------------------------
// Calculation ayanamsa (degrees)
//-----------------------------------------------------
function calcayan(jd) {
   t = (jd - 2415020) / 36525;
   // avg node len moon
   om = 259.183275 - 1934.142008333206 * t + 0.0020777778 * t * t + 0.0000022222222 * t * t * t;
   // avg len sun
   ls = 279.696678 + 36000.76892 * t + 0.0003025 * t * t;
   aya = 17.23 * Math.sin(d2r * om) + 1.27 * Math.sin(d2r * ls * 2) - (5025.64 + 1.11 * t) * t;
   aya = (aya - 80861.27) / 3600.0; // 84038.27 = Fagan-Bradley, 80861.27 = Lahiri

   return aya;
}

//------------------------------------------------------------------------------------------
// cal date by number of date mon and year
//------------------------------------------------------------------------------------------
function mdy2julian(m, d, y) {
   with (Math) {
      im = 12 * (y + 4800) + m - 3;
      j = (2 * (im - floor(im / 12) * 12) + 7 + 365 * im) / 12;
      j = floor(j) + d + floor(im / 48) - 32083;
      if (j > 2299171) j += floor(im / 4800) - floor(im / 1200) + 38;
      j -= 0.5;
   }
   return j;
}

function dTime(jd) {

   var efdt = [124, 85, 62, 48, 37, 26, 16, 10, 9, 10, 11, 11, 12, 13, 15, 16, 17, 17, 13.7, 12.5, 12, 7.5, 5.7, 7.1, 7.9, 1.6, -5.4, -5.9, -2.7, 10.5, 21.2, 24, 24.3, 29.2, 33.2, 40.2, 50.5, 56.9, 65.7, 75.5];
   s = calData(jd);
   dgod = kyear + (kmon - 1) / 12 + (kday - 1) / 365.25;
   t = (jd - 2378497) / 36525; // IN centuries rejection of 1800 bc
   //t = (jd - 2415020)/36525; // in cent rejection of 1900 bc
   if (dgod >= 1620 && dgod < 2010) {
      i1 = Math.floor((dgod - 1620) / 10);
      di = dgod - (1620 + i1 * 10);
      dt = (efdt[i1] + ((efdt[i1 + 1] - efdt[i1]) * di) / 10);
   }
   else {
      if (dgod >= 2010) dt = 25.5 * t * t - 39;
      //if (dgod >= 2010) dt = 29.949 * t * t - 56.796;
      //if (dgod < 1620) dt = 5 + 24.349 + 72.3165 * t + 29.949 * t * t;
      if (dgod >= 948 && dgod < 1620) dt = 25.5 * t * t;
      if (dgod < 948) dt = 1361.7 + 320 * t + 44.3 * t * t;
   }
   dt /= 3600;
   return dt;
}

//------------------------------------------------------------------------------------------
// cal date on calendar date
//------------------------------------------------------------------------------------------
function calData(jd) {
   with (Math) {
      z1 = jd + 0.5;
      z2 = floor(z1);
      f = z1 - z2;

      if (z2 < 2299161) a = z2;
      else {
         alf = floor((z2 - 1867216.25) / 36524.25);
         a = z2 + 1 + alf - floor(alf / 4);
      }

      b = a + 1524;
      c = floor((b - 122.1) / 365.25);
      d = floor(365.25 * c);
      e = floor((b - d) / 30.6001);

      days = b - d - floor(30.6001 * e) + f;
      kday = floor(days);

      if (e < 13.5) kmon = e - 1;
      else kmon = e - 13;

      if (kmon > 2.5) kyear = c - 4716;
      if (kmon < 2.5) kyear = c - 4715;

      hh1 = (days - kday) * 24;
      khr = floor(hh1);
      kmin = hh1 - khr;
      ksek = kmin * 60;
      kmin = floor(ksek);
      ksek = floor((ksek - kmin) * 60);
      s = new Date(kyear, kmon - 1, kday, khr, kmin, ksek, 0);
   }
   return s;
}

//------------------------------------------------------------------------------------------
// transalation deg logitudinal in degrees,min and sec
//------------------------------------------------------------------------------------------
function lon2dmsz(x) {
   with (Math) {
      var d, m, s;
      x = abs(x);
      z = floor(x / 30);
      d = floor(x);
      ss0 = round((x - d) * 3600);
      m = floor(ss0 / 60);
      s = (ss0 % 60) % 60;
      d %= 30;
      str = d + " " + m + "'" + s + "\" ";
   }
   return str;
}


//------------------------------------------------------------------------------------------
// translation degrees in deg, min and sec
//------------------------------------------------------------------------------------------
function lon2dms(x) {
   with (Math) {
      var d, m, s;
      x = abs(x);
      d = floor(x);
      ss0 = round((x - d) * 3600);
      m = floor(ss0 / 60);
      s = (ss0 % 60) % 60;
      str = d + " " + m + "'" + s + "\"";
   }
   return str;
}

//------------------------------------------------------------------------------------------
// fixing the angle within 360 degrees
//------------------------------------------------------------------------------------------
function fix360(v) {
   while (v < 0.0) v += 360.0;
   while (v > 360.0) v -= 360.0;
   return v;
}

//------------------------------------------------------------------------------------------
// Day of the Week
//------------------------------------------------------------------------------------------
function weekDay(jd) {
   // Julian date for the begin of the day
   jd0 = Math.floor(jd) + 0.5;
   if (jd < jd0) jd0 -= 1;

   // day
   jdn = jd0 + 1.5;
   dn1 = Math.floor(jdn / 7) * 7;


   wday = Math.floor(jdn - dn1);

   return wday;
}
// var tith = ["Padyami", "Vidhiya", "Thadiya", "Chavithi", "Panchami", "Shasti", "Sapthami", "Ashtami", "Navami", "Dasami", "Ekadasi", "Dvadasi", "Trayodasi", "Chaturdasi", "Punnami", "Padyami", "Vidhiya", "Thadiya", "Chaviti", "Panchami", "Shasti", "Sapthami", "Ashtami", "Navami", "Dasami", "Ekadasi", "Dvadasi", "Trayodasi", "Chaturdasi", "Amavasya"];
//var masa = ["Chaitra", "Vaishaakha", "Jyaishtha", "Aashaadha", "Shraavana", "Bhaadrapada", "Aashvayuja", "Kaartika", "Maargashiirsha", "Pushya", "Maagha", "Phaalguna"];
function getCurrentMonth (_dt) {
   for(var i = 0; i < panchang.months.length; i++) {
      if (_dt >= panchang.months[i].start && _dt <= panchang.months[i].end) {
         return panchang.months[i]
      }
   }
   return panchang.months[11];
}
panchang = {
   Day: {},
   Tithi: {},
   Nakshatra: {},
   Karna: {},
   Yoga: {},
   Ayanamsa: {},
   Raasi: {},
   Maasa : {},
   paksham : {},
   samvatsara : {},
   months : [],
   version: "0.2",
   calculate: function (d, cb) {
      var n_wday, n_tithi = 1, n_naksh = 1, n_karana, n_yoga;
      this.months = [];
      with (Math) {
         var day = d.getDate();
         var mon = d.getMonth() + 1;
         var year = d.getFullYear();
         var hr = d.getHours();
         hr += d.getMinutes() / 60;
         var tzone = d.getTimezoneOffset() / 60 * (-1);
      }
      inpmin = Math.floor(d.getMinutes());
      if (inpmin < 10) inpmin = "0" + inpmin;


      // Julian date in local p. LT:
      dayhr = day + hr / 24;
      jdlt = mdy2julian(mon, dayhr, year);

      // day:
      n_wday = weekDay(jdlt);
      this.Day.name = wd[n_wday];

      // julian day at the begining of the day
      jd0 = mdy2julian(mon, day, year);
      jdut = jd0 + (hr - tzone) / 24;

      dt = dTime(jdut);
      jd = jdut + dt / 24;

      //ayyanamsa
      ayanamsa = calcayan(jd);

      // length Moon
      Lmoon = moon(jd);

      // Logitudinal Sun
      Lsun = sun(jd);

      // yoga:
      dmoonYoga = (LmoonYoga + ayanamsa - 491143.07698973856);
      dsunYoga = (LsunYoga + ayanamsa - 36976.91240579201);
      zyoga = dmoonYoga + dsunYoga;
      n_yoga = zyoga * 6 / 80;
      while (n_yoga < 0) n_yoga += 27;
      while (n_yoga > 27) n_yoga -= 27;
      n3 = n_yoga;
      n_yoga = Math.floor(n_yoga);
      s_yoga = yoga(jd, zyoga, tzone);

      // Nakstra
      Lmoon0 = fix360(Lmoon + ayanamsa);
      n_naksh = Math.floor(Lmoon0 * 6 / 80);
      var s_naksh = nakshatra(jd, n_naksh, tzone);

      // tithi
      Lmoon0 = Lmoon;
      Lsun0 = Lsun;
      if (Lmoon0 < Lsun0) Lmoon0 += 360;
      n_tithi = Math.floor((Lmoon0 - Lsun0) / 12);
      var s_tithi = tithi(jd, n_tithi, tzone, 12);

      // Karana
      Lmoon0 = Lmoon;
      Lsun0 = Lsun;
      if (Lmoon0 < Lsun0) Lmoon0 += 360;
      nk = Math.floor((Lmoon0 - Lsun0) / 6);
      if (nk == 0) n_karana = 10;
      if (nk >= 57) n_karana = nk - 50;
      if (nk > 0 && nk < 57) n_karana = (nk - 1) - (Math.floor((nk - 1) / 7)) * 7;
      var s_karana = tithi(jd, nk, tzone, 6);

      with (Math) {
         var z = floor(abs(fix360(Lmoon + ayanamsa)) / 30);
      }

      this.Ayanamsa.name = lon2dms(ayanamsa);
      this.Raasi.name = zn[z];
      this.Nakshatra.name = naks[n_naksh];
      this.Nakshatra.start = s_naksh.start;
      this.Nakshatra.end = s_naksh.end;

      this.Karna.name = kar[n_karana];
      this.Karna.start = s_karana.start;
      this.Karna.end = s_karana.end;

      this.Yoga.name = yog[n_yoga];
      this.Yoga.start = s_yoga.start;
      this.Yoga.end = s_yoga.end;


      this.Tithi.name = tith[n_tithi];
      this.Tithi.start = s_tithi.start;
      this.Tithi.end = s_tithi.end;

      var cDate = new Date();
      if (cDate.getTime() >= this.Tithi.start.getTime() && cDate.getTime() <= this.Tithi.end.getTime()) {
         cDate = s_tithi.start;
      }
      var yearMod = (cDate.getFullYear() % 1927);
      if (yearMod > 60) {
         yearMod %= 60;
      } else {
         yearMod -= 60;
      }

      this.samvatsara.name = samvatsara[yearMod] ;
      this.paksham.name = n_tithi < 16 ? "శుక్ల" : "కృష్ణ" ;

      var d1 = new Date();
      d1.setDate(21); d1.setMonth(2);d1.setFullYear(2023);
      var m1 = null;
      var m2 = null;
      if (d1.getMonth() == 0 && d1.getMonth() == 1) { // jan & few ugadi will begin in previous year
         m1 = new Date(d1);
         m1.setFullYear(m1.getFullYear() - 1);
         m1.setMonth(2);
         m1.setDate(1);

         m1 = getNextNewMoon(m1);
      } else if (d1.getMonth() == 2 && d1.getDate() >= 15) { // mar and day >= 15 then new year will begin in this year
         // m1 = new Date(d1);
         // m1.setMonth(2);
        // m1.setDate(1);
         m1 = getNextNewMoon(d1);
      } else if (d1.getMonth() == 2 && d1.getDate() < 15) {
         m1 = new Date(d1);
         m1.setFullYear(m1.getFullYear() - 1);
         m1.setMonth(2);
         m1.setDate(1);

         m1 = getNextNewMoon(m1);
      } else {
         m1 = new Date(d1);
         m1.setMonth(2);
         m1.setDate(1);
         m1 = getNextNewMoon(m1);
      }
      
      m1.setDate(m1.getDate() + 1);
      console.log("first month starts at", m1);

      m2 = new Date(m1);
      m2.setDate(m2.getDate() + 1);
      m2 = getNextNewMoon(m2);

      console.log("first month ends at", m2)

      this.months.push({
         name:masa[0],
         start: m1,
         end: m2
      });

      for (var i = 1; i < 12; i++) {
         m1 = new Date(m2);
         m1.setDate(m1.getDate() + 1);

         m2 = new Date(m1);
         m2.setDate(m2.getDate() + 1);
         m2 = getNextNewMoon(m2);
         this.months.push({
            name:masa[i],
            start: m1,
            end: m2
         });
   
      }
//      console.log("current month", d1, getCurrentMonth(d1))
      if (cb) { cb(); }
   },
   init : function() {
      const fs = require('fs');
      let beginYear = 399;
      let years = [];
      let monthCount = -1;
      let d1 = new Date();
      var endYear = beginYear + 2023
      var adhikaMasaFlag = false;
      for (var y = beginYear; y <= (d1.getFullYear() + beginYear); y++) {
        adhikaMasaFlag = false;
        var curMonths = [];
         var fileName = "";
         
         if (y === beginYear) {
            d1.setFullYear(y);
            d1.setDate(01);
            d1.setMonth(02);
         }


         let s1 = getNextNewMoon(d1, false);
         if (monthCount == -1) {
            monthCount = s1.getMonth() + 1;
         }
         let s2 = new Date(s1);
         s2.setDate(s2.getDate() + 10);
         s2 = getNextNewMoon(s2)
         fileName = s1.getFullYear();
         curMonths.push({
            name : masa[0],
            start:s1,
            end:s2
         })
         monthCount++
         if (monthCount % 33 === 0) {
            adhikaMasaFlag = true;
            let astart = new Date(s2);
            astart.setDate(astart.getDate() + 1);
            astart = getNextNewMoon(astart);

            let aend = new Date(astart);
            aend.setDate(aend.getDate() + 1);
            aend = getNextNewMoon(aend);
            var adhikaMasa = {
               name : masa[0],
               start : astart,
               end : aend
            } 
            monthCount++
            curMonths.push(adhikaMasa);
            //console.log("1111 Extra month", y, adhikaMasa);
            s2 = new Date(aend);
            s2.setDate(s2.getDate() + 1)
         }

         for (var i = 1; i < 12; i++) {
            s1 = new Date(s2);
            s1.setDate(s1.getDate() + 1);
            s1 = getNextNewMoon(s1);
            s2 = new Date(s1);
            s2.setDate(s2.getDate() + 10);
            s2 = getNextNewMoon(s2);
            curMonths.push({
               name : masa[i],
               start:s1,
               end:s2
            })
            monthCount++
            if (monthCount % 33 === 0) {
               adhikaMasaFlag = true;
               let astart = new Date(s2);
               astart.setDate(astart.getDate() + 1);
               astart = getNextNewMoon(astart);

               let aend = new Date(astart);
               aend.setDate(aend.getDate() + 1);
               aend = getNextNewMoon(aend);
               var adhikaMasa = {
                  name : masa[i],
                  start : astart,
                  end : aend
               }
               //console.log("22222 Extra month", y, adhikaMasa);
               s2 = new Date(aend);
               s2.setDate(s2.getDate() + 1)
               monthCount++
               curMonths.push(adhikaMasa);             
            }
   
            if (i === 11) {
               d1 = new Date(s2);
               d1.setDate(d1.getDate() + 1)
            }
         }

        //  years.push({
        //     year : y,
        //     months : curMonths
        //  })
        let year = {
            year: s1.getFullYear(),
            months : curMonths
        };
        if (adhikaMasaFlag) {
            fileName += "-am"
        }
        console.log("processing ", fileName)
        
            let jsonContent = JSON.stringify(year);
            //console.log(jsonContent)
            fs.writeFileSync(`./data/${fileName}.json`, jsonContent, 'utf8', function (err) {
                if (err) {
                    console.log("An error occured while writing JSON Object to File.");
                    return console.log(err);
                }
                console.log("JSON file has been saved.");
            });
        
      }
      //console.log(years);
   },
   buildAdhikaMasa : function() {
    // var tith = ["Padyami", "Vidhiya", "Thadiya", "Chavithi", "Panchami", "Shasti", "Sapthami", "Ashtami", "Navami", "Dasami", "Ekadasi", "Dvadasi", "Trayodasi", "Chaturdasi", "Punnami", "Padyami", "Vidhiya", "Thadiya", "Chaviti", "Panchami", "Shasti", "Sapthami", "Ashtami", "Navami", "Dasami", "Ekadasi", "Dvadasi", "Trayodasi", "Chaturdasi", "Amavasya"];
    // var masa = ["Chaitra", "Vaishaakha", "Jyaishtha", "Aashaadha", "Shraavana", "Bhaadrapada", "Aashvayuja", "Kaartika", "Maargashiirsha", "Pushya", "Maagha", "Phaalguna"];
    
    var dt = new Date();
    dt.setFullYear(2023);
    dt.setMonth(02);
    dt.setDate(1);
    //console.log("next new moon", getNextNewMoon(dt, true))
    var s1 = null;
    var s2 = null;
    for (var i = 0; i < masa_en.length; i++) {
        if (s1 === null) {
            s1 = getNextNewMoon(dt,false);
        } else {
            s1 = panchang.months[panchang.months.length - 1].end;
        }
        s1.setDate(s1.getDate() + 1);
        s2 = new Date(s1);
        s2.setDate(s2.getDate() + 1);
        s2 = getNextNewMoon(s2, false);
        panchang.months.push({
            name : masa[i],
            start : s1,
            end : s2
        })
        dt = new Date(s2);
        dt.setDate(dt.getDate() + 1);
        console.log(panchang.months[panchang.months.length - 1]);
    }
   }
   
}


//console.log(panchang)
//panchang.init();
panchang.buildAdhikaMasa();