爱气象，爱气象家园！ 
·设为首页 开启辅助访问 新帖子新回复热门贴精华帖

气象家园

 找回密码
 立即注册

QQ登录

只需一步，快速开始

新浪微博登陆

只需一步, 快速开始

搜索
气象家园»论坛 › 编程作图 › GrADS › GrADS画时间序列图
返回列表 发新帖
查看: 6082|回复: 2

[秀图] GrADS画时间序列图

[复制链接]
气象家园蒙面人
气象家园蒙面人  发表于 2014-12-28 12:31:25 |阅读模式

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看，没有帐号？立即注册 新浪微博登陆

x
本帖最后由 匿名 于 2014-12-28 14:52 编辑

    大家好，我是刚刚接触GrADs软件，遇到了画图的问题，请大家指教帮忙，谢谢啦！遇到的问题是这样的：现在有预报结果postvar，一个月资料，现在想和NECP资料做对比分析，画出变量h,t,u,v,q2的时间序列和廓线图，我现在有老师给我的绘图软件包，但是里面只有h,t，u，v的绘图程序。想加上q2的程序，现在照着其他变量的程序写了一下q2的，但是一直没画出图，也找不到哪里的错误，请大家帮帮忙，可以加我的QQ：279995024，一起讨论，谢谢大家了！

程序有点长，是这样的：
ancor.sh_q:
cat > ancor_q.gs << EOF
'reinit'
txt1='${4}_${11}'
txt2='${5}_${11}'
txt3='${6}_${11}'
txt4='${7}_${11}'
if(${11}='xa')
  'open fnl.ctl'
  'open xa1.ctl'
  'open xa2.ctl'
  'open xa3.ctl'
  'open xa4.ctl'
endif
if(${11}='xb')
  'open fnl.ctl'
  'open xa5.ctl'
  'open xa6.ctl'
  'open xa7.ctl'
  'open xa8.ctl'
endif
if(${11}='xaxb')
  'open fnl.ctl'
  'open xa1.ctl'
  'open xa2.ctl'
  'open xa3.ctl'
  'open xa4.ctl'
  'open xa5.ctl'
  'open xa6.ctl'
  'open xa7.ctl'
  'open xa8.ctl'
  txt1='${4}_xa'
  txt2='${5}_xa'
  txt3='${6}_xa'
  txt4='${7}_xa'
  txt5='${4}_xb'
  txt6='${5}_xb'
  txt7='${6}_xb'
  txt8='${7}_xb'
endif
'enable print output/ana_grapes_${10}_${1}_q.gmf'
*------------------------------------------
cc.1=1
cc.2=6
cc.3=4
cc.4=3
ms.1=0
ms.2=0
ms.3=0
ms.4=0
mn.1=0
mn.2=0
mn.3=0
mn.4=0
lev.1=1000
lev.2=925
lev.3=850
lev.4=700
lev.5=500
lev.6=300
lev.7=250
lev.8=200
lev.9=150
lev.10=100
lev.11=70
lev.12=50
lev.13=30
lev.14=20
lev.15=10
*----------------------------------
lat0=20
lat1=90
'set t 1 $2 '
'set lat 0 '
'set lon 0'
i=$8
while(i<=$9)
'set grid off'
'set vrange -1.0 1.0'
if (i>12)
'set vrange -2.0 2.0'
endif
'set xlopts 1 5 0.15'
'set ylopts 1 5 0.15'
'set lat 0 '
'set lon 0'
'set lev 'lev.i' '
lat0=20
lat1=90
* plotting for SH
**bias
iks=1
ike=1
ltype.1=1
ltype.2=2
if(${11}='xaxb')
ike=2
endif
while(iks<=ike)
  j=1
  while(j<=$3)
    k=j+1
    if(iks=2)
      k=5+j
    endif
    'set grads off'
    'set grid off'
    'set ccolor 'cc.j''
    'set cmark 'ms.j''
    'set cstyle 'ltype.iks''
    'set cthick 8'
*   say k
    'd tloop(aave(q2.'k'-q.1,lon=0,lon=360,lat=-'lat1',lat=-'lat0'))'
    j=j+1
  endwhile
  iks=iks+1
endwhile
'set ccolor 4 '
'set cstyle 5 '
'set cmark 0'
'set cthick 8'
'd tloop(0.0)'
if (${11}!='xaxb')
  if ($3=1)
  'cbar_line_thick -x 2.5 -y 7 -c 'cc.1' -m 'ms.1' -l 'ltype.1' -t "'txt1'-S" '
  endif
  if ($3=2)
  'cbar_line_thick -x 2.5 -y 7 -c 'cc.1' 'cc.2' -m 'ms.1' 'ms.2' -l 'ltype.1' 'ltype.1' -t "'txt1'-S" "'txt2'-S" '
  endif
  if ($3=3)
  'cbar_line_thick -x 2.5 -y 7 -c 'cc.1' 'cc.2' 'cc.3' -m 'ms.1' 'ms.2' 'ms.3' -l 'ltype.1' 'ltype.1' 'ltype.1' -t "'txt1'-S" "'txt2'-S" "'txt3'-S" '
  endif
  if ($3=4)
  'cbar_line_thick -x 2.5 -y 7 -c 'cc.1' 'cc.2' 'cc.3' 'cc.4' -m 'ms.1' 'ms.2' 'ms.3' 'ms.4' -l 'ltype.1' 'ltype.1' 'ltype.1' 'ltype.1' -t "'txt1'-S" "'txt2'-S" "'txt3'-S" "'txt4'-S" '
  endif
else
  if ($3=1)
  'cbar_line_thick -x 2.5 -y 7 -c 'cc.1' 'cc.1' -m 'ms.1' 'ms.1' -l 'ltype.1' 'ltype.2' -t "'txt1'-S" "'txt5'-S"'
  endif
  if ($3=2)
  'cbar_line_thick -x 2.5 -y 7 -c 'cc.1' 'cc.2' 'cc.1' 'cc.2' -m 'ms.1' 'ms.2' 'ms.1' 'ms.2' -l 'ltype.1' 'ltype.1' 'ltype.2' 'ltype.2' -t "'txt1'-S" "'txt2'-S" "'txt5'-S" "'txt6'-S"'
  endif
  if ($3=3)
  'cbar_line_thick -x 2.5 -y 7 -c 'cc.1' 'cc.2' 'cc.3' 'cc.1' 'cc.2' 'cc.3' -m 'ms.1' 'ms.2' 'ms.3' 'ms.1' 'ms.2' 'ms.3' -l 'ltype.1' 'ltype.1' 'ltype.1' 'ltype.2' 'ltype.2' 'ltype.2' -t "'txt1'-S" "'txt2'-S" "'txt3'-S" "'txt5'-S" "'txt6'-S" "'txt7'-S"'
  endif
  if ($3=4)
  'cbar_line_thick -x 2.5 -y 7 -c 'cc.1' 'cc.2' 'cc.3' 'cc.4' 'cc.1' 'cc.2' 'cc.3' 'cc.4' -m 'ms.1' 'ms.2' 'ms.3' 'ms.4' 'ms.1' 'ms.2' 'ms.3' 'ms.4' -l 'ltype.1' 'ltype.1' 'ltype.1' 'ltype.1' 'ltype.2' 'ltype.2' 'ltype.2' 'ltype.2' -t "'txt1'-S" "'txt2'-S" "'txt3'-S" "'txt4'-S" "'txt5'-S" "'txt6'-S" "'txt7'-S" "'txt8'-S"'
  endif
