#! /usr/bin/python3

import WXP
import sys

#domstr = "ne"
domstr = sys.argv[1]
nf = sys.argv[2]

WXP.Prog.readResrcFile( "sfcmap" )

back = WXP.ImagePlot( "file=us_relief.gif" )

sat = WXP.SatPlot( "if=sat_ir_g16,ct=sat.clr,cof=sat_ir16.cfl,cutoff=-1:1,mod_mask" )
sat.setModelMask( "ruc,init:2m_ag:temp:C,-10" )

rad = WXP.RadPlot( "if=rad_wxpmos,ct=radar1.clr,cof=radmos2a.cfl,bar,sync" )

grd = WXP.GridPlot( "mo=rap,va=init:sl:presm:mb,pl=ln,in=4,ct=wxp.clr,coln=blue:w2" )

frt = WXP.FrontPlot( "shadow=3" )
wat = WXP.WatchPlot( "shadow=3" )
wat.setColor("yellow:wi=2:st=dsh,lred:wi=2:st=dsh" )
warn = WXP.WarnPlot("shadow=3,nolabel" )
#lgt = WXP.LightPlot( "mark=light" )

sfc = WXP.SfcPlot( "va=all,co=lcyan,sc=.65,pr=+1" )

data = WXP.Datum()
dplot = WXP.DatumPlot()

domain = WXP.Domain()
WXP.DomainTool.decode( domain, domstr )
#
#  Set up map
#
map = WXP.MapPlot()
map.setMapFiles( "fi:cnty.mpl,mk:39.24:-94.58:co=white:mk=+::we=2" )
#
#  Open new plot
#
plot = 0
if plot:
   plot = WXP.Plot()
   plot.open( "de=d,ge=1024x768" )
   plot.setDomain( domain )
   back.draw( plot )
   sat.draw( plot )
   rad.draw( plot )
   map.draw( plot )
   sfc.draw( plot )
   grd.draw( plot )
   wat.draw( plot )
   frt.draw( plot )
   data.set( 40.05491, -75.51085, WXP.Datum.MARK, "+", "co=white" )
   dplot.draw( plot, data )
   data.set( 39.91289, -75.57531, WXP.Datum.MARK, "+", "co=white" )
   dplot.draw( plot, data )
   plot.drawBorder()
   plot.drawLabels()
   plot.manage()
else:
   figure = WXP.Figure()
   figure.open( "de=d,ge=1050x850,nf="+nf+",refresh" )
   figure.startPlot( domain )
   figure.addDrawing( back )
   figure.addDrawing( sat )
   figure.addDrawing( rad )
   figure.addDrawing( map )
   figure.addDrawing( sfc )
   figure.addDrawing( grd )
#   figure.addDrawing( lgt )
   figure.addDrawing( wat )
   figure.addDrawing( warn )
   figure.addDrawing( frt )
   figure.endPlot()
   figure.manage()