! комментарии начинаются с восклицательного знака !-------------------------------------------------- ! ! Это список производных карт, которые вычисляются автоматически ! на основе существующих карт (см. файл map_List.txt) ! !-------------------------------- ! Расчёт продольного поля методом центров тяжести ! из положений крыльев профилей Стокса I+V и I-V ! на середине высоты профиля ! верхнюю точку берём на уровне 0.95 от уровня непрерывного спектра ! BIV1=mul(sub(LIVp195,LIVm195),348.23) ! напряжённость продольного поля для линии 6301 BIV2=mul(sub(LIVp295,LIVm295),233.6195) ! напряжённость продольного поля для линии 6302 BIV=mid(BIV2,BIV1) ! средняя напряжённость для двух линий aBIV=abs(BIV) ! средняя напряжённость без учёта полярности aBIV21=sub(abs(BIV2),abs(BIV1)) ! абсолютная величина градиента напряжённости dB=mul(sub(BIV1,BIV2),sig(BIV)) ! градиент напряжённости с учётом знака !-------------------------------- ! лучевые скорости км/с ! из положения середин крыльев профиля интенсивности LWb LWr - синего и красного крыла ! VW1=mul(add(mul(add(mid(LWb1,LWr1),-56.5),0.021549),+0.5792),47.57477079) ! скорость в линии 6301 VW2=mul(add(mul(add(mid(LWb2,LWr2),-56.5),0.021549),-0.4132),47.5672796) ! скорость в линии 6302 VW=mid(VW1,VW2) ! средняя по 2-м линиям !-------------------------------- ! нормированные (заданные относительно среднего по карте) лучевые скорости км/с ! - скорости после вычитания средних по карте скоростей $VW1 и $VW1. ! (Средние взяты только для немагнитных элементов) ! VW_1=sub(VW1,$VW1) ! нормированная скорость для 6301 VW_2=sub(VW2,$VW2) ! нормированная скорость для 6302 VW_=mid(VW_1,VW_2) ! градиент скорости (от верхнего уровня к нижнему) VW_21=sub(sub(VW2,VW1),0.16) ! градиент скорости от нижнему слоя к верхнему с учётом систематической разницы для двух линий 160 м/с !-------------------------------- ! для анализа скоростей по компонентам V и QU профилей ! из значений V,Q,U профилей переходим к абсолютным значениям поляризаций ! находим положения центров тяжести этих профилей в пикселях спектрографа ! переводим значения в длины волн, а затем в км/с ! VQUa1=mul(add(mul(add(LQUa195,-56.5),0.021549),+0.5792),47.57477079) ! скорость центра тяжести QU 6301 VQUa2=mul(add(mul(add(LQUa295,-56.5),0.021549),-0.4132),47.5672796) ! скорость центра тяжести QU 6302 VVa1=mul(add(mul(add(LVMa195,-56.5),0.021549),+0.5792),47.57477079) ! скорость центра тяжести V 6301 VVa2=mul(add(mul(add(LVMa295,-56.5),0.021549),-0.4132),47.5672796) ! скорость центра тяжести V 6302 VQUc1=mul(add(mul(add(LQUc195,-56.5),0.021549),+0.5792),47.57477079) ! скорость ц.т. центрального лепестка QU 6301 VQUc2=mul(add(mul(add(LQUc295,-56.5),0.021549),-0.4132),47.5672796) ! скорость ц.т. центрального лепестка QU 6302 VVa=mid(VVa1,VVa2) ! средняя скорость V профиля для 2-х линий VQUa=mid(VQUa1,VQUa2) ! средняя скорость Q+U профиля WQ95=add(add(abs(wQp195),abs(wQn195)),add(abs(wQp295),abs(wQn295))) WU95=add(add(abs(wUp195),abs(wUn195)),add(abs(wUp295),abs(wUn295))) WV95=add(add(abs(wVp195),abs(wVn195)),add(abs(wVp295),abs(wVn295))) !Va=rat(sub(add(wVp195,wVp295),add(wVn195,wVn295)),WV95) !Va2=rat(sub(wVp295,wVn295),add(abs(wVp295),abs(wVn295))) !!!!!!!WQU95=add(WQ95,WU95) !!!!!!!WVQU=add(WV95,WQU95) W2QU=add(W2QU195,W2QU295) W2QUV=add(W2QUV295,W2QUV195) QsV=sub(W2QU,WV95) wV95=add(add(wVp195,wVn195),add(wVp295,wVn295)) ! знак wVp и wWn одинаковый для антисимметричного профиля! sNCP=add(sub(wVp195,wVn195),sub(wVp295,wVn295)) CB0=rat(CONT,$CB0) VW_12=sub(VW_1,VW_2) rR=rat(sqrt(add(sq(X_Coordinate),sq(Y_Coordinate))),$SOLAR_RA) !FWHM2=sub(LWr2,LWb2) !FWHM1=sub(LWr1,LWb1) !FWHM=add(FWHM1,FWHM2) !!d0_2=rat(d0Full2,CONT) !!d0_1=rat(d0Full1,CONT) !!d0_12=sub(d0_1,d0_2) !!rRTh=add(neg(sqrt(add(neg(mul(rR_LM,rR_LM)),1))),1) ! = 1-SQRT(1-rR^2) !======= В Р Е М Е Н Н О ========>>> aBIV1=abs(BIV1) aBIV2=abs(BIV2) kBIV12=rat(aBIV1,aBIV2) Blos=mul(Field_Strength,cos(Field_Inclination)) Btang=mul(Field_Strength,sin(Field_Inclination)) aBlos=abs(Blos) kMEIV=rat(aBlos,aBIV) aInc=add(mul(abs(sub(Field_Inclination,90)),-1),90) B_los=mul(Blos,Stray_Light_Fill_Factor) aB_los=abs(B_los) B_tang=mul(Btang,Stray_Light_Fill_Factor) Bflux=mul(Field_Strength,Stray_Light_Fill_Factor) !k_MEIV=rat(aB_los,aBIV) !lg_MEIV=log(k_MEIV) !B_QsV=mul(sub(abs(Btang),aBlos),Stray_Light_Fill_Factor) ! !lgMEIV=log(kMEIV) !lgBIV12=log(kBIV12) !laBIV=log(aBIV) !laBIV1=log(aBIV1) !laBIV2=log(aBIV2) !======= В Р Е М Е Н Н О ========<<< !asimV2=rat(sub(abs(wVp295),abs(wVn295)),add(wVp295,wVn295)) !AsimV2=abs(asimV2) W+W=add(WI1,WI2) ! сумма экв.ширин kVmi=rat(WI2,WI1) BMV=mid(BMV1,BMV2) ! среднее продольное поле BML=mid(BML1,BML2) ! среднее поперечное поле Gm=mid(Gm1,Gm2) ! среднее значение угла наклона поля к лучу зрения aGm1=add(mul(abs(sub(Gm1,90)),-1),90) aGm2=add(mul(abs(sub(Gm2,90)),-1),90) daGm=sub(aGm1,aGm2) aGm=mid(aGm1,aGm2) dGm=sub(Gm1,Gm2) Xi=mid(Xi1,Xi2) ! средний азимут поля dXi=circsub(Xi1,Xi2,180) ! поворот поля lkBMV=log(abs(rat(BMV1,BMV2))) ! log отношения продольных полей lkBML=log(rat(BML1,BML2)) ! log отношения поперечных полей dBML=sub(BML1,BML2) dAzXi=circsub(Field_Azimuth,Xi,180) BM=sqrt(add(mul(BMV,BMV),mul(BML,BML))) BM1=sqrt(add(mul(BMV1,BMV1),mul(BML1,BML1))) BM2=sqrt(add(mul(BMV2,BMV2),mul(BML2,BML2))) dBM12=sub(BM1,BM2) ddBML12=rat(BML1,BML2) lBML12=log(rat(BML1,BML2))