分享一个画污染物浓度场的脚本

js250768097

这是一个画污染物浓度场的gs脚本，分享出来供大家参考
'reinit'
'open F:/ewaq/grads/data/camx/camxd04.2016052516.ctl'
'q time'
string=subwrd(result,3)
hour1=substr(string,1,2)
hour=hour1-16
day1=substr(string,4,2)
day=day1+1
day2=day+4
engmonth=substr(string,6,3)
year=substr(string,9,4)
tt=1
if(engmonth=JAN)
   month=1
   else
   if(engmonth=FEB)
   month=2
   else
   if(engmonth=Mar)
   month=3
   else
   if(engmonth=APR)
   month=4
   else
   if(engmonth=MAY)
   month=5
   else
   if(engmonth=JUN)
   month=6
   else
   if(engmonth=JUL)
   month=7
   else
   if(engmonth=AUG)
   month=8
   else
   if(engmonth=SEP)
   month=9
   else
   if(engmonth=OCT)
   month=10
   else
   if(engmonth=NOV)
   month=11
   else
   if(engmonth=DEC)
   month=12
   endif
   endif
   endif
   endif
   endif
   endif
   endif
   endif
   endif
   endif
   endif
   endif
while(day<day2)
while(hour<24)
   dd = day
   if(day < 10)
       dd=0%day
   endif
        if(hour<10)
          aaa=year%0%month%dd%0%hour
        else
          aaa=year%0%month%dd%hour
        endif
     'set grads off'
'set grid off'
     'set mpdraw off'
     'set lon 113.76 114.62'
     'set lat 22.45 22.87'
     'set t 'tt''
     'set gxout shaded'
'128colors.gs'
     'set clevs  0.1 1 2 6 15 25 30 40 60 80 100'
     'd so2'
     'set gxout vector'
     'set ccolor 1'
     'd skip(u,4,4);skip(v,4,4)'
     'draw shp CHN_adm3.shp'
     'cbar_matlab.gs'
     'draw title SO2_CAS-CAMx_d04_'%aaa%''
     'printim  F:\ewaq\grads\images\camx\2016050700\SO2_CAS-CAMX_'%aaa%'.png white'
     'c'
    'set grads off'
     'set gxout shaded'
     '128colors.gs'
    'set clevs  0.1 1 2 6 15 25 30 40 60 80 100'
     'd NO2'
     'set gxout vector'
     'set ccolor 1'
     'd skip(u,4,4);skip(v,4,4)'
     ' draw shp CHN_adm3.shp'
    'cbar_matlab.gs'
    'draw title NO2_CAS-CAMx_d04_'%aaa%''
    'printim  F:\ewaq\grads\images\camx\2016050700\NO2_CAS-CAMX_'%aaa%'.png x800 y600 white'
     'c'
     'set grads off'
     'set gxout shaded'
     '128colors.gs'
     'set clevs 0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 0.6 '
     'd CO'
     'set gxout vector'
     'set ccolor 1'
     'd skip(u,4,4);skip(v,4,4)'
     ' draw shp CHN_adm3.shp'
     'cbar_matlab.gs'
     'draw title CO_CAS-CAMx_d04_'%aaa%''
     'printim  F:\ewaq\grads\images\camx\2016050700\CO_CAS-CAMX_'%aaa%'.png x800 y600 white'
     'c'
     'set grads off'
     'set gxout shaded'
     '128colors.gs'
    'set clevs  5 15 30 45 60 75 90 105 125 135 145'
     'd O3'
     'set gxout vector'
     'set ccolor 1'
     'd skip(u,4,4);skip(v,4,4)'
     'draw shp CHN_adm3.shp'
     'cbar_matlab.gs'
     'draw title O3_CAS-CAMx_d04_'%aaa%''
     'printim  F:\ewaq\grads\images\camx\2016050700\O3_CAS-CAMX_'%aaa%'.png x800 y600 white'
     'c'
     'set grads off'
     'set gxout shaded'
     '128colors.gs'
    'set clevs  0.1 1 2 6 15 25 30 40 60 80 100'
     'd TPM25'
     'set gxout vector'
     'set ccolor 1'
     'd skip(u,4,4);skip(v,4,4)'
     ' draw shp CHN_adm3.shp'
    'cbar_matlab.gs'
    'draw title PM2.5_CAS-CAMx_d04_'%aaa%''
    'printim  F:\ewaq\grads\images\camx\2016050700\PM2.5_CAS-CAMX_'%aaa%'.png x800 y600 white'
     'c'
     'set grads off'
     'set gxout shaded'
     '128colors.gs'
    'set clevs  0.1 1 2 6 15 25 30 40 60 80 100'
     'd TPM10'
     'set gxout vector'
     'set ccolor 1'
     'd skip(u,4,4);skip(v,4,4)'
     ' draw shp CHN_adm3.shp'
     'cbar_matlab.gs'
     'draw title PM10_CAS-CAMx_d04__'%aaa%''
     'printim  F:\ewaq\grads\images\camx\2016050700\PM10_CAS-CAMX_'%aaa%'.png x800 y600 white'
     'c'
     hour=hour+1
     tt=tt+1
endwhile
hour=0
day=day+1
endwhile

'reinit'

林夕

谢谢楼主分享，不知楼主可否顺带分享一张楼主画的美图呀~

郭小侠V

楼主能不能说明什么数据，来源，并且图

qxtjls

请说明所用数据呀？