#!/usr/bin/perl

$cr = $ARGV[0];
$s = $ARGV[1];
$f = $ARGV[2];

$SYNOPDIR = "/synoptic/carrot/M/$cr";
system("mkdir $SYNOPDIR");

$MAGDIR = "/synoptic/carrot/fd_M_96m_01d";
$REMAPDIR = "/synoptic/carrot/fd_M_remap_96m_01d";

$ENV{'mdi'} = 'wd:/synoptic/carrot/{series}/{series}.{#%06d#series}';

$s = sprintf("%d", $s); $ss = $s - 4;
$f = sprintf("%d", $f); $ff = $f + 4;
foreach $d ($ss .. $ff) {
#    if (-d "$REMAPDIR/fd_M_remap_96m_01d.00$d") { next; }
    system("mkdir $REMAPDIR/fd_M_remap_96m_01d.00$d");
    system("rcp -r sonar:/mag/fd_M_96m_01d.00$d $MAGDIR");
    system("uncompress $MAGDIR/fd_M_96m_01d.00$d/*Z");
    system("v2helio z=1 MOFFSET=1 MCORLEV=1 VCORLEV=0 MAPRMAX=0.994 MAPLGMAX=90 MAPLGMIN=-90 MAPBMAX=90 LGSHIFT=2 MAPMMAX=1800 SINBDIVS=540 in=prog:mdi,level:lev1.8,series:fd_M_96m_01d\[$d\] out=prog:mdi,level:lev1.8,series:fd_M_remap_96m_01d\[$d\]");
    system("rm -r $MAGDIR/fd_M_96m_01d.00$d");
}

system("/home/arta/Dev/STAGING/bin/_sgi4/magsynop -f CR=$cr N=7 in=prog:mdi,level:lev1.8,series:fd_M_remap_96m_01d\[$s-$f\] out=$SYNOPDIR MAPMMAX=1800 MAPLGMAX=90 MAPLGMIN=-90 checkqual=1 qualmask=0x402c05f2 center=0 nsig=2.5 rejmin=100 2>1 > $SYNOPDIR/synop_M.$cr.log");
system("mv $SYNOPDIR/synop_M.$cr.fits $SYNOPDIR/synop_M_0.$cr.fits");
system("mv $SYNOPDIR/wt_M.$cr.fits $SYNOPDIR/wt_M_0.$cr.fits");

$ss = $s - 1; $ff = $f - 1;
system("/home/arta/Dev/STAGING/bin/_sgi4/magsynop -f CR=$cr N=7 in=prog:mdi,level:lev1.8,series:fd_M_remap_96m_01d\[$ss-$ff\] out=$SYNOPDIR MAPMMAX=1800 MAPLGMAX=90 MAPLGMIN=-90 checkqual=1 qualmask=0x402c05f2 center=-15 nsig=2.5 rejmin=100 2>1 > $SYNOPDIR/synop_M_15E.$cr.log");
system("mv $SYNOPDIR/synop_M.$cr.fits $SYNOPDIR/synop_M_15E.$cr.fits");
system("mv $SYNOPDIR/wt_M.$cr.fits $SYNOPDIR/wt_M_15E.$cr.fits");

$ss = $s - 2; $ff = $f - 2;
system("/home/arta/Dev/STAGING/bin/_sgi4/magsynop -f CR=$cr N=7 in=prog:mdi,level:lev1.8,series:fd_M_remap_96m_01d\[$ss-$ff\] out=$SYNOPDIR MAPMMAX=1800 MAPLGMAX=90 MAPLGMIN=-90 checkqual=1 qualmask=0x402c05f2 center=-30 nsig=2.5 rejmin=100 2>1 > $SYNOPDIR/synop_M_30E.$cr.log");
system("mv $SYNOPDIR/synop_M.$cr.fits $SYNOPDIR/synop_M_30E.$cr.fits");
system("mv $SYNOPDIR/wt_M.$cr.fits $SYNOPDIR/wt_M_30E.$cr.fits");

$ss = $s - 3; $ff = $f - 3;
system("/home/arta/Dev/STAGING/bin/_sgi4/magsynop -f CR=$cr N=7 in=prog:mdi,level:lev1.8,series:fd_M_remap_96m_01d\[$ss-$ff\] out=$SYNOPDIR MAPMMAX=1800 MAPLGMAX=90 MAPLGMIN=-90 checkqual=1 qualmask=0x402c05f2 center=-45 nsig=2.5 rejmin=100 2>1 > $SYNOPDIR/synop_M_45E.$cr.log");
system("mv $SYNOPDIR/synop_M.$cr.fits $SYNOPDIR/synop_M_45E.$cr.fits");
system("mv $SYNOPDIR/wt_M.$cr.fits $SYNOPDIR/wt_M_45E.$cr.fits");

