#转自论坛日志# grad用fwrite生成的二进制文件大小和ncl用fbinwrite生成的不一样,相差有8个字节,参考了以下网页,问题出在顺序存取和直接存取的差异. 以下是我的理解,欢迎指正:   Fortran里面,二进制文件的存储方式有sequential(顺序存取)和direct(直接存取)两种.其中,"顺序存取"不需要每次存取的数据(也称"记录“)的长度一致,而是在每笔记录两端加上了记录该记录所占长度的标识符.而"直接存取"要求每一笔记录的长度…

GrADS画特定经线方法 set clevs 23.5; d lat set clevs 120;d lon GrADS 查看多个打开文件的ctl q ctlinfo 1 q ctlinfo 2 ... 2个不同的缺测值 fwrite更改默认缺测值 GrADS用fwrite输出时,默认的缺测值是-9.9900000E+08,如果输入源的缺测值与之不一致,则会出现两种缺测值的情况.解决的办法,可以在fwrite 声明后面设置默认缺测值设置.比如 'reinit' 'open J:/ec79-15/…

经常有身边的学友问到用什么Python发行版比较好? 其实目前比较流行的Python科学计算发行版,主要有这么几个: Python(x,y) GUI基于PyQt,曾经是功能最全也是最强大的,而且是Windows系统中科学免费Python发行版的不二选择.不过今时已不同往昔! PythonXY里面的许多包为了兼容性的问题,无法使用最新的程序包.尤其是令人气愤的是MinGW到现在还是古董级的4.5版本,而TDM-GCC现在都4.8.1-3了.不过这个包在你安装了之后,除了占用较大的磁盘空间之外,基本…

作为一个测绘GIS专业的学生,会有很多绘制地图的需求,一般情况使用ArcGIS.QGIS就可以解决,但是在绘制如论文插图时需要使用更加专业可定制化程度更高的工具,专业传统一般使用GMT,几经比较,最终选择学习basemap,所以有了一下这篇博文. basemap与其他工具比较 工具简介 Basemap: Basemap是 Python 可视化库 Matplotlib 下的一个工具包,主要功能是绘制二维地图,是Python空间数据可视化的重要手段. GMT:即Generic Mapping Too…

一些朋友总是咨询关于二进制文件的读写和转化.这里就我自己的理解说一说. 一).一般问题 二进制文件与我们通常使用的文本文件储存方式有根本的不同.这样的不同很难用言语表达,自己亲自看一看,理解起来会容易得多.因此,我推荐学习二进制文件读写的朋友安装一款十六进制编辑器.这样的编辑器有很多,在我们的 CVF 附带的集成开发环境下就可以(将二进制文件拖动到 IDE 窗口后松开).Visual Studio 2005 也是可以的.(不过需要在 File 菜单下 Open,File) 另外推荐一款使用较多的…

clear,clc nt=735;ny=73;    %2.5*2.5格点的nx=144;    %2.5*2.5格点的f=netcdf('air.mon.mean.nc','nowrite');tt=f{'air'}(:);close(f)fid=fopen('a5.grd','w');  for it=1:nt     for j=1:ny        for i=1:nx         count= fwrite(fid,tt(it,j,i),'float');%根据例子这里用原始的数…

采用fortran进行数值计算的朋友们都应该听说过大名鼎鼎的lapack库,我就不多做介绍了,在此,我只是介绍一个编译好的lapack二进制包ACML(AMD Core Math Library),并且我强烈推荐,有几个原因: 1. ACML是预编译好的二进制文件,省去了许多重新编译的麻烦: 2. AMD对ACML中的许多函数针对x86构架进行了优化,执行效率很高: 3. ACML中有专门针对SMP(不要告诉我你不知道SMP是什么东东)的库,可以实现并行多处理,如果你是多核多线程处理器的话可以最…

使用gfortran将数据写成Grads格式的代码示例: !-----'Fortran4Grads.f90' program Fortran4Grads implicit none integer,parameter::xn=32 !-----经度(Longitude)网格数 integer,parameter::yn=18 !-----纬度(Latitude)网格数 integer,parameter::zn=7 !-----高度层数 integer,parameter::tn=5 !----…

网 上说要分c为主程序和fortran为主程序两种情况讨论,其实我觉得不用,只要你了解生成可执行文件的机制.这个机制就是:不论是单一语言模块之间的 链接还是不同语言之间的混合链接,本质目的都是要链接器能找到定义于其他模块中的符号,如果全部找到,则链接成功,生成可执行的二进制文件. 下面的内容比较基础,看烦了就跳过. 比如简单的一个c程序: /* main.c */ void FOO (); /* 其实按c语言的规定,调用程序之前不需声明, * c编译器会猜测函数的原型.需要在调用之前 * 定义或…

由于项目需要,需要对二进制文件进行读写.转换. 文件说明:由其他程序得到的二进制文件,文件内容为:包含23543个三角形.13270个顶点的三角网所对应的721组流速矢量(u.v)文件,通俗些说,一条数据包含两个双精度型的数值,每组数组包含23543条数据,如果以一个双精度数值为单位,则总共有23543 * 721 * 2 =33,949,006条数据.由Fortran程序以每 8 Byte存储一个数值的二进制文件存储,最终文件大小为下图所示: 测试:从该文件读出数据之后,转换为十进制,存储到另…