! jnl file that plots the station location of Argo floats used to make the grid
! requires input arguments for time range, depth range and lat/lon range

! to run the program in ferret, you will need to type the following line with the various values filled in:

!  go fer_temp_station_loc time_start time_end depth_min depth_max lat_min lat_max lon_min lon_max iland 

! If you choose not to give values for the various variables, default ones will be substituted for you.
! If you only want to choose some values, you can, but you cannot skip any values.  This means you can 
! fill in the first 2 and no others, or if you want to specify the lat/lon ranges, you must fill in values
! for the first 4 variables as well as the lat/lon ranges.

! depth values currently correspond to the surface region 0-5db
! time values are entered in corresponding to monthly values starting with Jan, 2004 

! iland  0=hig res, 1=low res go fland
! if no values are entered, the default region is 50S to 25N and 130E to 250E at the surface in Jan 2006
cancel/all data
cancel/all var
cancel/all sym
cancel region
set window/clear
go set_pixel_size 1024 768
set mem/size=80

! assigning variables to arguments passed in the go command as well as default values
let lat_min=`$5%-50%`
let lat_max=`$6%25%`
let lon_min=`$7%130%`
let lon_max=`$8%250%`
let depth_min=`$3%5%`
let depth_max=`$4%5%`
let time_start=`$1%49%`
let time_end=`$2%53%`
let iland=`$9%2%`
let dia=0.2
set mode interpolate
let lonmax=if lon_max gt 360 then 1 else 0 

!let rcheck=`lon_max`-`lon_min`

define axis/x=20.5E:379.5E:1/units=degrees lon
set axis/modulo lon
define axis/y= 89.5S: 89.5N:1/units=degrees lat
!define axis/t=15-jan-2004:15-dec-2009:1/units=months time
define grid/x=lon/y=lat bath


file/order=xy/format="(10f8.3)"/var=count/col=360/skip=0/grid=bath "station_density"
set variable/bad=0 count
shade/set/nokey/nolab/hlim=`lon_min`:`lon_max`/vlim=`lat_min`:`lat_max`/levels="(0,1000,1000)" count+5000
ppl ylab
ppl xlab
ppl shade

define axis/x/x=1:100000:1 longax
define grid/x=longax biggrid

file/var=fl_month/skip=0 "argo_list_map_sum_temp_2004"
file/var=xlon,xlat/grid=biggrid/skip=0 "argo_list_map_floats_temp_2004"
file/var=fl_month/skip=0 "argo_list_map_sum_temp_2005"
file/var=xlon,xlat/grid=biggrid/skip=0 "argo_list_map_floats_temp_2005"
file/var=fl_month/skip=0 "argo_list_map_sum_temp_2006"
file/var=xlon,xlat/grid=biggrid/skip=0 "argo_list_map_floats_temp_2006"
file/var=fl_month/skip=0 "argo_list_map_sum_temp_2007"
file/var=xlon,xlat/grid=biggrid/skip=0 "argo_list_map_floats_temp_2007"
file/var=fl_month/skip=0 "argo_list_map_sum_temp_2008"
file/var=xlon,xlat/grid=biggrid/skip=0 "argo_list_map_floats_temp_2008"
file/var=fl_month/skip=0 "argo_list_map_sum_temp_2009"
file/var=xlon,xlat/grid=biggrid/skip=0 "argo_list_map_floats_temp_2009"
file/var=fl_month/skip=0 "argo_list_map_sum_temp_2010"
file/var=xlon,xlat/grid=biggrid/skip=0 "argo_list_map_floats_temp_2010"
file/var=fl_month/skip=0 "argo_list_map_sum_temp_2011"
file/var=xlon,xlat/grid=biggrid/skip=0 "argo_list_map_floats_temp_2011"
file/var=fl_month/skip=0 "argo_list_map_sum_psal_2012"
file/var=xlon,xlat/grid=biggrid/skip=0 "argo_list_map_floats_psal_2012"

! defining constants to help make plots

