求助：sst数据

lee309

对于不同维数的数据进行四则运算时会出错
sst为全球海平面温度，数据为所有经纬度的数据
(179,345)    94.20001
(179,346)    94.20001
(179,347)    94.20001
(179,348)    94.20001
(179,349)    94.20001
(179,350)    94.20001
(179,351)    94.20001
(179,352)    94.20001
(179,353)    94.20001
(179,354)    94.20001
(179,355)    94.20001
(179,356)    94.20001
(179,357)    94.20001
(179,358)    94.20001
(179,359)    94.20001
t为空气温度，数据为网格数据
(41,29)    301.7791
(41,30)    301.8303
(41,31)    301.8813
(41,32)    301.9309
(41,33)    301.9777
(41,34)    302.0218
(41,35)    302.066
(41,36)    302.1104
(41,37)    302.1546
(41,38)    302.1984
(41,39)    302.2416
(41,40)    302.2839
(41,41)    302.3249
因为sst和t维数不一致对于n = sst - t 会报错 Dimension size, for dimension number 0, of operands does not match, can't continue
这种问题怎么解决？
另外sst里的time也很奇怪，不知道对应的年月日
Variable: time
Type: double
Total Size: 3168 bytes
            396 values
Number of Dimensions: 1
Dimensions and sizes:    [time | 396]
Coordinates:
            time: [66443..78466]
Number Of Attributes: 9
  bounds :    time_bnds
  axis :    T
  standard_name :    time
  prev_avg_period :    0000-00-07 00:00:00
  avg_period :    0000-01-00 00:00:00
  delta_t :    0000-01-00 00:00:00
  actual_range :    ( 66443, 78466 )
  long_name :    Time
  units :    days since 1800-1-1 00:00:00
(0)    66443
(1)    66474
(2)    66505
(3)    66533
(4)    66564
(5)    66594
(6)    66625
(7)    66655
(8)    66686
(9)    66717
(10)    66747
(11)    66778
(12)    66808
(13)    66839
(14)    66870
(15)    66898
(16)    66929
(17)    66959
(18)    66990
(19)    67020
(20)    67051
(21)    67082
(22)    67112
(23)    67143
(24)    67173
(25)    67204
(26)    67235
(27)    67264
(28)    67295
(29)    67325
(30)    67356
(31)    67386
(32)    67417
(33)    67448
(34)    67478
(35)    67509
(36)    67539
(37)    67570
(38)    67601
(39)    67629
(40)    67660
(41)    67690
(42)    67721
(43)    67751
(44)    67782
(45)    67813
(46)    67843
(47)    67874
(48)    67904
(49)    67935
(50)    67966
(51)    67994
(52)    68025
(53)    68055
(54)    68086
(55)    68116
(56)    68147
(57)    68178
(58)    68208
(59)    68239
(60)    68269
(61)    68300
(62)    68331
(63)    68359
(64)    68390
(65)    68420
(66)    68451
(67)    68481
(68)    68512
(69)    68543
(70)    68573
(71)    68604
(72)    68634
(73)    68665
(74)    68696
(75)    68725
(76)    68756
(77)    68786
(78)    68817
(79)    68847
(80)    68878
(81)    68909
(82)    68939
(83)    68970
(84)    69000
(85)    69031
(86)    69062
(87)    69090
(88)    69121
(89)    69151
(90)    69182
(91)    69212
(92)    69243
(93)    69274
(94)    69304
(95)    69335
(96)    69365
(97)    69396
(98)    69427
(99)    69455
(100)    69486
(101)    69516
(102)    69547
(103)    69577
(104)    69608
(105)    69639
(106)    69669
(107)    69700
(108)    69730
(109)    69761
(110)    69792
(111)    69820
(112)    69851
(113)    69881
(114)    69912
(115)    69942
(116)    69973
(117)    70004
(118)    70034
(119)    70065
(120)    70095
(121)    70126
(122)    70157
(123)    70186
(124)    70217
(125)    70247
(126)    70278
(127)    70308
(128)    70339
(129)    70370
(130)    70400
(131)    70431
(132)    70461
(133)    70492
(134)    70523
(135)    70551
(136)    70582
(137)    70612
(138)    70643
(139)    70673
(140)    70704
(141)    70735
(142)    70765
(143)    70796
(144)    70826
(145)    70857
(146)    70888
(147)    70916
(148)    70947
(149)    70977
(150)    71008
(151)    71038
(152)    71069
(153)    71100
(154)    71130
(155)    71161
(156)    71191
(157)    71222
(158)    71253
(159)    71281
(160)    71312
(161)    71342
(162)    71373
(163)    71403
(164)    71434
(165)    71465
(166)    71495
(167)    71526
(168)    71556
(169)    71587
(170)    71618
(171)    71647
(172)    71678
(173)    71708
(174)    71739
(175)    71769
(176)    71800
(177)    71831
(178)    71861
(179)    71892
(180)    71922
(181)    71953
(182)    71984
(183)    72012
(184)    72043
(185)    72073
(186)    72104
(187)    72134
(188)    72165
(189)    72196
(190)    72226
(191)    72257
(192)    72287
(193)    72318
(194)    72349
(195)    72377
(196)    72408
(197)    72438
(198)    72469
(199)    72499
(200)    72530
(201)    72561
(202)    72591
(203)    72622
(204)    72652
(205)    72683
(206)    72714
(207)    72742
(208)    72773
(209)    72803
(210)    72834
(211)    72864
(212)    72895
(213)    72926
(214)    72956
(215)    72987
(216)    73017
(217)    73048
(218)    73079
(219)    73108
(220)    73139
(221)    73169
(222)    73200
(223)    73230
(224)    73261
(225)    73292
(226)    73322
(227)    73353
(228)    73383
(229)    73414
(230)    73445
(231)    73473
(232)    73504
(233)    73534
(234)    73565
(235)    73595
(236)    73626
(237)    73657
(238)    73687
(239)    73718
(240)    73748
(241)    73779
(242)    73810
(243)    73838
(244)    73869
(245)    73899
(246)    73930
(247)    73960
(248)    73991
(249)    74022
(250)    74052
(251)    74083
(252)    74113
(253)    74144
(254)    74175
(255)    74203
(256)    74234
(257)    74264
(258)    74295
(259)    74325
(260)    74356
(261)    74387
(262)    74417
(263)    74448
(264)    74478
(265)    74509
(266)    74540
(267)    74569
(268)    74600
(269)    74630
(270)    74661
(271)    74691
(272)    74722
(273)    74753
(274)    74783
(275)    74814
(276)    74844
(277)    74875
(278)    74906
(279)    74934
(280)    74965
(281)    74995
(282)    75026
(283)    75056
(284)    75087
(285)    75118
(286)    75148
(287)    75179
(288)    75209
(289)    75240
(290)    75271
(291)    75299
(292)    75330
(293)    75360
(294)    75391
(295)    75421
(296)    75452
(297)    75483
(298)    75513
(299)    75544
(300)    75574
(301)    75605
(302)    75636
(303)    75664
(304)    75695
(305)    75725
(306)    75756
(307)    75786
(308)    75817
(309)    75848
(310)    75878
(311)    75909
(312)    75939
(313)    75970
(314)    76001
(315)    76030
(316)    76061
(317)    76091
(318)    76122
(319)    76152
(320)    76183
(321)    76214
(322)    76244
(323)    76275
(324)    76305
(325)    76336
(326)    76367
(327)    76395
(328)    76426
(329)    76456
(330)    76487
(331)    76517
(332)    76548
(333)    76579
(334)    76609
(335)    76640
(336)    76670
(337)    76701
(338)    76732
(339)    76760
(340)    76791
(341)    76821
(342)    76852
(343)    76882
(344)    76913
(345)    76944
(346)    76974
(347)    77005
(348)    77035
(349)    77066
(350)    77097
(351)    77125
(352)    77156
(353)    77186
(354)    77217
(355)    77247
(356)    77278
(357)    77309
(358)    77339
(359)    77370
(360)    77400
(361)    77431
(362)    77462
(363)    77491
(364)    77522
(365)    77552
(366)    77583
(367)    77613
(368)    77644
(369)    77675
(370)    77705
(371)    77736
(372)    77766
(373)    77797
(374)    77828
(375)    77856
(376)    77887
(377)    77917
(378)    77948
(379)    77978
(380)    78009
(381)    78040
(382)    78070
(383)    78101
(384)    78131
(385)    78162
(386)    78193
(387)    78221
(388)    78252
(389)    78282
(390)    78313
(391)    78343
(392)    78374
(393)    78405
(394)    78435
(395)    78466