endif
'draw title Bias of q(grapes)-q(${10}) S. Hemis ('lev.i'hPa)'
'print'
'c'
**rms
'set grads off'
'set grid off'
'set vrange 0 5'
iks=1
ike=1
ltype.1=1
ltype.2=2
if(${11}='xaxb')
ike=2
endif
while(iks<=ike)
  j=1
  while(j<=$3)
    k=j+1
    if(iks=2)
      k=5+j
    endif
    'set grads off'
    'set grid off'
    'set ccolor 'cc.j''
    'set cmark 'ms.j''
    'set cstyle 'ltype.iks''
    'set cthick 8'
    'd tloop(sqrt(aave((q2.'k'-q.1)*(q2.'k'-q.1),lon=0,lon=360,lat=-'lat1',lat=-'lat0')))'
    j=j+1
  endwhile
  iks=iks+1
endwhile
if (${11}!='xaxb')
  if ($3=1)
  'cbar_line_thick -x 2.5 -y 7 -c 'cc.1' -m 'ms.1' -l 'ltype.1' -t "'txt1'-S" '
  endif
  if ($3=2)
  'cbar_line_thick -x 2.5 -y 7 -c 'cc.1' 'cc.2' -m 'ms.1' 'ms.2' -l 'ltype.1' 'ltype.1' -t "'txt1'-S" "'txt2'-S" '
  endif
  if ($3=3)
  'cbar_line_thick -x 2.5 -y 7 -c 'cc.1' 'cc.2' 'cc.3' -m 'ms.1' 'ms.2' 'ms.3' -l 'ltype.1' 'ltype.1' 'ltype.1' -t "'txt1'-S" "'txt2'-S" "'txt3'-S" '
  endif
  if ($3=4)
  'cbar_line_thick -x 2.5 -y 7 -c 'cc.1' 'cc.2' 'cc.3' 'cc.4' -m 'ms.1' 'ms.2' 'ms.3' 'ms.4' -l 'ltype.1' 'ltype.1' 'ltype.1' 'ltype.1' -t "'txt1'-S" "'txt2'-S" "'txt3'-S" "'txt4'-S" '
  endif
else
  if ($3=1)
  'cbar_line_thick -x 2.5 -y 7 -c 'cc.1' 'cc.1' -m 'ms.1' 'ms.1' -l 'ltype.1' 'ltype.2' -t "'txt1'-S" "'txt5'-S"'
  endif
  if ($3=2)
  'cbar_line_thick -x 2.5 -y 7 -c 'cc.1' 'cc.2' 'cc.1' 'cc.2' -m 'ms.1' 'ms.2' 'ms.1' 'ms.2' -l 'ltype.1' 'ltype.1' 'ltype.2' 'ltype.2' -t "'txt1'-S" "'txt2'-S" "'txt5'-S" "'txt6'-S"'
  endif
  if ($3=3)
  'cbar_line_thick -x 2.5 -y 7 -c 'cc.1' 'cc.2' 'cc.3' 'cc.1' 'cc.2' 'cc.3' -m 'ms.1' 'ms.2' 'ms.3' 'ms.1' 'ms.2' 'ms.3' -l 'ltype.1' 'ltype.1' 'ltype.1' 'ltype.2' 'ltype.2' 'ltype.2' -t "'txt1'-S" "'txt2'-S" "'txt3'-S" "'txt5'-S" "'txt6'-S" "'txt7'-S"'
  endif
  if ($3=4)
  'cbar_line_thick -x 2.5 -y 7 -c 'cc.1' 'cc.2' 'cc.3' 'cc.4' 'cc.1' 'cc.2' 'cc.3' 'cc.4' -m 'ms.1' 'ms.2' 'ms.3' 'ms.4' 'ms.1' 'ms.2' 'ms.3' 'ms.4' -l 'ltype.1' 'ltype.1' 'ltype.1' 'ltype.1' 'ltype.2' 'ltype.2' 'ltype.2' 'ltype.2' -t "'txt1'-S" "'txt2'-S" "'txt3'-S" "'txt4'-S" "'txt5'-S" "'txt6'-S" "'txt7'-S" "'txt8'-S"'
  endif
endif
'draw title RMS of q(grapes)-q(${10}) S. Hemis ('lev.i'hPa)'
'print'
'c'
* plotting for NH
'set grads off'
'set grid off'
'set vrange -1.0 1.0'
if (i>12)
'set vrange -2.0 2.0'
endif
'set xlopts 1 5 0.15'
'set ylopts 1 5 0.15'
iks=1
ike=1
ltype.1=1
ltype.2=2
if(${11}='xaxb')
ike=2
endif
while(iks<=ike)
  j=1
  while(j<=$3)
    k=j+1
    if(iks=2)
      k=5+j
    endif
    'set grads off'
    'set grid off'
    'set ccolor 'cc.j''
    'set cmark 'ms.j''
    'set cstyle 'ltype.iks''
    'set cthick 8'
    'd tloop(aave(q2.'k'-q.1,lon=0,lon=360,lat='lat0',lat='lat1'))'
    j=j+1
  endwhile
  iks=iks+1
endwhile
'set ccolor 4 '
'set cstyle 5 '
'set cmark 0'
'set cthick 8'
'd tloop(0.0)'
if (${11}!='xaxb')
  if ($3=1)
  'cbar_line_thick -x 2.5 -y 7 -c 'cc.1' -m 'mn.1' -l 'ltype.1' -t "'txt1'-N" '
  endif
  if ($3=2)
  'cbar_line_thick -x 2.5 -y 7 -c 'cc.1' 'cc.2' -m 'mn.1' 'mn.2' -l 'ltype.1' 'ltype.1' -t "'txt1'-N" "'txt2'-N" '
  endif
  if ($3=3)
  'cbar_line_thick -x 2.5 -y 7 -c 'cc.1' 'cc.2' 'cc.3' -m 'mn.1' 'mn.2' 'mn.3' -l 'ltype.1' 'ltype.1' 'ltype.1' -t "'txt1'-N" "'txt2'-N" "'txt3'-N" '
  endif
  if ($3=4)
  'cbar_line_thick -x 2.5 -y 7 -c 'cc.1' 'cc.2' 'cc.3' 'cc.4' -m 'mn.1' 'mn.2' 'mn.3' 'mn.4' -l 'ltype.1' 'ltype.1' 'ltype.1' 'ltype.1' -t "'txt1'-N" "'txt2'-N" "'txt3'-N" "'txt4'-N" '
  endif