let s2004=if `time_start` le 12 then 1 else 0
let e2004=if `time_end` le 12 then 1 else 0
let s2005=if `time_start` le 24 and `time_start` gt 12 then 1 else 0
let e2005=if `time_end` le 24 and `time_end` gt 12 then 1 else 0
let s2006=if `time_start` le 36 and `time_start` gt 24 then 1 else 0
let e2006=if `time_end` le 36 and `time_end` gt 24 then 1 else 0
let s2007=if `time_start` le 48 and `time_start` gt 36 then 1 else 0
let e2007=if `time_end` le 48 and `time_end` gt 36 then 1 else 0
let s2008=if `time_start` le 60 and `time_start` gt 48 then 1 else 0
let e2008=if `time_end` le 60 and `time_end` gt 48 then 1 else 0
let s2009=if `time_start` le 72 and `time_start` gt 60 then 1 else 0
let e2009=if `time_end` le 72 and `time_end` gt 60 then 1 else 0
let s2010=if `time_start` le 84 and `time_start` gt 72 then 1 else 0
let e2010=if `time_end` le 84 and `time_end` gt 72 then 1 else 0
let s2011=if `time_start` le 96 and `time_start` gt 84 then 1 else 0
let e2011=if `time_end` le 96 and `time_end` gt 84 then 1 else 0
let s2012=if `time_start` le 108 and `time_start` gt 96 then 1 else 0
let e2012=if `time_end` le 108 and `time_end` gt 96 then 1 else 0

if `s2004` then 
if `e2004` then
let sumstart=fl_month[d=argo_list_map_sum_temp_2004,i=`time_start`]
let time_end_s=`time_end`+1
let sumend=fl_month[d=argo_list_map_sum_temp_2004,i=`time_end_s`]
let lonx=xlon[d=3,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=3,i=`sumstart`:`sumend`]
say 2004
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia


elif `e2004 eq 0` then

let sumstart=fl_month[d=argo_list_map_sum_temp_2004,i=`time_start`]
let sumend=fl_month[d=argo_list_map_sum_temp_2004,i=13]
let lonx=xlon[d=argo_list_map_floats_temp_2004,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_temp_2004,i=`sumstart`:`sumend`]
say 2004plus
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia

endif  !e2004

endif ! s2004



if `s2005` then
! starts in 2005 and ends in 2005
if `e2005` then
let time_start_s=`time_start`-12
let sumstart=fl_month[d=argo_list_map_sum_temp_2005,i=`time_start_s`]
let time_end_s=`time_end`-11
let sumend=fl_month[d=argo_list_map_sum_temp_2005,i=`time_end_s`]
let lonx=xlon[d=argo_list_map_floats_temp_2005,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_temp_2005,i=`sumstart`:`sumend`]
say 2005
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia


elif `e2004 eq 0` then
! starts in 2005 and ends later

let time_start_s=`time_start`-12
let sumstart=fl_month[d=argo_list_map_sum_temp_2005,i=`time_start_s`]
let sumend=fl_month[d=argo_list_map_sum_temp_2005,i=13]
let lonx=xlon[d=argo_list_map_floats_temp_2005,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_temp_2005,i=`sumstart`:`sumend`]
say 2005plus
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia

endif  !e2005

! starts before 2005 and ends in 2005
elif `e2005` then
let sumstart=fl_month[d=argo_list_map_sum_temp_2005,i=1]
let time_end_s=`time_end`-11
let sumend=fl_month[d=argo_list_map_sum_temp_2005,i=`time_end_s`]
let lonx=xlon[d=argo_list_map_floats_temp_2005,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_temp_2005,i=`sumstart`:`sumend`]
say 2005end
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia

endif  !s2005
 
! starts before 2005 and ends after 2005

let t5= if (`time_start` lt 12 and `time_end` gt 24) then 1 else 0
if `t5` then
let sumstart=fl_month[d=argo_list_map_sum_temp_2005,i=1]
let sumend=fl_month[d=argo_list_map_sum_temp_2005,i=13]
let lonx=xlon[d=argo_list_map_floats_temp_2005,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_temp_2005,i=`sumstart`:`sumend`]
say all2005
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia
endif !t5



