▍1. 一维可压N-S方程求解 c代码及fortran代码
利用二阶迎风型TVB差分格式求解一维可压缩黏性流动问题
利用二阶迎风型TVB差分格式求解一维可压缩黏性流动问题
利用气象常用程序fortron编写的小波变换程序。利用气象常用程序fortron编写的小波变换程序。利用气象常用程序fortron编写的小波变换程序。利用气象常用程序fortron编写的小波变换程序。利用气象常用程序fortron编写的小波变换程序。
40中差分格式的Fortran源代码40中差分格式的Fortran源代码40中差分格式的Fortran源代码40中差分格式的Fortran源代码40中差分格式的Fortran源代码40中差分格式的Fortran源代码40中差分格式的Fortran源代码40中差分格式的Fortran源代码40中差分格式的Fortran源代码40中差分格式的Fortran源代码40中差分格式的Fortran源代码40中差分格式的Fortran源代码
该压缩包包含求解线性波方程、浅水波方程和欧拉方程的一维、二维、三维的源代码,是计算流体力学知名学者Toro的教学用源码。
该程序为非饱和均匀沙一维数学计算模型.程序包含以下模块:1)定义主函数、变量、参数、数组;2)新建(打开)文件、读取数据;3)定义初始边界条件;4)天数循环(包括水流计算、泥沙计算分别的断面循环);5)定义函数、子程序。
通过黄平编写的润滑数值计算方法,总结一些关于滑动轴承的动压效应、挤压效应、能量方程,弹流方程等等运算规律,有兴趣的可以借鉴吸收,互相帮助学习。主要介绍面接触滑块数值计算方法
已知单层板正轴刚度求层合板刚度ABD,已知单层板正轴刚度求层合板刚度ABD,已知单层板正轴刚度求层合板刚度ABD,已知单层板正轴刚度求层合板刚度ABD,已知单层板正轴刚度求层合板刚度ABD,已知单层板正轴刚度求层合板刚度ABD.
卫星精密定轨是一项非常复杂的任务,由于受到各种外在因素的影响,为了保持卫星在预定的轨道运行,必须进行轨道机动调整,程序主要为卫星机动模块,采用fortran语言编写。
采用间断有限元方法求解 线性方程u_t+u_x=0, x in [0,2*pi],初始条件u(x,0)=sin(x),
该程序为利用追赶法原理求解三对角矩阵的程序示例,非常容易理解,便于各位同行根据需要进行扩展,尤其是对于电磁、动力学以及数理领域,对于时域有限差分法中的ADI-FDTD、HIE-FDTD非常有用
单相材料再结晶代码,deform二次开发再结晶代码例子。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。。
本程序可用做水文资料、气象资料自动分割,内含部分资料示例,使用者可参看数据格式替换掉数据文件进行操作,可以有效减少时间
sod激波管是经典的黎曼问题,计算流体力学基础问题之一。该例子使用一维结构网格,采用有限差分法,时间积分上应用三阶荣格库塔时间积分,求得激波管中密度,速度,压力的变化曲线。
利用fortran编写,采用压气机叶片优化设计常用的流线曲率法,对给定条件的压气机叶片进行优化设计。流道中沿径向共设置6个计算站,计算站方向沿径向,相邻计算站之间的距离为0.16米,可对沿不同叶高的数据进行测量。
根据6SV2.1,在原有基础上,添加HJ1B-CCD1的光谱响应函数,完成编译过程,生成sixsV2.1.EXE 所需文件 6s源码:http://6s.ltdri.org/pages/downloads.html 是2015年发布的V2.1版本。 所需卫星的光谱响应函数文件 HJ-1B卫星的CCD相机的光谱响应函数,可以在 http://www.cresda.com/site1/Downloads/gpxyhs/index.shtml 下载,这个网页还有GF、ZY等卫星的光谱响应函数文件。
Fortran 90是根据目前最新的Fortran 语法标准所建立,它比从前的Fortran 77的标 准语法有更多的优点,虽然对我们来说,或许需要花一点时间来习惯这种新的语法,但由 于它仍然是 Fortran,对以前的程序员而言,这并不是一件难事,更何况一种语言的更新 是必然的趋势 (目前已有人开始讨论Fortran 95该增加及删除那些语法) 。藉由本文之介 绍,相信各位读者对Fortran 90的新语法已有所认识,甚至您已能着手撰写Fortran 90的 程序了,让我们一起来熟悉它
用于计算不同阶数的第一类贝塞尔函数Fortran代码!第一类n阶BESSEL函数bessjfunction bessj(n,x)integer n,IACCreal(8) bessj,x,BIGNO,BIGNIPARAMETER (IACC=40,BIGNO=1.e10,BIGNI=1.e-10)INTEGER j,jsum,mreal(8) ax,bj,bjm,bjp,sum,tox,bessj0,bessj1if(nfloat(n))then tox=2./ax bjm=bessj0(ax) bj=bessj1(ax) do j=1,n-1 bjp=j*tox*bj-bjmbjm=bjbj=bjp end do bessj=bjelse tox=2./ax m=2*((n+int(sqrt(float(IACC*n))))/2) bessj=0. jsum=0 sum=0. bjp=0. bj=1. do j=m,1,-1 bjm=j*tox*bj-bjpbjp=bjbj=bjmif(abs(bj)>BIGNO)then bj=bj*BIGNI bjp=bjp*BIGNI bessj=bessj*BIGNI sum=sum*BIGNIendifif(jsum/=0) s