else
  if ($3=1)
  'cbar_line_thick -x 2.5 -y 7 -c 'cc.1' 'cc.1' -m 'mn.1' 'mn.1' -l 'ltype.1' 'ltype.2' -t "'txt1'-N" "'txt5'-N"'
  endif
  if ($3=2)
  'cbar_line_thick -x 2.5 -y 7 -c 'cc.1' 'cc.2' 'cc.1' 'cc.2' -m 'mn.1' 'mn.2' 'mn.1' 'mn.2' -l 'ltype.1' 'ltype.1' 'ltype.2' 'ltype.2' -t "'txt1'-N" "'txt2'-N" "'txt5'-N" "'txt6'-N"'
  endif
  if ($3=3)
  'cbar_line_thick -x 2.5 -y 7 -c 'cc.1' 'cc.2' 'cc.3' 'cc.1' 'cc.2' 'cc.3' -m 'mn.1' 'mn.2' 'mn.3' 'mn.1' 'mn.2' 'mn.3' -l 'ltype.1' 'ltype.1' 'ltype.1' 'ltype.2' 'ltype.2' 'ltype.2' -t "'txt1'-N" "'txt2'-N" "'txt3'-N" "'txt5'-N" "'txt6'-N" "'txt7'-N"'
  endif
  if ($3=4)
  'cbar_line_thick -x 2.5 -y 7 -c 'cc.1' 'cc.2' 'cc.3' 'cc.4' 'cc.1' 'cc.2' 'cc.3' 'cc.4' -m 'mn.1' 'mn.2' 'mn.3' 'mn.4' 'mn.1' 'mn.2' 'mn.3' 'mn.4' -l 'ltype.1' 'ltype.1' 'ltype.1' 'ltype.1' 'ltype.2' 'ltype.2' 'ltype.2' 'ltype.2' -t "'txt1'-N" "'txt2'-N" "'txt3'-N" "'txt4'-N" "'txt5'-N" "'txt6'-N" "'txt7'-N" "'txt8'-N"'
  endif
endif
'draw title Bias of q(grapes)-q(${10}) N. Hemis ('lev.i'hPa)'
'print'
'c'
**rms
'set grads off'
'set grid off'
'set vrange 0 4'
iks=1
ike=1
ltype.1=1
ltype.2=2
if(${11}='xaxb')
ike=2
endif
while(iks<=ike)
  j=1
  while(j<=$3)
    k=j+1
    if(iks=2)
      k=5+j
    endif
    'set grads off'
    'set grid off'
    'set ccolor 'cc.j''
    'set cmark 'mn.j''
    'set cstyle 'ltype.iks''
    'set cthick 8'
    'd tloop(sqrt(aave((q2.'k'-q.1)*(q2.'k'-q.1),lon=0,lon=360,lat='lat0',lat='lat1')))'
    j=j+1
  endwhile
  iks=iks+1
endwhile
if (${11}!='xaxb')
  if ($3=1)
  'cbar_line_thick -x 2.5 -y 7 -c 'cc.1' -m 'mn.1' -l 'ltype.1' -t "'txt1'-N" '
  endif
  if ($3=2)
  'cbar_line_thick -x 2.5 -y 7 -c 'cc.1' 'cc.2' -m 'mn.1' 'mn.2' -l 'ltype.1' 'ltype.1' -t "'txt1'-N" "'txt2'-N" '
  endif
  if ($3=3)
  'cbar_line_thick -x 2.5 -y 7 -c 'cc.1' 'cc.2' 'cc.3' -m 'mn.1' 'mn.2' 'mn.3' -l 'ltype.1' 'ltype.1' 'ltype.1' -t "'txt1'-N" "'txt2'-N" "'txt3'-N" '
  endif
  if ($3=4)
  'cbar_line_thick -x 2.5 -y 7 -c 'cc.1' 'cc.2' 'cc.3' 'cc.4' -m 'mn.1' 'mn.2' 'mn.3' 'mn.4' -l 'ltype.1' 'ltype.1' 'ltype.1' 'ltype.1' -t "'txt1'-N" "'txt2'-N" "'txt3'-N" "'txt4'-N" '
  endif
else
  if ($3=1)
  'cbar_line_thick -x 2.5 -y 7 -c 'cc.1' 'cc.1' -m 'mn.1' 'mn.1' -l 'ltype.1' 'ltype.2' -t "'txt1'-N" "'txt5'-N"'
  endif
  if ($3=2)
  'cbar_line_thick -x 2.5 -y 7 -c 'cc.1' 'cc.2' 'cc.1' 'cc.2' -m 'mn.1' 'mn.2' 'mn.1' 'mn.2' -l 'ltype.1' 'ltype.1' 'ltype.2' 'ltype.2' -t "'txt1'-N" "'txt2'-N" "'txt5'-N" "'txt6'-N"'
  endif
  if ($3=3)
  'cbar_line_thick -x 2.5 -y 7 -c 'cc.1' 'cc.2' 'cc.3' 'cc.1' 'cc.2' 'cc.3' -m 'mn.1' 'mn.2' 'mn.3' 'mn.1' 'mn.2' 'mn.3' -l 'ltype.1' 'ltype.1' 'ltype.1' 'ltype.2' 'ltype.2' 'ltype.2' -t "'txt1'-N" "'txt2'-N" "'txt3'-N" "'txt5'-N" "'txt6'-N" "'txt7'-N"'
  endif
  if ($3=4)
  'cbar_line_thick -x 2.5 -y 7 -c 'cc.1' 'cc.2' 'cc.3' 'cc.4' 'cc.1' 'cc.2' 'cc.3' 'cc.4' -m 'mn.1' 'mn.2' 'mn.3' 'mn.4' 'mn.1' 'mn.2' 'mn.3' 'mn.4' -l 'ltype.1' 'ltype.1' 'ltype.1' 'ltype.1' 'ltype.2' 'ltype.2' 'ltype.2' 'ltype.2' -t "'txt1'-N" "'txt2'-N" "'txt3'-N" "'txt4'-N" "'txt5'-N" "'txt6'-N" "'txt7'-N" "'txt8'-N"'
  endif
endif
'draw title RMS of q(grapes)-q(${10}) N. Hemis ('lev.i'hPa)'
'print'
'c'

* plotting for Tro.
**bias
'set grads off'
'set grid off'
'set vrange -1.0 1.0'
if (i>12)
'set vrange -2.0 2.0'
endif
iks=1
ike=1
ltype.1=1
ltype.2=2
if(${11}='xaxb')
ike=2
endif
while(iks<=ike)
  j=1
  while(j<=$3)
    k=j+1
    if(iks=2)
      k=5+j
    endif
    'set grads off'
    'set grid off'
    'set ccolor 'cc.j''
    'set cmark 'ms.j''
    'set cstyle 'ltype.iks''
    'set cthick 8'
    'd tloop(aave(q2.'k'-q.1,lon=0,lon=360,lat=-'lat0',lat='lat0'))'
    j=j+1
  endwhile
  iks=iks+1
endwhile
'set ccolor 4 '
'set cstyle 5 '
'set cmark 0'
'set cthick 8'
'd tloop(0.0)'
if (${11}!='xaxb')
  if ($3=1)
  'cbar_line_thick -x 2.5 -y 7 -c 'cc.1' -m 'mn.1' -l 'ltype.1' -t "'txt1'-Tro." '
  endif
  if ($3=2)
  'cbar_line_thick -x 2.5 -y 7 -c 'cc.1' 'cc.2' -m 'mn.1' 'mn.2' -l 'ltype.1' 'ltype.1' -t "'txt1'-Tro." "'txt2'-Tro." '
  endif
  if ($3=3)
  'cbar_line_thick -x 2.5 -y 7 -c 'cc.1' 'cc.2' 'cc.3' -m 'mn.1' 'mn.2' 'mn.3' -l 'ltype.1' 'ltype.1' 'ltype.1' -t "'txt1'-Tro." "'txt2'-Tro." "'txt3'-Tro." '
  endif
  if ($3=4)
  'cbar_line_thick -x 2.5 -y 7 -c 'cc.1' 'cc.2' 'cc.3' 'cc.4' -m 'mn.1' 'mn.2' 'mn.3' 'mn.4' -l 'ltype.1' 'ltype.1' 'ltype.1' 'ltype.1' -t "'txt1'-Tro." "'txt2'-Tro." "'txt3'-Tro." "'txt4'-Tro." '
  endif