if `s2006` then
! starts and ends in 2006
if `e2006` then
let time_start_s=`time_start`-24
let sumstart=fl_month[d=argo_list_map_sum_temp_2006,i=`time_start_s`]
let time_end_s=`time_end`-23
let sumend=fl_month[d=argo_list_map_sum_temp_2006,i=`time_end_s`]
let lonx=xlon[d=argo_list_map_floats_temp_2006,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_temp_2006,i=`sumstart`:`sumend`]
say 2006
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia


elif `e2006 eq 0` then

let time_start_s=`time_start`-24
let sumstart=fl_month[d=argo_list_map_sum_temp_2006,i=`time_start_s`]
let sumend=fl_month[d=argo_list_map_sum_temp_2006,i=13]
let lonx=xlon[d=argo_list_map_floats_temp_2006,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_temp_2006,i=`sumstart`:`sumend`]
say 2006plus
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia

endif ! e2006

! starts before 2006 and ends in 2006
elif `e2006` then

let sumstart=fl_month[d=argo_list_map_sum_temp_2006,i=1]
let time_end_s=`time_end`-23
let sumend=fl_month[d=argo_list_map_sum_temp_2006,i=`time_end_s`]
let lonx=xlon[d=argo_list_map_floats_temp_2006,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_temp_2006,i=`sumstart`:`sumend`]
say 2006end
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia

endif !s2006

! starts before 2006 and ends after 2006

let t6= if (`time_start` lt 25 and `time_end` gt 36) then 1 else 0
if `t6` then
let sumstart=fl_month[d=argo_list_map_sum_temp_2006,i=1]
let sumend=fl_month[d=argo_list_map_sum_temp_2006,i=13]
let lonx=xlon[d=argo_list_map_floats_temp_2006,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_temp_2006,i=`sumstart`:`sumend`]
say all2006
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia
endif !t6


if `s2007` then
! start and end in 2007
if `e2007` then
let time_start_s=`time_start`-36
let sumstart=fl_month[d=argo_list_map_sum_temp_2007,i=`time_start_s`]
let time_end_s=`time_end`-35
let sumend=fl_month[d=argo_list_map_sum_temp_2007,i=`time_end_s`]
let lonx=xlon[d=argo_list_map_floats_temp_2007,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_temp_2007,i=`sumstart`:`sumend`]
say 2007
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia

! starts in 2007 and ends later
elif `e2007 eq 0` then

let time_start_s=`time_start`-36
let sumstart=fl_month[d=argo_list_map_sum_temp_2007,i=`time_start_s`]
let sumend=fl_month[d=argo_list_map_sum_temp_2007,i=13]
let lonx=xlon[d=argo_list_map_floats_temp_2007,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_temp_2007,i=`sumstart`:`sumend`]
say 2007plus
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia

endif ! e2007

! starts before 2007 and ends in 2007
elif `e2007` then

let sumstart=fl_month[d=argo_list_map_sum_temp_2007,i=1]
let time_end_s=`time_end`-35
let sumend=fl_month[d=argo_list_map_sum_temp_2007,i=`time_end_s`]
let lonx=xlon[d=argo_list_map_floats_temp_2007,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_temp_2007,i=`sumstart`:`sumend`]
say 2007end
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia

endif !s2007

! starts before 2007 and ends after 2007

let t7= if (`time_start` lt 37 and `time_end` gt 48) then 1 else 0
if `t7` then
let sumstart=fl_month[d=argo_list_map_sum_temp_2007,i=1]
let sumend=fl_month[d=argo_list_map_sum_temp_2007,i=13]
let lonx=xlon[d=argo_list_map_floats_temp_2007,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_temp_2007,i=`sumstart`:`sumend`]
say all2007
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia
endif !t7



