chuki

引数に指定したユリウス日の直前の中気を求めます。

構文
chuki( JD )
引数
JD
ユリウス日
戻値
中気と太陽黄経を格納した配列(0 : 中気, 1 : 太陽黄経)

プログラム

////////////////////////////////////////////////// // 【引数】 // JD : ユリウス日 // 【戻値】 // 中気と太陽黄経を格納した配列(0 : 中気, 1 : 太陽黄経) ////////////////////////////////////////////////// FUNCTION chuki(JD) JD = JD - 9/24 DIM t = (JD + 0.5 - 2451545) / 36525 DIM λsun = longitudeSun(t) DIM λsun0 = 30 * INT(λsun/30) REPEAT t = (JD + 0.5 - 2451545) / 36525 λsun = longitudeSun(t) DIM Δλ = λsun - λsun0 SELECT TRUE CASE Δλ > 180 Δλ = Δλ - 360 CASE Δλ < -180 Δλ = Δλ + 180 SELEND DIM Δt = Δλ * 365.2 / 360 JD = JD - Δt UNTIL Δt <= 1/86400 JD = JD + 9/24 DIM arr[1] = JD, λsun RESULT = SLICE(arr) FEND ////////////////////////////////////////////////// // 【引数】 // deg : 度数法 // 【戻値】 // 度数法から弧度法に変換した値 ////////////////////////////////////////////////// FUNCTION degToRad(deg) DIM pi = 3.14159265358979 RESULT = deg * pi / 180 FEND ////////////////////////////////////////////////// // 【引数】 // JC : ユリウス世紀 // 【戻値】 // 太陽黄経(λsun) ////////////////////////////////////////////////// FUNCTION longitudeSun(JC) DIM A[18] = 0.0004, 0.0004, 0.0005, 0.0005, 0.0006, 0.0007, 0.0007, 0.0007, 0.0013, 0.0015, 0.0018, 0.0018, 0.0020, 0.0200, -0.0048*JC, 1.9147, 36000.7695*JC, 280.4659 DIM k[18] = 31557.0, 29930.0, 2281.0, 155.0, 33718.0, 9038.0, 3035.0, 65929.0, 22519.0, 45038.0, 445267.0, 19.0, 32964.0, 71998.1, 35999.05, 35999.05, 0, 0 DIM θ0[18] = 161.0, 48.0, 221.0, 118.0, 316.0, 64.0, 110.0, 45.0, 352.0, 254.0, 208.0, 159.0, 158.0, 265.1, 267.52, 267.52, 0, 0 DIM λsun[18] FOR n = 0 TO 15 DIM ang = normalizeAngle(k[n] * JC + θ0[n]) λsun[n] = A[n] * COS(degToRad(ang)) NEXT λsun[16] = normalizeAngle(A[16]) λsun[17] = normalizeAngle(A[17]) RESULT = normalizeAngle(CALCARRAY(λsun, CALC_ADD)) FEND ////////////////////////////////////////////////// // 【引数】 // deg : 度数法 // 【戻値】 // 度単位の角度を0~360度に正規化 ////////////////////////////////////////////////// FUNCTION normalizeAngle(deg) SELECT TRUE CASE deg >= 360 deg = deg - INT(deg / 360) * 360 CASE deg < 0 deg = deg + INT(ABS(deg / 360) + 1) * 360 SELEND RESULT = deg FEND

中気

二十四節気のうち正節を除く12の節気を中気ちゅうきと言います。雨水(330)・春分(0)・穀雨(30)・小満(60)・夏至(90)・大暑(120)・処暑(150)・秋分(180)・霜降(210)・小雪(240)・冬至(270)・大寒(300)がそれにあたります。( )内は太陽黄経の角度。

旧暦(天保暦)には、原則として1月から12月までの各月に節気が入る。また調整のため閏月が入ることもある。

プログラム実行例

今年の穀雨の日付を求める

穀雨は太陽黄経30度なので、翌年1月1日から直前の中気を太陽黄経が30度になるまで繰り返す。

GETTIME() DIM year = G_TIME_YY + 1 DIM month = 1 DIM day = 1 REPEAT DIM chuki = chuki(YMDToJD(year, month, day)) DIM d = JDToYMD(chuki[0]) year = d[0] month = d[1] day = d[2] UNTIL ROUND(chuki[1]) = 30 PRINT d[0] + "/" + d[1] + "/" + d[2] ////////////////////////////////////////////////// // 【引数】 // JD : ユリウス日 // 【戻値】 // 中気と太陽黄経を格納した配列(0 : 中気, 1 : 太陽黄経) ////////////////////////////////////////////////// FUNCTION chuki(JD) JD = JD - 9/24 DIM t = (JD + 0.5 - 2451545) / 36525 DIM λsun = longitudeSun(t) DIM λsun0 = 30 * INT(λsun/30) REPEAT t = (JD + 0.5 - 2451545) / 36525 λsun = longitudeSun(t) DIM Δλ = λsun - λsun0 SELECT TRUE CASE Δλ > 180 Δλ = Δλ - 360 CASE Δλ < -180 Δλ = Δλ + 180 SELEND DIM Δt = Δλ * 365.2 / 360 JD = JD - Δt UNTIL Δt <= 1/86400 JD = JD + 9/24 DIM arr[1] = JD, λsun RESULT = SLICE(arr) FEND ////////////////////////////////////////////////// // 【引数】 // deg : 度数法 // 【戻値】 // 度数法から弧度法に変換した値 ////////////////////////////////////////////////// FUNCTION degToRad(deg) DIM pi = 3.14159265358979 RESULT = deg * pi / 180 FEND ////////////////////////////////////////////////// // 【引数】 // JD : ユリウス日 // 【戻値】 // グレゴリオ暦を格納した配列(0 : 年, 1 : 月, 2 : 日, 3 : 時, 4 : 分, 5 : 秒) ////////////////////////////////////////////////// FUNCTION JDToYMD(JD) DIM x0 = INT(JD + 68570) DIM x1 = INT(x0 / 36524.25) DIM x2 = x0 - INT(36524.25 * x1 + 0.75) DIM x3 = INT((x2 + 1) / 365.2425) DIM x4 = x2 - INT(365.25 * x3) + 31 DIM x5 = INT(INT(x4) / 30.59) DIM x6 = INT(INT(x5) / 11) DIM t2 = x4 - INT(30.59 * x5) DIM t1 = x5 - 12 * x6 + 2 DIM t0 = 100 * (x1 - 49) + x3 + x6 IFB t1 = 2 AND t2 > 28 THEN SELECT TRUE CASE t0 MOD 100 = 0 AND t0 MOD 400 = 0 t2 = 29 CASE t0 MOD 4 = 0 t2 = 29 DEFAULT t2 = 28 SELEND ENDIF DIM tm = 86400 * (JD - INT(JD)) DIM t3 = INT(tm / 3600) DIM t4 = INT((tm - 3600 * t3) / 60) DIM t5 = INT(tm - 3600 * t3 - 60 * t4) DIM t[] = t0, t1, t2, t3, t4, t5 RESULT = SLICE(t) FEND ////////////////////////////////////////////////// // 【引数】 // JC : ユリウス世紀 // 【戻値】 // 太陽黄経(λsun) ////////////////////////////////////////////////// FUNCTION longitudeSun(JC) DIM A[18] = 0.0004, 0.0004, 0.0005, 0.0005, 0.0006, 0.0007, 0.0007, 0.0007, 0.0013, 0.0015, 0.0018, 0.0018, 0.0020, 0.0200, -0.0048*JC, 1.9147, 36000.7695*JC, 280.4659 DIM k[18] = 31557.0, 29930.0, 2281.0, 155.0, 33718.0, 9038.0, 3035.0, 65929.0, 22519.0, 45038.0, 445267.0, 19.0, 32964.0, 71998.1, 35999.05, 35999.05, 0, 0 DIM θ0[18] = 161.0, 48.0, 221.0, 118.0, 316.0, 64.0, 110.0, 45.0, 352.0, 254.0, 208.0, 159.0, 158.0, 265.1, 267.52, 267.52, 0, 0 DIM λsun[18] FOR n = 0 TO 15 DIM ang = normalizeAngle(k[n] * JC + θ0[n]) λsun[n] = A[n] * COS(degToRad(ang)) NEXT λsun[16] = normalizeAngle(A[16]) λsun[17] = normalizeAngle(A[17]) RESULT = normalizeAngle(CALCARRAY(λsun, CALC_ADD)) FEND ////////////////////////////////////////////////// // 【引数】 // deg : 度数法 // 【戻値】 // 度単位の角度を0~360度に正規化 ////////////////////////////////////////////////// FUNCTION normalizeAngle(deg) SELECT TRUE CASE deg >= 360 deg = deg - INT(deg / 360) * 360 CASE deg < 0 deg = deg + INT(ABS(deg / 360) + 1) * 360 SELEND RESULT = deg FEND ////////////////////////////////////////////////// // 【引数】 // year : 年 // month : 月 // day : 日 // hour : 時 // minute : 分 // second : 秒 // 【戻値】 // ユリウス日 ////////////////////////////////////////////////// FUNCTION YMDToJD(year, month, day, hour = 0, minute = 0, second = 0) IFB month < 3 THEN year = year - 1 month = month + 12 ENDIF DIM JD = INT(year * 365.25) JD = JD + INT(year / 400) JD = JD - INT(year / 100) JD = JD + INT((month - 2) * 30.59) JD = JD + 1721088 JD = JD + day DIM t = second / 3600 t = t + minute / 60 t = t + hour t = t / 24 JD = JD + t RESULT = JD FEND
結果
2020/4/19

