GrADS同时打开两个ctl，拥有一样的变量名，写命令的时候如何区分变量从哪个文件来

zyzhi

278803532 发表于 2014-6-27 17:28
不知道你那个是不是从文件复制下来了，主要中文逗号和英文逗号的区别

没有，软件中我是英文逗号，我还怀疑过是内存不够，重启了电脑不开任何程序直接grads，还是提示一样的错误，画不出speed2的ave，dell的工作站4G内存应该不是内存不够的问题。还有其他的办法算跨文件的平均吗？比如合并文件？如何合并和时次都是1-124合并的时候如何处理？非常感谢不厌其烦的回答我的问题！

278803532

你打开CTL看看，应该是范围不一样的，我看你的那个错误，可能不是在GS上，有可能是因为CTL的范围不一样，或者数据的范围比CTL的要小

zyzhi

278803532 发表于 2014-6-27 18:35
你打开CTL看看，应该是范围不一样的，我看你的那个错误，可能不是在GS上，有可能是因为CTL的范围不一样，或 ...

同时打开两个ctl截了个图，可以看到除了时间和时次（因为1月31天，4月30天）不一样，其他都是一样的

zyzhi

ctl和dat文件都是wrf跑出来后ARWpost转的，应该不存在数据和ctl范围不一致的问题

278803532

那你把你的CTL什么的都发上来，GS，什么全部发出来给大神看吧。。。

苏涛

数据有问题吗，（u.1 ；v.1） （u.2； v.2）  第一 第二个文件中的风场是没有问题的，是不是你数据维数或者范围与平均函数冲突

风儿飘飘

以变量名.1与变量名.2这样应用区分来自不同文件的变量名

chanch0103

實用的經驗分享

limoxi

不知楼主的问题解决了没，我用两个ctl文件求距平的时候显示如下图的错误，不知道哪里出了问题

1. 'reinit''open f:\Rain2014\Rain2014day.ctl''open f:\RainData\Rain1948-2012mon.ctl''set grads off''set grid off'*'define_colors''set lon 73 135''set lat 18 55''set lev 500''set mpdset cnworld''draw map''set dfile1''pre=0't=152while (t<=182)'set t 't'pre=pre+r.1't=t+1endwhile'set lon 73 135''set lat 18 55''set z 1''set dfile2''mean6=0't=390while (t<=738)'set t 't'mean6=mean6+R.2't=t+1endwhile'set gxout contour''define a=(pre-mean6)''cnbasemap a''draw title sum precipitation in summer 2014'*'gxprint d:\shuju\tu\sumrain.png white'*'reinit';

'reinit''open f:\Rain2014\Rain2014day.ctl''open f:\RainData\Rain1948-2012mon.ctl''set grads off''set grid off'*'define_colors''set lon 73 135''set lat 18 55''set lev 500''set mpdset cnworld''draw map''set dfile1''pre=0't=152while (t<=182)'set t 't'pre=pre+r.1't=t+1endwhile'set lon 73 135''set lat 18 55''set z 1''set dfile2''mean6=0't=390while (t<=738)'set t 't'mean6=mean6+R.2't=t+1endwhile'set gxout contour''define a=(pre-mean6)''cnbasemap a''draw title sum precipitation in summer 2014'*'gxprint d:\shuju\tu\sumrain.png white'*'reinit';

limoxi

不好意思啊，我重新发一下gs文件

1. 'reinit'
   
2. 'open f:\Rain2014\Rain2014day.ctl'
   
3. 'open f:\RainData\Rain1948-2012mon.ctl'
   
4. 'set grads off'
   
5. 'set grid off'
   
6. *'define_colors'
   
7. 'set lon 73 135'
   
8. 'set lat 18 55'
   
9. 'set lev 500'
   
10. 'set mpdset cnworld'
    
11. 'draw map'
    
12. 'set dfile1'
    
13. 'pre=0'
    
14. t=152
    
15. while (t<=182)
    
16. 'set t 't
    
17. 'pre=pre+r.1'
    
18. t=t+1
    
19. endwhile
    
20. 'set lon 73 135'
    
21. 'set lat 18 55'
    
22. 'set z 1'
    
23. 'set dfile2'
    
24. 'mean6=0'
    
25. t=390
    
26. while (t<=738)
    
27. 'set t 't
    
28. 'mean6=mean6+R.2'
    
29. t=t+1
    
30. endwhile
    
31. 'set gxout contour'
    
32. 'define a=(pre-mean6)'
    
33. 'cnbasemap a'
    
34. 'draw title sum precipitation in summer 2014'
    
35. *'gxprint d:\shuju\tu\sumrain.png white'
    
36. *'reinit'
    
37. ;

复制代码