if `s2008` then
! start and end in 2008
if `e2008` then
let time_start_s=`time_start`-48
let sumstart=fl_month[d=argo_list_map_sum_temp_2008,i=`time_start_s`]
let time_end_s=`time_end`-47
let sumend=fl_month[d=argo_list_map_sum_temp_2008,i=`time_end_s`]
let lonx=xlon[d=argo_list_map_floats_temp_2008,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx 
let latx=xlat[d=argo_list_map_floats_temp_2008,i=`sumstart`:`sumend`]
say 2008
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia

! starts in 2008 and ends later
elif `e2008 eq 0` then

let time_start_s=`time_start`-48
let sumstart=fl_month[d=argo_list_map_sum_temp_2008,i=`time_start_s`]
let sumend=fl_month[d=argo_list_map_sum_temp_2008,i=13]
let lonx=xlon[d=argo_list_map_floats_temp_2008,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_temp_2008,i=`sumstart`:`sumend`]
say 2008plus
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia

endif ! e2008

! starts before 2008 and ends in 2008
elif `e2008` then

let sumstart=fl_month[d=argo_list_map_sum_temp_2008,i=1]
let time_end_s=`time_end`-47
let sumend=fl_month[d=argo_list_map_sum_temp_2008,i=`time_end_s`]
let lonx=xlon[d=argo_list_map_floats_temp_2008,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_temp_2008,i=`sumstart`:`sumend`]
say 2008end
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia
endif !s2008

! starts before 2008 and ends after 2008

let t8= if (`time_start` lt 49 and `time_end` gt 60) then 1 else 0
if `t8` then
let sumstart=fl_month[d=argo_list_map_sum_temp_2008,i=1]
let sumend=fl_month[d=argo_list_map_sum_temp_2008,i=13]
let lonx=xlon[d=argo_list_map_floats_temp_2008,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_temp_2008,i=`sumstart`:`sumend`]
say all2008
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia
endif !t8



if `s2009` then
! start and end in 2009
if `e2009` then
let time_start_s=`time_start`-60
let sumstart=fl_month[d=argo_list_map_sum_temp_2009,i=`time_start_s`]
let time_end_s=`time_end`-59
let sumend=fl_month[d=argo_list_map_sum_temp_2009,i=`time_end_s`]
let lonx=xlon[d=argo_list_map_floats_temp_2009,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_temp_2009,i=`sumstart`:`sumend`]
say 2009
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia

! starts in 2009 and ends later
elif `e2009 eq 0` then

let time_start_s=`time_start`-60
let sumstart=fl_month[d=argo_list_map_sum_temp_2009,i=`time_start_s`]
let sumend=fl_month[d=argo_list_map_sum_temp_2009,i=13]
let lonx=xlon[d=argo_list_map_floats_temp_2009,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_temp_2009,i=`sumstart`:`sumend`]
say 2009plus
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia

endif ! e2009

! starts before 2009 and ends in 2009
elif `e2009` then

let sumstart=fl_month[d=argo_list_map_sum_temp_2009,i=1]
let time_end_s=`time_end`-59
let sumend=fl_month[d=argo_list_map_sum_temp_2009,i=`time_end_s`]
let lonx=xlon[d=argo_list_map_floats_temp_2009,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_temp_2009,i=`sumstart`:`sumend`]
say 2009end
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia
endif !s2009

! starts before 2009 and ends after 2009

let t9=if (`time_start lt 61` and `time_end gt 72`) then 1 else 0
if `t9` then 
let sumstart=fl_month[d=argo_list_map_sum_temp_2009,i=1]
let sumend=fl_month[d=argo_list_map_sum_temp_2009,i=13]
let lonx=xlon[d=argo_list_map_floats_temp_2009,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_temp_2009,i=`sumstart`:`sumend`]
say all2009
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia
endif !t9


if `s2010` then
! start and end in 2010
if `e2010` then
let time_start_s=`time_start`-72
let sumstart=fl_month[d=argo_list_map_sum_temp_2010,i=`time_start_s`]
let time_end_s=`time_end`-71
let sumend=fl_month[d=argo_list_map_sum_temp_2010,i=`time_end_s`]
let lonx=xlon[d=argo_list_map_floats_temp_2010,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_temp_2010,i=`sumstart`:`sumend`]
say 2010
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia

! starts in 2010 and ends later
elif `e2010 eq 0` then

let time_start_s=`time_start`-72
let sumstart=fl_month[d=argo_list_map_sum_temp_2010,i=`time_start_s`]
let sumend=fl_month[d=argo_list_map_sum_temp_2010,i=13]
let lonx=xlon[d=argo_list_map_floats_temp_2010,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_temp_2010,i=`sumstart`:`sumend`]
say 2010plus
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia

endif ! e2010

! starts before 2010 and ends in 2010
elif `e2010` then

let sumstart=fl_month[d=argo_list_map_sum_temp_2010,i=1]
let time_end_s=`time_end`-71
let sumend=fl_month[d=argo_list_map_sum_temp_2010,i=`time_end_s`]
let lonx=xlon[d=argo_list_map_floats_temp_2010,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_temp_2010,i=`sumstart`:`sumend`]
say 2010end
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia
endif !s2010

! starts before 2010 and ends after 2010

let t10=if (`time_start lt 73` and `time_end gt 84`) then 1 else 0
if `t10` then 
let sumstart=fl_month[d=argo_list_map_sum_temp_2010,i=1]
let sumend=fl_month[d=argo_list_map_sum_temp_2010,i=13]
let lonx=xlon[d=argo_list_map_floats_temp_2010,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_temp_2010,i=`sumstart`:`sumend`]
say all2010
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia
endif !t10


if `s2011` then
! start and end in 2011
if `e2011` then
let time_start_s=`time_start`-84
let sumstart=fl_month[d=argo_list_map_sum_temp_2011,i=`time_start_s`]
let time_end_s=`time_end`-83
let sumend=fl_month[d=argo_list_map_sum_temp_2011,i=`time_end_s`]
let lonx=xlon[d=argo_list_map_floats_temp_2011,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_temp_2011,i=`sumstart`:`sumend`]
say 2011
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia

! starts in 2011 and ends later
elif `e2011 eq 0` then

let time_start_s=`time_start`-84
let sumstart=fl_month[d=argo_list_map_sum_temp_2011,i=`time_start_s`]
let sumend=fl_month[d=argo_list_map_sum_temp_2011,i=13]
let lonx=xlon[d=argo_list_map_floats_temp_2011,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_temp_2011,i=`sumstart`:`sumend`]
say 2011plus
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia

endif ! e2011

! starts before 2011 and ends in 2011
elif `e2011` then

let sumstart=fl_month[d=argo_list_map_sum_temp_2011,i=1]
let time_end_s=`time_end`-83
let sumend=fl_month[d=argo_list_map_sum_temp_2011,i=`time_end_s`]
let lonx=xlon[d=argo_list_map_floats_temp_2011,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_temp_2011,i=`sumstart`:`sumend`]
say 2011end
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia
endif !s2011

! starts before 2011 and ends after 2011

let t11=if (`time_start lt 85` and `time_end gt 96`) then 1 else 0
if `t11` then 
let sumstart=fl_month[d=argo_list_map_sum_temp_2011,i=1]
let sumend=fl_month[d=argo_list_map_sum_temp_2011,i=13]
let lonx=xlon[d=argo_list_map_floats_temp_2011,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_temp_2011,i=`sumstart`:`sumend`]
say all2011
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia
endif !t11

if `s2012` then
! start and end in 2012
if `e2012` then
let time_start_s=`time_start`-96
let sumstart=fl_month[d=argo_list_map_sum_temp_2012,i=`time_start_s`]
let time_end_s=`time_end`-95
let sumend=fl_month[d=argo_list_map_sum_temp_2012,i=`time_end_s`]
let lonx=xlon[d=argo_list_map_floats_temp_2012,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_temp_2012,i=`sumstart`:`sumend`]
say 2012
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia

! starts in 2012 and ends later
elif `e2012 eq 0` then

let time_start_s=`time_start`-96
let sumstart=fl_month[d=argo_list_map_sum_temp_2012,i=`time_start_s`]
let sumend=fl_month[d=argo_list_map_sum_temp_2012,i=13]
let lonx=xlon[d=argo_list_map_floats_temp_2012,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_temp_2012,i=`sumstart`:`sumend`]
say 2012plus
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia

endif ! e2012

! starts before 2012 and ends in 2012
elif `e2012` then

let sumstart=fl_month[d=argo_list_map_sum_temp_2012,i=1]
let time_end_s=`time_end`-95
let sumend=fl_month[d=argo_list_map_sum_temp_2012,i=`time_end_s`]
let lonx=xlon[d=argo_list_map_floats_temp_2012,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_temp_2012,i=`sumstart`:`sumend`]
say 2012end
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia
endif !s2012

! starts before 2012 and ends after 2012

let t12=if (`time_start lt 97` and `time_end gt 108`) then 1 else 0
if `t12` then 
let sumstart=fl_month[d=argo_list_map_sum_temp_2012,i=1]
let sumend=fl_month[d=argo_list_map_sum_temp_2012,i=13]
let lonx=xlon[d=argo_list_map_floats_temp_2012,i=`sumstart`:`sumend`]
let lonxx=if lonmax eq 1 and lonx lt `lon_min` then lonx+360 else lonx
let latx=xlat[d=argo_list_map_floats_temp_2012,i=`sumstart`:`sumend`]
say all2012
go polymark poly/over/nolab/noaxes/nokey/palette=blue lonxx,latx, , circle, dia
endif !t12

! now use etopo6 to add bathymetry
! check zonal range to see if using high or low land resoluion
let rcheck=`lon_max`-`lon_min`
if `rcheck gt 120` then
if `iland eq 0` then
use etopo12_grid.nc
shade/over/nokey/nolab/pal=black/hlim=`lon_min`:`lon_max`/vlim=`lat_min`:`lat_max`/lev=(0,10000,10000) btdata12[y=`lat_min`:`lat_max`,x=`lon_min`:`lon_max`]
shade/over/nokey/nolab/palette=lightgreyscale/over/hlim=`lon_min`:`lon_max`/vlim=`lat_min`:`lat_max`/levels=(-1000,0,1000) btdata12[y=`lat_min`:`lat_max`,x=`lon_min`:`lon_max`]
shade/over/nokey/nolab/pal=grey/over/hlim=`lon_min`:`lon_max`/vlim=`lat_min`:`lat_max`/lev=(-200,0,200) btdata12[y=`lat_min`:`lat_max`,x=`lon_min`:`lon_max`]
endif
if `iland eq 1` then
use etopo12_grid.nc
fill/over/nokey/nolab/palette=black/hlim=`lon_min`:`lon_max`/vlim=`lat_min`:`lat_max`/levels=(0,10000,10000) btdata12[y=`lat_min`:`lat_max`,x=`lon_min`:`lon_max`]
endif
endif
if `rcheck le 120` then
if `iland eq 0` then
use etopo6_grid.nc
fill/over/nokey/nolab/palette=black/hlim=`lon_min`:`lon_max`/vlim=`lat_min`:`lat_max`/levels=(0,10000,10000) btdata[y=`lat_min`:`lat_max`,x=`lon_min`:`lon_max`]
fill/over/nokey/nolab/palette=lightgreyscale/over/hlim=`lon_min`:`lon_max`/vlim=`lat_min`:`lat_max`/levels=(-1000,0,1000) btdata[y=`lat_min`:`lat_max`,x=`lon_min`:`lon_max`]
fill/over/nokey/nolab/palette=grey/over/hlim=`lon_min`:`lon_max`/vlim=`lat_min`:`lat_max`/levels=(-200,0,200) btdata[y=`lat_min`:`lat_max`,x=`lon_min`:`lon_max`]
else
use etopo6_grid.nc
fill/over/nokey/nolab/palette=black/hlim=`lon_min`:`lon_max`/vlim=`lat_min`:`lat_max`/levels=(0,10000,10000) btdata[y=`lat_min`:`lat_max`,x=`lon_min`:`lon_max`]
endif

endif




say Plot is finished
say To adjust parameters, use the command line to pass arguments
say For this file, the form to pass arguments is
say go fer_station_loc time_start time_end depth_min depth_max lat_min lat_max lon_min lon_max