无意中在网上看到篇关于Wafa Sultan的评论，我才知道，有个如此彪悍的人，说过如此彪悍的话。后来，找来了半岛电视台的现场录像和文字记录，震惊不已。

Wafa Sultan大妈，早年毕业于叙利亚阿勒颇大学，后赴美获得心理学博士学位，留美任教。Wafa曾是虔诚的穆斯林，1979年，她目睹了伊斯兰激进分子冲进大学教室，当众残忍杀害正在讲课的教授的一幕，从此，对自己的宗教产生了怀疑。经过多年的痛苦思考，有了05年6月在半岛台大战伊斯兰宗教学士的一幕。

她的观点与亨廷顿有别：1，文明不冲突，只竞争。2，伊斯兰和西方的冲突不是文明冲突，而是时代冲突。

牛既不A也不C。偶强烈赞同她的第二个观点。我们一提伊斯兰，立马就冒出个一夫多妻。男人们窃笑的同时，也大叹：落后啊。。。    

落后在此指跟不上时代。伊斯兰教创立时期，偶们不得不说它是个超前的先进的宗教：他改变了阿拉伯的社会结构，发扬部落传承，抵制因物质水平提高而带来的不良思想，当时，那是一片“和谐”，“现代化”啊。而后，他也经历了辉煌，此时，应是伊斯兰与时代刚好适应的时候，既不超前也不滞后。而再往后，在西方基督教世界进步的同时，它开始相对以负加速度滞后。因此现在，伊斯兰世界的时间坐标停留在中世纪。

为什么会这样？我想，一个直接原因是没有大范围的宗教改革。西方出现过加尔文和路德，有教廷的多次分裂，这在客观上利于基督教内部产生竞争，优胜劣汰。而伊斯兰教，稍有宗教改革便被强大的保守力量扼杀：至今印度的一支激进改革教派仍被世界伊斯兰大会认定为叛教者，非洲的某些与当地宗教融合的较为温和教派也被视为异教徒。然而，为什么会这样？以前我一直以为，是伊斯兰教的教长权威过大所致，可这显然不是事实，现在各大小派别泛滥和争斗，伊斯兰世界乱作一团，这应该是典型的缺乏权威领袖所致。

究其根本，一切源自《古兰经》。穆斯林相信：《古兰经》是天启的最后一部经书，前面的全部作废（讨拉，圣经等），往后不再会有新的；穆罕默德是最后的使者，前面的都失败了（摩西和耶稣等），往后不再会有新的。这一基本信条成功的将犹太教和基督教列为异端，并且抑制新宗教的出现。按照这样的逻辑，世界在《古兰经》降世的这一刻，就处于了一种停止状态，一种最优的状态，并且这种状态将永远持续，不可终结。我想，穆罕默德当初之所以树立这样的信条，更多是出于政治原因保护新兴宗教。然而，当《古兰经》成为指导人们思维和行动的唯一依据的时候，这种一味否定新事物的态度无疑将是个大灾难。由于这一信条，人们对任何改良都持敌视态度，《古兰经》成了超越任何时代的唯一标准。任何人越是保守，越是原教旨，那就越得势。最终，伊斯兰沦为保守势力的堡垒。物质社会的时间轴在高速向前发展，而精神世界的时间轴则拒绝更新，矛盾必然出现。

至于文明冲突与否，偶想，那还是会有的。我们可以说是中世纪的伊斯兰文明和21实际的基督教文明在冲突。当然，更多的的确应该是竞争，就像华夏文明和基督教文明，偶们处于劣势哦，这个就不多说了。

在师弟推荐下，昨看了个感人的片《我的邻居山田君》，很值。

整个影片都在平淡中进行：琐碎的家庭生活，淡雅的电影配乐，简洁的Q版写意画风。一部不折不扣的小制作。讲的是一个普通家庭的故事，由高中生山田，幼稚园的妹妹野野子，面临强大工作压力的爸爸，职业家庭主妇妈妈以及倔强的小老太太奶奶组成，它是千万个日本家庭中抽象出来的一个典型。

影片中描述了这个可爱的家庭里发生的一些平常的不能再平常的事：一日三餐，一年四季。叙述的方式也非常特别，很喜欢那种比喻的手法，生动而幽默。最令人震撼的是大雨中爸爸在半路上遇到家人打伞来接他，一家五口人走在雨中又说又笑的场景：幽深的淡墨色的深巷里，几个温馨的背影渐渐远去。

其实，每个家庭都有这种温馨的时刻，只是他被淹没在了平日的琐碎里了而已。


P.S:  主题曲《ひとりぼっちはやめた(Quit Being Alone)》超好听，正好表达了这种悠悠的平淡的美。

实在是很恼人，想不到原来物性参数这样的复杂。

