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