\u25C7\u4E1C\u7ECF:'+jing+'.'+jings+'\u5EA6 \u25C7\u5317\u7EAC:'+wei+'.'+weis+'\u5EA6
'; } } //function hrsmin(hours) { // var hrs, h, m, dum; // hrs = Math.floor(hours * 60 + 0.5)/ 60.0; // h = Math.floor(hrs); // m = Math.floor(60 * (hrs - h) + 0.5); // if(h<10)h="0"+h; // if(m<10)m="0"+m; // dum = h +"\u65F6"+ m+"\u5206"; // // if (dum < 1000) dum = "0" + dum; // if (dum <100) dum = "0" + dum; // if (dum < 10) dum = "0" + dum; // return dum; // } //function ipart(x) { // var a; // if (x> 0) { // a = Math.floor(x); // } // else { // a = Math.ceil(x); // } // return a; // } // // //function frac(x) { // var a; // a = x - Math.floor(x); // if (a < 0) a += 1; // return a; // } function round(num, dp) { return Math.round (num * Math.pow(10, dp)) / Math.pow(10, dp); } function range(x) { var a, b; b = x / 360; a = 360 * (b - ipart(b)); if (a < 0 ) { a = a + 360 } return a } function mjd(day, month, year, hour) { var a, b; if (month <= 2) { month = month + 12; year = year - 1; } a = 10000.0 * year + 100.0 * month + day; if (a <= 15821004.1) { b = -2 * Math.floor((year + 4716)/4) - 1179; } else { b = Math.floor(year/400) - Math.floor(year/100) + Math.floor(year/4); } a = 365.0 * year - 679004.0; return (a + b + Math.floor(30.6001 * (month + 1)) + day + hour/24.0); } function quad(ym, yz, yp) { var nz, a, b, c, dis, dx, xe, ye, z1, z2, nz; var quadout = new Array; nz = 0; a = 0.5 * (ym + yp) - yz; b = 0.5 * (yp - ym); c = yz; xe = -b / (2 * a); ye = (a * xe + b) * xe + c; dis = b * b - 4.0 * a * c; if (dis > 0) { dx = 0.5 * Math.sqrt(dis) / Math.abs(a); z1 = xe - dx; z2 = xe + dx; if (Math.abs(z1) <= 1.0) nz += 1; if (Math.abs(z2) <= 1.0) nz += 1; if (z1 < -1.0) z1 = z2; } quadout[0] = nz; quadout[1] = z1; quadout[2] = z2; quadout[3] = xe; quadout[4] = ye; return quadout; } function lmst(mjd, glong) { var lst, t, d; d = mjd - 51544.5 t = d / 36525.0; lst = range(280.46061837 + 360.98564736629 * d + 0.000387933 *t*t - t*t*t / 38710000); return (lst/15.0 + glong/15); } function minisun(t) { var p2 = 6.283185307, coseps = 0.91748, sineps = 0.39778; var L, M, DL, SL, X, Y, Z, RHO, ra, dec; var suneq = new Array; M = p2 * frac(0.993133 + 99.997361 * t); DL = 6893.0 * Math.sin(M) + 72.0 * Math.sin(2 * M); L = p2 * frac(0.7859453 + M / p2 + (6191.2 * t + DL)/1296000); SL = Math.sin(L); X = Math.cos(L); Y = coseps * SL; Z = sineps * SL; RHO = Math.sqrt(1 - Z * Z); dec = (360.0 / p2) * Math.atan(Z / RHO); ra = (48.0 / p2) * Math.atan(Y / (X + RHO)); if (ra <0 ) ra += 24; suneq[1] = dec; suneq[2] = ra; return suneq; } function sin_alt(iobj, mjd0, hour, glong, cglat, sglat) { var mjd, t, ra, dec, tau, salt, rads = 0.0174532925; var objpos = new Array; mjd = mjd0 + hour/24.0; t = (mjd - 51544.5) / 36525.0; if (iobj == 1) { objpos = minimoon(t); } else { objpos = minisun(t); } ra = objpos[2]; dec = objpos[1]; tau = 15.0 * (lmst(mjd, glong) - ra); salt = sglat * Math.sin(rads*dec) + cglat * Math.cos(rads*dec) * Math.cos(rads*tau); return salt; } function getzttime(mjd, tz, glong) { var sglong, sglat, date, ym, yz, utrise, utset, j; var yp, nz, hour, z1, z2, iobj, rads = 0.0174532925; var quadout = new Array; sinho = Math.sin(rads * -0.833); date = mjd - tz/24; hour = 1.0; ym = sin_alt(2, date, hour - 1.0, glong, 1, 0) - sinho; while(hour < 25) { yz = sin_alt(2, date, hour, glong, 1, 0) - sinho; yp = sin_alt(2, date, hour + 1.0, glong, 1, 0) - sinho; quadout = quad(ym, yz, yp); nz = quadout[0]; z1 = quadout[1]; z2 = quadout[2]; xe = quadout[3]; ye = quadout[4]; if (nz == 1) { if (ym < 0.0) utrise = hour + z1; else utset = hour + z1; } if (nz == 2) { if (ye < 0.0) { utrise = hour + z2; utset = hour + z1; } else { utrise = hour + z1; utset = hour + z2; } } ym = yp; hour += 2.0; } var zt=(utrise+utset)/2; if(zt
\u25CE\u65E5\u51FA:' + hrsmin(utrise); // else outstring += "\u25CE\u65E5\u51FA:\u65E0\u65E5\u51FA"+" "; // if (sett == true) outstring += " \u25CE\u65E5\u843D:" + hrsmin(utset)+"
"; // else outstring += " \u25CE\u65E5\u843D:\u65E0\u65E5\u843D
"; // var zt=getzttime(mjd, tz, glong) // outstring+= " \u25CE\u65E5\u4E2D:" + hrsmin(zt); // outstring += " \u25CE\u663C\u957F:"+hrsmin(utset-utrise); // // } // else { // if (above == true){ // outstring += always_up; // var zt=getzttime(mjd, tz, glong); // outstring+="
"+"\u65E5\u4E2D:"+hrsmin(zt); // outstring += "
\u663C\u957F:"+hrsmin(24); // } // else outstring += always_down; // // } return outstring; } function initial() { dStyle = document.getElementById("detail").style; ccc(); document.getElementById("CLD").SY.selectedIndex=tY-1900; document.getElementById("CLD").SM.selectedIndex=tM; drawCld(tY,tM); //tick(); } function getString(isString,n){ var isResult=isString.substr(isString.length-n,1); return isResult; } function deltaT(year){return (58 + 0.8 * (year - 1990));} function MoonLong(T) { var kaku,i,ans; T /= 36525.0; ans = 0.0; for (i = 62; i >= 1 ; i--) { kaku = (mlb[i] * T + mlc[i]) * d2r; ans += mla[i] * Math.cos(kaku); } kaku = (mlb[0] * T + mlc[0] * d2r); ans += mla[0] * T * Math.cos(kaku); ans = ans - Math.floor(ans / 360.0) * 360.0; return ans; } function MoonLat(T) { var kaku,i,ans; T /= 36525.0; ans = 0.0; for (i = 44; i >= 0 ; i--) { kaku = (mab[i] * T + mac[i]) * d2r; ans += maa[i] * Math.cos(kaku); } return ans; } function MoonHP(T) { var kaku,i,ans; T /= 36525.0; ans = 0.0; for (i = 43; i >= 0 ; i--) { kaku = (mhb[i] * T + mhc[i]) * d2r; ans += mha[i] * Math.cos(kaku); } return ans; } function GetEpci(T) { var Epci; T /= 36525.0; Epci = 0.00256 * Math.cos((1934 * T + 235) * d2r) + 0.00015 * Math.cos((72002 * T + 201) * d2r); Epci += 23.43928 - 0.01301 * T; return Epci * d2r; } function GetGMT(UT,T) { var gmt; gmt = UT + 6.69736 + 24.000513 * T / 365.25; gmt += 0.00029 * Math.sin((1934 * T / 36525 + 235) * d2r); gmt = gmt / 24.0; gmt = (gmt - Math.floor(gmt)) * 24.0; return gmt; } function POL2XYZ(p,x) { with(Math) { x[2] = p[0] * sin(p[2]); x[1] = p[0] * cos(p[2]); x[0] = x[1] * cos(p[1]); x[1] *= sin(p[1]); } } function XYZ2POL(x,p) { var r; r = x[0] * x[0] + x[1] * x[1]; with(Math) { p[0] = sqrt(x[2] * x[2] + r); r = sqrt(r); p[2] = atan2(x[2],r); p[1] = atan2(x[1],x[0]); } } function ROTXYZ(x,p,base) { var nx = new Array(3); var e1,e2; if (base == 0) {e1 = 1; e2 = 2;} else if (base == 1) {e1 = 0; e2 = 2;} else {e1 = 0; e2 = 1;} with(Math) { nx[base] = x[base]; nx[e1] = cos(p) * x[e1] - sin(p) * x[e2]; nx[e2] = sin(p) * x[e1] + cos(p) * x[e2]; } for (e1 in nx) x[e1] = nx[e1]; } function Koudou2Sekidou(Kou,Sek,T,FLAG) { var e; var xyz = new Array(3); e = GetEpci(T); if (FLAG != 1) e = -e; POL2XYZ(Kou,xyz); ROTXYZ(xyz,e,0); XYZ2POL(xyz,Sek); } var Now_dcd; function GetEL(T,ra,dc,lng,lat) { var hangl,el,latd,hour; hour = T - 0.5; hour = (hour - Math.floor(hour)) * 24.0; hangl = GetGMT(hour,T) * 15.0 * d2r; hangl = ra - (hangl + lng); Now_dcd = 90 * d2r - dc; latd = 90 * d2r - lat; with (Math) { el = cos(latd) * cos(Now_dcd) + sin(latd) * sin(Now_dcd) * cos(hangl); el = atan2(sqrt(1.0 - el * el),el); } return Math.PI/2.0 - el; } function addZero(n) { if (n < 10) return '0' + n; else return '' + n; } function RiseSB(ts,dtime,lng,lat,objectflag) { var dtimemin,TD; var firstflag; var dt,t0,t1; var el0,el1; var elbase; var Kou = new Array(3); var Sek = new Array(3); firstflag = 1; dt = deltaT(yy) / 86400.0; dtimemin = 0.5 / 1440; t1 = ts; el0 = 0; while (dtime >= dtimemin) { TD = t1 + dt; if (objectflag == 1) { Kou[0] = 1.0; Kou[1] = SunLong(TD) * d2r; Kou[2] = 0.0; elbase = 0.902 * d2r; } else { elbase = MoonHP(TD); Kou[0] = 1.0 / Math.sin(elbase); Kou[1] = MoonLong(TD) * d2r; Kou[2] = MoonLat(TD) * d2r; elbase = (0.586 - elbase) * d2r; } Koudou2Sekidou(Kou,Sek,TD,1); el1 = GetEL(t1,Sek[1],Sek[2],lng,lat) + elbase; if (firstflag > 0) { el0 = el1; firstflag = 0; } if (el0 * el1 >= 0) { t0 = t1; el0 = el1; } else { dtime /= 2.0; } t1 = t0 + dtime; } RiseFlag = el1; return (t0 + dtime); } function MoonRise(v,Te){ var T,tmp,moont,moond,XZ,sitelong,sitelat,jwnum; var sObj=document.getElementById('SD'+ v); var d=sObj.innerHTML-1; if(sObj.innerHTML!='') { yy=cld[d].sYear; mm=cld[d].sMonth; jwnum=document.getElementById("CLD").per_des; if(jwnum==undefined){ sitelong = 113.07; sitelat = 23.02; }else{ XZ=getString(jwnum,5); if(XZ==0)XZ=''; sitelong=(XZ+''+getString(jwnum,4)+''+getString(jwnum,3)+'.'+getString(jwnum,2)+''+getString(jwnum,1))*1; sitelat=(getString(jwnum,9)+''+getString(jwnum,8)+'.'+getString(jwnum,7)+''+getString(jwnum,6))*1; } if(Te==7){dd =tD;}else{dd =cld[d].sDay;}} RiseTime = -999; SetTime = -999; T = Ymd2Jd(yy,mm,dd) - Ymd2Jd(2000,1,1.5) - 9/24.0; tmp = (RiseSB(T,(8.0/24.0),(sitelong * d2r),(sitelat * d2r),2) - T) * 24.0; if (RiseFlag < 0) { tmp = (RiseSB((T + (tmp + 1)/24),(8.0/24.0),(sitelong * d2r),(sitelat * d2r),2) - T) * 24.0; } RiseTime = tmp; tmp = (RiseSB(T,(8.0/24.0),(sitelong * d2r),(sitelat * d2r),2) - T) * 24.0; if (RiseFlag >= 0) { tmp = (RiseSB((T + (tmp + 0.5)/24),(8.0/24.0),(sitelong * d2r),(sitelat * d2r),2) - T) * 24.0; } SetTime = tmp; RiseTime += 0.5 / 1440.0; SetTime += 0.5 / 1440.0; /*if ((RiseTime >= 0) && (RiseTime < 24.0)) { hour = Math.floor(RiseTime); min = Math.floor((RiseTime - hour) * 60.0); moont='\u25CE\u6708\u51FA:'+addZero(hour) + '\u65F6' + addZero(min) + '\u5206
'; } else { moont='\u25CE\u6708\u51FA:00\u65F600\u5206
';} if ((SetTime >= 0) && (SetTime < 24.0)) { hour = Math.floor(SetTime); min = Math.floor((SetTime - hour) * 60.0); moond='\u25CE\u6708\u843D:'+addZero(hour) + '\u65F6' + addZero(min) + '\u5206'; } else { moond='\u25CE\u6708\u843D:00\u65F600\u5206';} document.getElementById("rcrl2").innerHTML ='
'+moont + moond;*/ }