爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 9049|回复: 11

[求助] 用fortran计算得到的数据所画的图形中不该出现的斜线

[复制链接]

新浪微博达人勋

发表于 2012-11-16 22:03:18 | 显示全部楼层 |阅读模式

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册 新浪微博登陆

x
我用arwpost后处理的wrf数据作为初始输入数据,用fortran计算散度风,用所得到的结果画图,发现在500hpa以上图形是正常的,但是850hpa的图形中有一些多余的线,原始数据的缺测值是1e30,因为fortran中进行迭代的数据是原始数据的0.1倍,所以程序中设置的缺测值是1e29,请大家帮我分析一下出现这种现象的原因,是因为缺测值的问题还是ctl文件有误?非常感谢!
图中有一些斜线
描述文件如下
  1. dset e:\div-wind\850\vwx.grd
  2. undef 1.e29
  3. title  OUTPUT FROM WRF V3.3 MODEL
  4. xdef  259 levels
  5.    100.9291   
  6.    101.0459   
  7.    101.1628   
  8.    101.2796   
  9.    101.3964   
  10.    101.5133   
  11.    101.6301   
  12.    101.7469   
  13.    101.8637   
  14.    101.9806   
  15.    102.0974   
  16.    102.2142   
  17.    102.3311   
  18.    102.4479   
  19.    102.5647   
  20.    102.6815   
  21.    102.7984   
  22.    102.9152   
  23.    103.0320   
  24.    103.1489   
  25.    103.2657   
  26.    103.3825   
  27.    103.4993   
  28.    103.6162   
  29.    103.7330   
  30.    103.8498   
  31.    103.9667   
  32.    104.0835   
  33.    104.2003   
  34.    104.3171   
  35.    104.4340   
  36.    104.5508   
  37.    104.6676   
  38.    104.7845   
  39.    104.9013   
  40.    105.0181   
  41.    105.1349   
  42.    105.2518   
  43.    105.3686   
  44.    105.4854   
  45.    105.6023   
  46.    105.7191   
  47.    105.8359   
  48.    105.9527   
  49.    106.0696   
  50.    106.1864   
  51.    106.3032   
  52.    106.4201   
  53.    106.5369   
  54.    106.6537   
  55.    106.7705   
  56.    106.8874   
  57.    107.0042   
  58.    107.1210   
  59.    107.2379   
  60.    107.3547   
  61.    107.4715   
  62.    107.5883   
  63.    107.7052   
  64.    107.8220   
  65.    107.9388   
  66.    108.0557   
  67.    108.1725   
  68.    108.2893   
  69.    108.4061   
  70.    108.5230   
  71.    108.6398   
  72.    108.7566   
  73.    108.8735   
  74.    108.9903   
  75.    109.1071   
  76.    109.2239   
  77.    109.3408   
  78.    109.4576   
  79.    109.5744   
  80.    109.6913   
  81.    109.8081   
  82.    109.9249   
  83.    110.0417   
  84.    110.1586   
  85.    110.2754   
  86.    110.3922   
  87.    110.5091   
  88.    110.6259   
  89.    110.7427   
  90.    110.8595   
  91.    110.9764   
  92.    111.0932   
  93.    111.2100   
  94.    111.3269   
  95.    111.4437   
  96.    111.5605   
  97.    111.6773   
  98.    111.7942   
  99.    111.9110   
  100.    112.0278   
  101.    112.1447   
  102.    112.2615   
  103.    112.3783   
  104.    112.4951   
  105.    112.6120   
  106.    112.7288   
  107.    112.8456   
  108.    112.9625   
  109.    113.0793   
  110.    113.1961   
  111.    113.3129   
  112.    113.4298   
  113.    113.5466   
  114.    113.6634   
  115.    113.7803   
  116.    113.8971   
  117.    114.0139   
  118.    114.1307   
  119.    114.2476   
  120.    114.3644   
  121.    114.4812   
  122.    114.5981   
  123.    114.7149   
  124.    114.8317   
  125.    114.9485   
  126.    115.0654   
  127.    115.1822   
  128.    115.2990   
  129.    115.4159   
  130.    115.5327   
  131.    115.6495   
  132.    115.7663   
  133.    115.8832   
  134.    116.0000   
  135.    116.1168   
  136.    116.2337   
  137.    116.3505   
  138.    116.4673   
  139.    116.5841   
  140.    116.7010   
  141.    116.8178   
  142.    116.9346   
  143.    117.0515   
  144.    117.1683   
  145.    117.2851   
  146.    117.4019   
  147.    117.5188   
  148.    117.6356   
  149.    117.7524   
  150.    117.8693   
  151.    117.9861   
  152.    118.1029   
  153.    118.2197   
  154.    118.3366   
  155.    118.4534   
  156.    118.5702   
  157.    118.6871   
  158.    118.8039   
  159.    118.9207   
  160.    119.0375   
  161.    119.1544   
  162.    119.2712   
  163.    119.3880   
  164.    119.5049   
  165.    119.6217   
  166.    119.7385   
  167.    119.8553   
  168.    119.9722   
  169.    120.0890   
  170.    120.2058   
  171.    120.3227   
  172.    120.4395   
  173.    120.5563   
  174.    120.6731   
  175.    120.7900   
  176.    120.9068   
  177.    121.0236   
  178.    121.1405   
  179.    121.2573   
  180.    121.3741   
  181.    121.4909   
  182.    121.6078   
  183.    121.7246   
  184.    121.8414   
  185.    121.9583   
  186.    122.0751   
  187.    122.1919   
  188.    122.3087   
  189.    122.4256   
  190.    122.5424   
  191.    122.6592   
  192.    122.7761   
  193.    122.8929   
  194.    123.0097   
  195.    123.1265   
  196.    123.2434   
  197.    123.3602   
  198.    123.4770   
  199.    123.5939   
  200.    123.7107   
  201.    123.8275   
  202.    123.9443   
  203.    124.0612   
  204.    124.1780   
  205.    124.2948   
  206.    124.4117   
  207.    124.5285   
  208.    124.6453   
  209.    124.7621   
  210.    124.8790   
  211.    124.9958   
  212.    125.1126   
  213.    125.2295   
  214.    125.3463   
  215.    125.4631   
  216.    125.5799   
  217.    125.6968   
  218.    125.8136   
  219.    125.9304   
  220.    126.0473   
  221.    126.1641   
  222.    126.2809   
  223.    126.3977   
  224.    126.5146   
  225.    126.6314   
  226.    126.7482   
  227.    126.8651   
  228.    126.9819   
  229.    127.0987   
  230.    127.2155   
  231.    127.3324   
  232.    127.4492   
  233.    127.5660   
  234.    127.6829   
  235.    127.7997   
  236.    127.9165   
  237.    128.0333   
  238.    128.1502   
  239.    128.2670   
  240.    128.3838   
  241.    128.5007   
  242.    128.6175   
  243.    128.7343   
  244.    128.8512   
  245.    128.9680   
  246.    129.0848   
  247.    129.2016   
  248.    129.3185   
  249.    129.4353   
  250.    129.5521   
  251.    129.6689   
  252.    129.7858   
  253.    129.9026   
  254.    130.0194   
  255.    130.1363   
  256.    130.2531   
  257.    130.3699   
  258.    130.4867   
  259.    130.6036   
  260.    130.7204   
  261.    130.8372   
  262.    130.9541   
  263.    131.0709   
  264. ydef  254 levels   
  265.    7.758774   
  266.    7.874512   
  267.    7.990227   
  268.    8.105911   
  269.    8.221542   
  270.    8.337158   
  271.    8.452728   
  272.    8.568275   
  273.    8.683784   
  274.    8.799255   
  275.    8.914688   
  276.    9.030090   
  277.    9.145447   
  278.    9.260780   
  279.    9.376060   
  280.    9.491318   
  281.    9.606522   
  282.    9.721695   
  283.    9.836823   
  284.    9.951920   
  285.    10.06696   
  286.    10.18197   
  287.    10.29694   
  288.    10.41187   
  289.    10.52676   
  290.    10.64159   
  291.    10.75639   
  292.    10.87115   
  293.    10.98585   
  294.    11.10052   
  295.    11.21514   
  296.    11.32972   
  297.    11.44424   
  298.    11.55872   
  299.    11.67316   
  300.    11.78755   
  301.    11.90189   
  302.    12.01618   
  303.    12.13043   
  304.    12.24463   
  305.    12.35877   
  306.    12.47287   
  307.    12.58691   
  308.    12.70091   
  309.    12.81485   
  310.    12.92875   
  311.    13.04259   
  312.    13.15637   
  313.    13.27011   
  314.    13.38380   
  315.    13.49742   
  316.    13.61099   
  317.    13.72452   
  318.    13.83798   
  319.    13.95139   
  320.    14.06474   
  321.    14.17804   
  322.    14.29129   
  323.    14.40447   
  324.    14.51760   
  325.    14.63067   
  326.    14.74368   
  327.    14.85664   
  328.    14.96952   
  329.    15.08236   
  330.    15.19513   
  331.    15.30785   
  332.    15.42049   
  333.    15.53309   
  334.    15.64561   
  335.    15.75809   
  336.    15.87049   
  337.    15.98284   
  338.    16.09512   
  339.    16.20734   
  340.    16.31949   
  341.    16.43159   
  342.    16.54360   
  343.    16.65557   
  344.    16.76746   
  345.    16.87930   
  346.    16.99105   
  347.    17.10275   
  348.    17.21438   
  349.    17.32594   
  350.    17.43743   
  351.    17.54886   
  352.    17.66022   
  353.    17.77150   
  354.    17.88272   
  355.    17.99387   
  356.    18.10495   
  357.    18.21596   
  358.    18.32690   
  359.    18.43777   
  360.    18.54855   
  361.    18.65929   
  362.    18.76994   
  363.    18.88052   
  364.    18.99102   
  365.    19.10146   
  366.    19.21181   
  367.    19.32210   
  368.    19.43230   
  369.    19.54244   
  370.    19.65250   
  371.    19.76250   
  372.    19.87241   
  373.    19.98223   
  374.    20.09200   
  375.    20.20168   
  376.    20.31128   
  377.    20.42080   
  378.    20.53025   
  379.    20.63962   
  380.    20.74891   
  381.    20.85812   
  382.    20.96725   
  383.    21.07631   
  384.    21.18528   
  385.    21.29417   
  386.    21.40298   
  387.    21.51171   
  388.    21.62037   
  389.    21.72894   
  390.    21.83743   
  391.    21.94582   
  392.    22.05415   
  393.    22.16238   
  394.    22.27054   
  395.    22.37862   
  396.    22.48660   
  397.    22.59450   
  398.    22.70232   
  399.    22.81006   
  400.    22.91771   
  401.    23.02527   
  402.    23.13274   
  403.    23.24014   
  404.    23.34745   
  405.    23.45466   
  406.    23.56179   
  407.    23.66884   
  408.    23.77580   
  409.    23.88267   
  410.    23.98945   
  411.    24.09614   
  412.    24.20274   
  413.    24.30926   
  414.    24.41568   
  415.    24.52202   
  416.    24.62827   
  417.    24.73443   
  418.    24.84048   
  419.    24.94646   
  420.    25.05235   
  421.    25.15813   
  422.    25.26383   
  423.    25.36945   
  424.    25.47495   
  425.    25.58038   
  426.    25.68572   
  427.    25.79096   
  428.    25.89610   
  429.    26.00115   
  430.    26.10611   
  431.    26.21096   
  432.    26.31573   
  433.    26.42042   
  434.    26.52499   
  435.    26.62948   
  436.    26.73387   
  437.    26.83815   
  438.    26.94235   
  439.    27.04645   
  440.    27.15046   
  441.    27.25436   
  442.    27.35816   
  443.    27.46188   
  444.    27.56550   
  445.    27.66901   
  446.    27.77244   
  447.    27.87576   
  448.    27.97899   
  449.    28.08211   
  450.    28.18513   
  451.    28.28805   
  452.    28.39088   
  453.    28.49361   
  454.    28.59624   
  455.    28.69876   
  456.    28.80119   
  457.    28.90352   
  458.    29.00574   
  459.    29.10786   
  460.    29.20989   
  461.    29.31181   
  462.    29.41363   
  463.    29.51535   
  464.    29.61696   
  465.    29.71848   
  466.    29.81989   
  467.    29.92120   
  468.    30.02240   
  469.    30.12350   
  470.    30.22450   
  471.    30.32540   
  472.    30.42619   
  473.    30.52688   
  474.    30.62746   
  475.    30.72794   
  476.    30.82831   
  477.    30.92857   
  478.    31.02875   
  479.    31.12881   
  480.    31.22876   
  481.    31.32861   
  482.    31.42835   
  483.    31.52798   
  484.    31.62751   
  485.    31.72694   
  486.    31.82626   
  487.    31.92545   
  488.    32.02456   
  489.    32.12356   
  490.    32.22246   
  491.    32.32124   
  492.    32.41991   
  493.    32.51846   
  494.    32.61693   
  495.    32.71529   
  496.    32.81351   
  497.    32.91165   
  498.    33.00968   
  499.    33.10760   
  500.    33.20541   
  501.    33.30310   
  502.    33.40068   
  503.    33.49816   
  504.    33.59553   
  505.    33.69279   
  506.    33.78995   
  507.    33.88698   
  508.    33.98391   
  509.    34.08073   
  510.    34.17744   
  511.    34.27402   
  512.    34.37051   
  513.    34.46690   
  514.    34.56316   
  515.    34.65931   
  516.    34.75535   
  517.    34.85129   
  518.    34.94711        
  519. zdef   1 levels  850.00000
  520. tdef   33 linear 06Z18JUN2009     180MN      
  521. VARS    1
  522. U   1  0  u divergence wind
  523. ENDVARS

