Допустим, мы хотим видеть карты, которые отображают 
  продольную   Blos    и
  поперечную   Btang   составляющие магнитного поля.

У нас есть карты, которые рассчитаны с помощью модели фотосферы Милна - Эддингтона 
  Field_Strength    - абсолютная величина поля и
  Field_Inclination - угол между вектором поля и лучом зрения

Тогда мы можем рассчитать требуемые карты с помощью процедур преобразования.
На вход процедуры преобразования надо подать строки вида:
  Blos =mul(Field_Strength,cos(Field_Inclination))
  Btang=mul(Field_Strength,sin(Field_Inclination))

Процедура преобразования рекурсивно парсит строки
  она находит два операнда, разделённых запятыми
  когда видит в одном операнде вложенные скобки, то продолжает парсить
  и видит унарную операцию. 
К каждой точке карт процедура применяет sin или cos к Field_Inclination 
  затем делает mul (умножение) значений Field_Strength и sin/cos(Field_Inclination)
Результат записывается в структуры, которые содержат карты. 
Этим вновь созданным структурам (=картам) присваиваются имена   Blos   и   Btang,
Теперь эти карты можно использовать анлогично картам, прочитаным из FITS файла.