【目前用下来,听书声音最全最好用的App,集成4大语音合成引擎,超100种音色,更是支持离线朗读的换源神器,huanyuanapp.org 换源App】
此时此刻。
五十多块屏幕前都坐着一位表情正式的计算机操作员,有男有女,年纪都在二十多岁上下。
他们是国内最早一批...也就是58级的计算机专业毕业生,毕业院校正是中科大。
如果徐云能够飞到这里和其中某位短的女同志摆个把子,他还能成为后世科大数学系某位教授的义父......
“小林,小王,小陈,来帮忙一下数据!”
带着批数据进入机房后。
吴几康快步走到左前方的几位组长面前,分起了各自的数据,并且对众人嘱咐道:
“大家拿到文稿后就可以开始导入数据了,你们的任务事关重大,切记要务必谨慎,千万不能犯错!都明白了吗?”
台下很快响起了一致的回答:
“明白了!”
吴几康则依旧沉着脸,看起来很严肃。
其实吧。
考虑到这次录入的数据很多,个别数据要是出了错误其实不会太过影响最终结果,因此上头还给了一个3%的容错率。
不过吴几康并没有把这事儿告诉几位操作员,怕他们会因为放松而出现意外。
几分钟后。
每位操作员的手上都拿到了相应的数据,并且开始进行了导入。
早先提及过。
老郭他们原先仪器采集到的数据有五百多组,但后来经过叶笃正的计算现这些数据完全不够用。
想要完成气象模拟,必须要多出十万倍才可以。
也就是需要的数据足足有五千多万组。
不过这只是计算需要的数据,导入的工作量并没有这么多。
因为在经过静电分析和X双波段这两个关卡后,这些数据会被归类到一个更规范的区间:
这个区间大概一小时刷新一次,其中雷同的数值会被叠加,通过光电转换元件的电信号直接进行录入。
因此传到吴几康他们手中、需要纯手工导入的数据远远没有五千万组那么多。
以兔子们第一台气象多普勒雷达为例。
国内第一台气象多普勒雷达型号叫做714SD,它单次采集的数据多达1400万组,周期十五分钟一次。
不过在真正计算的时候,电转换元件过滤了98%的雷同数值。
所以真正需要人工导入的数据也就小几十万组罢了——注意,导入和计算是两码事,实际计算的数据依旧是1400万组。
眼下5000万组数据虽然看起来多,但真正需要手工录入的大概也就百万组上下。
如今机房内有五十多个人,每位操作员实际上要录入的数据差不多是三万组左右吧。
真·日更三万。
接着很快。
众多操作员便开始输入起了数据。
“气温26.7℃.....”
“湿度27%.....”
“经向风18.5/s.....”
“散射辐射220W/^2.......”
很早以前提及过。
不同于后世计算天气时的RF-CMAQ等一众模型。
此时的气象学依旧处于一个非常原始的阶段。
因此即便有茫茫多的数据,也仍然只能通过解大气方程的方法进行预测。
也就是此前说过的将地球表面格点化,然后将各气象要素的值进行插值运算,得到未来某一个时刻的值。
因此实际的模型推导还需要叶笃正他们那边进行,中科院的这台104机主要是把各种数据拟合成更精确的参数而已。
用后世的话来说就是......
加工成半成品?
当然了。
这个年代虽然科技比较落后,但终究脱离了近代史早期。
加之计算机所也好歹是兔子们最前端的一处技术中枢,因此此时104机倒也没原始到特别凄惨的地步。
目前104机虽然还是脱离不了解方程的局限,但它使用的算法...或者说思路其实还挺特别的:
时间分辨率为1小时,覆盖了6层垂直高度,dcodr只有两级。
更关键的是....
它的数学本底机理,出自华罗庚老爷子带领的数学研究所。
没错。
华罗庚。
因此虽然依旧是暴力破解计算。
但在数据充分的情况下,104机模拟的还算很顺利。
就这样。
每间隔一段时间,便有一批数据从221厂传送至计算机所。
计算机所的操作员则对数据导入,104机用老爷车般的算力缓缓产出着结果。
.......
八个小时后。
一份厚达五厘米的数据结果正式出炉。
原本零散的大气粒子数据被归类成了各个时域的参数、风场、气压场、密度场等众多更详细的模型。
紧接着。
这些数据被通过加密方式传回了221厂,最终来到了叶笃正和陶诗言的手里。
至此。
整个气象推演也到了最后一个环节。
那就是气象中心的......
人力修正与推导。
“诸位同志。”
一处临时搭建的大型帐篷内。
叶笃正与陶诗言看了眼面前的众多同事,彼此对视一眼,只听叶笃正说道:
“一个星期之前,我们的计算结果很遗憾的失败了,原因是数据的量级不够,我们低估了大自然的神秘。”
“但在今天,我们拥有了更多量级的数据,经过计算机处理,汇总成了我面前的这份参数文件。”
“原本都的竺可桢老先生提出过将这份文件交给都气象局处理的想法,但被我拒绝了——这一次是我们221厂气象中心雪耻的机会,谁都不能和咱们抢!”
“为此我还在厂领导面前立下了军令状,要是这次再计算失误,包括我在内,所有参与数据推导的人全都立马去畜牧副业队报道!”
说罢。
叶笃正环视了现场一圈,语气略微一缓,继续道:
“当然了,如果对此有异议的同志,可以现在退出接下来的计算任务。”
“退出后的成员可以继续在气象中心工作,毕竟基地不可能完全放弃这么一个部门,所以......”
“现在想退出的同志,你们可以离开了。”
此话说完。
现场瞬间变得毫无声息,仿佛呼吸都凝固了。
但在这片死寂之中。
却没有一人从位置上站起离去。
现场每个人的心中,都有一把火在燃烧。
这把火不是怒火,更不是委屈。
而是单纯的想要证明自己的价值的冲动。
在过去这些天里。
虽然厂里头没有处罚他们,但气象中心的成员们多少都听到了一些风言风语。
没办法。
毕竟整个221厂有上万员工呢,总是会有一些碎嘴或者彼此看不顺眼的人存在。
要是有本事的还好说。
像老郭陆光达那种专家都受人尊敬,没人会诋毁他们。
但如果是那种既不干活但又没成果的人,必然就会受到一些非议了。
更关键的是.....
现场的这些人都是气象领域的从业者,在相关方面从业了数年、十数年甚至数十年。
从他们的‘道心’出,也绝不愿意见到自己所学的知识到头来没有任何用处。
因此这是一场自我雪耻的战斗,也是一次证明自己的机会。
如果惧怕担责任,他们压根连基地的大门都不会踏进来。
众多同事沉默但坚定的选择很快也感染了台上的叶笃正,只见他勐然一拍桌子:
“好!很好!都是好样的!”
“既然大家都不愿意离开这里,那么今晚咱们就好好的算上一场!”
说着。
叶笃正又拍了拍面前的文件,快速说道:
“各位同志,大家应该都知道,气象数据存在很强的时效性。”
“所以现在我们长话短说,请小组长们上前认领各自需要计算的文档吧。”
“先是第一组......”
随后在陶诗言和几位组长的协助下。
这些文件被快速下到了每个人的手中。
接着短短五分钟不到。
帐篷内便响起了刷刷刷的笔算和噼里啪啦的算盘声。
作为项目的总负责人,叶笃正这次同样也亲自下场进行了计算,并且负责的是最复杂的涡度场。
涡度。
这是和散度是差不多用处的概念。
不过散度描的是述气流的离散程度,一般正值为气流辐散,负值气流辐合。
而涡度有绝对涡度和相对涡度之分。
它们的关系可以通过【绝对涡度=相对涡度+2Ω】(其中Ω为地球自转角速度)来计算。
这部分计算是叶笃正主动申请下来的,毕竟.....
在之前的计算过程中,他就曾经在三维空间流体方面栽了个跟头。
当时他将笛卡尔坐标系转化为曲面坐标,将连续方程拆分成水平和垂直两个方向分别计算。
同时在痕量物质方面依据雷诺分解,把瞬时浓度分解为了均值项和湍流项。
但后来实际情况证明他的思路是错误的,他低估了垂直梯度的实际变动量。
换而言之.....
他必须要重新设计出一个模型。
想到这里。
叶笃正先在算纸上写下了一个方程:
Du/Dt=??(p/ρ)+v?2u
这是很有名的纳维-斯托克斯方程,提出于一百多年前,属于一个描述流体情况的方程组。
其中的斯托克斯想必有些同学会感觉眼熟——没错,这个斯托克斯就是1850副本中徐云的便宜导师......
它关于u的边界条件是 u=0。
接着叶笃正很快又写道:
δT=(?T/?t)δt+(?T/?x)δx+…
δx=uxδt,进而
DT/Dt=?T?t+ux?T?x+uy?T?y+uz?T?z=?T?t+(u??)T.......
DA/Dt=?A?t+(u??)A.....
所以oks方程可以改写为:
Du/Dt=?u?t+(u??)u=??(p/ρ)+v?2u。
写到这里。
叶笃正不由笔尖一顿。
上头这部分推导是他在前些天想出来的优化形式,弥补了自己原先思路的不足。
但是.....
到了变式后的这一步。
叶笃正就不知道该如何继续了。
没错。
不是计算或者推导不出哪个数值。
而是不知道该怎么推导了。
为此他还请教过都的竺可桢老先生,但即便是竺老也没什么办法。
竺老只是给出了一个考虑非线性项的想法,但叶笃正总感觉这样算有点问题。
而就在叶笃正一筹莫展之际。
不知道为什么,叶笃正忽然感觉自己的耳边有些异样。
怎么说呢.......
仿佛有人在对着自己的耳边吹风?
于是叶笃正下意识转过头,结果正正好对上了一张裹的跟木乃尹似的脸。
叶笃正:
“?!”
......
注:
要改的都改好了,明天开始继续日万,冲鸭!
!