#!/bin/csh
set cr = 2065
set rhofile = rho012.hdf
set rhofileper = rhoper.hdf

set vf="newkirk.hdf"
mkper -z $rhofile $rhofileper
cp /Users/pete/common/templates/newkirk.hdf .
set nlon = 72
@ dlon = 360 / $nlon
set longend = 360

# compute the pb files

#set vf="fs3d09.hdf"
#set option="-rocc 1. -power 3.5 -fac 1. -vf $vf"
set option="-rocc 1.02 -power 0.0 -fac 1. -vf $vf"

@ long = 0
@ num = 0

while ( $long <= $longend )
  @ num = $num + 1
  @ seqn = 1000 + $num
  set seqn=`echo $seqn | sed 's/^1//'`
  set outfile=pb$seqn'.hdf'
  set b0 = `carr2sunpos -cr $cr -long $long`
  echo b0 $b0
  getpb -long $long -b0 $b0 $option $rhofileper $outfile
  echo getpb -long $long -b0 $b0 $option $rhofileper $outfile
stoc -fmt hdf -nx 512 -nz 512 -x0 -3 -x1 3 -z0 -3 -z1 3 $outfile xy_$outfile
movie -scale -min 0 -max 2E-9 -gif -seq $seqn $seqn xy_pb xy_pb_$seqn -palette /Users/pete/common/templates/pals/Grayscale.pal -d 512 512
# oversaturate to look at stalks
#stoc -fmt hdf -nx 512 -nz 512 -x0 -5 -x1 5 -z0 -5 -z1 5 $outfile xy_$outfile
#movie -scale -min 0 -max 5E-10 -gif -seq $seqn $seqn xy_pb xy_pb_$seqn -palette /Users/pete/common/templates/pals/Grayscale.pal -d 512 512

@ long = $long + $dlon

end

#movie -scale -min 0 -max 5E-10 -gif -seq 1 73 xy_pb xy_pb_movie -palette /Users/pete/common/templates/pals/Grayscale.pal -d 512 512
movie -scale -min 0 -max 2E-9 -gif -seq 1 73 xy_pb xy_pb_movie -palette /Users/pete/common/templates/pals/Grayscale.pal -d 512 512

exit 0