关于grads批处理的问题

Sue.Chang

我在批处理的时候所产生的最后结果的DAT的内容都是一样的
运行整个gs的时候并没有报错
这个是我的ctl之一 其他的类似

1. 
   
2. <P> dset H:\q11\%y4.dat
   
3. options template
   
4. undef -9.99E+08
   
5. title mean daily NMC reanalysis<1983>
   
6. xdef  144 linear    0  2.5
   
7. ydef   73 linear  -90  2.5
   
8. zdef   12 levels 1000 925 850 700 600 500 400 300 250 200 150 100
   
9. tdef  23376 linear 00Z01JAN1948 1dy
   
10. vars    1
    
11. q11     12  99 ** local change                :J/s.kg
    
12. endvars</P>

复制代码

这个是我的gs

1. 
   
2. <P>'reinit'
   
3. 'OPEN E:\q1\rsj\q11.ctl'
   
4. 'OPEN E:\q1\rsj\q12.ctl'
   
5. 'OPEN E:\q1\rsj\q13.ctl'
   
6. 'open E:\q1\rsj\pres.CTL'
   
7. yy=1948
   
8. j=0
   
9. id=1
   
10. while(yy<=2011)
    
11. 'set gxout fwrite'
    
12. 'set fwrite H:\q1\'yy'.grd'
    
13. 'set x 1 144'
    
14. 'set y 1 73'
    
15. if ((math_mod(yy,4)=0&math_mod(yy,100)!=0)|math_mod(yy,400)=0)</P>
    
16. <P>  while(id<=366+j)
    
17.     'set t 'id''</P>
    
18. <P>      'd vint(pres.4/100,q11.1+q12.2+(q13.3+q14.3+q15.3)/3,100)/9.8'</P>
    
19. <P>        id=id+1
    
20.   endwhile
    
21.   j=j+366
    
22. else
    
23.      
    
24.   while(id<=365+j)
    
25.     'set t 'id''</P>
    
26. <P>      'd vint(pres.4/100,q11.1+q12.2+(q13.3+q14.3+q15.3)/3,100)/9.8'</P>
    
27. <P>        id=id+1
    
28.         endwhile
    
29.         j=j+365
    
30.         endif
    
31.         'close 1'
    
32.         yy=yy+1
    
33. 'disable fwrite'
    
34. endwhile
    
35. 'print'
    
36. 'reinit'
    
37. ;</P>
    
38. <P> </P>
    
39. <P> </P>
    
40. <P> </P>

复制代码

希望大神们抽空看看 帮我解决一下