1，对于FROST前辈的那篇关于焓流电弧模型中的SF6转移特性的文章，仍旧需要加深理解。尽管老先生是离散的测出了1，2，4，8，16个大气压下的各种参数随温度变化的特征，但是我们应当认识到，这些物性参数应该是温度和压强两个变量的函数才是。因此，一系列参数表现为phi=f(P,T)形式。

2，文中提到的各种物性的表征参数以及参数的单位和PHOENICS中的有很大差异。电导率sigma的单位文中为mho/cm，实际上即为siemens/cm，使用中需要换算为siemens/m，并且这里的电导率呈现一维的表述，它是否为单位面积上的电导与长度的比值？对于粘度，有许多种不同的表示方法。文中的ita为动力粘度(dynamic viscosity)，单位poise，即g/cm*sec，而phoenics中默认输入的则为更为国际化的运动粘度(kinematic viscosity，v)，单位stoke，即cm2/sec，两者转换方法为v＝ita/roh，roh为密度。

3，对于SF6的静辐射系数(net emission coefficient)，从AUBRECHT前辈的文章看来是仅与温度和电弧半径有关，在实际中，认为在低温(<10000)时不存在，在电焊兄的文中，判断为电导率sig<10e-4，即sig=0时不存在。

4，终于明白了电焊兄的程序中的几个奇怪的10×50数组是用来干啥的，原来对应着P和T，插值求各物性参数。唉，开始还想着用拟合，一看那些怪怪的曲线，我投降了，还是插值吧，可怜的，得自己编程序在GROUND里……

天气是我每天都关注的。今天一大早，天色灰白，阵风四起。我断定，今天必然会有好东西。

接着，我们眼看着天空变得越来越黑暗。正午时分，写字楼里的人们不得不打开所有的照明设备，向外望去，一种错觉：似是午夜。高空中黑雾蒸腾，压低，再压低，最后滚滚的罩在我们头上。伴随着强风撕心裂肺的嘶吼，暴雨如期而至：硕大的雨点斜打在玻璃窗上，每一个都足有酒瓶底那么大。

雨声渐渐紧密，四周电闪雷鸣。停电了，四周一片漆黑，除了雷电像照相机的闪光灯一样，不断对周围的事物进行曝光。过来些时候，雷电远去了，只留下失去靠山的雨，安静的下了一整下午。

爽，这个暴雨同学应该表扬……

前些日子INFORM老是出问题，电势边界条件怎么也加不上。现在，俺可以大声的说：让INFORM见鬼去吧！

INFORM虽然语句简单，但其并不像PLANT一样最终将语句写入GROUND里，而是写在Q1EAR里。这个Q1EAR是一个与Q1和GROUND同样具有GROUP对应的FOR文件。

实际上，通过GROUND可以做到任何事情，电势边界条件也不例外：

1，在区域边界处加入无厚度的NewObject，并将其类型设定为UserDefined。

2，在此形体Attribute中有Patch和Coval设置，通过其对电势FI的Coefficient和Value进行设定，即可确定电势边界条件。

此处，首先需要对Patch进行命名，并对其类型进行设定（诸如EAST，NORTH，HIGH之类）。之后对下方感兴趣的参数进行设定，如电势FI的Coefficient为Fixval，Value则为电势的值。

3，以往需要在PLANT中写入的电场强度计算公式，可以直接在GROUND中写入，其规律可由PLANT自动生成代码获得。

但是，在此步骤会遇到最后一行列网格由于缺少后续的FI值而出现极大或极小错误数据的情况，需要在GROUND的写入时注意，对于最后的行列专门用语句描述。

这里，描述语句将写入GROUP 19. Special calls to GROUND from EARTH下的SECTION 6.Finish of iz slab，如下：

C   * ------------------- SECTION 6 ---- Finish of iz slab.
C      Special calls name: SC0601    
       IF(ISTEP.GE.1       .AND.ISTEP.LE.LSTEP   ) THEN                         
       IF(IZ.GE.1       .AND.IZ.LE.NZ      ) THEN                               
       LFEX   =L0F(INAME('EX    '))                                            
       LFFI  =L0F(INAME('FI    '))                                             
       LFDXG2=L0F(DXG2D )                                                      
       DO 19601 IX=1       ,59                                                 
        IADD=NY*(IX-1)                                                         
       DO 19601 IY=1       ,NY                                                
        I=IY+IADD                                                              
       L0EX  =LFEX  +I                                                         
       L0FI  =LFFI  +I                                                         
       L0DXG2=LFDXG2+I                                                         
19601 F(L0EX    )=(F(L0FI+NY)-F(L0FI))/F(L0DXG2)                               
      ENDIF                                                                    
      ENDIF                                                     