今年の春彼岸の日付を求める

翌年1月1日の直前の太陽黄経0度の日付(春分の日)を求め、その日を中日とする前後3を合わせた7日間を求める。

GETTIME() DIM year = G_TIME_YY + 1 DIM month = 1 DIM day = 1 REPEAT DIM chuki = chuki(YMDToJD(year, month, day)) DIM d = JDToYMD(chuki[0]) year = d[0] month = d[1] day = d[2] UNTIL normalizeAngle(ROUND(chuki[1])) = 0 RESIZE(d, 2) arrayMap("TEXT(%val%, <#DBL>00<#DBL>)", d) date = JOIN(d, "/") PRINT "彼岸入り<#TAB>" + dateAdd("d", -3, date) PRINT "春分の日<#TAB>" + date PRINT "彼岸明け<#TAB>" + dateAdd("d", 3, date) ////////////////////////////////////////////////// // 【引数】 // callback : 配列の各要素に適用するコールバック関数 // array : コールバック関数を適用する配列 // 【戻値】 // callback関数を適用した後、 適用後の要素を含む配列 ////////////////////////////////////////////////// PROCEDURE arrayMap(callback, var array[]) DIM tmp[RESIZE(array)] DIM n = 0 FOR %val% IN array tmp[n] = EVAL(callback) n = n + 1 NEXT RESIZE(array, RESIZE(tmp)) FOR n = 0 TO RESIZE(tmp) array[n] = tmp[n] NEXT FEND ////////////////////////////////////////////////// // 【引数】 // JD : ユリウス日 // 【戻値】 // 中気と太陽黄経を格納した配列(0 : 中気, 1 : 太陽黄経) ////////////////////////////////////////////////// FUNCTION chuki(JD) JD = JD - 9/24 DIM t = (JD + 0.5 - 2451545) / 36525 DIM λsun = longitudeSun(t) DIM λsun0 = 30 * INT(λsun/30) REPEAT t = (JD + 0.5 - 2451545) / 36525 λsun = longitudeSun(t) DIM Δλ = λsun - λsun0 SELECT TRUE CASE Δλ > 180 Δλ = Δλ - 360 CASE Δλ < -180 Δλ = Δλ + 180 SELEND DIM Δt = Δλ * 365.2 / 360 JD = JD - Δt UNTIL Δt <= 1/86400 JD = JD + 9/24 DIM arr[1] = JD, λsun RESULT = SLICE(arr) FEND ////////////////////////////////////////////////// // 【引数】 // interval : 加算する時間間隔を表す文字列式(yyyy:年、m:月、d:日、ww:週、h:時、n:分、s:秒) // num : dateに加算する値。未来は正、過去は負で指定 // date : 時間間隔を加算する日付 // 【戻値】 // 日時(date)に、指定した単位(interval)の時間(num)を加算して返します ////////////////////////////////////////////////// FUNCTION dateAdd(interval, num, date) GETTIME(0, date) SELECT interval CASE "yyyy" d = (G_TIME_YY + num) + "/" + G_TIME_MM2 + "/" + G_TIME_DD2 CASE "m" IFB num > 0 THEN year = G_TIME_YY + INT((G_TIME_MM + num) / 12) month = REPLACE(FORMAT(((G_TIME_MM + num) MOD 12), 2), " ", "0") ELSE year = G_TIME_YY + CEIL((G_TIME_MM + num) / 12 - 1) month = REPLACE(FORMAT(G_TIME_MM - (ABS(num) MOD 12), 2), " ", "0") ENDIF IF month = "00" THEN month = 12 day = G_TIME_DD2 d = "" + year + month + day IFB !isDate(d) THEN d = year + "/" + month + "/" + "01" d = getEndOfMonth(d) ELSE d = year + "/" + month + "/" + day ENDIF CASE "d" GETTIME(num, date) d = G_TIME_YY4 + "/" + G_TIME_MM2 + "/" + G_TIME_DD2 CASE "ww" GETTIME(num * 7, date) d = G_TIME_YY4 + "/" + G_TIME_MM2 + "/" + G_TIME_DD2 CASE "h" GETTIME(num / 24, date) d = G_TIME_YY4 + "/" + G_TIME_MM2 + "/" + G_TIME_DD2 + " " + G_TIME_HH2 + ":" + G_TIME_NN2 + ":" + G_TIME_SS2 CASE "n" GETTIME(num / 1440, date) d = G_TIME_YY4 + "/" + G_TIME_MM2 + "/" + G_TIME_DD2 + " " + G_TIME_HH2 + ":" + G_TIME_NN2 + ":" + G_TIME_SS2 CASE "s" GETTIME(num / 86400, date) d = G_TIME_YY4 + "/" + G_TIME_MM2 + "/" + G_TIME_DD2 + " " + G_TIME_HH2 + ":" + G_TIME_NN2 + ":" + G_TIME_SS2 SELEND result = d FEND ////////////////////////////////////////////////// // 【引数】 // interval : 時間単位(yyyy︰年、m︰月、d︰日) // date1 : 日時1 // date2 : 日時2 // 【戻値】 // date2からdate1を引いた日数を求めます。 ////////////////////////////////////////////////// FUNCTION dateDiff(interval, date1, date2) SELECT interval CASE "yyyy" GETTIME(0, date1) y1 = G_TIME_YY GETTIME(0, date2) y2 = G_TIME_YY d = y2 - y1 CASE "m" GETTIME(0, date1) y1 = G_TIME_YY m1 = G_TIME_MM GETTIME(0, date2) y2 = G_TIME_YY m2 = G_TIME_MM d = (y2 - y1) * 12 + m2 - m1 CASE "d" d1 = GETTIME(0, date1) d2 = GETTIME(0, date2) d = (d2 - d1) / 86400 SELEND RESULT = d FEND ////////////////////////////////////////////////// // 【引数】 // deg : 度数法 // 【戻値】 // 度数法から弧度法に変換した値 ////////////////////////////////////////////////// FUNCTION degToRad(deg) DIM pi = 3.14159265358979 RESULT = deg * pi / 180 FEND ////////////////////////////////////////////////// // 【引数】 // arr : 最大公約数を求める数値を格納した配列 // 【戻値】 // 最大公約数 ////////////////////////////////////////////////// FUNCTION GCD(arr[]) c = LENGTH(arr) rem = arr[c-1] MOD arr[c-2] IFB rem = 0 THEN IFB LENGTH(arr) = 2 THEN RESULT = arr[c-2] EXIT ENDIF RESIZE(arr, c-2) RESULT = GCD(arr) EXIT ENDIF arr[c-1] = arr[c-2] arr[c-2] = rem RESULT = GCD(arr) FEND ////////////////////////////////////////////////// // 【引数】 // date : 日付(”YYYYMMDD” or “YYYY/MM/DD” or “YYYY-MM-DD” or “YYYYMMDDHHNNSS” or “YYYY/MM/DD HH:NN:SS”) // m : 第一引数の指定日からプラスマイナス m 月とする(デフォルト=0) // 【戻値】 // date から m 月後の月末の日付 ////////////////////////////////////////////////// FUNCTION getEndOfMonth(date, m = 0) GETTIME(0, date) IFB m >= 0 THEN year = G_TIME_YY + INT((G_TIME_MM + m + 1) / 12) month = REPLACE(FORMAT(((G_TIME_MM + m + 1) MOD 12), 2), " ", "0") ELSE year = G_TIME_YY + CEIL((G_TIME_MM + m + 1) / 12 - 1) month = REPLACE(FORMAT(G_TIME_MM - (ABS(m + 1) MOD 12), 2), " ", "0") IF EVAL(month) < 0 THEN month = REPLACE(FORMAT(12 + EVAL(month), 2), " ", "0") ENDIF IF month = "00" THEN month = "12" day = "01" d = "" + year + month + day GETTIME(-1, d) RESULT = G_TIME_YY4 + "/" + G_TIME_MM2 + "/" + G_TIME_DD2 FEND ////////////////////////////////////////////////// // 【引数】 // 数値 : GETTIMEで取得した曜日の番号(G_TIME_WW=0:日曜….6:土曜) // フォーマット : 以下のアルファベットは数値で指定した対応する曜日に置き換えられます。(aaa : 日〜土、aaaa : 日曜日〜土曜日、ddd : Sun〜Sat、dddd : Sunday〜Saturday) // 【戻値】 // 数値をフォーマットした文字列を返します。 ////////////////////////////////////////////////// FUNCTION getWeekdayName(num, format = "aaa") DIM re = CREATEOLEOBJ("VBScript.RegExp") DIM aaa[] = "日", "月", "火", "水", "木", "金", "土"; DIM aaaa[] = "日曜日", "月曜日", "火曜日", "水曜日", "木曜日", "金曜日", "土曜日"; DIM ddd[] = "Sun", "Mon", "Tue", "Wed", "Thu", "Fri", "Sat"; DIM dddd[] = "Sunday", "Monday", "Tuesday", "Wednesday", "Thursday", "Friday", "Saturday"; match = reExecute(format, "(a|d)+") type = match.Item(0).Value RESULT = reReplace(format, EVAL(type + "[" + num + "]"), type) FEND ////////////////////////////////////////////////// // 【引数】 // シリアル値 : 時間を表すシリアル値を指定 // 【戻値】 // 例)0…00、0.5…12、1.0…24 ////////////////////////////////////////////////// FUNCTION Hour(serial) RESULT = REPLACE(FORMAT(INT(serial * 24), 2), " ", "0") FEND ////////////////////////////////////////////////// // 【引数】 // 日付 : YYYYMMDD // 【戻値】 // TRUE : 日付として認識できる、FALSE : 日付として認識できない ////////////////////////////////////////////////// FUNCTION isDate(date) Pattern = "^(?!([02468][1235679]|[13579][01345789])000229)(([0-9]{4}(01|03|05|07|08|10|12)(0[1-9]|[12][0-9]|3[01]))|([0-9]{4}(04|06|09|11)(0[1-9]|[12][0-9]|30))|([0-9]{4}02(0[1-9]|1[0-9]|2[0-8]))|([0-9]{2}([02468][048]|[13579][26])0229))$" match = reExecute(date, Pattern) SELECT match.Count CASE 1 res = TRUE CASE 0 res = FALSE SELEND RESULT = res FEND ////////////////////////////////////////////////// // 【引数】 // JD : ユリウス日 // 【戻値】 // グレゴリオ暦を格納した配列(0 : 年, 1 : 月, 2 : 日, 3 : 時, 4 : 分, 5 : 秒) ////////////////////////////////////////////////// FUNCTION JDToYMD(JD) DIM x0 = INT(JD + 68570) DIM x1 = INT(x0 / 36524.25) DIM x2 = x0 - INT(36524.25 * x1 + 0.75) DIM x3 = INT((x2 + 1) / 365.2425) DIM x4 = x2 - INT(365.25 * x3) + 31 DIM x5 = INT(INT(x4) / 30.59) DIM x6 = INT(INT(x5) / 11) DIM t2 = x4 - INT(30.59 * x5) DIM t1 = x5 - 12 * x6 + 2 DIM t0 = 100 * (x1 - 49) + x3 + x6 IFB t1 = 2 AND t2 > 28 THEN SELECT TRUE CASE t0 MOD 100 = 0 AND t0 MOD 400 = 0 t2 = 29 CASE t0 MOD 4 = 0 t2 = 29 DEFAULT t2 = 28 SELEND ENDIF DIM tm = 86400 * (JD - INT(JD)) DIM t3 = INT(tm / 3600) DIM t4 = INT((tm - 3600 * t3) / 60) DIM t5 = INT(tm - 3600 * t3 - 60 * t4) DIM t[] = t0, t1, t2, t3, t4, t5 RESULT = SLICE(t) FEND ////////////////////////////////////////////////// // 【引数】 // JC : ユリウス世紀 // 【戻値】 // 太陽黄経(λsun) ////////////////////////////////////////////////// FUNCTION longitudeSun(JC) DIM A[18] = 0.0004, 0.0004, 0.0005, 0.0005, 0.0006, 0.0007, 0.0007, 0.0007, 0.0013, 0.0015, 0.0018, 0.0018, 0.0020, 0.0200, -0.0048*JC, 1.9147, 36000.7695*JC, 280.4659 DIM k[18] = 31557.0, 29930.0, 2281.0, 155.0, 33718.0, 9038.0, 3035.0, 65929.0, 22519.0, 45038.0, 445267.0, 19.0, 32964.0, 71998.1, 35999.05, 35999.05, 0, 0 DIM θ0[18] = 161.0, 48.0, 221.0, 118.0, 316.0, 64.0, 110.0, 45.0, 352.0, 254.0, 208.0, 159.0, 158.0, 265.1, 267.52, 267.52, 0, 0 DIM λsun[18] FOR n = 0 TO 15 DIM ang = normalizeAngle(k[n] * JC + θ0[n]) λsun[n] = A[n] * COS(degToRad(ang)) NEXT λsun[16] = normalizeAngle(A[16]) λsun[17] = normalizeAngle(A[17]) RESULT = normalizeAngle(CALCARRAY(λsun, CALC_ADD)) FEND ////////////////////////////////////////////////// // 【引数】 // deg : 度数法 // 【戻値】 // 度単位の角度を0~360度に正規化 ////////////////////////////////////////////////// FUNCTION normalizeAngle(deg) SELECT TRUE CASE deg >= 360 deg = deg - INT(deg / 360) * 360 CASE deg < 0 deg = deg + INT(ABS(deg / 360) + 1) * 360 SELEND RESULT = deg FEND ////////////////////////////////////////////////// // 【引数】 // str : 正規表現による検索の対象となる文字列 // Pattern : 正規表現で使用するパターンを設定 // IgnoreCase : 大文字・小文字を区別しない場合はTrue、区別する場合はFalse // Global : 文字列全体を検索する場合はTrue、しない場合はFalse // 【戻値】 // 正規表現で検索した結果をMatchesコレクションとして返します。 ////////////////////////////////////////////////// FUNCTION reExecute(str, Pattern, IgnoreCase = TRUE, Global = TRUE) DIM re = CREATEOLEOBJ("VBScript.RegExp") re.Pattern = Pattern re.IgnoreCase = IgnoreCase re.Global = Global RESULT = re.Execute(str) FEND ////////////////////////////////////////////////// // 【引数】 // str1 : 置換される文字列 // str2 : 置換後の文字列 // Pattern : 置換する文字列のパターン // IgnoreCase : 大文字・小文字を区別しない場合はTrue、区別する場合はFalse // Global : 文字列全体を検索する場合はTrue、しない場合はFalse // 【戻値】 // 正規表現置換後の文字列 ////////////////////////////////////////////////// FUNCTION reReplace(str1, str2, Pattern, IgnoreCase = TRUE, Global = TRUE) DIM re = CREATEOLEOBJ("VBScript.RegExp") re.Pattern = Pattern re.IgnoreCase = IgnoreCase re.Global = Global RESULT = re.Replace(str1, str2) FEND ////////////////////////////////////////////////// // 【引数】 // str : 正規表現による検索の対象となる文字列 // Pattern : 正規表現で使用するパターンを設定 // IgnoreCase : 大文字・小文字を区別しない場合はTrue、区別する場合はFalse // Global : 文字列全体を検索する場合はTrue、しない場合はFalse // 【戻値】 // 正規表現にマッチするかどうか ////////////////////////////////////////////////// FUNCTION reTest(str, Pattern, IgnoreCase = TRUE, Global = TRUE) DIM re = CREATEOLEOBJ("VBScript.RegExp") re.Pattern = Pattern re.IgnoreCase = IgnoreCase re.Global = Global RESULT = re.Test(str) FEND ////////////////////////////////////////////////// // 【引数】 // inputs : 繰り返す文字列 // multiplier : inputsを繰り返す回数 // 【戻値】 // inputsをmultiplier回を繰り返した文字列を返します ////////////////////////////////////////////////// FUNCTION strRepeat(inputs, multiplier) DIM res = "" FOR n = 1 TO multiplier res = res + inputs NEXT RESULT = res FEND ////////////////////////////////////////////////// // 【引数】 // serial : シリアル値 // format : フォーマット // 【戻値】 // 数値を表示書式に基づいて変換した文字列 ////////////////////////////////////////////////// FUNCTION TEXT(serial, format) SELECT TRUE CASE reTest(format, "\[h+\]") Matches = reExecute(format, "\[(h+)\]") RESULT = TEXT(Hour(serial), strRepeat("0", LENGTH(Matches.Item(0).SubMatches(0)))) CASE reTest(format, "^h+$") Matches = reExecute(format, "^(h+)$") RESULT = TEXT(Hour(serial) MOD 24, strRepeat("0", LENGTH(Matches.Item(0).SubMatches(0)))) CASE reTest(format, "\[m+\]") Matches = reExecute(format, "\[(m+)\]") RESULT = TEXT(serial * 1440, strRepeat("0", LENGTH(Matches.Item(0).SubMatches(0)))) CASE format = "m" GETTIME(serial, "1899/12/30") RESULT = TEXT(G_TIME_MM, "0") CASE format = "mm" GETTIME(serial, "1899/12/30") RESULT = G_TIME_MM2 CASE format = "s" GETTIME(serial, "1899/12/30") RESULT = TEXT(G_TIME_SS, "0") CASE format = "ss" GETTIME(serial, "1899/12/30") RESULT = G_TIME_SS2 CASE format = "yyyy" GETTIME(serial, "1899/12/30") RESULT = G_TIME_YY4 CASE format = "yy" GETTIME(serial, "1899/12/30") RESULT = COPY(G_TIME_YY4, 3, 2) CASE format = "e" SELECT TRUE // 令和 CASE dateDiff("d", "2019/05/01", TEXT(serial, "yyyy/mm/dd")) >= 0 RESULT = TEXT(serial, "yyyy") - 2018 // 平成 CASE dateDiff("d", "1989/01/08", TEXT(serial, "yyyy/mm/dd")) >= 0 RESULT = TEXT(serial, "yyyy") - 1988 // 昭和 CASE dateDiff("d", "1926/12/25", TEXT(serial, "yyyy/mm/dd")) >= 0 RESULT = TEXT(serial, "yyyy") - 1925 // 大正 CASE dateDiff("d", "1912/07/30", TEXT(serial, "yyyy/mm/dd")) >= 0 RESULT = TEXT(serial, "yyyy") - 1911 // 明治 CASE dateDiff("d", "1868/01/25", TEXT(serial, "yyyy/mm/dd")) >= 0 RESULT = TEXT(serial, "yyyy") - 1867 SELEND CASE format = "g" SELECT TRUE // 令和 CASE dateDiff("d", "2019/05/01", TEXT(serial, "yyyy/mm/dd")) >= 0 RESULT = "R" // 平成 CASE dateDiff("d", "1989/01/08", TEXT(serial, "yyyy/mm/dd")) >= 0 RESULT = "H" // 昭和 CASE dateDiff("d", "1926/12/25", TEXT(serial, "yyyy/mm/dd")) >= 0 RESULT = "S" // 大正 CASE dateDiff("d", "1912/07/30", TEXT(serial, "yyyy/mm/dd")) >= 0 RESULT = "T" // 明治 CASE dateDiff("d", "1868/01/25", TEXT(serial, "yyyy/mm/dd")) >= 0 RESULT = "M" SELEND CASE format = "gg" RESULT = COPY(TEXT(serial, "ggg"), 1, 1) CASE format = "ggg" SELECT TRUE // 令和 CASE dateDiff("d", "2019/05/01", TEXT(serial, "yyyy/mm/dd")) >= 0 RESULT = "令和" // 平成 CASE dateDiff("d", "1989/01/08", TEXT(serial, "yyyy/mm/dd")) >= 0 RESULT = "平成" // 昭和 CASE dateDiff("d", "1926/12/25", TEXT(serial, "yyyy/mm/dd")) >= 0 RESULT = "昭和" // 大正 CASE dateDiff("d", "1912/07/30", TEXT(serial, "yyyy/mm/dd")) >= 0 RESULT = "大正" // 明治 CASE dateDiff("d", "1868/01/25", TEXT(serial, "yyyy/mm/dd")) >= 0 RESULT = "明治" SELEND CASE format = "mmmmm" RESULT = COPY(TEXT(serial, "mmmm"), 1, 1) CASE format = "mmmm" SELECT TEXT(serial, "m") CASE 1; RESULT = "January" CASE 2; RESULT = "February" CASE 3; RESULT = "March" CASE 4; RESULT = "April" CASE 5; RESULT = "May" CASE 6; RESULT = "June" CASE 7; RESULT = "July" CASE 8; RESULT = "August" CASE 9; RESULT = "September" CASE 10; RESULT = "October" CASE 11; RESULT = "November" CASE 12; RESULT = "December" SELEND CASE format = "mmm" RESULT = COPY(TEXT(serial, "mmmm"), 1, 3) CASE format = "dd" GETTIME(serial, "1899/12/30") RESULT = TEXT(G_TIME_DD2, "00") CASE format = "d" GETTIME(serial, "1899/12/30") RESULT = TEXT(G_TIME_DD, "0") CASE format = "aaaa" GETTIME(serial, "1899/12/30") RESULT = getWeekdayName(G_TIME_WW, "aaaa") CASE format = "aaa" GETTIME(serial, "1899/12/30") RESULT = getWeekdayName(G_TIME_WW, "aaa") CASE format = "dddd" GETTIME(serial, "1899/12/30") RESULT = getWeekdayName(G_TIME_WW, "dddd") CASE format = "ddd" GETTIME(serial, "1899/12/30") RESULT = getWeekdayName(G_TIME_WW, "ddd") CASE reTest(format, "(0+\.?0+)?%") Matches = reExecute(format, "(0+\.?0+)?%") RESULT = TEXT(serial * 100, Matches.Item(0).SubMatches(0)) + "%" CASE reTest(format, "^(.*?)\/(.*?)$") // 分母、分子 DIM demon, numer // 分数 DIM frac[1] DIM Pattern = "(\d+\.\d*?(\d+)\2+)" IFB reTest(serial, Pattern) THEN Matches = reExecute(serial, Pattern) // 循環節の桁数 digit = LENGTH(Matches.Item(0).SubMatches(1)) denom = POWER(10, digit) - 1 numer = serial * POWER(10, digit) - serial // 分子が小数なら整数に通分する num = POWER(10, LENGTH(numer) - POS(".", numer)) denom = denom * num numer = numer * num ELSE denom = POWER(10, LENGTH(serial) - POS(".", serial)) numer = serial * denom ENDIF frac[0] = denom frac[1] = numer // 最大公約数を求める n = GCD(frac) denom = frac[0] numer = frac[1] // 既約分数にする denom = denom / n numer = numer / n RESULT = numer + "/" + denom CASE reTest(format, "((#+),?)+") RESULT = reReplace(serial, "$1,", "(\d)(?=(\d{3})+(?!\d))") CASE reTest(format, "(0+)\.?(0+)?") Matches = reExecute(format, "(0+)\.?(0+)?") len1 = LENGTH(Matches.Item(0).SubMatches(0)) len2 = LENGTH(Matches.Item(0).SubMatches(1)) DIM arr[] = LENGTH(INT(serial)), len1 IFB POS(".", format) THEN RESULT = REPLACE(FORMAT(serial, CALCARRAY(arr, CALC_MAX) + len2 + 1, len2), " ", "0") ELSE RESULT = REPLACE(FORMAT(serial, CALCARRAY(arr, CALC_MAX)), " ", "0") ENDIF CASE reTest(format, "(.+)?((?:AM?|am?)\/(?:PM?|pm?))(.+)?") DIM arr[2] Matches = reExecute(format, "(.+)?((?:AM?|am?)\/(?:PM?|pm?))(.+)?") res = "" FOR n = 0 TO Matches.Item(0).SubMatches.Count - 1 arr[n] = Matches.Item(0).SubMatches(n) NEXT DIM res[2] FOR n = 0 TO 2 STEP 2 Matches = reExecute(arr[n], "([^yegmdahs]+)?(([yegmdahs])\3{0,})([^yegmdahs]+)?") FOR m = 0 TO Matches.Count - 1 SELECT TRUE // 「m」を含む CASE reTest(Matches.Item(m).SubMatches(1), "m+") // アルファベット部分の前後で「m」が月か分か判断 SELECT TRUE // 「m」が最初に出現するアルファベットなら CASE m = 0 // 次に出現するアルファベットが「d」なら「m」は月 IFB reTest(Matches.Item(m+1).SubMatches(1), "d+") THEN res[n] = res[n] + TEXT(serial, Matches.Item(m).SubMatches(1)) ENDIF // 次に出現するアルファベットが「s」なら「m」は分 IFB reTest(Matches.Item(m+1).SubMatches(1), "s+") THEN res[n] = res[n] + TEXT(INT(serial * 1440) MOD 60, strRepeat("0", LENGTH(Matches.Item(m).SubMatches(1)))) ENDIF // 「m」が2番目(2-1=1)以降に出現するアルファベットなら CASE m >= 1 // 前に出現したアルファベットが「y」なら「m」は月 IFB reTest(Matches.Item(m-1).SubMatches(1), "y+") THEN res[n] = res[n] + TEXT(serial, Matches.Item(m).SubMatches(1)) ENDIF // 前に出現したアルファベットが「h」なら「m」は分 IFB reTest(Matches.Item(m-1).SubMatches(1), "h+") THEN res[n] = res[n] + TEXT(INT(serial * 1440) MOD 60, strRepeat("0", LENGTH(Matches.Item(m).SubMatches(1)))) ENDIF SELEND CASE reTest(Matches.Item(m).SubMatches(1), "^h+$") IFB TEXT(serial, Matches.Item(m).SubMatches(1)) < 12 THEN res[1] = COPY(arr[1], 1, POS("/", arr[1]) - 1) // 午前 ELSE res[1] = COPY(arr[1], POS("/", arr[1]) + 1) // 午後 ENDIF res[n] = res[n] + TEXT(serial, Matches.Item(m).SubMatches(1)) // それ以外ならば DEFAULT res[n] = res[n] + TEXT(serial, Matches.Item(m).SubMatches(1)) SELEND res[n] = Matches.Item(m).SubMatches(0) + res[n] + Matches.Item(m).SubMatches(3) NEXT NEXT RESULT = JOIN(res, "") CASE reTest(format, "(([yegmdahs])\2{0,})([^yegmdahs]+)?") res = "" Matches = reExecute(format, "(([yegmdahs])\2{0,})([^yegmdahs]+)?") FOR n = 0 TO Matches.Count - 1 // mの前がyもしくはmの後がdならば→mは月 SELECT TRUE // 「m」を含む CASE reTest(Matches.Item(n).SubMatches(0), "m+") // アルファベット部分の前後で「m」が月か分か判断 SELECT TRUE // 「m」が最初に出現するアルファベットなら CASE n = 0 // 次に出現するアルファベットが「d」なら「m」は月 IFB reTest(Matches.Item(n+1).SubMatches(0), "d+") THEN res = res + TEXT(serial, Matches.Item(n).SubMatches(0)) ENDIF // 次に出現するアルファベットが「s」なら「m」は分 IFB reTest(Matches.Item(n+1).SubMatches(0), "s+") THEN res = res + TEXT(INT(serial * 1440) MOD 60, strRepeat("0", LENGTH(Matches.Item(n).SubMatches(0)))) ENDIF // 「m」が2番目(2-1=1)以降に出現するアルファベットなら CASE n >= 1 // 前に出現したアルファベットが「y」なら「m」は月 IFB reTest(Matches.Item(n-1).SubMatches(0), "y+") THEN res = res + TEXT(serial, Matches.Item(n).SubMatches(0)) ENDIF // 前に出現したアルファベットが「h」なら「m」は分 IFB reTest(Matches.Item(n-1).SubMatches(0), "h+") THEN res = res + TEXT(INT(serial * 1440) MOD 60, strRepeat("0", LENGTH(Matches.Item(n).SubMatches(0)))) ENDIF SELEND // それ以外ならば DEFAULT res = res + TEXT(serial, Matches.Item(n).SubMatches(0)) SELEND res = res + Matches.Item(n).SubMatches(2) NEXT RESULT = res SELEND FEND ////////////////////////////////////////////////// // 【引数】 // year : 年 // month : 月 // day : 日 // hour : 時 // minute : 分 // second : 秒 // 【戻値】 // ユリウス日 ////////////////////////////////////////////////// FUNCTION YMDToJD(year, month, day, hour = 0, minute = 0, second = 0) IFB month < 3 THEN year = year - 1 month = month + 12 ENDIF DIM JD = INT(year * 365.25) JD = JD + INT(year / 400) JD = JD - INT(year / 100) JD = JD + INT((month - 2) * 30.59) JD = JD + 1721088 JD = JD + day DIM t = second / 3600 t = t + minute / 60 t = t + hour t = t / 24 JD = JD + t RESULT = JD FEND
結果
彼岸入り	2020/03/17
春分の日	2020/03/20
彼岸明け	2020/03/23