else
  if ($3=1)
  'cbar_line_thick -x 2.5 -y 7 -c 'cc.1' 'cc.1' -m 'mn.1' 'mn.1' -l 'ltype.1' 'ltype.2' -t "'txt1'-Tro." "'txt5'-Tro."'
  endif
  if ($3=2)
  'cbar_line_thick -x 2.5 -y 7 -c 'cc.1' 'cc.2' 'cc.1' 'cc.2' -m 'mn.1' 'mn.2' 'mn.1' 'mn.2' -l 'ltype.1' 'ltype.1' 'ltype.2' 'ltype.2' -t "'txt1'-Tro." "'txt2'-Tro." "'txt5'-Tro." "'txt6'-Tro."'
  endif
  if ($3=3)
  'cbar_line_thick -x 2.5 -y 7 -c 'cc.1' 'cc.2' 'cc.3' 'cc.1' 'cc.2' 'cc.3' -m 'mn.1' 'mn.2' 'mn.3' 'mn.1' 'mn.2' 'mn.3' -l 'ltype.1' 'ltype.1' 'ltype.1' 'ltype.2' 'ltype.2' 'ltype.2' -t "'txt1'-Tro." "'txt2'-Tro." "'txt3'-Tro." "'txt5'-Tro." "'txt6'-Tro." "'txt7'-Tro."'
  endif
  if ($3=4)
  'cbar_line_thick -x 2.5 -y 7 -c 'cc.1' 'cc.2' 'cc.3' 'cc.4' 'cc.1' 'cc.2' 'cc.3' 'cc.4' -m 'mn.1' 'mn.2' 'mn.3' 'mn.4' 'mn.1' 'mn.2' 'mn.3' 'mn.4' -l 'ltype.1' 'ltype.1' 'ltype.1' 'ltype.1' 'ltype.2' 'ltype.2' 'ltype.2' 'ltype.2' -t "'txt1'-Tro." "'txt2'-Tro." "'txt3'-Tro." "'txt4'-Tro." "'txt5'-Tro." "'txt6'-Tro." "'txt7'-Tro." "'txt8'-Tro."'
  endif
endif
'draw title Bias of q(grapes)-q(${10}) Tropics ('lev.i'hPa)'
'print'
'c'
**rms
'set grads off'
'set grid off'
'set vrange 0 4'
iks=1
ike=1
ltype.1=1
ltype.2=2
if(${11}='xaxb')
ike=2
endif
while(iks<=ike)
  j=1
  while(j<=$3)
    k=j+1
    if(iks=2)
      k=5+j
    endif
    'set grads off'
    'set grid off'
    'set ccolor 'cc.j''
    'set cmark 'mn.j''
    'set cstyle 'ltype.iks''
    'set cthick 8'
    'd tloop(sqrt(aave((q2.'k'-q.1)*(q2.'k'-q.1),lon=0,lon=360,lat=-'lat0',lat='lat0')))'
    j=j+1
  endwhile
  iks=iks+1
endwhile
if (${11}!='xaxb')
  if ($3=1)
  'cbar_line_thick -x 2.5 -y 7 -c 'cc.1' -m 'mn.1' -l 'ltype.1' -t "'txt1'-Tro." '
  endif
  if ($3=2)
  'cbar_line_thick -x 2.5 -y 7 -c 'cc.1' 'cc.2' -m 'mn.1' 'mn.2' -l 'ltype.1' 'ltype.1' -t "'txt1'-Tro." "'txt2'-Tro." '
  endif
  if ($3=3)
  'cbar_line_thick -x 2.5 -y 7 -c 'cc.1' 'cc.2' 'cc.3' -m 'mn.1' 'mn.2' 'mn.3' -l 'ltype.1' 'ltype.1' 'ltype.1' -t "'txt1'-Tro." "'txt2'-Tro." "'txt3'-Tro." '
  endif
  if ($3=4)
  'cbar_line_thick -x 2.5 -y 7 -c 'cc.1' 'cc.2' 'cc.3' 'cc.4' -m 'mn.1' 'mn.2' 'mn.3' 'mn.4' -l 'ltype.1' 'ltype.1' 'ltype.1' 'ltype.1' -t "'txt1'-Tro." "'txt2'-Tro." "'txt3'-Tro." "'txt4'-Tro." '
  endif
else
  if ($3=1)
  'cbar_line_thick -x 2.5 -y 7 -c 'cc.1' 'cc.1' -m 'mn.1' 'mn.1' -l 'ltype.1' 'ltype.2' -t "'txt1'-Tro." "'txt5'-Tro."'
  endif
  if ($3=2)
  'cbar_line_thick -x 2.5 -y 7 -c 'cc.1' 'cc.2' 'cc.1' 'cc.2' -m 'mn.1' 'mn.2' 'mn.1' 'mn.2' -l 'ltype.1' 'ltype.1' 'ltype.2' 'ltype.2' -t "'txt1'-Tro." "'txt2'-Tro." "'txt5'-Tro." "'txt6'-Tro."'
  endif
  if ($3=3)
  'cbar_line_thick -x 2.5 -y 7 -c 'cc.1' 'cc.2' 'cc.3' 'cc.1' 'cc.2' 'cc.3' -m 'mn.1' 'mn.2' 'mn.3' 'mn.1' 'mn.2' 'mn.3' -l 'ltype.1' 'ltype.1' 'ltype.1' 'ltype.2' 'ltype.2' 'ltype.2' -t "'txt1'-Tro." "'txt2'-Tro." "'txt3'-Tro." "'txt5'-Tro." "'txt6'-Tro." "'txt7'-Tro."'
  endif
  if ($3=4)
  'cbar_line_thick -x 2.5 -y 7 -c 'cc.1' 'cc.2' 'cc.3' 'cc.4' 'cc.1' 'cc.2' 'cc.3' 'cc.4' -m 'mn.1' 'mn.2' 'mn.3' 'mn.4' 'mn.1' 'mn.2' 'mn.3' 'mn.4' -l 'ltype.1' 'ltype.1' 'ltype.1' 'ltype.1' 'ltype.2' 'ltype.2' 'ltype.2' 'ltype.2' -t "'txt1'-Tro." "'txt2'-Tro." "'txt3'-Tro." "'txt4'-Tro." "'txt5'-Tro." "'txt6'-Tro." "'txt7'-Tro." "'txt8'-Tro."'
  endif