vwx.ctl (8.76 KB, 下载次数: 12)
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2012-11-16 22:11:41 | 显示全部楼层
850hpa的散度风如下附件,有很多斜线

div-wind.gmf

5.63 MB, 下载次数: 110, 下载积分: 金钱 -5

850hpa散度风

密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2012-11-16 22:55:03 | 显示全部楼层
直接贴出来比较好吧,别人下载还得花钱呢
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2012-11-17 08:32:18 | 显示全部楼层
就是!!!
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2012-11-17 09:35:33 | 显示全部楼层
我也想直接贴上来,可是ctl太长了,图片搞了半天也弄不上来,只好这样了,辛苦大家了!
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2012-11-17 21:50:24 | 显示全部楼层

我看ctl文件里面的xdef和ydef不都是下面这样的一个形式吗
xdef  687 linear  118.49510   0.00450450
ydef  586 linear   30.16783   0.00450450
你那个表示方式没用过,不知道是不是那个问题呢,改一下试试吧
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2012-11-17 22:22:32 | 显示全部楼层
river 发表于 2012-11-17 21:50
我看ctl文件里面的xdef和ydef不都是下面这样的一个形式吗
xdef  687 linear  118.49510   0.00450450
y ...

我用的是墨卡托投影,必须是这种形式的xdef,ydef不然会定位不准确
密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2012-11-17 22:44:07 | 显示全部楼层
Ivy 发表于 2012-11-17 22:22
我用的是墨卡托投影,必须是这种形式的xdef,ydef不然会定位不准确

