对多时次大区域（多省份）maskout的小小改进

mofangbao · 发表于 2013-4-17 23:24:50

平流层的萝卜发表于 2013-4-17 22:45

登录/注册后可看大图

关于传说提出的修改mask描述ctl中的时间设置，我还是不懂。因为mask文件的大小对应的就是一个时次的，如果 ...

如果对于每一个时次maskout要屏蔽的部分都是一样的话，加一个括号指定t=1就行了，比如1,10个时次都是要屏蔽中国以外的区域，你 define a=maskout(var,g.2(t=1)) 类似这样就可以了。
如果每一个时次需要屏蔽的内容不一样，就需要每一个时次配一个了，因为maskout的本质是将第二个参数为0时候的值设置为缺测。
另外，fortran生成时候插值的道理也一样，oacres函数本质上用的是你grid文件的格点信息，如果你每个时次要插值的格点是相同的，那么指定一个时间即可，无需每个时次都配，如果不同的话，也是要另外处理的。
MeteoInfo制作的应该是一个时次的，你想变为多时次的，还需要自己生成。

平流层的萝卜 · 发表于 2013-4-17 23:45:52

mofangbao 发表于 2013-4-17 23:24

登录/注册后可看大图

如果对于每一个时次maskout要屏蔽的部分都是一样的话，加一个括号指定t=1就行了，比如1,10个时次都是要屏 ...

嗯，说的对，谢啦

刘振 · 发表于 2013-4-20 14:20:41

这样的话请问mask的ctl需要跟第一个时次的资料时间一致吗，还有个问题就是，mask资料设置的格点必须和数据资料额格点设置相匹配吗，我画区域平均总是提醒有错误。

平流层的萝卜 · 发表于 2013-4-20 15:16:21

刘振发表于 2013-4-20 14:20

登录/注册后可看大图

这样的话请问mask的ctl需要跟第一个时次的资料时间一致吗，还有个问题就是，mask资料设置的格点必须和数据资 ...

第一个时间我设置的是一致的，格点不涉及到一致的问题吧，因为maskout底图是直接生成的，站点资料是通过程序筛选的，所说的一致是什么意思呢？

刘振 · 发表于 2013-4-20 15:45:05

平流层的萝卜发表于 2013-4-20 15:16

登录/注册后可看大图

第一个时间我设置的是一致的，格点不涉及到一致的问题吧，因为maskout底图是直接生成的，站点资料是通过程 ...

谢谢你的回复，我说的格点的意思是区域资料格点的经纬度要和地图ctl设置的格点的经纬度一样，我做的是只要陆地的降水区域平均，一致的话才可以判断每个格点是否是水体的坐标和资料的坐标联系起来，不然就可能有交错，试了下好像是这样。额，我描述的自己都感觉不清楚，你有联系方式吗，我自己都断定不了求的平均到底对不对，我q286909655，谢谢啦

中华L · 发表于 2013-4-22 20:58:36

感谢楼主分享！学习了！

刘振 · 发表于 2013-4-23 09:29:16

mofangbao 发表于 2013-4-17 23:24

登录/注册后可看大图

如果对于每一个时次maskout要屏蔽的部分都是一样的话，加一个括号指定t=1就行了，比如1,10个时次都是要屏 ...

想问下管理员，如果是三维的呢，我可以更改maskout中的层数吗，为了计算下垫面不同的三维结构，因为每次屏蔽的区域还是一样的。

平流层的萝卜 · 发表于 2013-4-23 10:03:09

刘振发表于 2013-4-20 15:45

登录/注册后可看大图

谢谢你的回复，我说的格点的意思是区域资料格点的经纬度要和地图ctl设置的格点的经纬度一样，我做的是只要 ...

额为什么这条留言没看到呢？

平流层的萝卜 · 发表于 2013-4-23 10:04:15

刘振发表于 2013-4-23 09:29

登录/注册后可看大图

想问下管理员，如果是三维的呢，我可以更改maskout中的层数吗，为了计算下垫面不同的三维结构，因为每次屏 ...

你看ctl里是zdef和tdef都有，那不妨mask（t=1,z=1)试下，不知道可行不？

mofangbao · 发表于 2013-4-23 10:05:02

刘振发表于 2013-4-23 09:29

登录/注册后可看大图

想问下管理员，如果是三维的呢，我可以更改maskout中的层数吗，为了计算下垫面不同的三维结构，因为每次屏 ...

你试试就知道了

		自动登录	找回密码
密码			立即注册

对多时次大区域（多省份）maskout的小小改进

RE: 对多时次大区域（多省份）maskout的小小改进