endif
'draw title RMS of q(grapes)-q(${10}) Tropics ('lev.i'hPa)'
'print'
'c'
i=i+1
endwhile
'disable print'
* for averaged profile
lat0=20
lat1=90
'set lat 0'
'set lon 0'
'set fwrite output/prof_${1}_q.dat'
'set gxout fwrite'
tt=1
while(tt<=$2)
iks=1
ike=1
if(${11}='xaxb')
ike=2
endif
while(iks<=ike)
j=1                                                                        
while(j<=$3)                                                                 
  i=1
  while(i<=$9)
     'set lev 'lev.i' '
     'set t 'tt' '
      k=j+1
      if(iks=2)
        k=5+j
      endif
      'define aa=aave(q2.'k'-q.1,lon=0,lon=360,lat=-'lat1',lat=-'lat0')'
      'd aa'
    i=i+1
  endwhile
  i=1                                                                        
  while(i<=$9)                                                               
     'set lev 'lev.i' '                                                      
     'set t 'tt' '                                                            
      k=j+1
      if(iks=2)
        k=5+j
      endif
      'define aa=sqrt(aave((q2.'k'-q.1)*(q2.'k'-q.1),lon=0,lon=360,lat=-'lat1',lat=-'lat0'))'
      'd aa'                                                                  
    i=i+1
  endwhile
  i=1
  while(i<=$9)
     'set lev 'lev.i' '
     'set t 'tt' '
      k=j+1
      if(iks=2)
        k=5+j
      endif
      'define bb=aave(q2.'k'-q.1,lon=0,lon=360,lat='lat0',lat='lat1')'
      'd bb'
    i=i+1
  endwhile
  i=1
  while(i<=$9)                                                               
     'set lev 'lev.i' '                                                      
     'set t 'tt' '                                                            
      k=j+1
      if(iks=2)
        k=5+j
      endif
      'define bb=sqrt(aave((q2.'k'-q.1)*(q2.'k'-q.1),lon=0,lon=360,lat='lat0',lat='lat1'))'
      'd bb'                                                                  
    i=i+1                                                                     
  endwhile
  i=1
  while(i<=$9)                                                                  
     'set lev 'lev.i' '                                                         
     'set t 'tt' '                                                               
      k=j+1
      if(iks=2)
        k=5+j
      endif
      'define cc=aave(q2.'k'-q.1,lon=0,lon=360,lat=-'lat0',lat='lat0')'            
      'd cc'                                                                     
    i=i+1
  endwhile
  i=1
  while(i<=$9)                                                                  
     'set lev 'lev.i' '                                                         
     'set t 'tt' '                                                               
      k=j+1
      if(iks=2)
        k=5+j
      endif
      'define cc=sqrt(aave((q2.'k'-q.1)*(q2.'k'-q.1),lon=0,lon=360,lat=-'lat0',lat='lat0'))'
      'd cc'                                                                     
    i=i+1                                                                        
  endwhile
  j=j+1
endwhile
iks=iks+1
endwhile
tt=tt+1
endwhile
'disable fwrite'
quit
EOF

pltprof.sh_q:
cat > pltprof_q.gs << EOF
'reinit'
'open output/prof_${1}_q.ctl'
'set zlog on'
'set xlopts 1 6 0.15'
'set ylopts 1 6 0.15'
t1='${4}_${9}'
t2='${5}_${9}'
t3='${6}_${9}'
t4='${7}_${9}'
iks=1
ike=1
if(${9}='xaxb')
  t1='${4}_xa'
  t2='${5}_xa'
  t3='${6}_xa'
  t4='${7}_xa'
  t5='${4}_xb'
  t6='${5}_xb'
  t7='${6}_xb'
  t8='${7}_xb'
  ike=2
endif
rmsn.1=rms_n1
rmsn.2=rms_n2
rmsn.3=rms_n3
rmsn.4=rms_n4
rmsn.5=rms_n1_f
rmsn.6=rms_n2_f
rmsn.7=rms_n3_f
rmsn.8=rms_n4_f
rmss.1=rms_s1
rmss.2=rms_s2
rmss.3=rms_s3
rmss.4=rms_s4
rmss.5=rms_s1_f
rmss.6=rms_s2_f
rmss.7=rms_s3_f
rmss.8=rms_s4_f
rmseq.1=rms_eq1
rmseq.2=rms_eq2
rmseq.3=rms_eq3
rmseq.4=rms_eq4
rmseq.5=rms_eq1_f
rmseq.6=rms_eq2_f
rmseq.7=rms_eq3_f
rmseq.8=rms_eq4_f
biasn.1=bias_n1
biasn.2=bias_n2
biasn.3=bias_n3
biasn.4=bias_n4
biasn.5=bias_n1_f
biasn.6=bias_n2_f
biasn.7=bias_n3_f
biasn.8=bias_n4_f
biass.1=bias_s1
biass.2=bias_s2
biass.3=bias_s3
biass.4=bias_s4
biass.5=bias_s1_f
biass.6=bias_s2_f
biass.7=bias_s3_f
biass.8=bias_s4_f
biaseq.1=bias_eq1
biaseq.2=bias_eq2
biaseq.3=bias_eq3
biaseq.4=bias_eq4
biaseq.5=bias_eq1_f
biaseq.6=bias_eq2_f
biaseq.7=bias_eq3_f
biaseq.8=bias_eq4_f
time1=${2}
cc.1=1
cc.2=6
cc.3=4
cc.4=3
ltype.1=1
ltype.2=2
'enable print output/prof_${8}_${1}_q.gmf'
while(iks<=ike)
j=1
while(j<=${3})
  k=j
  if(iks=2)
    k=4+j
  endif
  'set vrange 5 50'
  'set ylevs 1000 925 850 700 500 300 250 200 150 100 70 50 30 20 10'
  'set cthick 8'
  'set ccolor 'cc.j' '
  'set cmark 0'
  'set cstyle 'ltype.iks''
  'set grads off'
  'd ave('rmsn.k',t=1,t='time1')'
  j=j+1
endwhile
iks=iks+1
endwhile
if(${9}!='xaxb')
  if(${3}=1)
    './cbar_line_thick.gs -x 5.2 -y 2.8 -c 'cc.1' -m 0 -l 'ltype.1' -t "'t1'-N" '
  endif
  if(${3}=2)
    './cbar_line_thick.gs -x 5.2 -y 2.8 -c 'cc.1' 'cc.2' -m 0 0 -l 'ltype.1' 'ltype.1' -t "'t1'-N" "'t2'-N" '
  endif
  if(${3}=3)
    './cbar_line_thick.gs -x 5.2 -y 2.8 -c 'cc.1' 'cc.2' 'cc.3' -m 0 0 0 -l 'ltype.1' 'ltype.1' 'ltype.1' -t "'t1'-N" "'t2'-N" "'t3'-N"'
  endif
  if(${3}=4)
    './cbar_line_thick.gs -x 5.2 -y 2.8 -c 'cc.1' 'cc.2' 'cc.3' 'cc.4' -m 0 0 0 0 -l 'ltype.1' 'ltype.1' 'ltype.1' 'ltype.1' -t "'t1'-N" "'t2'-N" "'t3'-N" "'t4'-N"'
  endif
