#!/bin/bash
# Purpose: converts SEL coordinates to lunar centric latitude,longitude,altitude, br,btheta, and bphi
#usage: ./llabrth_moon.bash

for approach in `cat file_list`
do

cat >llabrth_moon.f <<EOF
*       
*     llabrth_moon.f 
*********************************************************************
*   From Level 1b daily data calculates Br, Bth, and Bph
*      and selongraphic latitude, longitude, and altitude
*      using a mean radius of 1737.1 km
*******************************************************************
* 
      Character*60 mag_file_name,mag_output
      Character*21 pt
      real*8 dt
      Character*230 full_line
      Data dg /57.295779510/
      Data radius /1737.1/
      itotal = 0
C
c     Print *,'Enter input file name: '
c     Read(5,37) mag_file_name 
 37   Format(A60)
C
 240  Open(11,file='$approach',status='old',
     &        access='sequential',form='formatted')
      Print *,'Enter output file name: '
      Open(40,file='$approach.lla',status='new',
     &        access='sequential',form='formatted')
C
  55  Read(11,444,end=99) pt,dt,bx,by,bz,xpos,ypos,zpos,isun
 444  Format(A21,F12.6,F9.3,2(F8.3),34X,F10.2,2F9.2,28X,I3)
      itotal=itotal+1
      If (mod(itotal,1000).eq.0) print *,itotal, ' obs processed'
      alt = sqrt(xpos**2+ypos**2+zpos**2)-radius
      tr = atan2(sqrt(xpos**2+ypos**2),zpos)
      td = tr*dg
      pr = atan2(ypos,xpos)
      pd = pr*dg
      qlat = 90. - td
C
      rbr = sin(tr)*cos(pr)*bx + sin(tr)*sin(pr)*by + cos(tr)*bz
      rbt = cos(tr)*cos(pr)*bx + cos(tr)*sin(pr)*by - sin(tr)*bz
      rbp =        -sin(pr)*bx +         cos(pr)*by
C
      Write(40,445) qlat,pd,alt,rbr,rbt,rbp,isun,dt,pt
 445  Format(6F10.3,I3,F12.6,2x,A21)
      Go to 55
C
  99  Print *,itotal,' observation locations processed'
      End
EOF
g77 -O -o llabrth_moon llabrth_moon.f
./llabrth_moon

done
rm llabrth_moon llabrth_moon.f