急！利用fnl通过grads求位温时出现values<=zero????

7795216223

'reinit'
'set gxout fwrite'
'set fwrite d:\shuju\2010\thse.dat'
'open d:\shuju\2010\data.ctl'
i=1
while(i<=20)
'set lon 50 160'
'set lat 10 70'
'set t 'i''
'set z 1'
'define prs=100000'  
'define qv=SPFHprs'
'define e=prs*qv/(0.62197+qv)+1e-10'
'define tlcl=55.0+2840.0/(3.5*log(TMPprs)-log(e)-4.805)'
'undefine e'
'define theta=TMPprs*pow((100000/prs),(0.2854*(1.0-0.28*qv)))'
'define thse=theta*exp(((3376./tlcl)-2.54)*qv*(1.0+0.81*qv))'
'd thse'
'set z 2'
'define prs=97500'
'define qv=SPFHprs'
'define e=prs*qv/(0.62197+qv)+1e-10'
'define tlcl=55.0+2840.0/(3.5*log(TMPprs)-log(e)-4.775)'
'undefine e'
'define theta=TMPprs*pow((100000/prs),(0.2854*(1.0-0.28*qv)))'
'define thse=theta*exp(((3376./tlcl)-2.54)*qv*(1.0+0.81*qv))'
'd thse'
'set z 3'
'define prs=95000'
'define qv=SPFHprs'
'define e=prs*qv/(0.62197+qv)+1e-10'
'define tlcl=55.0+2840.0/(3.5*log(TMPprs)-log(e)-4.775)'
'undefine e'
'define theta=TMPprs*pow((100000/prs),(0.2854*(1.0-0.28*qv)))'
'define thse=theta*exp(((3376./tlcl)-2.54)*qv*(1.0+0.81*qv))'
'd thse'
'set z 4'
'define prs=92500'  
'define qv=SPFHprs'
'define e=prs*qv/(0.62197+qv)+1e-10'
'define tlcl=55.0+2840.0/(3.5*log(TMPprs)-log(e)-4.775)'
'undefine e'
'define theta=TMPprs*pow((100000/prs),(0.2854*(1.0-0.28*qv)))'
'define thse=theta*exp(((3376./tlcl)-2.54)*qv*(1.0+0.81*qv))'
'd thse'
'set z 5'
'define prs=90000'
'define qv=SPFHprs'
'define e=prs*qv/(0.62197+qv)+1e-10'
'define tlcl=55.0+2840.0/(3.5*log(TMPprs)-log(e)-4.775)'
'undefine e'
'define theta=TMPprs*pow((100000/prs),(0.2854*(1.0-0.28*qv)))'
'define thse=theta*exp(((3376./tlcl)-2.54)*qv*(1.0+0.81*qv))'
'd thse'
'set z 6'
'define prs=87500'  
'define qv=SPFHprs'
'define e=prs*qv/(0.62197+qv)+1e-10'
'define tlcl=55.0+2840.0/(3.5*log(TMPprs)-log(e)-4.775)'
'undefine e'
'define theta=TMPprs*pow((100000/prs),(0.2854*(1.0-0.28*qv)))'
'define thse=theta*exp(((3376./tlcl)-2.54)*qv*(1.0+0.81*qv))'
'd thse'
'set z 7'
'define prs=85000'
'define qv=SPFHprs'
'define e=prs*qv/(0.62197+qv)+1e-10'
'define tlcl=55.0+2840.0/(3.5*log(TMPprs)-log(e)-4.775)'
'undefine e'
'define theta=TMPprs*pow((100000/prs),(0.2854*(1.0-0.28*qv)))'
'define thse=theta*exp(((3376./tlcl)-2.54)*qv*(1.0+0.81*qv))'
'd thse'
'set z 8'
'define prs=82500' '
'define qv=SPFHprs'
'define e=prs*qv/(0.62197+qv)+1e-10'
'define tlcl=55.0+2840.0/(3.5*log(TMPprs)-log(e)-4.775)'
'undefine e'
'define theta=TMPprs*pow((100000/prs),(0.2854*(1.0-0.28*qv)))'
'define thse=theta*exp(((3376./tlcl)-2.54)*qv*(1.0+0.81*qv))'
'd thse'
'set z 9'
'define prs=80000'
'define qv=SPFHprs'
'define e=prs*qv/(0.62197+qv)+1e-10'
'define tlcl=55.0+2840.0/(3.5*log(TMPprs)-log(e)-4.775)'
'undefine e'
'define theta=TMPprs*pow((100000/prs),(0.2854*(1.0-0.28*qv)))'
'define thse=theta*exp(((3376./tlcl)-2.54)*qv*(1.0+0.81*qv))'
'd thse'
'set z 10'
'define prs=77500'
'define qv=SPFHprs'
'define e=prs*qv/(0.62197+qv)+1e-10'
'define tlcl=55.0+2840.0/(3.5*log(TMPprs)-log(e)-4.775)'
'undefine e'
'define theta=TMPprs*pow((100000/prs),(0.2854*(1.0-0.28*qv)))'
'define thse=theta*exp(((3376./tlcl)-2.54)*qv*(1.0+0.81*qv))'
'd thse'
'set z 11'
'define prs=75000'
'define qv=SPFHprs'
'define e=prs*qv/(0.62197+qv)+1e-10'
'define tlcl=55.0+2840.0/(3.5*log(TMPprs)-log(e)-4.775)'
'undefine e'
'define theta=TMPprs*pow((100000/prs),(0.2854*(1.0-0.28*qv)))'
'define thse=theta*exp(((3376./tlcl)-2.54)*qv*(1.0+0.81*qv))'
'd thse'
'set z 12'
'define prs=70000'  
'define qv=SPFHprs'
'define e=prs*qv/(0.62197+qv)+1e-10'
'define tlcl=55.0+2840.0/(3.5*log(TMPprs)-log(e)-4.775)'
'undefine e'
'define theta=TMPprs*pow((100000/prs),(0.2854*(1.0-0.28*qv)))'
'define thse=theta*exp(((3376./tlcl)-2.54)*qv*(1.0+0.81*qv))'
'd thse'
'set z 13'
'define prs=65000''
'define qv=SPFHprs'
'define e=prs*qv/(0.62197+qv)+1e-10'
'define tlcl=55.0+2840.0/(3.5*log(TMPprs)-log(e)-4.775)'
'undefine e'
'define theta=TMPprs*pow((100000/prs),(0.2854*(1.0-0.28*qv)))'
'define thse=theta*exp(((3376./tlcl)-2.54)*qv*(1.0+0.81*qv))'
'd thse'
'set z 14'
'define prs=60000''
'define qv=SPFHprs'
'define e=prs*qv/(0.62197+qv)+1e-10'
'define tlcl=55.0+2840.0/(3.5*log(TMPprs)-log(e)-4.775)'
'undefine e'
'define theta=TMPprs*pow((100000/prs),(0.2854*(1.0-0.28*qv)))'
'define thse=theta*exp(((3376./tlcl)-2.54)*qv*(1.0+0.81*qv))'
'd thse'
'set z 15'
'define prs=55000'
'define qv=SPFHprs'
'define e=prs*qv/(0.62197+qv)+1e-10'
'define tlcl=55.0+2840.0/(3.5*log(TMPprs)-log(e)-4.775)'
'undefine e'
'define theta=TMPprs*pow((100000/prs),(0.2854*(1.0-0.28*qv)))'
'define thse=theta*exp(((3376./tlcl)-2.54)*qv*(1.0+0.81*qv))'
'd thse'
'set z 16'
'define prs=50000'
'define qv=SPFHprs'
'define e=prs*qv/(0.62197+qv)+1e-10'
'define tlcl=55.0+2840.0/(3.5*log(TMPprs)-log(e)-4.775)'
'undefine e'
'define theta=TMPprs*pow((100000/prs),(0.2854*(1.0-0.28*qv)))'
'define thse=theta*exp(((3376./tlcl)-2.54)*qv*(1.0+0.81*qv))'
'd thse'
'set z 17'
'define prs=45000'
'define qv=SPFHprs'
'define e=prs*qv/(0.62197+qv)+1e-10'
'define tlcl=55.0+2840.0/(3.5*log(TMPprs)-log(e)-4.775)'
'undefine e'
'define theta=TMPprs*pow((100000/prs),(0.2854*(1.0-0.28*qv)))'
'define thse=theta*exp(((3376./tlcl)-2.54)*qv*(1.0+0.81*qv))'
'd thse'
'set z 18'
'define prs=40000'
'define qv=SPFHprs'
'define e=prs*qv/(0.62197+qv)+1e-10'
'define tlcl=55.0+2840.0/(3.5*log(TMPprs)-log(e)-4.775)'
'undefine e'
'define theta=TMPprs*pow((100000/prs),(0.2854*(1.0-0.28*qv)))'
'define thse=theta*exp(((3376./tlcl)-2.54)*qv*(1.0+0.81*qv))'
'd thse'
'set z 19'
'define prs=35000'
'define qv=SPFHprs'
'define e=prs*qv/(0.62197+qv)+1e-10'
'define tlcl=55.0+2840.0/(3.5*log(TMPprs)-log(e)-4.775)'
'undefine e'
'define theta=TMPprs*pow((100000/prs),(0.2854*(1.0-0.28*qv)))'
'define thse=theta*exp(((3376./tlcl)-2.54)*qv*(1.0+0.81*qv))'
'd thse'
'set z 20'
'define prs=30000'  
'define qv=SPFHprs'
'define e=prs*qv/(0.62197+qv)+1e-10'
'define tlcl=55.0+2840.0/(3.5*log(TMPprs)-log(e)-4.775)'
'undefine e'
'define theta=TMPprs*pow((100000/prs),(0.2854*(1.0-0.28*qv)))'
'define thse=theta*exp(((3376./tlcl)-2.54)*qv*(1.0+0.81*qv))'
'd thse'
'set z 21'
'define prs=25000' '
'define qv=SPFHprs'
'define e=prs*qv/(0.62197+qv)+1e-10'
'define tlcl=55.0+2840.0/(3.5*log(TMPprs)-log(e)-4.775)'
'undefine e'
'define theta=TMPprs*pow((100000/prs),(0.2854*(1.0-0.28*qv)))'
'define thse=theta*exp(((3376./tlcl)-2.54)*qv*(1.0+0.81*qv))'
'd thse'
'set z 22'
'define prs=22500'  
'define qv=SPFHprs'
'define e=prs*qv/(0.62197+qv)+1e-10'
'define tlcl=55.0+2840.0/(3.5*log(TMPprs)-log(e)-4.775)'
'undefine e'
'define theta=TMPprs*pow((100000/prs),(0.2854*(1.0-0.28*qv)))'
'define thse=theta*exp(((3376./tlcl)-2.54)*qv*(1.0+0.81*qv))'
'd thse'
'set z 23'
'define prs=20000'  
'define qv=SPFHprs'
'define e=prs*qv/(0.62197+qv)+1e-10'
'define tlcl=55.0+2840.0/(3.5*log(TMPprs)-log(e)-4.775)'
'undefine e'
'define theta=TMPprs*pow((100000/prs),(0.2854*(1.0-0.28*qv)))'
'define thse=theta*exp(((3376./tlcl)-2.54)*qv*(1.0+0.81*qv))'
'd thse'
'set z 24'
'define prs=17500'
'define qv=SPFHprs'
'define e=prs*qv/(0.62197+qv)+1e-10'
'define tlcl=55.0+2840.0/(3.5*log(TMPprs)-log(e)-4.775)'
'undefine e'
'define theta=TMPprs*pow((100000/prs),(0.2854*(1.0-0.28*qv)))'
'define thse=theta*exp(((3376./tlcl)-2.54)*qv*(1.0+0.81*qv))'
'd thse'
'set z 25'
'define prs=15000'  
'define qv=SPFHprs'
'define e=prs*qv/(0.62197+qv)+1e-10'
'define tlcl=55.0+2840.0/(3.5*log(TMPprs)-log(e)-4.775)'
'undefine e'
'define theta=TMPprs*pow((100000/prs),(0.2854*(1.0-0.28*qv)))'
'define thse=theta*exp(((3376./tlcl)-2.54)*qv*(1.0+0.81*qv))'
'd thse'
'set z 26'
'define prs=12500'
'define qv=SPFHprs'
'define e=prs*qv/(0.62197+qv)+1e-10'
'define tlcl=55.0+2840.0/(3.5*log(TMPprs)-log(e)-4.775)'
'undefine e'
'define theta=TMPprs*pow((100000/prs),(0.2854*(1.0-0.28*qv)))'
'define thse=theta*exp(((3376./tlcl)-2.54)*qv*(1.0+0.81*qv))'
'd thse'
'set z 27'
'define prs=10000'
'define qv=SPFHprs'
'define e=prs*qv/(0.62197+qv)+1e-10'
'define tlcl=55.0+2840.0/(3.5*log(TMPprs)-log(e)-4.775)'
'undefine e'
'define theta=TMPprs*pow((100000/prs),(0.2854*(1.0-0.28*qv)))'
'define thse=theta*exp(((3376./tlcl)-2.54)*qv*(1.0+0.81*qv))'
'd thse'
   j=1
while(j<=27)
   'set lon 50 160'
   'set lat 10 70'
   'set t 'i''
   'set z 'j''
   'd ugrdprs'
    j=j+1
endwhile
    j=1
while(j<=27)
   'set lon 50 160'
   'set lat 10 70'
   'set t 'i''
   'set z 'j''
   'd vgrdprs'
    j=j+1
endwhile
i=i+1
endwhile
'disable fwrite'
'reinit'

lqouc

http://bbs.06climate.com/forum.php?mod=viewthread&tid=4571
去看看吧。