C      Special calls name: SC0602     
       IF(ISTEP.GE.1       .AND.ISTEP.LE.LSTEP   ) THEN                         
       IF(IZ.GE.1       .AND.IZ.LE.NZ      ) THEN                               
       LFEY   =L0F(INAME('EY    '))                                            
       LFFI  =L0F(INAME('FI    '))                                             
       LFDYG2=L0F(DYG2D )                                                      
       DO 19602 IX=1       ,NX                                                 
        IADD=NY*(IX-1)                                                         
       DO 19602 IY=1       ,59                                                 
        I=IY+IADD                                                              
       L0EY  =LFEY  +I                                                         
       L0FI  =LFFI  +I                                                         
       L0DYG2=LFDYG2+I                                                         
19602 F(L0EY    )=(F(L0FI+1)-F(L0FI))/F(L0DYG2)                                
      ENDIF                                                                    
      ENDIF             
C      Special calls name: SC0604      
       IF(ISTEP.GE.1       .AND.ISTEP.LE.LSTEP   ) THEN                         
       IF(IZ.GE.1       .AND.IZ.LE.NZ      ) THEN                               
       LFEY   =L0F(INAME('EY    '))                                            
       LFFI  =L0F(INAME('FI    '))                                             
       LFDYG2=L0F(DYG2D )                                                      
       DO 19604 IX=1       ,NX                                                 
        IADD=NY*(IX-1)                                                         
       DO 19604 IY=60      ,NY                                                 
        I=IY+IADD                                                              
       L0EY  =LFEY  +I                                                         
       L0FI  =LFFI  +I                                                         
       L0DYG2=LFDYG2+I                                                         
19604 F(L0EY    )=(F(L0FI)-F(L0FI-1))/F(L0DYG2)                                
      ENDIF                                                                    
      ENDIF                       
C      Special calls name: SC0605 
       IF(ISTEP.GE.1       .AND.ISTEP.LE.LSTEP   ) THEN                         
       IF(IZ.GE.1       .AND.IZ.LE.NZ      ) THEN                               
       LFEX   =L0F(INAME('EX    '))                                            
       LFFI  =L0F(INAME('FI    '))                                             
       LFDXG2=L0F(DXG2D )                                                      
       DO 19605 IX=60       ,NX                                                 
        IADD=NY*(IX-1)                                                         
       DO 19605 IY=1       ,NY                                                 
        I=IY+IADD                                                              
       L0EX  =LFEX  +I                                                         
       L0FI  =LFFI  +I                                                         
       L0DXG2=LFDXG2+I                                                         
19605 F(L0EX    )=(F(L0FI)-F(L0FI-NY))/F(L0DXG2)                               
      ENDIF                                                                    
      ENDIF                    
C      Special calls name: SC0600   
       IF(ISTEP.GE.1       .AND.ISTEP.LE.LSTEP   ) THEN                         
       IF(IZ.GE.1       .AND.IZ.LE.NZ      ) THEN                               
       LFE    =L0F(INAME('E     '))                                            
       LFEX  =L0F(INAME('EX    '))                                             
       LFEY  =L0F(INAME('EY    '))                                             
       DO 19600 IX=1       ,NX                                                 
        IADD=NY*(IX-1)                                                         
       DO 19600 IY=1       ,NY                                                 
        I=IY+IADD                                                              
       L0E   =LFE   +I                                                         
       L0EX  =LFEX  +I                                                         
       L0EY  =LFEY  +I                                                         
19600 F(L0E     )=SQRT(F(L0EX)**2+F(L0EY)**2)                                  
      ENDIF                                                                    
      ENDIF                     

可见，在对一个x×y为60×60的区域进行计算时，需要先对1－59行的Ex进行计算，在对1－59的Ey进行计算，之后，在分别对最后一行列的Ex和Ey进行计算，而E则始终在整个计算域。

在GROUND中写入的语句可靠性好，可读性强，条例清晰，便于检查和修正。但是，似乎只可以对无厚度边界进行定义，即意味着对于形体复杂的物体这种方法可能相当繁复。为了获得无厚边界我们也许不得不用BFC网格，这样工作量就更大了。看来，若要对矩形结构网格中的复杂形体进行电势边界加载INFORM的确是有它的好处的。

时间过得真快，转眼间，g结束都有两周了。又有好几天都没来更新了，看来每天一更真的不那么简单。

今天又下了雨，让人清新了一番。雨一停，甚么又都回去了：课题看不到进展，尽管还有两周就要答辩；文也发不了；T也报不上；找工作的事也没眉目。看看周围，怎么都有一种乱七八糟的感觉，不对路。

暑假快来了，以前种种令人兴奋的安排估计都不会成为现实。其实很想和虫去西藏。

希望虫子的过敏快点好起来。هى ثًقفس   

然后，猛做课题去……