24 subroutine bwplotcv(ele,fname,bw,be,bs,bn,jm,b1,b2,maxplots,
25 *olddtm,newdtm,region,elecap,ij,igridsec,fn)
47 character*3 ele,elecap
49 character*10 olddtm,newdtm,region
50 real*8 bw(maxplots),be(maxplots),bs(maxplots),bn(maxplots)
52 real*4 b1(maxplots),b2(maxplots)
53 character*10 fn(maxplots)
69 if(fname(3:3).eq.
't')
then 71 elseif(fname(3:3).eq.
'd')
then 73 elseif(fname(3:3).eq.
'a')
then 75 elseif(fname(3:3).eq.
'r')
then 81 1
format(
'FATAL in bwplotcv. Bad character in spot 3: ',a)
83 if(fname(1:2).ne.
'cv')
then 87 2
format(
'FATAL in bwplotcv. Bad character in spots 1-2: ',a)
89 if(.not.(fname(4:5).eq.
'cd' .or. fname(4:5).eq.
'dd'))
then 93 3
format(
'FATAL in bwplotcv. Bad character in spots 4-5: ',a)
95 if(fname(6:8).ne.ele)
then 96 write(6,4)trim(fname),ele
99 4
format(
'FATAL in bwplotcv. Bad match of fname / ele: ',a,1x,a)
103 if(fname(1:3).eq.
'cvacd')
then 107 if(igridsec.le.0)
then 110 write(gridnote,10)igridsec
112 10
format(
'(',i0,
' sec)')
127 write(99,991)ele,extra,trim(region),trim(fn(ij)),
128 *ele,extra,trim(region),trim(fn(ij))
131 write(99,904)trim(fname),bw(ij),be(ij),bs(ij),bn(ij),
132 *jm(ij),b1(ij),b2(ij),trim(newdtm),trim(olddtm),elecap,
133 *extra,trim(gridnote),trim(region),trim(fn(ij))
142 write(99,910)trim(fname),trim(fn(ij))
150 *
'# -----------------------------------------------------',/,
151 *
'# coverage in ',a,a,1x,a,1x,a,/,
152 *
'# -----------------------------------------------------',/,
153 *
'echo ...coverage in ',a,a,1x,a,1x,a)
158 904
format(
'psxy ',a,
' -Xc -R',f0.6,
'/',f0.6,
'/',sp,f0.6,
'/',f0.6,
159 * ss,
' -JM',f3.1,
'i -B',f0.6,
'/',f0.6,
':."',
160 *
'NADCON v5.0 ',a,
' minus ',a,
' ',a3,a5,a,
163 *
'-Gblack -K > plot.ps')
189 905
format(
'ps2raster plot.ps -Tj -P -A ')
193 910
format(
'mv -f plot.jpg ',a,
'.',a,
'.jpg',/,
197 1910
format(
'mv -f plot.jpg cv',a3,
'.',a,
'.',a,
'.',a,
198 *
'.',a,
'.jpg',/,
'rm -f plot.ps')
subroutine plotcoast(region, ifnum)
Subroutine to write GMT-based commands to create a shoreline Write GMT-based commands to create a sho...
subroutine bwplotcv(ele, fname, bw, be, bs, bn, jm, b1, b2, maxplots, olddtm, newdtm, region, elecap, ij, igridsec, fn)
Subroutine to make GMT calls to do a B/W coverage plot.