65 subroutine coplotwcv(ele,fname,bw,be,bs,bn,jm,b1,b2,maxplots,
66 *olddtm,newdtm,region,elecap,ij,cptlo,cpthi,cptin6,suffixused,
134 implicit real*8(a-h,o-z)
140 character*3 ele,elecap
141 character*200 fname,cvfname
142 character*10 olddtm,newdtm,region
143 real*8 bw(maxplots),be(maxplots),bs(maxplots),bn(maxplots)
145 real*4 b1(maxplots),b2(maxplots)
146 character*10 fn(maxplots)
148 character*200 suffixused
150 character*20 gridnote
151 character*20 gridnote2
156 if(.not.(fname(1:3).eq.
'vmt'.or.fname(1:3).eq.
'vst' .or.
157 * fname(1:3).eq.
'vmr'.or.fname(1:3).eq.
'vsr'))
then 158 write(6,2)trim(fname)
161 2
format(
'FATAL in coplotwcv. Bad character in spots 1-3: ',a)
163 if(.not.(fname(4:5).eq.
'cd' .or. fname(4:5).eq.
'dd'))
then 164 write(6,3)trim(fname)
167 3
format(
'FATAL in coplotwcv. Bad character in spots 4-5: ',a)
169 if(fname(6:8).ne.ele)
then 170 write(6,4)trim(fname),ele
173 4
format(
'FATAL in coplotwcv. Bad match of fname / ele: ',a,1x,a)
175 if(fname(2:2).eq.
'm')
then 177 elseif(fname(2:2).eq.
's')
then 180 write(6,5)trim(fname)
183 5
format(
'FATAL in coplotwcv. Bad units in name: ',a)
185 if(igridsec.le.0)
then 188 write(gridnote,10)igridsec
189 write(gridnote2,11)igridsec
191 6
format(
'FATAL in coplotwcv. Bad igridsec: ',i0)
192 10
format(
'(',i0,
' sec)')
220 write(99,991)ele,trim(units),trim(region),trim(fn(ij)),
222 *ele,trim(units),trim(region),trim(fn(ij)),
231 write(99,901)cptlo,cpthi,cptin
252 *
'gmtset D_FORMAT %.12G',/,
253 *
'makecpt -Crainbow -T',f0.10,
'/',f0.10,
'/',f0.10,
255 *
'gmtset D_FORMAT %.3G')
266 *
'grdgradient ',a,
'.',a,
267 *
'.grd -N4.8 -A90 -M -Gtempi.grd')
274 write(99,903)fname(1:8),trim(suffixused),
275 *bw(ij),be(ij),bs(ij),bn(ij),
276 *jm(ij),b1(ij),b2(ij),
277 *trim(newdtm),trim(olddtm),elecap,igridsec,
278 *trim(region),trim(fn(ij)),
'method noise' 313 *
'grdimage ',a,
'.',a,
'.grd',
315 *
' -R',f0.6,
'/',f0.6,
'/',f0.6,
'/',f0.6,
316 *
' -JM',f3.1,
'i -B',f0.6,
'/',f0.6,
317 *
':."NADCON v5.0 ',a,
' minus ',a,
' ',a3,
320 *
':" -Ctemp.cpt -K > plot.ps')
331 write(99,904)qqa,trim(units)
334 *
'psscale -D',f3.1,
'i/-0.4i/4.0i/0.2ih',
335 *
' -Ctemp.cpt -I0.4 -B/:',a,
336 *
': -O -K >> plot.ps')
346 ll = len(trim(suffixused))
347 if(suffixused(ll-1:ll).eq.
'd3')
goto 807
348 if(suffixused(ll-1:ll).eq.
'09')
goto 807
349 if(fname(3:5).eq.
'ete')
goto 807
352 write(99,908)fname(1:8),trim(suffixused)
365 *
'grdcontour ',a,
'.',a,
'.grd',
369 *
' -A- -O -K >> plot.ps')
396 if(suffixused(ll-1:ll).eq.
'd3' .and. ij.eq.1)
goto 808
397 if(suffixused(ll-1:ll).eq.
'09' .and. ij.eq.1)
goto 808
398 if(fname(3:5).eq.
'ete' .and. ij.eq.1)
goto 808
401 write(99,909)trim(cvfname)
404 *
'psxy ',a,
' -R -J -Sc0.02i -Gblack -O -K >> plot.ps')
416 *
'ps2raster plot.ps -Tj -P -A')
419 write(99,906)fname(2:8),trim(suffixused),trim(fn(ij))
420 906
format(
'mv plot.jpg c',a,
'.',a,
'.',a,
'.jpg')
437 *
'# -----------------------------------------------------',/,
438 *
'# color plots with coverage in ',a,1x,a,1x,a,1x,a,
440 *
'# -----------------------------------------------------',/,
441 *
'echo ...color plots with coverage in ',a,1x,a,1x,a,1x,a,
subroutine coplotwcv(ele, fname, bw, be, bs, bn, jm, b1, b2, maxplots, olddtm, newdtm, region, elecap, ij, cptlo, cpthi, cptin6, suffixused, igridsec, fn, cvfname)
Subroutine to make GMT calls to do a color raster rendering of gridded data, with coverage overlaid...
subroutine plotcoast(region, ifnum)
Subroutine to write GMT-based commands to create a shoreline Write GMT-based commands to create a sho...