else
  if(${3}=1)
    './cbar_line_thick.gs -x 5.2 -y 2.8 -c 'cc.1' 'cc.1'  -m 0 0 -l 'ltype.1' 'ltype.2' -t "'t1'-N" "'t5'-N" '
  endif
  if(${3}=2)
    './cbar_line_thick.gs -x 5.2 -y 2.8 -c 'cc.1' 'cc.2' 'cc.1' 'cc.2' -m 0 0 0 0 -l 'ltype.1' 'ltype.1' 'ltype.2' 'ltype.2' -t "'t1'-N" "'t2'-N" "'t5'-N" "'t6'-N"'
  endif
  if(${3}=3)
    './cbar_line_thick.gs -x 5.2 -y 2.8 -c 'cc.1' 'cc.2' 'cc.3' 'cc.1' 'cc.2' 'cc.3' -m 0 0 0 0 0 0 -l 'ltype.1' 'ltype.1' 'ltype.1' 'ltype.2' 'ltype.2' 'ltype.2' -t "'t1'-N" "'t2'-N" "'t3'-N" "'t5'-N" "'t6'-N" "'t7'-N"'
  endif
  if(${3}=4)
    './cbar_line_thick.gs -x 5.2 -y 2.8 -c 'cc.1' 'cc.2' 'cc.3' 'cc.4' 'cc.1' 'cc.2' 'cc.3' 'cc.4' -m 0 0 0 0 0 0 0 0 -l 'ltype.1' 'ltype.1' 'ltype.1' 'ltype.1' 'ltype.2' 'ltype.2' 'ltype.2' 'ltype.2' -t "'t1'-N" "'t2'-N" "'t3'-N" "'t4'-N" "'t5'-N" "'t6'-N" "'t7'-N" "'t8'-N"'
  endif
  
endif
'draw title time-averaged q(grapes)-q(${8}) RMS of N.Hemis'
'draw ylab Pressure (unit: hPa) '
'draw xlab RMS'
'print'
'c'
iks=1
while(iks<=ike)
j=1
while(j<=${3})
  k=j
  if(iks=2)
    k=4+j
  endif
  'set vrange -15 15'
  'set ylevs 1000 925 850 700 500 300 250 200 150 100 70 50 30 20 10'
  'set cthick 8'
  'set ccolor 'cc.j' '
  'set cmark 0'
  'set cstyle 'ltype.iks''
*  say 'biasn.k ' biasn.k
  'set grads off'
  'd ave('biasn.k',t=1,t='time1')'
  j=j+1
endwhile
iks=iks+1
endwhile
'set cstyle 5'
'set ccolor 1'
'set cmark 0'
'set cthick 8'
'd 'biasn.k'-'biasn.k''
if(${9}!='xaxb')
  if(${3}=1)
    './cbar_line_thick.gs -x 5.2 -y 2.8 -c 'cc.1' -m 0 -l 'ltype.1' -t "'t1'-N" '
  endif
  if(${3}=2)
    './cbar_line_thick.gs -x 5.2 -y 2.8 -c 'cc.1' 'cc.2' -m 0 0 -l 'ltype.1' 'ltype.1' -t "'t1'-N" "'t2'-N" '
  endif
  if(${3}=3)
    './cbar_line_thick.gs -x 5.2 -y 2.8 -c 'cc.1' 'cc.2' 'cc.3' -m 0 0 0 -l 'ltype.1' 'ltype.1' 'ltype.1' -t "'t1'-N" "'t2'-N" "'t3'-N"'
  endif
  if(${3}=4)
    './cbar_line_thick.gs -x 5.2 -y 2.8 -c 'cc.1' 'cc.2' 'cc.3' 'cc.4' -m 0 0 0 0 -l 'ltype.1' 'ltype.1' 'ltype.1' 'ltype.1' -t "'t1'-N" "'t2'-N" "'t3'-N" "'t4'-N"'
  endif
else
  if(${3}=1)
    './cbar_line_thick.gs -x 5.2 -y 2.8 -c 'cc.1' 'cc.1'  -m 0 0 -l 'ltype.1' 'ltype.2' -t "'t1'-N" "'t5'-N" '
  endif
  if(${3}=2)
    './cbar_line_thick.gs -x 5.2 -y 2.8 -c 'cc.1' 'cc.2' 'cc.1' 'cc.2' -m 0 0 0 0 -l 'ltype.1' 'ltype.1' 'ltype.2' 'ltype.2' -t "'t1'-N" "'t2'-N" "'t5'-N" "'t6'-N"'
  endif
  if(${3}=3)
    './cbar_line_thick.gs -x 5.2 -y 2.8 -c 'cc.1' 'cc.2' 'cc.3' 'cc.1' 'cc.2' 'cc.3' -m 0 0 0 0 0 0 -l 'ltype.1' 'ltype.1' 'ltype.1' 'ltype.2' 'ltype.2' 'ltype.2' -t "'t1'-N" "'t2'-N" "'t3'-N" "'t5'-N" "'t6'-N" "'t7'-N"'
  endif
  if(${3}=4)
    './cbar_line_thick.gs -x 5.2 -y 2.8 -c 'cc.1' 'cc.2' 'cc.3' 'cc.4' 'cc.1' 'cc.2' 'cc.3' 'cc.4' -m 0 0 0 0 0 0 0 0 -l 'ltype.1' 'ltype.1' 'ltype.1' 'ltype.1' 'ltype.2' 'ltype.2' 'ltype.2' 'ltype.2' -t "'t1'-N" "'t2'-N" "'t3'-N" "'t4'-N" "'t5'-N" "'t6'-N" "'t7'-N" "'t8'-N"'
  endif
  
endif
'draw title time-averaged q(grapes)-q(${8}) Bias of N.Hemis'
'draw ylab Pressure (unit: hPa) '
'draw xlab Bias'
'print'
'c'
iks=1
while(iks<=ike)
j=1
while(j<=${3})
  k=j
  if(iks=2)
    k=4+j  
  endif
  'set vrange 5 80'
  'set ylevs 1000 925 850 700 500 300 250 200 150 100 70 50 30 20 10'
  'set cthick 8'
  'set ccolor 'cc.j' '
  'set cmark 0'
  'set cstyle 'ltype.iks''
  'set grads off'
  'd ave('rmss.k',t=1,t='time1')'
  j=j+1
endwhile
iks=iks+1
endwhile
if(${9}!='xaxb')
  if(${3}=1)
    './cbar_line_thick.gs -x 5.2 -y 2.8 -c 'cc.1' -m 0 -l 'ltype.1' -t "'t1'-S" '
  endif
  if(${3}=2)
    './cbar_line_thick.gs -x 5.2 -y 2.8 -c 'cc.1' 'cc.2' -m 0 0 -l 'ltype.1' 'ltype.1' -t "'t1'-S" "'t2'-S" '
  endif
  if(${3}=3)
    './cbar_line_thick.gs -x 5.2 -y 2.8 -c 'cc.1' 'cc.2' 'cc.3' -m 0 0 0 -l 'ltype.1' 'ltype.1' 'ltype.1' -t "'t1'-S" "'t2'-S" "'t3'-S"'
  endif
  if(${3}=4)
    './cbar_line_thick.gs -x 5.2 -y 2.8 -c 'cc.1' 'cc.2' 'cc.3' 'cc.4' -m 0 0 0 0 -l 'ltype.1' 'ltype.1' 'ltype.1' 'ltype.1' -t "'t1'-S" "'t2'-S" "'t3'-S" "'t4'-S"'
  endif