今年の秋彼岸の日付を求める

翌年1月1日の直前の太陽黄経180度の日付(秋分の日)を求め、その日を中日とする前後3日を合わせた7日間を求める。

GETTIME() DIM year = G_TIME_YY + 1 DIM month = 1 DIM day = 1 REPEAT DIM chuki = chuki(YMDToJD(year, month, day)) DIM d = JDToYMD(chuki[0]) year = d[0] month = d[1] day = d[2] UNTIL normalizeAngle(ROUND(chuki[1])) = 180 RESIZE(d, 2) arrayMap("TEXT(%val%, <#DBL>00<#DBL>)", d) date = JOIN(d, "/") PRINT "彼岸入り<#TAB>" + dateAdd("d", -3, date) PRINT "秋分の日<#TAB>" + date PRINT "彼岸明け<#TAB>" + dateAdd("d", 3, date) ////////////////////////////////////////////////// // 【引数】 // callback : 配列の各要素に適用するコールバック関数 // array : コールバック関数を適用する配列 // 【戻値】 // callback関数を適用した後、 適用後の要素を含む配列 ////////////////////////////////////////////////// PROCEDURE arrayMap(callback, var array[]) DIM tmp[RESIZE(array)] DIM n = 0 FOR %val% IN array tmp[n] = EVAL(callback) n = n + 1 NEXT RESIZE(array, RESIZE(tmp)) FOR n = 0 TO RESIZE(tmp) array[n] = tmp[n] NEXT FEND ////////////////////////////////////////////////// // 【引数】 // JD : ユリウス日 // 【戻値】 // 中気と太陽黄経を格納した配列(0 : 中気, 1 : 太陽黄経) ////////////////////////////////////////////////// FUNCTION chuki(JD) JD = JD - 9/24 DIM t = (JD + 0.5 - 2451545) / 36525 DIM λsun = longitudeSun(t) DIM λsun0 = 30 * INT(λsun/30) REPEAT t = (JD + 0.5 - 2451545) / 36525 λsun = longitudeSun(t) DIM Δλ = λsun - λsun0 SELECT TRUE CASE Δλ > 180 Δλ = Δλ - 360 CASE Δλ < -180 Δλ = Δλ + 180 SELEND DIM Δt = Δλ * 365.2 / 360 JD = JD - Δt UNTIL Δt <= 1/86400 JD = JD + 9/24 DIM arr[1] = JD, λsun RESULT = SLICE(arr) FEND ////////////////////////////////////////////////// // 【引数】 // interval : 加算する時間間隔を表す文字列式(yyyy:年、m:月、d:日、ww:週、h:時、n:分、s:秒) // num : dateに加算する値。未来は正、過去は負で指定 // date : 時間間隔を加算する日付 // 【戻値】 // 日時(date)に、指定した単位(interval)の時間(num)を加算して返します ////////////////////////////////////////////////// FUNCTION dateAdd(interval, num, date) GETTIME(0, date) SELECT interval CASE "yyyy" d = (G_TIME_YY + num) + "/" + G_TIME_MM2 + "/" + G_TIME_DD2 CASE "m" IFB num > 0 THEN year = G_TIME_YY + INT((G_TIME_MM + num) / 12) month = REPLACE(FORMAT(((G_TIME_MM + num) MOD 12), 2), " ", "0") ELSE year = G_TIME_YY + CEIL((G_TIME_MM + num) / 12 - 1) month = REPLACE(FORMAT(G_TIME_MM - (ABS(num) MOD 12), 2), " ", "0") ENDIF IF month = "00" THEN month = 12 day = G_TIME_DD2 d = "" + year + month + day IFB !isDate(d) THEN d = year + "/" + month + "/" + "01" d = getEndOfMonth(d) ELSE d = year + "/" + month + "/" + day ENDIF CASE "d" GETTIME(num, date) d = G_TIME_YY4 + "/" + G_TIME_MM2 + "/" + G_TIME_DD2 CASE "ww" GETTIME(num * 7, date) d = G_TIME_YY4 + "/" + G_TIME_MM2 + "/" + G_TIME_DD2 CASE "h" GETTIME(num / 24, date) d = G_TIME_YY4 + "/" + G_TIME_MM2 + "/" + G_TIME_DD2 + " " + G_TIME_HH2 + ":" + G_TIME_NN2 + ":" + G_TIME_SS2 CASE "n" GETTIME(num / 1440, date) d = G_TIME_YY4 + "/" + G_TIME_MM2 + "/" + G_TIME_DD2 + " " + G_TIME_HH2 + ":" + G_TIME_NN2 + ":" + G_TIME_SS2 CASE "s" GETTIME(num / 86400, date) d = G_TIME_YY4 + "/" + G_TIME_MM2 + "/" + G_TIME_DD2 + " " + G_TIME_HH2 + ":" + G_TIME_NN2 + ":" + G_TIME_SS2 SELEND result = d FEND ////////////////////////////////////////////////// // 【引数】 // interval : 時間単位(yyyy︰年、m︰月、d︰日) // date1 : 日時1 // date2 : 日時2 // 【戻値】 // date2からdate1を引いた日数を求めます。 ////////////////////////////////////////////////// FUNCTION dateDiff(interval, date1, date2) SELECT interval CASE "yyyy" GETTIME(0, date1) y1 = G_TIME_YY GETTIME(0, date2) y2 = G_TIME_YY d = y2 - y1 CASE "m" GETTIME(0, date1) y1 = G_TIME_YY m1 = G_TIME_MM GETTIME(0, date2) y2 = G_TIME_YY m2 = G_TIME_MM d = (y2 - y1) * 12 + m2 - m1 CASE "d" d1 = GETTIME(0, date1) d2 = GETTIME(0, date2) d = (d2 - d1) / 86400 SELEND RESULT = d FEND ////////////////////////////////////////////////// // 【引数】 // deg : 度数法 // 【戻値】 // 度数法から弧度法に変換した値 ////////////////////////////////////////////////// FUNCTION degToRad(deg) DIM pi = 3.14159265358979 RESULT = deg * pi / 180 FEND ////////////////////////////////////////////////// // 【引数】 // arr : 最大公約数を求める数値を格納した配列 // 【戻値】 // 最大公約数 ////////////////////////////////////////////////// FUNCTION GCD(arr[]) c = LENGTH(arr) rem = arr[c-1] MOD arr[c-2] IFB rem = 0 THEN IFB LENGTH(arr) = 2 THEN RESULT = arr[c-2] EXIT ENDIF RESIZE(arr, c-2) RESULT = GCD(arr) EXIT ENDIF arr[c-1] = arr[c-2] arr[c-2] = rem RESULT = GCD(arr) FEND ////////////////////////////////////////////////// // 【引数】 // date : 日付(”YYYYMMDD” or “YYYY/MM/DD” or “YYYY-MM-DD” or “YYYYMMDDHHNNSS” or “YYYY/MM/DD HH:NN:SS”) // m : 第一引数の指定日からプラスマイナス m 月とする(デフォルト=0) // 【戻値】 // date から m 月後の月末の日付 ////////////////////////////////////////////////// FUNCTION getEndOfMonth(date, m = 0) GETTIME(0, date) IFB m >= 0 THEN year = G_TIME_YY + INT((G_TIME_MM + m + 1) / 12) month = REPLACE(FORMAT(((G_TIME_MM + m + 1) MOD 12), 2), " ", "0") ELSE year = G_TIME_YY + CEIL((G_TIME_MM + m + 1) / 12 - 1) month = REPLACE(FORMAT(G_TIME_MM - (ABS(m + 1) MOD 12), 2), " ", "0") IF EVAL(month) < 0 THEN month = REPLACE(FORMAT(12 + EVAL(month), 2), " ", "0") ENDIF IF month = "00" THEN month = "12" day = "01" d = "" + year + month + day GETTIME(-1, d) RESULT = G_TIME_YY4 + "/" + G_TIME_MM2 + "/" + G_TIME_DD2 FEND ////////////////////////////////////////////////// // 【引数】 // 数値 : GETTIMEで取得した曜日の番号(G_TIME_WW=0:日曜….6:土曜) // フォーマット : 以下のアルファベットは数値で指定した対応する曜日に置き換えられます。(aaa : 日〜土、aaaa : 日曜日〜土曜日、ddd : Sun〜Sat、dddd : Sunday〜Saturday) // 【戻値】 // 数値をフォーマットした文字列を返します。 ////////////////////////////////////////////////// FUNCTION getWeekdayName(num, format = "aaa") DIM re = CREATEOLEOBJ("VBScript.RegExp") DIM aaa[] = "日", "月", "火", "水", "木", "金", "土"; DIM aaaa[] = "日曜日", "月曜日", "火曜日", "水曜日", "木曜日", "金曜日", "土曜日"; DIM ddd[] = "Sun", "Mon", "Tue", "Wed", "Thu", "Fri", "Sat"; DIM dddd[] = "Sunday", "Monday", "Tuesday", "Wednesday", "Thursday", "Friday", "Saturday"; match = reExecute(format, "(a|d)+") type = match.Item(0).Value RESULT = reReplace(format, EVAL(type + "[" + num + "]"), type) FEND ////////////////////////////////////////////////// // 【引数】 // シリアル値 : 時間を表すシリアル値を指定 // 【戻値】 // 例)0…00、0.5…12、1.0…24 ////////////////////////////////////////////////// FUNCTION Hour(serial) RESULT = REPLACE(FORMAT(INT(serial * 24), 2), " ", "0") FEND ////////////////////////////////////////////////// // 【引数】 // 日付 : YYYYMMDD // 【戻値】 // TRUE : 日付として認識できる、FALSE : 日付として認識できない ////////////////////////////////////////////////// FUNCTION isDate(date) Pattern = "^(?!([02468][1235679]|[13579][01345789])000229)(([0-9]{4}(01|03|05|07|08|10|12)(0[1-9]|[12][0-9]|3[01]))|([0-9]{4}(04|06|09|11)(0[1-9]|[12][0-9]|30))|([0-9]{4}02(0[1-9]|1[0-9]|2[0-8]))|([0-9]{2}([02468][048]|[13579][26])0229))$" match = reExecute(date, Pattern) SELECT match.Count CASE 1 res = TRUE CASE 0 res = FALSE SELEND RESULT = res FEND ////////////////////////////////////////////////// // 【引数】 // JD : ユリウス日 // 【戻値】 // グレゴリオ暦を格納した配列(0 : 年, 1 : 月, 2 : 日, 3 : 時, 4 : 分, 5 : 秒) ////////////////////////////////////////////////// FUNCTION JDToYMD(JD) DIM x0 = INT(JD + 68570) DIM x1 = INT(x0 / 36524.25) DIM x2 = x0 - INT(36524.25 * x1 + 0.75) DIM x3 = INT((x2 + 1) / 365.2425) DIM x4 = x2 - INT(365.25 * x3) + 31 DIM x5 = INT(INT(x4) / 30.59) DIM x6 = INT(INT(x5) / 11) DIM t2 = x4 - INT(30.59 * x5) DIM t1 = x5 - 12 * x6 + 2 DIM t0 = 100 * (x1 - 49) + x3 + x6 IFB t1 = 2 AND t2 > 28 THEN SELECT TRUE CASE t0 MOD 100 = 0 AND t0 MOD 400 = 0 t2 = 29 CASE t0 MOD 4 = 0 t2 = 29 DEFAULT t2 = 28 SELEND ENDIF DIM tm = 86400 * (JD - INT(JD)) DIM t3 = INT(tm / 3600) DIM t4 = INT((tm - 3600 * t3) / 60) DIM t5 = INT(tm - 3600 * t3 - 60 * t4) DIM t[] = t0, t1, t2, t3, t4, t5 RESULT = SLICE(t) FEND ////////////////////////////////////////////////// // 【引数】 // JC : ユリウス世紀 // 【戻値】 // 太陽黄経(λsun) ////////////////////////////////////////////////// FUNCTION longitudeSun(JC) DIM A[18] = 0.0004, 0.0004, 0.0005, 0.0005, 0.0006, 0.0007, 0.0007, 0.0007, 0.0013, 0.0015, 0.0018, 0.0018, 0.0020, 0.0200, -0.0048*JC, 1.9147, 36000.7695*JC, 280.4659 DIM k[18] = 31557.0, 29930.0, 2281.0, 155.0, 33718.0, 9038.0, 3035.0, 65929.0, 22519.0, 45038.0, 445267.0, 19.0, 32964.0, 71998.1, 35999.05, 35999.05, 0, 0 DIM θ0[18] = 161.0, 48.0, 221.0, 118.0, 316.0, 64.0, 110.0, 45.0, 352.0, 254.0, 208.0, 159.0, 158.0, 265.1, 267.52, 267.52, 0, 0 DIM λsun[18] FOR n = 0 TO 15 DIM ang = normalizeAngle(k[n] * JC + θ0[n]) λsun[n] = A[n] * COS(degToRad(ang)) NEXT λsun[16] = normalizeAngle(A[16]) λsun[17] = normalizeAngle(A[17]) RESULT = normalizeAngle(CALCARRAY(λsun, CALC_ADD)) FEND ////////////////////////////////////////////////// // 【引数】 // deg : 度数法 // 【戻値】 // 度単位の角度を0~360度に正規化 ////////////////////////////////////////////////// FUNCTION normalizeAngle(deg) SELECT TRUE CASE deg >= 360 deg = deg - INT(deg / 360) * 360 CASE deg < 0 deg = deg + INT(ABS(deg / 360) + 1) * 360 SELEND RESULT = deg FEND ////////////////////////////////////////////////// // 【引数】 // str : 正規表現による検索の対象となる文字列 // Pattern : 正規表現で使用するパターンを設定 // IgnoreCase : 大文字・小文字を区別しない場合はTrue、区別する場合はFalse // Global : 文字列全体を検索する場合はTrue、しない場合はFalse // 【戻値】 // 正規表現で検索した結果をMatchesコレクションとして返します。 ////////////////////////////////////////////////// FUNCTION reExecute(str, Pattern, IgnoreCase = TRUE, Global = TRUE) DIM re = CREATEOLEOBJ("VBScript.RegExp") re.Pattern = Pattern re.IgnoreCase = IgnoreCase re.Global = Global RESULT = re.Execute(str) FEND ////////////////////////////////////////////////// // 【引数】 // str1 : 置換される文字列 // str2 : 置換後の文字列 // Pattern : 置換する文字列のパターン // IgnoreCase : 大文字・小文字を区別しない場合はTrue、区別する場合はFalse // Global : 文字列全体を検索する場合はTrue、しない場合はFalse // 【戻値】 // 正規表現置換後の文字列 ////////////////////////////////////////////////// FUNCTION reReplace(str1, str2, Pattern, IgnoreCase = TRUE, Global = TRUE) DIM re = CREATEOLEOBJ("VBScript.RegExp") re.Pattern = Pattern re.IgnoreCase = IgnoreCase re.Global = Global RESULT = re.Replace(str1, str2) FEND ////////////////////////////////////////////////// // 【引数】 // str : 正規表現による検索の対象となる文字列 // Pattern : 正規表現で使用するパターンを設定 // IgnoreCase : 大文字・小文字を区別しない場合はTrue、区別する場合はFalse // Global : 文字列全体を検索する場合はTrue、しない場合はFalse // 【戻値】 // 正規表現にマッチするかどうか ////////////////////////////////////////////////// FUNCTION reTest(str, Pattern, IgnoreCase = TRUE, Global = TRUE) DIM re = CREATEOLEOBJ("VBScript.RegExp") re.Pattern = Pattern re.IgnoreCase = IgnoreCase re.Global = Global RESULT = re.Test(str) FEND ////////////////////////////////////////////////// // 【引数】 // inputs : 繰り返す文字列 // multiplier : inputsを繰り返す回数 // 【戻値】 // inputsをmultiplier回を繰り返した文字列を返します ////////////////////////////////////////////////// FUNCTION strRepeat(inputs, multiplier) DIM res = "" FOR n = 1 TO multiplier res = res + inputs NEXT RESULT = res FEND ////////////////////////////////////////////////// // 【引数】 // serial : シリアル値 // format : フォーマット // 【戻値】 // 数値を表示書式に基づいて変換した文字列 ////////////////////////////////////////////////// FUNCTION TEXT(serial, format) SELECT TRUE CASE reTest(format, "\[h+\]") Matches = reExecute(format, "\[(h+)\]") RESULT = TEXT(Hour(serial), strRepeat("0", LENGTH(Matches.Item(0).SubMatches(0)))) CASE reTest(format, "^h+$") Matches = reExecute(format, "^(h+)$") RESULT = TEXT(Hour(serial) MOD 24, strRepeat("0", LENGTH(Matches.Item(0).SubMatches(0)))) CASE reTest(format, "\[m+\]") Matches = reExecute(format, "\[(m+)\]") RESULT = TEXT(serial * 1440, strRepeat("0", LENGTH(Matches.Item(0).SubMatches(0)))) CASE format = "m" GETTIME(serial, "1899/12/30") RESULT = TEXT(G_TIME_MM, "0") CASE format = "mm" GETTIME(serial, "1899/12/30") RESULT = G_TIME_MM2 CASE format = "s" GETTIME(serial, "1899/12/30") RESULT = TEXT(G_TIME_SS, "0") CASE format = "ss" GETTIME(serial, "1899/12/30") RESULT = G_TIME_SS2 CASE format = "yyyy" GETTIME(serial, "1899/12/30") RESULT = G_TIME_YY4 CASE format = "yy" GETTIME(serial, "1899/12/30") RESULT = COPY(G_TIME_YY4, 3, 2) CASE format = "e" SELECT TRUE // 令和 CASE dateDiff("d", "2019/05/01", TEXT(serial, "yyyy/mm/dd")) >= 0 RESULT = TEXT(serial, "yyyy") - 2018 // 平成 CASE dateDiff("d", "1989/01/08", TEXT(serial, "yyyy/mm/dd")) >= 0 RESULT = TEXT(serial, "yyyy") - 1988 // 昭和 CASE dateDiff("d", "1926/12/25", TEXT(serial, "yyyy/mm/dd")) >= 0 RESULT = TEXT(serial, "yyyy") - 1925 // 大正 CASE dateDiff("d", "1912/07/30", TEXT(serial, "yyyy/mm/dd")) >= 0 RESULT = TEXT(serial, "yyyy") - 1911 // 明治 CASE dateDiff("d", "1868/01/25", TEXT(serial, "yyyy/mm/dd")) >= 0 RESULT = TEXT(serial, "yyyy") - 1867 SELEND CASE format = "g" SELECT TRUE // 令和 CASE dateDiff("d", "2019/05/01", TEXT(serial, "yyyy/mm/dd")) >= 0 RESULT = "R" // 平成 CASE dateDiff("d", "1989/01/08", TEXT(serial, "yyyy/mm/dd")) >= 0 RESULT = "H" // 昭和 CASE dateDiff("d", "1926/12/25", TEXT(serial, "yyyy/mm/dd")) >= 0 RESULT = "S" // 大正 CASE dateDiff("d", "1912/07/30", TEXT(serial, "yyyy/mm/dd")) >= 0 RESULT = "T" // 明治 CASE dateDiff("d", "1868/01/25", TEXT(serial, "yyyy/mm/dd")) >= 0 RESULT = "M" SELEND CASE format = "gg" RESULT = COPY(TEXT(serial, "ggg"), 1, 1) CASE format = "ggg" SELECT TRUE // 令和 CASE dateDiff("d", "2019/05/01", TEXT(serial, "yyyy/mm/dd")) >= 0 RESULT = "令和" // 平成 CASE dateDiff("d", "1989/01/08", TEXT(serial, "yyyy/mm/dd")) >= 0 RESULT = "平成" // 昭和 CASE dateDiff("d", "1926/12/25", TEXT(serial, "yyyy/mm/dd")) >= 0 RESULT = "昭和" // 大正 CASE dateDiff("d", "1912/07/30", TEXT(serial, "yyyy/mm/dd")) >= 0 RESULT = "大正" // 明治 CASE dateDiff("d", "1868/01/25", TEXT(serial, "yyyy/mm/dd")) >= 0 RESULT = "明治" SELEND CASE format = "mmmmm" RESULT = COPY(TEXT(serial, "mmmm"), 1, 1) CASE format = "mmmm" SELECT TEXT(serial, "m") CASE 1; RESULT = "January" CASE 2; RESULT = "February" CASE 3; RESULT = "March" CASE 4; RESULT = "April" CASE 5; RESULT = "May" CASE 6; RESULT = "June" CASE 7; RESULT = "July" CASE 8; RESULT = "August" CASE 9; RESULT = "September" CASE 10; RESULT = "October" CASE 11; RESULT = "November" CASE 12; RESULT = "December" SELEND CASE format = "mmm" RESULT = COPY(TEXT(serial, "mmmm"), 1, 3) CASE format = "dd" GETTIME(serial, "1899/12/30") RESULT = TEXT(G_TIME_DD2, "00") CASE format = "d" GETTIME(serial, "1899/12/30") RESULT = TEXT(G_TIME_DD, "0") CASE format = "aaaa" GETTIME(serial, "1899/12/30") RESULT = getWeekdayName(G_TIME_WW, "aaaa") CASE format = "aaa" GETTIME(serial, "1899/12/30") RESULT = getWeekdayName(G_TIME_WW, "aaa") CASE format = "dddd" GETTIME(serial, "1899/12/30") RESULT = getWeekdayName(G_TIME_WW, "dddd") CASE format = "ddd" GETTIME(serial, "1899/12/30") RESULT = getWeekdayName(G_TIME_WW, "ddd") CASE reTest(format, "(0+\.?0+)?%") Matches = reExecute(format, "(0+\.?0+)?%") RESULT = TEXT(serial * 100, Matches.Item(0).SubMatches(0)) + "%" CASE reTest(format, "^(.*?)\/(.*?)$") // 分母、分子 DIM demon, numer // 分数 DIM frac[1] DIM Pattern = "(\d+\.\d*?(\d+)\2+)" IFB reTest(serial, Pattern) THEN Matches = reExecute(serial, Pattern) // 循環節の桁数 digit = LENGTH(Matches.Item(0).SubMatches(1)) denom = POWER(10, digit) - 1 numer = serial * POWER(10, digit) - serial // 分子が小数なら整数に通分する num = POWER(10, LENGTH(numer) - POS(".", numer)) denom = denom * num numer = numer * num ELSE denom = POWER(10, LENGTH(serial) - POS(".", serial)) numer = serial * denom ENDIF frac[0] = denom frac[1] = numer // 最大公約数を求める n = GCD(frac) denom = frac[0] numer = frac[1] // 既約分数にする denom = denom / n numer = numer / n RESULT = numer + "/" + denom CASE reTest(format, "((#+),?)+") RESULT = reReplace(serial, "$1,", "(\d)(?=(\d{3})+(?!\d))") CASE reTest(format, "(0+)\.?(0+)?") Matches = reExecute(format, "(0+)\.?(0+)?") len1 = LENGTH(Matches.Item(0).SubMatches(0)) len2 = LENGTH(Matches.Item(0).SubMatches(1)) DIM arr[] = LENGTH(INT(serial)), len1 IFB POS(".", format) THEN RESULT = REPLACE(FORMAT(serial, CALCARRAY(arr, CALC_MAX) + len2 + 1, len2), " ", "0") ELSE RESULT = REPLACE(FORMAT(serial, CALCARRAY(arr, CALC_MAX)), " ", "0") ENDIF CASE reTest(format, "(.+)?((?:AM?|am?)\/(?:PM?|pm?))(.+)?") DIM arr[2] Matches = reExecute(format, "(.+)?((?:AM?|am?)\/(?:PM?|pm?))(.+)?") res = "" FOR n = 0 TO Matches.Item(0).SubMatches.Count - 1 arr[n] = Matches.Item(0).SubMatches(n) NEXT DIM res[2] FOR n = 0 TO 2 STEP 2 Matches = reExecute(arr[n], "([^yegmdahs]+)?(([yegmdahs])\3{0,})([^yegmdahs]+)?") FOR m = 0 TO Matches.Count - 1 SELECT TRUE // 「m」を含む CASE reTest(Matches.Item(m).SubMatches(1), "m+") // アルファベット部分の前後で「m」が月か分か判断 SELECT TRUE // 「m」が最初に出現するアルファベットなら CASE m = 0 // 次に出現するアルファベットが「d」なら「m」は月 IFB reTest(Matches.Item(m+1).SubMatches(1), "d+") THEN res[n] = res[n] + TEXT(serial, Matches.Item(m).SubMatches(1)) ENDIF // 次に出現するアルファベットが「s」なら「m」は分 IFB reTest(Matches.Item(m+1).SubMatches(1), "s+") THEN res[n] = res[n] + TEXT(INT(serial * 1440) MOD 60, strRepeat("0", LENGTH(Matches.Item(m).SubMatches(1)))) ENDIF // 「m」が2番目(2-1=1)以降に出現するアルファベットなら CASE m >= 1 // 前に出現したアルファベットが「y」なら「m」は月 IFB reTest(Matches.Item(m-1).SubMatches(1), "y+") THEN res[n] = res[n] + TEXT(serial, Matches.Item(m).SubMatches(1)) ENDIF // 前に出現したアルファベットが「h」なら「m」は分 IFB reTest(Matches.Item(m-1).SubMatches(1), "h+") THEN res[n] = res[n] + TEXT(INT(serial * 1440) MOD 60, strRepeat("0", LENGTH(Matches.Item(m).SubMatches(1)))) ENDIF SELEND CASE reTest(Matches.Item(m).SubMatches(1), "^h+$") IFB TEXT(serial, Matches.Item(m).SubMatches(1)) < 12 THEN res[1] = COPY(arr[1], 1, POS("/", arr[1]) - 1) // 午前 ELSE res[1] = COPY(arr[1], POS("/", arr[1]) + 1) // 午後 ENDIF res[n] = res[n] + TEXT(serial, Matches.Item(m).SubMatches(1)) // それ以外ならば DEFAULT res[n] = res[n] + TEXT(serial, Matches.Item(m).SubMatches(1)) SELEND res[n] = Matches.Item(m).SubMatches(0) + res[n] + Matches.Item(m).SubMatches(3) NEXT NEXT RESULT = JOIN(res, "") CASE reTest(format, "(([yegmdahs])\2{0,})([^yegmdahs]+)?") res = "" Matches = reExecute(format, "(([yegmdahs])\2{0,})([^yegmdahs]+)?") FOR n = 0 TO Matches.Count - 1 // mの前がyもしくはmの後がdならば→mは月 SELECT TRUE // 「m」を含む CASE reTest(Matches.Item(n).SubMatches(0), "m+") // アルファベット部分の前後で「m」が月か分か判断 SELECT TRUE // 「m」が最初に出現するアルファベットなら CASE n = 0 // 次に出現するアルファベットが「d」なら「m」は月 IFB reTest(Matches.Item(n+1).SubMatches(0), "d+") THEN res = res + TEXT(serial, Matches.Item(n).SubMatches(0)) ENDIF // 次に出現するアルファベットが「s」なら「m」は分 IFB reTest(Matches.Item(n+1).SubMatches(0), "s+") THEN res = res + TEXT(INT(serial * 1440) MOD 60, strRepeat("0", LENGTH(Matches.Item(n).SubMatches(0)))) ENDIF // 「m」が2番目(2-1=1)以降に出現するアルファベットなら CASE n >= 1 // 前に出現したアルファベットが「y」なら「m」は月 IFB reTest(Matches.Item(n-1).SubMatches(0), "y+") THEN res = res + TEXT(serial, Matches.Item(n).SubMatches(0)) ENDIF // 前に出現したアルファベットが「h」なら「m」は分 IFB reTest(Matches.Item(n-1).SubMatches(0), "h+") THEN res = res + TEXT(INT(serial * 1440) MOD 60, strRepeat("0", LENGTH(Matches.Item(n).SubMatches(0)))) ENDIF SELEND // それ以外ならば DEFAULT res = res + TEXT(serial, Matches.Item(n).SubMatches(0)) SELEND res = res + Matches.Item(n).SubMatches(2) NEXT RESULT = res SELEND FEND ////////////////////////////////////////////////// // 【引数】 // year : 年 // month : 月 // day : 日 // hour : 時 // minute : 分 // second : 秒 // 【戻値】 // ユリウス日 ////////////////////////////////////////////////// FUNCTION YMDToJD(year, month, day, hour = 0, minute = 0, second = 0) IFB month < 3 THEN year = year - 1 month = month + 12 ENDIF DIM JD = INT(year * 365.25) JD = JD + INT(year / 400) JD = JD - INT(year / 100) JD = JD + INT((month - 2) * 30.59) JD = JD + 1721088 JD = JD + day DIM t = second / 3600 t = t + minute / 60 t = t + hour t = t / 24 JD = JD + t RESULT = JD FEND
結果
彼岸入り	2020/09/19
秋分の日	2020/09/22
彼岸明け	2020/09/25

この記事は役に立ちましたか?
役に立った 役に立たなかった