grads阴影图等值线风矢量图叠加颜色报错

nuistwangqin

* 地转风、涡度和散度综合分析 - 四张图
* 固定色标范围：-2.0到2.0，间隔0.5，使用天气图常规颜色

'reinit'
'open D:\sx\vorticity_scheme1.ctl'
'set grid off'
'set lon 60 150'
'set lat 20 60'
'set t 5'
'set z 3'
'set gxout shaded'
* 设置固定色标范围和天气图常规颜色
'set clevs -2.0 -1.5 -1.0 -0.5 0.0 0.5 1.0 1.5 2.0'
*'set ccols 15 14 13 12 11 10 7 6 5 4'
'd vort * 1e4'
'run E:\OpenGrADS-2.2\Contents\Resources\Scripts\cbar.gs'
'open D:\sx\geowind.ctl'
'set dfile 2'
'set lon 60 150'
'set lat 20 60'
'set t 5'
'set z 3'
'set gxout vector'
'set arrscl 0.8  50'
'set ccolor 0'
'd ug; vg'
'open D:\sx\divergence.ctl'
'set dfile 3'
'set lon 60 150'
'set lat 20 60'
'set t 5'
'set z 3'
'set gxout contour'
'set ccolor 1'     *这里换其他颜色会报错
'd div * 1e4'
'draw title 25th day 00Z 850hPa Geostrophic Wind & Vorticity & Divergence'
'printim 850hPa_25day00Z_geo_vort_div.png white'

'c'

'reinit'
'open D:\sx\vorticity_scheme1.ctl'
'set grid off'
'set lon 60 150'
'set lat 20 60'
'set t 8'
'set z 3'
'set gxout shaded'
* 设置固定色标范围和天气图常规颜色
'set clevs -2.0 -1.5 -1.0 -0.5 0.0 0.5 1.0 1.5 2.0'
*'set ccols 4 5 6 7 10 11 12 13 14 15'
'd vort * 1e4'
'run E:\OpenGrADS-2.2\Contents\Resources\Scripts\cbar.gs'
'open D:\sx\geowind.ctl'
'set dfile 2'
'set lon 60 150'
'set lat 20 60'
'set t 8'
'set z 3'
'set gxout vector'
'set arrscl 0.8 50'
'set ccolor 0'
'd ug; vg'
'open D:\sx\divergence.ctl'
'set dfile 3'
'set lon 60 150'
'set lat 20 60'
'set t 8'
'set z 3'
'set gxout contour'
'set ccolor 1'   *这里换其他颜色会报错，其他地方类似
'd div * 1e4'
'draw title 26th day 12Z 850hPa Geostrophic Wind & Vorticity & Divergence'
'printim 850hPa_26day12Z_geo_vort_div.png white'

'c'

'reinit'
'open D:\sx\vorticity_scheme1.ctl'
'set grid off'
'set lon 60 150'
'set lat 20 60'
'set t 5'
'set z 5'
'set gxout shaded'
* 设置固定色标范围和天气图常规颜色
'set clevs -2.0 -1.5 -1.0 -0.5 0.0 0.5 1.0 1.5 2.0'
*'set ccols 4 5 6 7 10 11 12 13 14 15'
'd vort * 1e4'
'run E:\OpenGrADS-2.2\Contents\Resources\Scripts\cbar.gs'
'open D:\sx\geowind.ctl'
'set dfile 2'
'set lon 60 150'
'set lat 20 60'
'set t 5'
'set z 5'
'set gxout vector'
'set arrscl 0.8 50'
'set ccolor 0'
'd ug; vg'
'open D:\sx\divergence.ctl'
'set dfile 3'
'set lon 60 150'
'set lat 20 60'
'set t 5'
'set z 5'
'set gxout contour'
'set ccolor 1'
'd div * 1e4'
'draw title 25th day 00Z 500hPa Geostrophic Wind & Vorticity & Divergence'
'printim 500hPa_25day00Z_geo_vort_div.png white'

'c'

'reinit'
'open D:\sx\vorticity_scheme1.ctl'
'set grid off'
'set lon 60 150'
'set lat 20 60'
'set t 8'
'set z 5'
'set gxout shaded'
* 设置固定色标范围和天气图常规颜色
'set clevs -2.0 -1.5 -1.0 -0.5 0.0 0.5 1.0 1.5 2.0'
*'set ccols 4 5 6 7 10 11 12 13 14 15'
'd vort * 1e4'
'run E:\OpenGrADS-2.2\Contents\Resources\Scripts\cbar.gs'
'open D:\sx\geowind.ctl'
'set dfile 2'
'set lon 60 150'
'set lat 20 60'
'set t 8'
'set z 5'
'set gxout vector'
'set arrscl 0.8 50'
'set ccolor 0'
'd ug; vg'
'open D:\sx\divergence.ctl'
'set dfile 3'
'set lon 60 150'
'set lat 20 60'
'set t 8'
'set z 5'
'set gxout contour'
'set ccolor 1'
'd div * 1e4'
'draw title 26th day 12Z 500hPa Geostrophic Wind & Vorticity & Divergence'
'printim 500hPa_26day12Z_geo_vort_div.png white'

'say "四张图片已全部生成完成"'

nuistwangqin

E:\OpenGrADS-2.2\Contents\Resources\SampleDatasets\500hPa_25day00Z_geo_vort_div.png

nuistwangqin

http://bbs.06climate.com/forum.php?mod=attachment&aid=MTE2ODk0fDNlOTk3YjFjNzE0YjQ0YzY0YWRhY2Q5NzQ0YmI1NWNhfDE3NjQyNjE0NjM%3D&request=yes&_f=.png