!#INC default.dat # ! эта строка должна быть в начале каждого задания !======================================================================= ! ! РАЗМЕЩЕНИЕ ИСХОДНЫХ ТАБЛИЦ МОДЕЛЕЙ ФОТОСФЕРЫ ! #MOD_ROOT APP # ! размещение моделей APP/PROJ = общие модели/специфичные для данного проекта #MOD_DIR MODS # ! подкаталог моделей !----------------------------------------------------------------------- ! ! ДАННЫЕ ДЛЯ ПЕРЕРАСЧЕТА МОДЕЛИ ФОТОСФЕРЫ ! ! Задача =ION= может быть запущена отдельно с помощью #RUN ION # ! или в цепочке задачи =CONT= с помощью #RUN CONT # при параметре #ION = 1 ! !------- ОПЦИЯ ПРОГРАММЫ - перерасчет Pg и Pe из условия гидростатического равновесия #ION 0 # ! 1/0 = надо пересчитывать/не надо пересчитывать !======================================================================= ! ДАЛЕЕ СЛЕДУЕТ РЯД ПАРАМЕТРОВ, КОТОРЫЕ ЗАДЕЙСТВУЮТСЯ, КОГДА #ION = 1 ! !------- ускорение силы тяжести ----------------------------------------- #GRAVITY 2.738e4 # ! 2.738E4 - ускорение свободного падения на Солнце !------- таблица обилий ------------------------------------------------- #ABUND T HOLMU # ! обилия в системах HOLMU и ! T SW ! Stellmacher - Wiehr !-------- список атомов для расчета среднего атомного веса солнечного вещества ------ #ATP H HE LI BE C N O NA MG AL = SI S CA SC TI V CR MN FE CO = NI CU ZN Y ZR NE AR K P F = RB SR MO RU CD BA OS PB # !-------- список атомов для расчета ионизации солнечного вещества ------------------- #ATI H HE C O NA MG AL SI S FE = CA NI CR K # !-------- средн.ат.масса на один атом (ядро?) водорода ------------------------------ ! эта величина обычно вычисляется из таблицы обилий и ! списка атомов для расчета среднего атомного веса ! вручную заполняется, только если её надо задать принудительно !!! ! !#MJUPERH 1.435 # ! !------- m0 - вес столба атмосферы над верхней границей модели --------------------- #MAS0 5.0E-2 # ! возможные значения: 0, (1-7 e-5) !------- если #ION = 1 и мы будем делать перерасчет газового и эл.давлений, то требуется !------- исходная таблица, откуда берется температурное распределения для пересчета !------- и формула преобразования этой таблицы ! ! #PREMOD T KNEER =# #MOD T KNEER_Recalc =# ! вариант : !#PREMOD Q HOLMU 0.2 =# !#MOD I HMP02 =# ! таблица HMP02 получается после предварительного запуска ION и ручного переименования результата! !======================================================================= ! дальнейшие параметры требуются независимо от устанвки #ION !----------------------------------------------------------------------- ! ! собственно ТАБЛИЦА МОДЕЛИ ФОТОСФЕРЫ ! !----------------------------------------------------------------------- ! Модели фотосферы, а также параметры м.поля и скоростей в виде таблиц !----------------------------------------------------------------------- #MOD T HOLMU =# !======================================================================= ! УЧЕТ НЕПРЕРЫВНОГО ПОГЛОЩЕНИЯ !----------------------------------------------------------------------- #OPACITY LANDI # ! LANDI/ALLEN/ALLER = ф-лы LANDI / таблицы ALLEN / таблицы ALLER к-тов непрер.поглощ. (отриц.ионом водорода) !======================================================================= ! РАСЧЕТ ЗАТУХАНИЯ В ЛИНИИ !----------------------------------------------------------------------- #A_OPT 46_04 # ! EWAN/46_03/46_04/46+03/46+04 = с учетом EWAN/Г4+Г6 ; (T/10000)^0.3/0.4/-0.3/-0.4 #EWAN 2.0 # ! поправка к параметру затухания (используется, только когда #A_OPT = EWAN) !======================================================================= ! ГЛОБАЛЬНЫЕ ПАРАМЕТРЫ РАСЧЕТА !----------------------------------------------------------------------- #CTETA 1.00 # ! косинус гелиоцентрического угла #VMA 0.0 # ! макро-турбулент.скорости !======================================================================= ! ЗАДАНИЕ ПОЛЯ И СКОРОСТЕЙ !----------------------------------------------------------------------- #H 0.0 # ! напряженностьь поля #GM 0.0 # ! угол к лучу зрения #XI 0.0 # ! азимут поля #VD 0.0 # ! лучевые скорости #VMI 1.0 # ! микро-турбулент.скорости !======================================================================= ! ЛИНИЯ И ЕЁ ПАРАМЕТРЫ !----------------------------------------------------------------- #LIN FE6302 # ! идентификатор линии (ЭЛ-Т+дл.волны) !#EWAN == 3.5 # !#LGGF == 4.48 # #LGGF 12.0 # ! если #LGGF = 12 , то lg(gf) берется из таблицы. описывающей сп.линию #LAMB 0 # ! если #LAMB=0, то дл.волны нач.точки профиля = длине волны линии #LIN !----------------------------------------------------------------- ! Метод формирования массива длин волн ! ! ВАРИАНТ ! в этом случае симметрия профиля не учитывается ! #LAM_L 6148.900 # ! а в расчет входят все линии, какие есть в БД на данном отрезке ! #LAM_U 6149.500 # ! #LAMB 6149.000 # ! если LAN_L, LAM_U не заданы, то вводим данные о массиве длин волн расчета из переменных: #NLAM 61 # ! N : ЧИСЛО ТОЧЕК ПРОФИЛЯ #DL0 0. # ! DL0 : ДЕЛЬТА_ЛЯМБДА_0 #DLK 240. # ! DLK : ДЕЛЬТА_ЛЯМБДА_КОНЕЧНОЕ !----------------------------------------------------------- ! Управление форматами данных !----------------------------------------------------------- #FMOD 1 # ! номер фомата модели (таблицы глубин) для расчета в справочнике форматов -см.файл ..\once\rout.dat #FOUT 0 # ! номер сетки вывода -см.файл ..\once\rout.dat !----------------------------------------------------------- ! Управление точностью: !----------------------------------------------------------- #ERR 1E-6 # ! относительная ошибка 0.001 2E-6 #MAXSTEP 0.1 # ! MAX ШАГ ИНТЕГРИРОВАНИЯ #X0LO 1.5 # ! 1.8 -------поднимаем начало интегрирования до X0LO ------ #X0UP -3.0 # ! 0.8 ------ опускаем начало интегрирования до X0UP ------ !----------------------------------------------------------- ! Опция, определяющая учет аномальной дисперсии !----------------------------------------------------------- #KVOIG 6 # ! 1/2/4/8 = безАД / *"-" / *2 / *1/H(0,A) !#RUN CONT #