Corsica调研与工作进展报告人：邢哲导师：肖炳甲研究员ASIPP报告.PPT

Corsica调研与工作进展 报告人：邢哲 导师：肖炳甲 研究员 报告内容 1、Corsica简介 2、BASIS初步 3、Corsica构建平衡 4、Corsica应用调研 5、工作进展和计划 一、Corsica简介 Corsica是由美国劳伦斯利弗莫尔国家实验室（LLNL）开发的一个综合的托卡马克放电模拟程序。开发的初衷是通过Corsica实现对托卡马克中所有物理过程的综合模拟，分析放电试验中的物理现象，优化放电控制，也可以用于优化托卡马克设计，比如ITER Corsica代码耦合二维平衡计算和一维输运模拟，可以使用自由边界平衡解或者固定边界平衡解模拟放电平衡演化。目前的Corsica III 加入了二维的边界输运和湍流输运模型 Corsica包含的程序包可以数值计算等离子体平衡、模拟等离子体中的输运（如芯部输运、中心-边界输运、湍流输运等）、分析MHD稳定性、辅助加热等等。 Corsica的程序包还包括常用数学计算、多项式拟合、Bessel函数、快速傅里叶变换、图形绘制、BASIS开发环境等。 Corsica中还预留与Mdsplus和外部控制器的接口 二、BASIS初步 Corsica的一个优点是代码的适应性好，用户可以通过交互界面自由的定义变量，编译程序(使用体验类似Matlab)；甚至可以根据需要向Corsica加入编译好的程序包，方便国际合作。 这一优势源于Corsica基于BASIS开发环境，BASIS提供了交互式语言解释器，内部包括大量嵌入的函数、错误恢复、宏、存取变量、标准化I/O等。 BASIS语言是一种解释性语言，语法类似Fortran90语言。用户可以自定义变量和函数，之后保存为一个脚本文件.bas(类似Matlab中.m文件) 举例：test.bas # caltrans test.bas read tokamak.bas # Variables used by weqdsk... shotName=probname;shotTime=1 # Create a nominal direct-solve eq from scratch... ds(trim(probname)//.inp) # Converge to q(0)=1.1, q(a)=2.9, betaN=1... probid=DCON //trim(probname)// direct-solve nctot=3 vo=[ctroy,qsrf(1),qsrf(msrf)] vo0=[1,1.1,2.9] vi=[betaj,alfa(0),plcm] x0=[betaj,alfa(0),plcm] ihy=20; run # Save the direct-solve eq... saveq(trim(probname)//.sav) # Make some plots... layout;profiles # Write direct-solve EQDSK files... weqdsk(ag) # Standard a- and g-files quit Corsica中读写数据(1) 使用标准I/O设备读写文本文件 读入data1.dat： 写入data2.dat： integer n=100 character*1 tab=char(9) real big_array(2,n) io=basopen(data2.dat,w) integer io=basopen(data1.dat,r) io time current io big_array do $i=1,n call basclose(io) io format(time($i),10,3,1) \ real time(n)=big_array(1,) tab format(current($i),16,8,2) real current(n)=big_array(2,) enddo call basclose(io) 使用PFB设备读写二进制文件.pfb: 创