else
  if(${3}=1)
    './cbar_line_thick.gs -x 5.2 -y 2.8 -c 'cc.1' 'cc.1'  -m 0 0 -l 'ltype.1' 'ltype.2' -t "'t1'-S" "'t5'-S" '
  endif
  if(${3}=2)
    './cbar_line_thick.gs -x 5.2 -y 2.8 -c 'cc.1' 'cc.2' 'cc.1' 'cc.2' -m 0 0 0 0 -l 'ltype.1' 'ltype.1' 'ltype.2' 'ltype.2' -t "'t1'-S" "'t2'-S" "'t5'-S" "'t6'-S"'
  endif
  if(${3}=3)
    './cbar_line_thick.gs -x 5.2 -y 2.8 -c 'cc.1' 'cc.2' 'cc.3' 'cc.1' 'cc.2' 'cc.3' -m 0 0 0 0 0 0 -l 'ltype.1' 'ltype.1' 'ltype.1' 'ltype.2' 'ltype.2' 'ltype.2' -t "'t1'-S" "'t2'-S" "'t3'-S" "'t5'-S" "'t6'-S" "'t7'-S"'
  endif
  if(${3}=4)
    './cbar_line_thick.gs -x 5.2 -y 2.8 -c 'cc.1' 'cc.2' 'cc.3' 'cc.4' 'cc.1' 'cc.2' 'cc.3' 'cc.4' -m 0 0 0 0 0 0 0 0 -l 'ltype.1' 'ltype.1' 'ltype.1' 'ltype.1' 'ltype.2' 'ltype.2' 'ltype.2' 'ltype.2' -t "'t1'-S" "'t2'-S" "'t3'-S" "'t4'-S" "'t5'-S" "'t6'-S" "'t7'-S" "'t8'-S"'
  endif
  
endif
'draw title time-averaged q(grapes)-q(${8}) RMS of S.Hemis'
'draw ylab Pressure (unit: hPa) '
'draw xlab RMS'
'print'
'c'
iks=1
while(iks<=ike)
j=1
while(j<=${3})
  k=j
  if(iks=2)
    k=4+j  
  endif
  'set vrange -40 20'
  'set ylevs 1000 925 850 700 500 300 250 200 150 100 70 50 30 20 10'
  'set cthick 8'
  'set ccolor 'cc.j' '
  'set cmark 0'
  'set cstyle 'ltype.iks''
  'set grads off'
  'd ave('biass.k',t=1,t='time1')'
  j=j+1
endwhile
iks=iks+1
endwhile
'set cstyle 5'
'set ccolor 1'
'set cmark 0'
'set cthick 8'
'd 'biass.k'-'biass.k''
if(${9}!='xaxb')
  if(${3}=1)
    './cbar_line_thick.gs -x 5.2 -y 2.8 -c 'cc.1' -m 0 -l 'ltype.1' -t "'t1'-S" '
  endif
  if(${3}=2)
    './cbar_line_thick.gs -x 5.2 -y 2.8 -c 'cc.1' 'cc.2' -m 0 0 -l 'ltype.1' 'ltype.1' -t "'t1'-S" "'t2'-S" '
  endif
  if(${3}=3)
    './cbar_line_thick.gs -x 5.2 -y 2.8 -c 'cc.1' 'cc.2' 'cc.3' -m 0 0 0 -l 'ltype.1' 'ltype.1' 'ltype.1' -t "'t1'-S" "'t2'-S" "'t3'-S"'
  endif
  if(${3}=4)
    './cbar_line_thick.gs -x 5.2 -y 2.8 -c 'cc.1' 'cc.2' 'cc.3' 'cc.4' -m 0 0 0 0 -l 'ltype.1' 'ltype.1' 'ltype.1' 'ltype.1' -t "'t1'-S" "'t2'-S" "'t3'-S" "'t4'-S"'
  endif
else
  if(${3}=1)
    './cbar_line_thick.gs -x 5.2 -y 2.8 -c 'cc.1' 'cc.1'  -m 0 0 -l 'ltype.1' 'ltype.2' -t "'t1'-S" "'t5'-S" '
  endif
  if(${3}=2)
    './cbar_line_thick.gs -x 5.2 -y 2.8 -c 'cc.1' 'cc.2' 'cc.1' 'cc.2' -m 0 0 0 0 -l 'ltype.1' 'ltype.1' 'ltype.2' 'ltype.2' -t "'t1'-S" "'t2'-S" "'t5'-S" "'t6'-S"'
  endif
  if(${3}=3)
    './cbar_line_thick.gs -x 5.2 -y 2.8 -c 'cc.1' 'cc.2' 'cc.3' 'cc.1' 'cc.2' 'cc.3' -m 0 0 0 0 0 0 -l 'ltype.1' 'ltype.1' 'ltype.1' 'ltype.2' 'ltype.2' 'ltype.2' -t "'t1'-S" "'t2'-S" "'t3'-S" "'t5'-S" "'t6'-S" "'t7'-S"'
  endif
  if(${3}=4)
    './cbar_line_thick.gs -x 5.2 -y 2.8 -c 'cc.1' 'cc.2' 'cc.3' 'cc.4' 'cc.1' 'cc.2' 'cc.3' 'cc.4' -m 0 0 0 0 0 0 0 0 -l 'ltype.1' 'ltype.1' 'ltype.1' 'ltype.1' 'ltype.2' 'ltype.2' 'ltype.2' 'ltype.2' -t "'t1'-S" "'t2'-S" "'t3'-S" "'t4'-S" "'t5'-S" "'t6'-S" "'t7'-S" "'t8'-S"'
  endif
  
endif
'draw title time-averaged q(grapes)-q(${8}) Bias of S.Hemis'
'draw ylab Pressure (unit: hPa) '
'draw xlab Bias'
'print'
'c'
iks=1
while(iks<=ike)
j=1
while(j<=${3})
  k=j
  if(iks=2)
    k=4+j  
  endif
  'set vrange 0 45'
  'set ylevs 1000 925 850 700 500 300 250 200 150 100 70 50 30 20 10'
  'set cthick 8'
  'set ccolor 'cc.j' '
  'set cmark 0'
  'set cstyle 'ltype.iks''
  'set grads off'
  'd ave('rmseq.k',t=1,t='time1')'
  j=j+1
endwhile
iks=iks+1
endwhile
if(${9}!='xaxb')
  if(${3}=1)
    './cbar_line_thick.gs -x 5.2 -y 2.8 -c 'cc.1' -m 0 -l 'ltype.1' -t "'t1'-tro." '
  endif
  if(${3}=2)
    './cbar_line_thick.gs -x 5.2 -y 2.8 -c 'cc.1' 'cc.2' -m 0 0 -l 'ltype.1' 'ltype.1' -t "'t1'-tro." "'t2'-tro." '
  endif
  if(${3}=3)
    './cbar_line_thick.gs -x 5.2 -y 2.8 -c 'cc.1' 'cc.2' 'cc.3' -m 0 0 0 -l 'ltype.1' 'ltype.1' 'ltype.1' -t "'t1'-tro." "'t2'-tro." "'t3'-tro."'
  endif
  if(${3}=4)
    './cbar_line_thick.gs -x 5.2 -y 2.8 -c 'cc.1' 'cc.2' 'cc.3' 'cc.4' -m 0 0 0 0 -l 'ltype.1' 'ltype.1' 'ltype.1' 'ltype.1' -t "'t1'-tro." "'t2'-tro." "'t3'-tro." "'t4'-tro."'
  endif
