WRF初始场中尺度滤波试验出错求助！

schliezer

     大家好，我又来了。
     最近在做一个关于WRF初始场的滤波试验。主要过程是首先./real.exe生成三个wrfinput（我是三重双向嵌套）和一个wrfbdy，然后分别对三个input的高度场、风场进行barnes带通滤波，用滤波后中尺度场分别对d3，d2、d3，d1、d2、d3的input文件相应部分进行替换，最后驱动模式模拟。
     那么现在问题来了：仅替换d3以及替换d2、d3后./wrf.exe都能顺利转完，但同时替换d1、d2、d3后没转多久就出现segmentation fault ：core  dumped错误，网上有说是数组越界错误，我检查了d1高度场、风场的数组大小，比对了下也没发现前后有出入。不知道是不是与最外一层区域初始场要和边界场协调或是其他什么物理、动力原因相关，而非单纯的计算、内存错误？还请各位帮忙分析下，谢谢！
贴出log和嵌套设置图：


Sender: LSF System <lsfadmin@c0306>
Subject: Job 242867: <./wrf.exe> in cluster <hpc> Exited

Job <./wrf.exe> was submitted from host <ln01> by user <zuohch1> in cluster <hpc>.
Job was executed on host(s) <c0306>, in queue <blade>, as user <zuohch1> in cluster <hpc>.
</gpfshome/zuohch1> was used as the home directory.
</gpfshome/zuohch1/zhuy/usr/local/wrf/WRFV3/run> was used as the working directory.
Started at Wed Dec 17 15:17:20 2014
Results reported at Wed Dec 17 15:36:02 2014

Your job looked like:

------------------------------------------------------------
# LSBATCH: User input
./wrf.exe
------------------------------------------------------------

Exited with exit code 139.

Resource usage summary:

    CPU time :               1105.08 sec.
    Max Memory :             2761 MB
    Average Memory :         2325.79 MB
    Total Requested Memory : -
    Delta Memory :           -
    (Delta: the difference between total requested memory and actual max usage.)
    Max Swap :               2970 MB

    Max Processes :          3
    Max Threads :            4

The output (if any) follows:

Namelist logging not found in namelist.input. Using registry defaults for variables in logging.
/gpfshome/zuohch1/.lsbatch/1418798948.242867: line 8: 14646 Segmentation fault      (core dumped) ./wrf.exe

晓小

你好，我也正在学滤波，请问可以给我个滤波的程序吗

Mid_Farmer

在运行exe之前输入 ulimit -s unlimited试试

schliezer

Mid_Farmer 发表于 2014-12-18 12:05
在运行exe之前输入 ulimit -s unlimited试试

试过、不行

unice

可能是滤波得到的中尺度场动力不协调吧

schliezer

unice 发表于 2014-12-18 12:50
可能是滤波得到的中尺度场动力不协调吧

怎么个不协调法呢？为什么改了d3或同时改了d2、d3都至少可以正常转通，而d1也滤波就不行了，难道是d1和d2、d3动力不协调？再者，如果真是动力不协调，你觉得这个实验又该怎么设计？谢谢！

unice

schliezer 发表于 2014-12-18 12:54
怎么个不协调法呢？为什么改了d3或同时改了d2、d3都至少可以正常转通，而d1也滤波就不行了，难道是d1和d2 ...

这种试验我没做过，所以具体就不清楚了

freekiller

肯定不是数组越界，积分不稳定了。应该出现了NaN或者Inf

schliezer

freekiller 发表于 2014-12-19 08:32
肯定不是数组越界，积分不稳定了。应该出现了NaN或者Inf

肯定不是数组越界或积分不稳定？你是说wrfout里出现NaN或Inf吗？用ncdump看了三个wrfout的PHB（基态位势高度），都是正常的，但是三个wrfout的大小远远小于以前正常转出来的结果，另外每次出现上述错误时还会生成一个core.(几个数字，不知道啥意思)文件，近3000M大，不是标准的nc格式，ncdump无法查看。

freekiller

core文件就是出现Segmentation fault 错误的时候输出的，可以调试查看core文件，百度一下。但由于程序不是自己写的，所以基本解决不了你的问题。还是琢磨你修改的地方吧