关于查看domain的NCL问题

FrankieLJY

在网上找到一个查看domain的ncl的脚本，算是作为WRFDomainWizard画图出来的补充
但是由于本人对NCL的语法不是很了解，经常在运行时报如下错误，不知道是不是WRFDomainWizard中的参数哪里设置不对：
GKS ERROR NUMBER   51 ISSUED FROM SUBROUTINE GSVP  :
--RECTANGLE DEFINITION IS INVALID
GKS ERROR NUMBER   51 ISSUED FROM SUBROUTINE GSVP  :
--RECTANGLE DEFINITION IS INVALID
GKS ERROR NUMBER   51 ISSUED FROM SUBROUTINE GSVP  :
--RECTANGLE DEFINITION IS INVALID
warning:TickMarkSetValues: Zero span between tmXBDataLeftF and tmXBDataRightF, turning bottom tickmarks off:[errno=1104]
warning:TickMarkSetValues: Zero span between tmXTDataLeftF and tmXTDataRightF, turning top tickmarks off:[errno=1104]
warning:TickMarkSetValues: Zero span between tmYLDataBottomF and tmYLDataTopF, turning left tickmarks off:[errno=1104]
warning:TickMarkSetValues: Zero span between tmYRDataBottomF and tmYRDataTopF, turning right tickmarks off:[errno=1104]
GKS ERROR NUMBER   51 ISSUED FROM SUBROUTINE GSVP  :
--RECTANGLE DEFINITION IS INVALID
GKS ERROR NUMBER   51 ISSUED FROM SUBROUTINE GSVP  :
--RECTANGLE DEFINITION IS INVALID
warning:TickMarkSetValues: Zero span between tmXBDataLeftF and tmXBDataRightF, turning bottom tickmarks off:[errno=1104]
warning:TickMarkSetValues: Zero span between tmXTDataLeftF and tmXTDataRightF, turning top tickmarks off:[errno=1104]
warning:TickMarkSetValues: Zero span between tmYLDataBottomF and tmYLDataTopF, turning left tickmarks off:[errno=1104]
warning:TickMarkSetValues: Zero span between tmYRDataBottomF and tmYRDataTopF, turning right tickmarks off:[errno=1104]
GKS ERROR NUMBER   51 ISSUED FROM SUBROUTINE GSVP  :
--RECTANGLE DEFINITION IS INVALID
GKS ERROR NUMBER   51 ISSUED FROM SUBROUTINE GSVP  :
--RECTANGLE DEFINITION IS INVALID
GKS ERROR NUMBER   51 ISSUED FROM SUBROUTINE GSVP  :
--RECTANGLE DEFINITION IS INVALID
GKS ERROR NUMBER   51 ISSUED FROM SUBROUTINE GSVP  :
--RECTANGLE DEFINITION IS INVALID
GKS ERROR NUMBER   51 ISSUED FROM SUBROUTINE GSVP  :
--RECTANGLE DEFINITION IS INVALID
我的namelist的数据情况如下：
parent_id         =   1,       1,      2,3,
parent_grid_ratio =   1,       3,      3,3,
i_parent_start    =   1,       6,      16,27
j_parent_start    =   1,       9,     42,12
s_we              =   1,       1,      1,1
s_sn              =   1,       1,      1,1
e_we              =   100,     265,    700,283
e_sn              =   103,     259,    574,184
geog_data_res     = '10m', '2m', '30s','30s'
dx = 69201,
dy = 69201,
map_proj = 'lambert',
ref_lat   =  35.072,
ref_lon   =  100.15,
truelat1  =  35.072,
truelat2  =  100.15,
stand_lon =  100.15,

同时把我的ncl文件也附上 希望贴吧的朋友能帮我看看什么问题！

FrankieLJY

我了个去 我并没有设置下载要扣金钱的啊，这个NCL太长了 不太好粘贴出来啊

FrankieLJY

此问题已经解决，是自己太粗心，将truelat2的值设置成了经度值

华华工工石田立

楼主的ncl脚本收了，谢谢楼主