67 subroutine coplot(ele,fname,bw,be,bs,bn,jm,b1,b2,maxplots,
68 *olddtm,newdtm,region,elecap,ij,cptlo,cpthi,cptin6,suffixused,
140 implicit real*8(a-h,o-z)
146 character*3 ele,elecap
148 character*10 olddtm,newdtm,region
149 real*8 bw(maxplots),be(maxplots),bs(maxplots),bn(maxplots)
151 real*4 b1(maxplots),b2(maxplots)
152 character*10 fn(maxplots)
154 character*200 suffixused
156 character*20 gridnote
157 character*20 gridnote2
162 if(.not.(fname(1:3).eq.
'vmt'.or.fname(1:3).eq.
'vst' .or.
163 * fname(1:3).eq.
'vmr'.or.fname(1:3).eq.
'vsr' .or.
164 * fname(1:3).eq.
'vme'.or.fname(1:3).eq.
'vse'))
then 165 write(6,2)trim(fname)
168 2
format(
'FATAL in coplot. Bad character in spots 1-3: ',a)
170 if(.not.(fname(4:5).eq.
'cd' .or. fname(4:5).eq.
'dd' .or.
171 * fname(4:5).eq.
'te' ))
then 172 write(6,3)trim(fname)
175 3
format(
'FATAL in coplot. Bad character in spots 4-5: ',a)
177 if(fname(6:8).ne.ele)
then 178 write(6,4)trim(fname),ele
181 4
format(
'FATAL in coplot. Bad match of fname / ele: ',a,1x,a)
183 if(fname(2:2).eq.
'm')
then 185 elseif(fname(2:2).eq.
's')
then 188 write(6,5)trim(fname)
191 5
format(
'FATAL in coplot. Bad units in name: ',a)
193 if(igridsec.le.0)
then 196 write(gridnote,10)igridsec
197 write(gridnote2,11)igridsec
199 6
format(
'FATAL in coplot. Bad igridsec: ',i0)
200 10
format(
'(',i0,
' sec)')
226 write(99,991)ele,trim(units),trim(region),trim(fn(ij)),
228 *ele,trim(units),trim(region),trim(fn(ij)),
236 write(99,901)cptlo,cpthi,cptin
257 *
'gmtset D_FORMAT %.12G',/,
258 *
'makecpt -Crainbow -T',f0.10,
'/',f0.10,
'/',f0.10,
260 *
'gmtset D_FORMAT %.3G')
270 *
'grdgradient ',a,
'.',a,
271 *
'.grd -N4.8 -A90 -M -Gtempi.grd')
278 write(99,903)fname(1:8),trim(suffixused),
279 *bw(ij),be(ij),bs(ij),bn(ij),
280 *jm(ij),b1(ij),b2(ij),
281 *trim(newdtm),trim(olddtm),elecap,igridsec,
282 *trim(region),trim(fn(ij)),fname(2:5)
316 *
'grdimage ',a,
'.',a,
'.grd',
318 *
' -R',f0.6,
'/',f0.6,
'/',f0.6,
'/',f0.6,
319 *
' -JM',f3.1,
'i -B',f0.6,
'/',f0.6,
320 *
':."NADCON v5.0 ',a,
' minus ',a,
' ',a3,
323 *
':" -Ctemp.cpt -K > plot.ps')
335 write(99,904)qqa,trim(units)
338 *
'psscale -D',f3.1,
'i/-0.4i/4.0i/0.2ih',
339 *
' -Ctemp.cpt -I0.4 -B/:',a,
340 *
': -O -K >> plot.ps')
353 ll = len(trim(suffixused))
354 if(suffixused(ll-1:ll).eq.
'd3')
goto 807
355 if(suffixused(ll-1:ll).eq.
'09')
goto 807
356 if(fname(3:5).eq.
'ete')
goto 807
360 write(99,908)fname(1:8),trim(suffixused)
372 *
'grdcontour ',a,
'.',a,
'.grd',
376 *
' -A- -O -K >> plot.ps')
409 *
'ps2raster plot.ps -Tj -P -A')
412 write(99,906)fname(2:8),trim(suffixused),trim(fn(ij))
413 906
format(
'mv plot.jpg c',a,
'.',a,
'.',a,
'.jpg')
431 *
'# -----------------------------------------------------',/,
432 *
'# color plots in ',a,1x,a,1x,a,1x,a,1x,
' Grid:',a /,
433 *
'# -----------------------------------------------------',/,
434 *
'echo ...color plots in ',a,1x,a,1x,a,1x,a,1x,
' Grid:',a)
subroutine coplot(ele, fname, bw, be, bs, bn, jm, b1, b2, maxplots, olddtm, newdtm, region, elecap, ij, cptlo, cpthi, cptin6, suffixused, igridsec, fn)
Subroutine to make GMT calls to do Color Raster Rendering of Gridded Data.
subroutine plotcoast(region, ifnum)
Subroutine to write GMT-based commands to create a shoreline Write GMT-based commands to create a sho...