用wrfout绘制的月降雨有负值

风铃bling

模拟了2017年的降雨，采用的参数化方案为WSM6/BMJ/MYJ/MO/NOAH/CAM，计算了每个月的降雨，得到的月降雨均值竟然有负值，导致整体的降雨趋势很奇怪。附ncl脚本的核心部分：
-----------------------------------------------------------
    rainnc1 = wrf_user_getvar(a,"RAINNC",0)
    rainnc2 = wrf_user_getvar(a,"RAINNC",248)
    rainnc3 = wrf_user_getvar(a,"RAINNC",496)
    rainnc4 = wrf_user_getvar(a,"RAINNC",720)
    rainnc5 = wrf_user_getvar(a,"RAINNC",968)
    rainnc6 = wrf_user_getvar(a,"RAINNC",1208)

    rainc1 = wrf_user_getvar(a,"RAINC",0)
    rainc2 = wrf_user_getvar(a,"RAINC",248)
    rainc3 = wrf_user_getvar(a,"RAINC",496)
    rainc4 = wrf_user_getvar(a,"RAINC",720)
    rainc5 = wrf_user_getvar(a,"RAINC",968)
    rainc6 = wrf_user_getvar(a,"RAINC",1208)
    r1612=rainnc2+rainc2-rainnc1-rainc1
    r1701=rainnc3+rainc3-rainnc2-rainc2
    r1702=rainnc4+rainc4-rainnc3-rainc3
    r1703=rainnc5+rainc5-rainnc4-rainc4
    r1704=rainnc6+rainc6-rainnc5-rainc5
    r1705=rainnc7+rainc7-rainnc6-rainc6
   a0 = avg(r1612)
   a1 = avg(r1701)
   a2 = avg(r1702)
   a3 = avg(r1703)
   a4 = avg(r1704)
   a5 = avg(r1705)
   print(a0)
   print(a1)
   print(a2)
   print(a3)
   print(a4)
   print(a5)


-----------------------------------------------------------
以下为脚本运行的部分结果，可以看到一个负值，这让我怀疑别的数据的可靠性。不知道是否计算方法有误，还望指点。
-----------------------------------------------------------
Variable: a0
Type: float
Total Size: 4 bytes
            1 values
Number of Dimensions: 1
Dimensions and sizes: [1]
Coordinates:
(0) 57.04472

Variable: a1
Type: float
Total Size: 4 bytes
            1 values
Number of Dimensions: 1
Dimensions and sizes: [1]
Coordinates:
(0) 81.98467

Variable: a2
Type: float
Total Size: 4 bytes
            1 values
Number of Dimensions: 1
Dimensions and sizes: [1]
Coordinates:
(0) 101.5718

Variable: a3
Type: float
Total Size: 4 bytes
            1 values
Number of Dimensions: 1
Dimensions and sizes: [1]
Coordinates:
(0) 87.54868

Variable: a4
Type: float
Total Size: 4 bytes
            1 values
Number of Dimensions: 1
Dimensions and sizes: [1]
Coordinates:
(0) 91.13332

Variable: a5
Type: float
Total Size: 4 bytes
            1 values
Number of Dimensions: 1
Dimensions and sizes: [1]
Coordinates:
(0) -27.22223


------------------------------------------------
计算的部分大部分是重复的，绘图的部分问题不大，所以想请教一下各位有否遇过这种情况。在此先感谢了。
------------------------------------------------

freekiller

rainnc7+rainc7 哪里来的

风铃bling

freekiller 发表于 2019-6-14 22:24
rainnc7+rainc7 哪里来的

是这样的，因为转了一年的数据，代码重复的太多了，所以只截取了出现负值月份附近几个月的代码，疏忽之中没有保证前后统一，望谅解！