! 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_psal_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 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

define viewport/xlim=.0,.95/ylim=0,1 top
! setting map scales properly
let beq=if (`lon_min`) eq (`lon_max`) then 1 else 0
if `beq` then
let lon_max=`lon_max`+360
endif
let blm=if (`lon_min`) gt (`lon_max`) then 1 else 0
if `blm` then
let lon_max=`lon_max`+360
endif
define symbol xdeg `lon_max-lon_min`
define symbol ydeg `lat_max-lat_min`
define symbol alasp `($ydeg)/($xdeg)`
set win/asp=($alasp):axis
!set window/asp=0.5/siz=2.

PPL DFLTFNT DR
ppl conpre @P1@DR
ppl axlsze 0.15,.15
ppl axset 1,1,1,1
ppl pen 0,7
ppl pen 1,7
ppl labset .15,.15,.15


ppl pen 0,7
ppl cross,0


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 grid/x=lon/y=lat bath

file/order=xy/format="(10f8.3)"/var=count/col=360/skip=0/grid=bath "c:\marineatlas\ferret\TMAP\fer_dsets\data\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

let a1=if `time_start` eq 1 then 1 else 0
if `a1` then say 'Jan04'
if `a1` then
file/format=(2f8.3,25x,f1.0)/variables="xlat,xlon,date"/skip=0 "c:\marineatlas\ferret\TMAP\fer_dsets\data\station_location\station_psal_0000-0005_012004"
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon,xlat, , circle, dia
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon+360,xlat, , circle, dia
endif
let a2=if `time_start` eq 2 or `time_end` ge 2 and `time_start` le 2 then 1 else 0
if `a2` then say 'Feb04'
if `a2` then
file/format=(2f8.3,25x,f1.0)/variables="xlat,xlon,mode"/skip=0 "c:\marineatlas\ferret\TMAP\fer_dsets\data\station_location\station_psal_0000-0005_022004"
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon,xlat, , circle, dia
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon+360,xlat, , circle, dia
endif
let a3=if `time_start` eq 3 or `time_end` ge 3 and `time_start` le 3 then 1 else 0
if `a3` then say 'Mar04'
if `a3` then
file/format=(2f8.3,25x,f1.0)/variables="xlat,xlon,mode"/skip=0 "c:\marineatlas\ferret\TMAP\fer_dsets\data\station_location\station_psal_0000-0005_032004"
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon,xlat, , circle, dia
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon+360,xlat, , circle, dia
endif
let a4=if `time_start` eq 4 or `time_end` ge 4 and `time_start` le 4 then 1 else 0
if `a4` then say 'Apr04'
if `a4` then
file/format=(2f8.3,25x,f1.0)/variables="xlat,xlon,mode"/skip=0 "c:\marineatlas\ferret\TMAP\fer_dsets\data\station_location\station_psal_0000-0005_042004"
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon,xlat, , circle, dia
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon+360,xlat, , circle, dia
endif
let a5=if `time_start` eq 5 or `time_end` ge 5 and `time_start` le 5 then 1 else 0
if `a5` then say 'May04'
if `a5` then
file/format=(2f8.3,25x,f1.0)/variables="xlat,xlon,mode"/skip=0 "c:\marineatlas\ferret\TMAP\fer_dsets\data\station_location\station_psal_0000-0005_052004"
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon,xlat, , circle, dia
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon+360,xlat, , circle, dia
endif
let a6=if `time_start` eq 6 or `time_end` ge 6 and `time_start` le 6 then 1 else 0
if `a6` then say 'Jun04'
if `a6` then
file/format=(2f8.3,25x,f1.0)/variables="xlat,xlon,mode"/skip=0 "c:\marineatlas\ferret\TMAP\fer_dsets\data\station_location\station_psal_0000-0005_062004"
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon,xlat, , circle, dia
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon+360,xlat, , circle, dia
endif
let a7=if `time_start` eq 7 or `time_end` ge 7 and `time_start` le 7 then 1 else 0
if `a7` then say 'Jul04'
if `a7` then
file/format=(2f8.3,25x,f1.0)/variables="xlat,xlon,mode"/skip=0 "c:\marineatlas\ferret\TMAP\fer_dsets\data\station_location\station_psal_0000-0005_072004"
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon,xlat, , circle, dia
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon+360,xlat, , circle, dia
endif
let a8=if `time_start` eq 8 or `time_end` ge 8 and `time_start` le 8 then 1 else 0
if `a8` then say 'Aug04'
if `a8` then
file/format=(2f8.3,25x,f1.0)/variables="xlat,xlon,mode"/skip=0 "c:\marineatlas\ferret\TMAP\fer_dsets\data\station_location\station_psal_0000-0005_082004"
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon,xlat, , circle, dia
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon+360,xlat, , circle, dia
endif
let a9=if `time_start` eq 9 or `time_end` ge 9 and `time_start` le 9 then 1 else 0
if `a9` then say 'Sep04'
if `a9` then
file/format=(2f8.3,25x,f1.0)/variables="xlat,xlon,mode"/skip=0 "c:\marineatlas\ferret\TMAP\fer_dsets\data\station_location\station_psal_0000-0005_092004"
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon,xlat, , circle, dia
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon+360,xlat, , circle, dia
endif
let a10=if `time_start` eq 10 or `time_end` ge 10 and `time_start` le 10 then 1 else 0
if `a10` then say 'Oct04'
if `a10` then
file/format=(2f8.3,25x,f1.0)/variables="xlat,xlon,mode"/skip=0 "c:\marineatlas\ferret\TMAP\fer_dsets\data\station_location\station_psal_0000-0005_102004"
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon,xlat, , circle, dia
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon+360,xlat, , circle, dia
endif
let a11=if `time_start` eq 11 or `time_end` ge 11 and `time_start` le 11 then 1 else 0
if `a11` then say 'Nov04'
if `a11` then
file/format=(2f8.3,25x,f1.0)/variables="xlat,xlon,mode"/skip=0 "c:\marineatlas\ferret\TMAP\fer_dsets\data\station_location\station_psal_0000-0005_112004"
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon,xlat, , circle, dia
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon+360,xlat, , circle, dia
endif
let a12=if `time_start` eq 12 or `time_end` ge 12 and `time_start` le 12 then 1 else 0
if `a12` then say 'Dec04'
if `a12` then
file/format=(2f8.3,25x,f1.0)/variables="xlat,xlon,mode"/skip=0 "c:\marineatlas\ferret\TMAP\fer_dsets\data\station_location\station_psal_0000-0005_122004"
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon,xlat, , circle, dia
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon+360,xlat, , circle, dia
endif
let a13=if `time_start` eq 13 or `time_end` ge 13 and `time_start` le 13 then 1 else 0
if `a13` then say 'Jan05'
if `a13` then
file/format=(2f8.3,25x,f1.0)/variables="xlat,xlon,mode"/skip=0 "c:\marineatlas\ferret\TMAP\fer_dsets\data\station_location\station_psal_0000-0005_012005"
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon,xlat, , circle, dia
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon+360,xlat, , circle, dia
endif
let a14=if `time_start` eq 14 or `time_end` ge 14 and `time_start` le 14 then 1 else 0
if `a14` then say 'Feb05'
if `a14` then
file/format=(2f8.3,25x,f1.0)/variables="xlat,xlon,mode"/skip=0 "c:\marineatlas\ferret\TMAP\fer_dsets\data\station_location\station_psal_0000-0005_022005"
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon,xlat, , circle, dia
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon+360,xlat, , circle, dia
endif
let a15=if `time_start` eq 15 or `time_end` ge 15 and `time_start` le 15 then 1 else 0
if `a15` then say 'Mar05'
if `a15` then
file/format=(2f8.3,25x,f1.0)/variables="xlat,xlon,mode"/skip=0 "c:\marineatlas\ferret\TMAP\fer_dsets\data\station_location\station_psal_0000-0005_032005"
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon,xlat, , circle, dia
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon+360,xlat, , circle, dia
endif
let a16=if `time_start` eq 16 or `time_end` ge 16 and `time_start` le 16 then 1 else 0
if `a16` then say 'Apr05'
if `a16` then
file/format=(2f8.3,25x,f1.0)/variables="xlat,xlon,mode"/skip=0 "c:\marineatlas\ferret\TMAP\fer_dsets\data\station_location\station_psal_0000-0005_042005"
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon,xlat, , circle, dia
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon+360,xlat, , circle, dia
endif
let a17=if `time_start` eq 17 or `time_end` ge 17 and `time_start` le 17 then 1 else 0
if `a17` then say 'May05'
if `a17` then
file/format=(2f8.3,25x,f1.0)/variables="xlat,xlon,mode"/skip=0 "c:\marineatlas\ferret\TMAP\fer_dsets\data\station_location\station_psal_0000-0005_052005"
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon,xlat, , circle, dia
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon+360,xlat, , circle, dia
endif
let a18=if `time_start` eq 18 or `time_end` ge 18 and `time_start` le 18 then 1 else 0
if `a18` then say 'Jun05'
if `a18` then
file/format=(2f8.3,25x,f1.0)/variables="xlat,xlon,mode"/skip=0 "c:\marineatlas\ferret\TMAP\fer_dsets\data\station_location\station_psal_0000-0005_062005"
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon,xlat, , circle, dia
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon+360,xlat, , circle, dia
endif
let a19=if `time_start` eq 19 or `time_end` ge 19 and `time_start` le 19 then 1 else 0
if `a19` then say 'Jul05'
if `a19` then
file/format=(2f8.3,25x,f1.0)/variables="xlat,xlon,mode"/skip=0 "c:\marineatlas\ferret\TMAP\fer_dsets\data\station_location\station_psal_0000-0005_072005"
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon,xlat, , circle, dia
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon+360,xlat, , circle, dia
endif
let a20=if `time_start` eq 20 or `time_end` ge 20 and `time_start` le 20 then 1 else 0
if `a20` then say 'Aug05'
if `a20` then
file/format=(2f8.3,25x,f1.0)/variables="xlat,xlon,mode"/skip=0 "c:\marineatlas\ferret\TMAP\fer_dsets\data\station_location\station_psal_0000-0005_082005"
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon,xlat, , circle, dia
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon+360,xlat, , circle, dia
endif
let a21=if `time_start` eq 21 or `time_end` ge 21 and `time_start` le 21 then 1 else 0
if `a21` then say 'Sep05'
if `a21` then
file/format=(2f8.3,25x,f1.0)/variables="xlat,xlon,mode"/skip=0 "c:\marineatlas\ferret\TMAP\fer_dsets\data\station_location\station_psal_0000-0005_092005"
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon,xlat, , circle, dia
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon+360,xlat, , circle, dia
endif
let a22=if `time_start` eq 22 or `time_end` ge 22 and `time_start` le 22 then 1 else 0
if `a22` then say 'Oct05'
if `a22` then
file/format=(2f8.3,25x,f1.0)/variables="xlat,xlon,mode"/skip=0 "c:\marineatlas\ferret\TMAP\fer_dsets\data\station_location\station_psal_0000-0005_102005"
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon,xlat, , circle, dia
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon+360,xlat, , circle, dia
endif
let a23=if `time_start` eq 23 or `time_end` ge 23 and `time_start` le 23 then 1 else 0
if `a23` then say 'Nov05'
if `a23` then
file/format=(2f8.3,25x,f1.0)/variables="xlat,xlon,mode"/skip=0 "c:\marineatlas\ferret\TMAP\fer_dsets\data\station_location\station_psal_0000-0005_112005"
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon,xlat, , circle, dia
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon+360,xlat, , circle, dia
endif
let a24=if `time_start` eq 24 or `time_end` ge 24 and `time_start` le 24 then 1 else 0
if `a24` then say 'Dec05'
if `a24` then
file/format=(2f8.3,25x,f1.0)/variables="xlat,xlon,mode"/skip=0 "c:\marineatlas\ferret\TMAP\fer_dsets\data\station_location\station_psal_0000-0005_122005"
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon,xlat, , circle, dia
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon+360,xlat, , circle, dia
endif
let a25=if `time_start` eq 25 or `time_end` ge 25 and `time_start` le 25 then 1 else 0
if `a25` then say 'Jan06'
if `a25` then
file/format=(2f8.3,25x,f1.0)/variables="xlat,xlon,mode"/skip=0 "c:\marineatlas\ferret\TMAP\fer_dsets\data\station_location\station_psal_0000-0005_012006"
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon,xlat, , circle, dia
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon+360,xlat, , circle, dia
endif
let a26=if `time_start` eq 26 or `time_end` ge 26 and `time_start` le 26 then 1 else 0
if `a26` then say 'Feb06'
if `a26` then
file/format=(2f8.3,25x,f1.0)/variables="xlat,xlon,mode"/skip=0 "c:\marineatlas\ferret\TMAP\fer_dsets\data\station_location\station_psal_0000-0005_022006"
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon,xlat, , circle, dia
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon+360,xlat, , circle, dia
endif
let a27=if `time_start` eq 27 or `time_end` ge 27 and `time_start` le 27 then 1 else 0
if `a27` then say 'Mar06'
if `a27` then
file/format=(2f8.3,25x,f1.0)/variables="xlat,xlon,mode"/skip=0 "c:\marineatlas\ferret\TMAP\fer_dsets\data\station_location\station_psal_0000-0005_032006"
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon,xlat, , circle, dia
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon+360,xlat, , circle, dia
endif
let a28=if `time_start` eq 28 or `time_end` ge 28 and `time_start` le 28 then 1 else 0
if `a28` then say 'Apr06'
if `a28` then
file/format=(2f8.3,25x,f1.0)/variables="xlat,xlon,mode"/skip=0 "c:\marineatlas\ferret\TMAP\fer_dsets\data\station_location\station_psal_0000-0005_042006"
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon,xlat, , circle, dia
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon+360,xlat, , circle, dia
endif
let a29=if `time_start` eq 29 or `time_end` ge 29 and `time_start` le 29 then 1 else 0
if `a29` then say 'May06'
if `a29` then
file/format=(2f8.3,25x,f1.0)/variables="xlat,xlon,mode"/skip=0 "c:\marineatlas\ferret\TMAP\fer_dsets\data\station_location\station_psal_0000-0005_052006"
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon,xlat, , circle, dia
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon+360,xlat, , circle, dia
endif
let a30=if `time_start` eq 30 or `time_end` ge 30 and `time_start` le 30 then 1 else 0
if `a30` then say 'Jun06'
if `a30` then
file/format=(2f8.3,25x,f1.0)/variables="xlat,xlon,mode"/skip=0 "c:\marineatlas\ferret\TMAP\fer_dsets\data\station_location\station_psal_0000-0005_062006"
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon,xlat, , circle, dia
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon+360,xlat, , circle, dia
endif
let a31=if `time_start` eq 31 or `time_end` ge 31 and `time_start` le 31 then 1 else 0
if `a31` then say 'Jul06'
if `a31` then
file/format=(2f8.3,25x,f1.0)/variables="xlat,xlon,mode"/skip=0 "c:\marineatlas\ferret\TMAP\fer_dsets\data\station_location\station_psal_0000-0005_072006"
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon,xlat, , circle, dia
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon+360,xlat, , circle, dia
endif
let a32=if `time_start` eq 32 or `time_end` ge 32 and `time_start` le 32 then 1 else 0
if `a32` then say 'Aug06'
if `a32` then
file/format=(2f8.3,25x,f1.0)/variables="xlat,xlon,mode"/skip=0 "c:\marineatlas\ferret\TMAP\fer_dsets\data\station_location\station_psal_0000-0005_082006"
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon,xlat, , circle, dia
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon+360,xlat, , circle, dia
endif
let a33=if `time_start` eq 33 or `time_end` ge 33 and `time_start` le 33 then 1 else 0
if `a33` then say 'Sep06'
if `a33` then
file/format=(2f8.3,25x,f1.0)/variables="xlat,xlon,mode"/skip=0 "c:\marineatlas\ferret\TMAP\fer_dsets\data\station_location\station_psal_0000-0005_092006"
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon,xlat, , circle, dia
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon+360,xlat, , circle, dia
endif
let a34=if `time_start` eq 34 or `time_end` ge 34 and `time_start` le 34 then 1 else 0
if `a34` then say 'Oct06'
if `a34` then
file/format=(2f8.3,25x,f1.0)/variables="xlat,xlon,mode"/skip=0 "c:\marineatlas\ferret\TMAP\fer_dsets\data\station_location\station_psal_0000-0005_102006"
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon,xlat, , circle, dia
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon+360,xlat, , circle, dia
endif
let a35=if `time_start` eq 35 or `time_end` ge 35 and `time_start` le 35 then 1 else 0
if `a35` then say 'Nov06'
if `a35` then
file/format=(2f8.3,25x,f1.0)/variables="xlat,xlon,mode"/skip=0 "c:\marineatlas\ferret\TMAP\fer_dsets\data\station_location\station_psal_0000-0005_112006"
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon,xlat, , circle, dia
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon+360,xlat, , circle, dia
endif
let a36=if `time_start` eq 36 or `time_end` ge 36 and `time_start` le 36 then 1 else 0
if `a36` then say 'Dec06'
if `a36` then
file/format=(2f8.3,25x,f1.0)/variables="xlat,xlon,mode"/skip=0 "c:\marineatlas\ferret\TMAP\fer_dsets\data\station_location\station_psal_0000-0005_122006"
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon,xlat, , circle, dia
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon+360,xlat, , circle, dia
endif
let a37=if `time_start` eq 37 or `time_end` ge 37 and `time_start` le 37 then 1 else 0
if `a37` then say 'Jan07'
if `a37` then
file/format=(2f8.3,25x,f1.0)/variables="xlat,xlon,mode"/skip=0 "c:\marineatlas\ferret\TMAP\fer_dsets\data\station_location\station_psal_0000-0005_012007"
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon,xlat, , circle, dia
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon+360,xlat, , circle, dia
endif
let a38=if `time_start` eq 38 or `time_end` ge 38 and `time_start` le 38 then 1 else 0
if `a38` then say 'Feb07'
if `a38` then
file/format=(2f8.3,25x,f1.0)/variables="xlat,xlon,mode"/skip=0 "c:\marineatlas\ferret\TMAP\fer_dsets\data\station_location\station_psal_0000-0005_022007"
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon,xlat, , circle, dia
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon+360,xlat, , circle, dia
endif
let a39=if `time_start` eq 39 or `time_end` ge 39 and `time_start` le 39 then 1 else 0
if `a39` then say 'Mar07'
if `a39` then
file/format=(2f8.3,25x,f1.0)/variables="xlat,xlon,mode"/skip=0 "c:\marineatlas\ferret\TMAP\fer_dsets\data\station_location\station_psal_0000-0005_032007"
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon,xlat, , circle, dia
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon+360,xlat, , circle, dia
endif
let a40=if `time_start` eq 40 or `time_end` ge 40 and `time_start` le 40 then 1 else 0
if `a40` then say 'Apr07'
if `a40` then
file/format=(2f8.3,25x,f1.0)/variables="xlat,xlon,mode"/skip=0 "c:\marineatlas\ferret\TMAP\fer_dsets\data\station_location\station_psal_0000-0005_042007"
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon,xlat, , circle, dia
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon+360,xlat, , circle, dia
endif
let a41=if `time_start` eq 41 or `time_end` ge 41 and `time_start` le 41 then 1 else 0
if `a41` then say 'May07'
if `a41` then
file/format=(2f8.3,25x,f1.0)/variables="xlat,xlon,mode"/skip=0 "c:\marineatlas\ferret\TMAP\fer_dsets\data\station_location\station_psal_0000-0005_052007"
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon,xlat, , circle, dia
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon+360,xlat, , circle, dia
endif
let a42=if `time_start` eq 42 or `time_end` ge 42 and `time_start` le 42 then 1 else 0
if `a42` then say 'Jun07'
if `a42` then
file/format=(2f8.3,25x,f1.0)/variables="xlat,xlon,mode"/skip=0 "c:\marineatlas\ferret\TMAP\fer_dsets\data\station_location\station_psal_0000-0005_062007"
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon,xlat, , circle, dia
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon+360,xlat, , circle, dia
endif
let a43=if `time_start` eq 43 or `time_end` ge 43 and `time_start` le 43 then 1 else 0
if `a43` then say 'Jul07'
if `a43` then
file/format=(2f8.3,25x,f1.0)/variables="xlat,xlon,mode"/skip=0 "c:\marineatlas\ferret\TMAP\fer_dsets\data\station_location\station_psal_0000-0005_072007"
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon,xlat, , circle, dia
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon+360,xlat, , circle, dia
endif
let a44=if `time_start` eq 44 or `time_end` ge 44 and `time_start` le 44 then 1 else 0
if `a44` then say 'Aug07'
if `a44` then
file/format=(2f8.3,25x,f1.0)/variables="xlat,xlon,mode"/skip=0 "c:\marineatlas\ferret\TMAP\fer_dsets\data\station_location\station_psal_0000-0005_082007"
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon,xlat, , circle, dia
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon+360,xlat, , circle, dia
endif
let a45=if `time_start` eq 45 or `time_end` ge 45 and `time_start` le 45 then 1 else 0
if `a45` then say 'Sep07'
if `a45` then
file/format=(2f8.3,25x,f1.0)/variables="xlat,xlon,mode"/skip=0 "c:\marineatlas\ferret\TMAP\fer_dsets\data\station_location\station_psal_0000-0005_092007"
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon,xlat, , circle, dia
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon+360,xlat, , circle, dia
endif
let a46=if `time_start` eq 46 or `time_end` ge 46 and `time_start` le 46 then 1 else 0
if `a46` then say 'Oct07'
if `a46` then
file/format=(2f8.3,25x,f1.0)/variables="xlat,xlon,mode"/skip=0 "c:\marineatlas\ferret\TMAP\fer_dsets\data\station_location\station_psal_0000-0005_102007"
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon,xlat, , circle, dia
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon+360,xlat, , circle, dia
endif
let a47=if `time_start` eq 47 or `time_end` ge 47 and `time_start` le 47 then 1 else 0
if `a47` then say 'Nov07'
if `a47` then
file/format=(2f8.3,25x,f1.0)/variables="xlat,xlon,mode"/skip=0 "c:\marineatlas\ferret\TMAP\fer_dsets\data\station_location\station_psal_0000-0005_112007"
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon,xlat, , circle, dia
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon+360,xlat, , circle, dia
endif
let a48=if `time_start` eq 48 or `time_end` ge 48 and `time_start` le 48 then 1 else 0
if `a48` then say 'Dec07'
if `a48` then
file/format=(2f8.3,25x,f1.0)/variables="xlat,xlon,mode"/skip=0 "c:\marineatlas\ferret\TMAP\fer_dsets\data\station_location\station_psal_0000-0005_122007"
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon,xlat, , circle, dia
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon+360,xlat, , circle, dia
endif
let a49=if `time_start` eq 49 or `time_end` ge 49 and `time_start` le 49 then 1 else 0
if `a49` then say 'Jan08'
if `a49` then
file/format=(2f8.3,25x,f1.0)/variables="xlat,xlon,mode"/skip=0 "c:\marineatlas\ferret\TMAP\fer_dsets\data\station_location\station_psal_0000-0005_012008"
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon,xlat, , circle, dia
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon+360,xlat, , circle, dia
endif
let a50=if `time_start` eq 50 or `time_end` ge 50 and `time_start` le 50 then 1 else 0
if `a50` then say 'Feb08'
if `a50` then
file/format=(2f8.3,25x,f1.0)/variables="xlat,xlon,mode"/skip=0 "c:\marineatlas\ferret\TMAP\fer_dsets\data\station_location\station_psal_0000-0005_022008"
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon,xlat, , circle, dia
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon+360,xlat, , circle, dia
endif
let a51=if `time_start` eq 51 or `time_end` ge 51 and `time_start` le 51 then 1 else 0
if `a51` then say 'Mar08'
if `a51` then
file/format=(2f8.3,25x,f1.0)/variables="xlat,xlon,mode"/skip=0 "c:\marineatlas\ferret\TMAP\fer_dsets\data\station_location\station_psal_0000-0005_032008"
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon,xlat, , circle, dia
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon+360,xlat, , circle, dia
endif
let a52=if `time_start` eq 52 or `time_end` ge 52 and `time_start` le 52 then 1 else 0
if `a52` then say 'Apr08'
if `a52` then
file/format=(2f8.3,25x,f1.0)/variables="xlat,xlon,mode"/skip=0 "c:\marineatlas\ferret\TMAP\fer_dsets\data\station_location\station_psal_0000-0005_042008"
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon,xlat, , circle, dia
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon+360,xlat, , circle, dia
endif
let a53=if `time_start` eq 53 or `time_end` ge 53 and `time_start` le 53 then 1 else 0
if `a53` then say 'May08'
if `a53` then
file/format=(2f8.3,25x,f1.0)/variables="xlat,xlon,mode"/skip=0 "c:\marineatlas\ferret\TMAP\fer_dsets\data\station_location\station_psal_0000-0005_052008"
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon,xlat, , circle, dia
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon+360,xlat, , circle, dia
endif
let a54=if `time_start` eq 54 or `time_end` ge 54 and `time_start` le 54 then 1 else 0
if `a54` then say 'June08'
if `a54` then
file/format=(2f8.3,25x,f1.0)/variables="xlat,xlon,mode"/skip=0 "c:\marineatlas\ferret\TMAP\fer_dsets\data\station_location\station_psal_0000-0005_062008"
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon,xlat, , circle, dia
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon+360,xlat, , circle, dia
endif
let a55=if `time_start` eq 55 or `time_end` ge 55 and `time_start` le 55 then 1 else 0
if `a55` then say 'July08'
if `a55` then
file/format=(2f8.3,25x,f1.0)/variables="xlat,xlon,mode"/skip=0 "c:\marineatlas\ferret\TMAP\fer_dsets\data\station_location\station_psal_0000-0005_072008"
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon,xlat, , circle, dia
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon+360,xlat, , circle, dia
endif
let a56=if `time_start` eq 56 or `time_end` ge 56 and `time_start` le 56 then 1 else 0
if `a56` then say 'Aug08'
if `a56` then
file/format=(2f8.3,25x,f1.0)/variables="xlat,xlon,mode"/skip=0 "c:\marineatlas\ferret\TMAP\fer_dsets\data\station_location\station_psal_0000-0005_082008"
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon,xlat, , circle, dia
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon+360,xlat, , circle, dia
endif
let a57=if `time_start` eq 57 or `time_end` ge 57 and `time_start` le 57 then 1 else 0
if `a57` then say 'Sep08'
if `a57` then
file/format=(2f8.3,25x,f1.0)/variables="xlat,xlon,mode"/skip=0 "c:\marineatlas\ferret\TMAP\fer_dsets\data\station_location\station_psal_0000-0005_092008"
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon,xlat, , circle, dia
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon+360,xlat, , circle, dia
endif
let a58=if `time_start` eq 58 or `time_end` ge 58 and `time_start` le 58 then 1 else 0
if `a58` then say 'Oct08'
if `a58` then
file/format=(2f8.3,25x,f1.0)/variables="xlat,xlon,mode"/skip=0 "c:\marineatlas\ferret\TMAP\fer_dsets\data\station_location\station_psal_0000-0005_102008"
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon,xlat, , circle, dia
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon+360,xlat, , circle, dia
endif
let a59=if `time_start` eq 59 or `time_end` ge 59 and `time_start` le 59 then 1 else 0
if `a59` then say 'Nov08'
if `a59` then
file/format=(2f8.3,25x,f1.0)/variables="xlat,xlon,mode"/skip=0 "c:\marineatlas\ferret\TMAP\fer_dsets\data\station_location\station_psal_0000-0005_112008"
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon,xlat, , circle, dia
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon+360,xlat, , circle, dia
endif
let a60=if `time_start` eq 60 or `time_end` ge 60 and `time_start` le 60 then 1 else 0
if `a60` then say 'Dec08'
if `a60` then
file/format=(2f8.3,25x,f1.0)/variables="xlat,xlon,mode"/skip=0 "c:\marineatlas\ferret\TMAP\fer_dsets\data\station_location\station_psal_0000-0005_122008"
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon,xlat, , circle, dia
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon+360,xlat, , circle, dia
endif
let a61=if `time_start` eq 61 or `time_end` ge 61 and `time_start` le 61 then 1 else 0
if `a61` then say 'Jan09'
if `a61` then
file/format=(2f8.3,25x,f1.0)/variables="xlat,xlon,mode"/skip=0 "c:\marineatlas\ferret\TMAP\fer_dsets\data\station_location\station_psal_0000-0005_012009"
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon,xlat, , circle, dia
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon+360,xlat, , circle, dia
endif
let a62=if `time_start` eq 62 or `time_end` ge 62 and `time_start` le 62 then 1 else 0
if `a62` then say 'Feb09'
if `a62` then
file/format=(2f8.3,25x,f1.0)/variables="xlat,xlon,mode"/skip=0 "c:\marineatlas\ferret\TMAP\fer_dsets\data\station_location\station_psal_0000-0005_022009"
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon,xlat, , circle, dia
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon+360,xlat, , circle, dia
endif
let a63=if `time_start` eq 63 or `time_end` ge 63 and `time_start` le 63 then 1 else 0
if `a63` then say 'Mar09'
if `a63` then
file/format=(2f8.3,25x,f1.0)/variables="xlat,xlon,mode"/skip=0 "c:\marineatlas\ferret\TMAP\fer_dsets\data\station_location\station_psal_0000-0005_032009"
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon,xlat, , circle, dia
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon+360,xlat, , circle, dia
endif
let a64=if `time_start` eq 64 or `time_end` ge 64 and `time_start` le 64 then 1 else 0
if `a64` then say 'Apr09'
if `a64` then
file/format=(2f8.3,25x,f1.0)/variables="xlat,xlon,mode"/skip=0 "c:\marineatlas\ferret\TMAP\fer_dsets\data\station_location\station_psal_0000-0005_042009"
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon,xlat, , circle, dia
go polymark poly/over/nolab/noaxes/nokey/palette=blue xlon+360,xlat, , circle, dia
endif


! 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_psal_station_loc time_start time_end depth_min depth_max lat_min lat_max lon_min lon_max