走进不科学 第1043章

作者:新手钓鱼人

  如今机房内有五十多个人,每位操作员实际上要录入的数据差不多是三万组左右吧。

  真·日更三万。

  接着很快。

  众多操作员便开始输入起了数据。

  “气温26.7℃……”

  “湿度27%……”

  “经向风18.5m/s……”

  “散射辐射220W/m^2……”

  很早以前提及过。

  不同于后世计算天气时的WRF、WRF/Chem、WRF-CMAQ等一众模型。

  此时的气象学依旧处于一个非常原始的阶段。

  因此即便有茫茫多的数据,也仍然只能通过解大气方程的方法进行预测。

  也就是此前说过的将地球表面格点化,然后将各气象要素的值进行插值运算,得到未来某一个时刻的值。

  因此实际的模型推导还需要叶笃正他们那边进行,中科院的这台104机主要是把各种数据拟合成更精确的参数而已。

  用后世的话来说就是……

  加工成半成品?

  当然了。

  这个年代虽然科技比较落后,但终究脱离了近代史早期。

  加之计算机所也好歹是兔子们最前端的一处技术中枢,因此此时104机倒也没原始到特别凄惨的地步。

  目前104机虽然还是脱离不了解方程的局限,但它使用的算法……或者说思路其实还挺特别的:

  时间分辨率为1小时,覆盖了6层垂直高度,decoder只有两级。

  更关键的是……

  它的数学本底机理,出自华罗庚老爷子带领的数学研究所。

  没错。

  华罗庚。

  因此虽然依旧是暴力破解计算。

  但在数据充分的情况下,104机模拟的还算很顺利。

  就这样。

  每间隔一段时间,便有一批数据从221厂传送至计算机所。

  计算机所的操作员则对数据导入,104机用老爷车般的算力缓缓产出着结果。

  ……

  八个小时后。

  一份厚达五厘米的数据结果正式出炉。

  原本零散的大气粒子数据被归类成了各个时域的参数、风场、气压场、密度场等众多更详细的模型。

  紧接着。

  这些数据被通过加密方式传回了221厂,最终来到了叶笃正和陶诗言的手里。

  至此。

  整个气象推演也到了最后一个环节。

  那就是气象中心的……

  人力修正与推导。

  “诸位同志。”

  一处临时搭建的大型帐篷内。

  叶笃正与陶诗言看了眼面前的众多同事,彼此对视一眼,只听叶笃正说道:

  “一个星期之前,我们的计算结果很遗憾的失败了,原因是数据的量级不够,我们低估了大自然的神秘。”

  “但在今天,我们拥有了更多量级的数据,经过计算机处理,汇总成了我面前的这份参数文件。”

  “原本首都的竺可桢老先生提出过将这份文件交给首都气象局处理的想法,但被我拒绝了——这一次是我们221厂气象中心雪耻的机会,谁都不能和咱们抢!”

  “为此我还在厂领导面前立下了军令状,要是这次再计算失误,包括我在内,所有参与数据推导的人全都立马去畜牧副业队报道!”

  说罢。

  叶笃正环视了现场一圈,语气略微一缓,继续道:

  “当然了,如果对此有异议的同志,可以现在退出接下来的计算任务。”

  “退出后的成员可以继续在气象中心工作,毕竟基地不可能完全放弃这么一个部门,所以……”

  “现在想退出的同志,你们可以离开了。”

  此话说完。

  现场瞬间变得毫无声息,仿佛呼吸都凝固了。

  但在这片死寂之中。

  却没有一人从位置上站起离去。

  现场每个人的心中,都有一把火在燃烧。

  这把火不是怒火,更不是委屈。

  而是单纯的想要证明自己的价值的冲动。

  在过去这些天里。

  虽然厂里头没有处罚他们,但气象中心的成员们多少都听到了一些风言风语。

  没办法。

  毕竟整个221厂有上万员工呢,总是会有一些碎嘴或者彼此看不顺眼的人存在。

  要是有本事的还好说。

  像老郭陆光达那种专家都受人尊敬,没人会诋毁他们。

  但如果是那种既不干活但又没成果的人,必然就会受到一些非议了。

  更关键的是……

  现场的这些人都是气象领域的从业者,在相关方面从业了数年、十数年甚至数十年。

  从他们的‘道心’出发,也绝不愿意见到自己所学的知识到头来没有任何用处。

  因此这是一场自我雪耻的战斗,也是一次证明自己的机会。

  如果惧怕担责任,他们压根连基地的大门都不会踏进来。

  众多同事沉默但坚定的选择很快也感染了台上的叶笃正,只见他猛然一拍桌子:

  “好!很好!都是好样的!”

  “既然大家都不愿意离开这里,那么今晚咱们就好好的算上一场!”

  说着。

  叶笃正又拍了拍面前的文件,快速说道:

  “各位同志,大家应该都知道,气象数据存在很强的时效性。”

  “所以现在我们长话短说,请小组长们上前认领各自需要计算的文档吧。”

  “首先是第一组……”

  随后在陶诗言和几位组长的协助下。

  这些文件被快速下发到了每个人的手中。

  接着短短五分钟不到。

  帐篷内便响起了刷刷刷的笔算和噼里啪啦的算盘声。

  作为项目的总负责人,叶笃正这次同样也亲自下场进行了计算,并且负责的是最复杂的涡度场。

  涡度。

  这是和散度是差不多用处的概念。

  不过散度描的是述气流的离散程度,一般正值为气流辐散,负值气流辐合。

  而涡度有绝对涡度和相对涡度之分。

  它们的关系可以通过【绝对涡度=相对涡度+2Ω】(其中Ω为地球自转角速度)来计算。

  这部分计算是叶笃正主动申请下来的,毕竟……

  在之前的计算过程中,他就曾经在三维空间流体方面栽了个跟头。

  当时他将笛卡尔坐标系转化为曲面坐标,将连续方程拆分成水平和垂直两个方向分别计算。

  同时在痕量物质方面依据雷诺分解,把瞬时浓度分解为了均值项和湍流项。

  但后来实际情况证明他的思路是错误的,他低估了垂直梯度的实际变动量。

  换而言之……

  他必须要重新设计出一个模型。

  想到这里。

  叶笃正先在算纸上写下了一个方程:

  Du/Dt=-▽(p/ρ)+v▽^2u

  这是很有名的纳维-斯托克斯方程,提出于一百多年前,属于一个描述流体情况的方程组。