求质量流函数并画图

kba541

'reinit'
'sdfopen E:\newdata\vwnd.mon.mean.nc'
'run E:\newdata\season_mean(summer)_mass_Stream_function_Hadley\18colors.gs'
lev.1=1000 ; lev.2=925 ; lev.3=850 ; lev.4=700  ; lev.5=600  ; lev.6=500
lev.7=400  ; lev.8=300 ; lev.9=250 ; lev.10=200 ; lev.11=150 ; lev.12=100
lev.13=70  ; lev.14=50 ; lev.15=30 ; lev.16=20  ; lev.17=10
  'set t 846 848'
   'set x 1'   
    'set y 1 73'
    ceng=1
    while(ceng<=16)
          'set z ' ceng
        temp= ceng +1
        dp.ceng =(lev.temp-lev.ceng )*100
    ceng=ceng+1
  endwhile
  dp.17=(lev.17-lev.16)*100   

*time=1
*while(time<=848)
  *  'set t' time
     *ceng=1
  *  while(ceng<=17)   
        'set t 846 848'
        'set x 1'   
        'set y 1 73'
        'set z  1 17'
        'define vba=ave(vwnd.1,x=1,x=144,-b)'                           
       * ceng=ceng+1
    *endwhile

    ceng=17
    'set z'ceng
     'define psi'ceng'=vba*0'

       ceng=16
    while(ceng>=1)   
        'set x 1'  
        'set t 846 848'
        'set y 1 73'
        'set z'ceng
        temp1= ceng +1
          'define aa=2*3.141592654*6370000*cos(lat*3.141592654/180)/9.8'  
          'define psi=psi'temp1'-aa*vba*'dp.ceng''
     ceng=ceng-1
    endwhile

    *ceng = 1
   * while(ceng<=17)  
        *'d psi'ceng
     *endwhile

        *time=time+1
*endwhile


        'set x 1'   
        'set t 1'
        'set y 1 73'
        'set z 1 17'
        'define psijun=ave(psi,t=846,t=848,1)'                  

  'set font 4'
'set parea 1.3 10 1.2 7.8'
'set t 1'
'set lat -90 90'
'set lon 0 '
'set z 1 13'
'set csmooth on'
'set cterp on'
'set gxout shaded'
'set grads off'
'set grid off'
'set cint 2'
'set xlopts 1 5 0.2'
'set ylopts 1 5 0.2'   
'set zlog on'
'run E:\newdata\season_mean(summer)_mass_Stream_function_Hadley\rainbow.gs'
'd psijun/1e10'

'set gxout contour'
'set grads off'
'set grid off'
'set clab masked'
'set ccolor 1'   
'set cthick 6'
'set csmooth on'
'set cterp on'
'set xlpos -10'   
'set ylpos -10'
'set cint 2'
'set clopts 1 5 0.15'
'set zlog on'
'd psijun/1e10'

'draw xlab Latitude(`3.`1)'
'draw ylab Pressure(hpa)'
'draw title Zonal mean meridional mass streamfunction'

'cbar_interp 1 1 0'
'gxprint  E:\newdata\season_mean(summer)_mass_Stream_function_Hadley\season_mean(summer)_mass_Stream_function_Hadley.png white'
'gxprint E:\newdata\season_mean(summer)_mass_Stream_function_Hadley\season_mean(summer)_mass_Stream_function_Hadley.eps white'
;