yuxiao19910

握爪~！！ 我正在搞这个 WRF+SST， 以后还要多请教下 楼主

lee309

有什么办法把sst数据也变成网格数据么？t的数据是wrfout文件中的数据，两者怎么转化？

lee309

哪位用过sst数据的能解释下么，这方面小白

letsdo8

这个最好是读wrflowinp_d01（netCDF格式的），这文件中应该有WRF格点上的SST。

namelist.input里得有sst_update=1

lee309

letsdo8 发表于 2014-12-16 10:38
这个最好是读wrflowinp_d01（netCDF格式的），这文件中应该有WRF格点上的SST。

namelist.input里得有sst ...

请问你说的wrflowinp_d01文件是什么文件，我怎么没见过？
这个sst_update=1我在手册里见过，但不知道什么意思，能解释下么？谢谢！！！！！！

lee309

lee309 发表于 2014-12-16 10:54
请问你说的wrflowinp_d01文件是什么文件，我怎么没见过？
这个sst_update=1我在手册里见过，但不知道什 ...

namelist.input中sst_update=1,我添加了以下语句
io_form_auxinput4                   = 2
auxinput4_inname                    = “wrflowinp_d<domain>”
auxinput4_interval                  = 360, 360, 360,
sst_update                              =1
一直报错
FATAL CALLED FROM FILE:  <stdin>  LINE:    9884
ERRORS while reading one or more namelists from namelist.input.
请问这是怎么回事？除了添加语句还需要额外的sst数据文件吗？

lee309

lee309 发表于 2014-12-16 10:54
请问你说的wrflowinp_d01文件是什么文件，我怎么没见过？
这个sst_update=1我在手册里见过，但不知道什 ...

顺便问下，这个SST是海表面空气温度，还是海表面海水温度？

lee309

letsdo8 发表于 2014-12-16 10:38
这个最好是读wrflowinp_d01（netCDF格式的），这文件中应该有WRF格点上的SST。

namelist.input里得有sst ...

顺便问下，这个SST是海表面空气温度，还是海表面海水温度？

letsdo8

lee309 发表于 2014-12-16 13:29
顺便问下，这个SST是海表面空气温度，还是海表面海水温度？

海表面海水温度