$ss = $s - 4; $ff = $f - 4;
system("/home/arta/Dev/STAGING/bin/_sgi4/magsynop -f CR=$cr N=7 in=prog:mdi,level:lev1.8,series:fd_M_remap_96m_01d\[$ss-$ff\] out=$SYNOPDIR MAPMMAX=1800 MAPLGMAX=90 MAPLGMIN=-90 checkqual=1 qualmask=0x402c05f2 center=-60 nsig=2.5 rejmin=100 2>1 > $SYNOPDIR/synop_M_60E.$cr.log");
system("mv $SYNOPDIR/synop_M.$cr.fits $SYNOPDIR/synop_M_60E.$cr.fits");
system("mv $SYNOPDIR/wt_M.$cr.fits $SYNOPDIR/wt_M_60E.$cr.fits");

$ss = $s + 1; $ff = $f + 1;
system("/home/arta/Dev/STAGING/bin/_sgi4/magsynop -f CR=$cr N=7 in=prog:mdi,level:lev1.8,series:fd_M_remap_96m_01d\[$ss-$ff\] out=$SYNOPDIR MAPMMAX=1800 MAPLGMAX=90 MAPLGMIN=-90 checkqual=1 qualmask=0x402c05f2 center=15 nsig=2.5 rejmin=100 2>1 > $SYNOPDIR/synop_M_15W.$cr.log");
system("mv $SYNOPDIR/synop_M.$cr.fits $SYNOPDIR/synop_M_15W.$cr.fits");
system("mv $SYNOPDIR/wt_M.$cr.fits $SYNOPDIR/wt_M_15W.$cr.fits");

$ss = $s + 2; $ff = $f + 2;
system("/home/arta/Dev/STAGING/bin/_sgi4/magsynop -f CR=$cr N=7 in=prog:mdi,level:lev1.8,series:fd_M_remap_96m_01d\[$ss-$ff\] out=$SYNOPDIR MAPMMAX=1800 MAPLGMAX=90 MAPLGMIN=-90 checkqual=1 qualmask=0x402c05f2 center=30 nsig=2.5 rejmin=100 2>1 > $SYNOPDIR/synop_M_30W.$cr.log");
system("mv $SYNOPDIR/synop_M.$cr.fits $SYNOPDIR/synop_M_30W.$cr.fits");
system("mv $SYNOPDIR/wt_M.$cr.fits $SYNOPDIR/wt_M_30W.$cr.fits");

$ss = $s + 3; $ff = $f + 3;
system("/home/arta/Dev/STAGING/bin/_sgi4/magsynop -f CR=$cr N=7 in=prog:mdi,level:lev1.8,series:fd_M_remap_96m_01d\[$ss-$ff\] out=$SYNOPDIR MAPMMAX=1800 MAPLGMAX=90 MAPLGMIN=-90 checkqual=1 qualmask=0x402c05f2 center=45 nsig=2.5 rejmin=100 2>1 > $SYNOPDIR/synop_M_45W.$cr.log");
system("mv $SYNOPDIR/synop_M.$cr.fits $SYNOPDIR/synop_M_45W.$cr.fits");
system("mv $SYNOPDIR/wt_M.$cr.fits $SYNOPDIR/wt_M_45W.$cr.fits");

$ss = $s + 4; $ff = $f + 4;
system("/home/arta/Dev/STAGING/bin/_sgi4/magsynop -f CR=$cr N=7 in=prog:mdi,level:lev1.8,series:fd_M_remap_96m_01d\[$ss-$ff\] out=$SYNOPDIR MAPMMAX=1800 MAPLGMAX=90 MAPLGMIN=-90 checkqual=1 qualmask=0x402c05f2 center=60 nsig=2.5 rejmin=100 2>1 > $SYNOPDIR/synop_M_60W.$cr.log");
system("mv $SYNOPDIR/synop_M.$cr.fits $SYNOPDIR/synop_M_60W.$cr.fits");
system("mv $SYNOPDIR/wt_M.$cr.fits $SYNOPDIR/wt_M_60W.$cr.fits");

system("gzip $SYNOPDIR/*.log");
system("rm -r $REMAPDIR/*");