哦,原来是这样啊。你把gs也贴上来吧,ctl看着没什么问题
密码修改失败请联系微信:mofangbao

新浪微博达人勋

 楼主| 发表于 2012-11-18 09:25:27 | 显示全部楼层
river 发表于 2012-11-17 22:44
哦,原来是这样啊。你把gs也贴上来吧,ctl看着没什么问题

谢谢你,不清楚是不是fortran缺测值设置的问题,昨天改了半天缺测值还是这样子
'reinit'
'open e:\div-wind\850\vwx.ctl'
'open e:\div-wind\850\vwy.ctl'
'enable print e:\div-wind\850\div-wind3.gmf'

'set grads off'
'set timelab off'
'set mpdset cnworld'
'set xlopts 1 6 0.18'
'set ylopts 1 6 0.18'

i=1
while(i<=33)
'set lev 850'
'set t 'i''

'set gxout vector'
'set cthick 6'
'set arrscl 0.5 5'
'd skip(u,5,5);skip(v.2,5,5)'
'q time'
xx=subwrd(result,3)
'draw title 0903 850hPa t='xx''
'print'
'c'
i=i+1
endwhile
'disable print'
;

密码修改失败请联系微信:mofangbao

新浪微博达人勋

发表于 2012-11-18 11:34:38 | 显示全部楼层
Ivy 发表于 2012-11-18 09:25
谢谢你,不清楚是不是fortran缺测值设置的问题,昨天改了半天缺测值还是这样子
'reinit'
'open e:\div- ...

这看着没什么问题啊,奇怪了。你把那几个修饰的语句*掉试试,直接d u;v.2也试一下。
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表