else
  if(${3}=1)
    './cbar_line_thick.gs -x 5.2 -y 2.8 -c 'cc.1' 'cc.1'  -m 0 0 -l 'ltype.1' 'ltype.2' -t "'t1'-tro." "'t5'-tro." '
  endif
  if(${3}=2)
    './cbar_line_thick.gs -x 5.2 -y 2.8 -c 'cc.1' 'cc.2' 'cc.1' 'cc.2' -m 0 0 0 0 -l 'ltype.1' 'ltype.1' 'ltype.2' 'ltype.2' -t "'t1'-tro." "'t2'-tro." "'t5'-tro." "'t6'-tro."'
  endif
  if(${3}=3)
    './cbar_line_thick.gs -x 5.2 -y 2.8 -c 'cc.1' 'cc.2' 'cc.3' 'cc.1' 'cc.2' 'cc.3' -m 0 0 0 0 0 0 -l 'ltype.1' 'ltype.1' 'ltype.1' 'ltype.2' 'ltype.2' 'ltype.2' -t "'t1'-tro." "'t2'-tro." "'t3'-tro." "'t5'-tro." "'t6'-tro." "'t7'-tro."'
  endif
  if(${3}=4)
    './cbar_line_thick.gs -x 5.2 -y 2.8 -c 'cc.1' 'cc.2' 'cc.3' 'cc.4' 'cc.1' 'cc.2' 'cc.3' 'cc.4' -m 0 0 0 0 0 0 0 0 -l 'ltype.1' 'ltype.1' 'ltype.1' 'ltype.1' 'ltype.2' 'ltype.2' 'ltype.2' 'ltype.2' -t "'t1'-tro." "'t2'-tro." "'t3'-tro." "'t4'-tro." "'t5'-tro." "'t6'-tro." "'t7'-tro." "'t8'-tro."'
  endif
  
endif
'draw title time-averaged q(grapes)-q(${8}) RMS of Tropics'
'draw ylab Pressure (unit: hPa) '
'draw xlab RMS'
'print'
'c'
iks=1
while(iks<=ike)
j=1
while(j<=${3})
  k=j
  if(iks=2)
    k=4+j  
  endif
  'set vrange -45 40'
  'set ylevs 1000 925 850 700 500 300 250 200 150 100 70 50 30 20 10'
  'set cthick 8'
  'set ccolor 'cc.j' '
  'set cmark 0'
  'set cstyle 'ltype.iks''
  'set grads off'
  'd ave('biaseq.k',t=1,t='time1')'
  j=j+1
endwhile
iks=iks+1
endwhile
'set cstyle 5'
'set ccolor 1'
'set cmark 0'
'set cthick 8'
'd 'biass.k'-'biass.k''
if(${9}!='xaxb')
  if(${3}=1)
    './cbar_line_thick.gs -x 5.2 -y 2.8 -c 'cc.1' -m 0 -l 'ltype.1' -t "'t1'-tro." '
  endif
  if(${3}=2)
    './cbar_line_thick.gs -x 5.2 -y 2.8 -c 'cc.1' 'cc.2' -m 0 0 -l 'ltype.1' 'ltype.1' -t "'t1'-tro." "'t2'-tro." '
  endif
  if(${3}=3)
    './cbar_line_thick.gs -x 5.2 -y 2.8 -c 'cc.1' 'cc.2' 'cc.3' -m 0 0 0 -l 'ltype.1' 'ltype.1' 'ltype.1' -t "'t1'-tro." "'t2'-tro." "'t3'-tro."'
  endif
  if(${3}=4)
    './cbar_line_thick.gs -x 5.2 -y 2.8 -c 'cc.1' 'cc.2' 'cc.3' 'cc.4' -m 0 0 0 0 -l 'ltype.1' 'ltype.1' 'ltype.1' 'ltype.1' -t "'t1'-tro." "'t2'-tro." "'t3'-tro." "'t4'-tro."'
  endif
else
  if(${3}=1)
    './cbar_line_thick.gs -x 5.2 -y 2.8 -c 'cc.1' 'cc.1'  -m 0 0 -l 'ltype.1' 'ltype.2' -t "'t1'-tro." "'t5'-tro." '
  endif
  if(${3}=2)
    './cbar_line_thick.gs -x 5.2 -y 2.8 -c 'cc.1' 'cc.2' 'cc.1' 'cc.2' -m 0 0 0 0 -l 'ltype.1' 'ltype.1' 'ltype.2' 'ltype.2' -t "'t1'-tro." "'t2'-tro." "'t5'-tro." "'t6'-tro."'
  endif
  if(${3}=3)
    './cbar_line_thick.gs -x 5.2 -y 2.8 -c 'cc.1' 'cc.2' 'cc.3' 'cc.1' 'cc.2' 'cc.3' -m 0 0 0 0 0 0 -l 'ltype.1' 'ltype.1' 'ltype.1' 'ltype.2' 'ltype.2' 'ltype.2' -t "'t1'-tro." "'t2'-tro." "'t3'-tro." "'t5'-tro." "'t6'-tro." "'t7'-tro."'
  endif
  if(${3}=4)
    './cbar_line_thick.gs -x 5.2 -y 2.8 -c 'cc.1' 'cc.2' 'cc.3' 'cc.4' 'cc.1' 'cc.2' 'cc.3' 'cc.4' -m 0 0 0 0 0 0 0 0 -l 'ltype.1' 'ltype.1' 'ltype.1' 'ltype.1' 'ltype.2' 'ltype.2' 'ltype.2' 'ltype.2' -t "'t1'-tro." "'t2'-tro." "'t3'-tro." "'t4'-tro." "'t5'-tro." "'t6'-tro." "'t7'-tro." "'t8'-tro."'
  endif
  
endif
'draw title time-averaged q(grapes)-q(${8}) Bias of Tropics'
'draw ylab Pressure (unit: hPa) '
'draw xlab Bias'
'print'
'c'
'disable print'
quit
EOF


ancor.sh_q

17.25 KB, 下载次数: 3, 下载积分: 金钱 -5

pltprof.sh_q

13.03 KB, 下载次数: 2, 下载积分: 金钱 -5

回复

使用道具 举报

xiaocaoqiqiao

新浪微博达人勋
发表于 2014-12-28 13:40:04 | 显示全部楼层
你应该把程序贴出来，想要帮你回答问题还要花10个金钱下载你的文件
密码修改失败请联系微信：mofangbao
回复 支持 反对

使用道具 举报

279995024

新浪微博达人勋
发表于 2014-12-28 14:44:39 | 显示全部楼层
xiaocaoqiqiao 发表于 2014-12-28 13:40
你应该把程序贴出来，想要帮你回答问题还要花10个金钱下载你的文件

好的，谢谢啦，我把程序粘上
密码修改失败请联系微信：mofangbao
回复 支持 反对

使用道具 举报

返回列表 发新帖
高级模式
B Color Image Link Quote Code Smilies
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

小黑屋|帮助中心|清除痕迹

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表，不代表气象家园立场，禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表