当前位置：文档之家 › FNL_1x1(grib1)资料认识与应用(GrADS)-兰溪整合版

FNL_1x1(grib1)资料认识与应用(GrADS)-兰溪整合版

  • 格式：pdf
  • 大小：26.17 MB
  • 文档页数：7

FNL 1X1 资料认识与应用兰溪之水整合版2011-12-16 14:34:05FNL 1.0X1.0数据下载地址：/datasets/ds083.2/需要注册通过方可下载。

FNL 1.0X1.0 资料包含的物理量及其缩写缩写参数名称No4LFTXsfc 近地表四层等压面的抬升指数No5WAVAprs 500 hPa等压面位势高度距平No5WAVHprs 500 hPa等压面位势高度ABSVprs 绝对涡度CAPE 对流有效位能CIN 对流抑制能CLWMRprs 云水CWATclm 气柱云水GPAprs 位势高度距平HGT 位势高度HPBLsfc 地表行星边界层高度ICECsfc 海冰密集度LANDsfc 陆地覆盖LFTXsfc 地表抬升指数O3MRprs 臭氧层混合比POTsig995 位温PRE 气压PWATclm 可降水量RH 相对湿度SOILW 土壤体积含水量SPFH 比湿TCDCcvl 对流云总云量TM 温度TOZNEclm 臭氧含量UGRD u分量VGRD v 分量VVEL 垂直速度VWSH 垂直风切变WEASDsfc 累积雪量GrADS处理FNL 1.0X1.0（grib1）数据处理前需要先将grib2ctl.exe放到GrADS安装文件夹GrADS19\win32（1.9版本）或OpenGrADS\Contents\Cygwin\Versions\2.0.a9.oga.1\i686（2.0版本）下，方便操作。

gribmap.exe（这个GrADS自带的）第一步：先生成一个描述文件ctl打开命令提示符，C:\Documents and Settings\Administrator>grib2ctl grib_file > grib_file.ctl如：或者进入GrADS：gs>!grib2ctl grib_file > grib_file.ctl （！表示调用外部的程序，注意路径用“/”）如：这样fnl_20101013_00_00_c.ctl描述文件就生成了！第二步：利用GrADS自带的gribmap.exe生成索引文件：C:\Documents and Settings\Administrator>gribmap –v –i grib_file.ctl > b.txt如：或进入GrADS：这样就会在H:\test路径下生成名为fnl_20101013_00_00_c.idx的索引文件，并会把整个映射过程写入到b.txt文件中，当然这里也可省略查看映射过程，即C:\Documents and Settings\Administrator>gribmap –i grib_file.ctl这样就可以开始用GrADS读取FNL文件画图了！不过我们处理FNL资料的时候一般都是不止一个文件，所以我们就需要用到批处理了！第一步：利用grib2ctl.exe生成初始时刻ncep数据的ctl文件；示例中生成的是fnl_20101013_00_00_c数据的文件。

第二步：将ctl中相应的信息进行修改：a、将dset的数据集文件名写成代换模块格式；如：dset H:/test/fnl_201010%d2_%h2_00_c％y2 代表两位数年％y4 代表四位数年％m1 代表一位或者两位数的月％m2 代表两位数月（用0补齐1位数）％mc 3个字符月份的缩写％d1 1或2位天％d2 两位天％h1 1或者2位时％h2 2位时b、在第一行后添加一行“options template”定义时间范围和增量；c、tdef number LINEAR start increment，将其中的number和increment（ncep时间间隔为6小时，可写成6hr）设成与所拥有的ncep资料一致；（修改后的信息如图所示,里面的路径需自己设置）；第三步：利用gribmap.exe生成已改过后的ctl相对应的idx文件；这样后，我们就可以用GrADS打开多个fnl文件进行数据处理。

############################################################################### ### 一种更简单的方法### ############################################################################## 第一步：第二步：############################################################################################################################################################## ### 还有一种方法### ############################################################################## 就是把多个文件合并成一个，然后再使用grib2ctl.exe 和gribmap.exe生成对应的ctl文件和idx文件：Windows下一般使用copy命令，linux下一般使用cat命令,copy /b file1 file2 > newfile, 文件比较多可以使用通配符：copy /b fnl* > newfile或cat file1 file2 > newfile, 使用通配符：cat fnl* > newfile然后对合并后的newfile用grib2ctl.exe 和gribmap.exe生成对应的ctl文件和idx文件，这个上面的重复，就不详细说明了！附：如果我们只需要fnl文件中的某些变量，我们可以编写gs文件提取所需要时刻的ncep 数据了。

gs如下：'reinit''open H:/test/fnl_201010Megi.ctl''set gxout fwrite''set fwrite H:/test/ Megi.dat''set lat 15 40''set lon 110 145't=1while(t<=60)'set t ' tz=1while(z<=26)'set z ' z'd HGTprs''d RHprs''d TMPprs''d UGRDprs''d VGRDprs''d VVELprs'z=z+1endwhilet=t+1endwhile'disable fwrite''reinit';然后编写描述文件：Megi.ctl还可以写成批处理.bat文件，如：@echo offrem 此程序用于处理NECP FNL资料，并生成.ctl和.idx文件，然后使用grads绘图rem 此程序使用的为Grads2.0echo 正在生成.ctl文件,Please wait ...rem --------------------------------------------------------------------------- pushd C:\cd C:\OpenGrADS\Contents\Cygwin\Versions\2.0.a9.oga.1\i686\grib2ctl H:\test\fnl_20101013_00_00_c > H:\test\fnl_20101013_00_00_c.ctlrem --------------------------------------------------------------------------- echo 正在生成.idx索引文件，请稍后...pushd C:\cd C:\OpenGrADS\Contents\Cygwin\Versions\2.0.a9.oga.1\i686\gribmap -v -i H:\test\fnl_20101013_00_00_c.ctlrem --------------------------------------------------------------------------- pushd H:\test\if exist H:\test\fnl_500hpa.gs (del H:\test\fnl_500hpa.gs) else (echo sorry!there is no file)rem ---------------------------------------------------------------------------rem 下面开始编写GS文件...set f=H:\test\set n=fnl_20101013_00_00_cset l=.ctlecho 'reinit' >> fnl_500hpa.gsecho 'open %f%%n%%l% ' >>H:\test\fnl_500hpa.gsecho 'set grads off' >>H:\test\fnl_500hpa.gsecho 'set grid on' >>H:\test\fnl_500hpa.gsecho 'set gxout contour' >>H:\test\fnl_500hpa.gsecho *自定义颜色>> fnl_500hpa.gsecho 'set rgb 16 51 255 0'>> fnl_500hpa.gsecho 'set rgb 17 255 153 0'>> fnl_500hpa.gsecho 'set rgb 18 255 0 0'>> fnl_500hpa.gsecho 'set rgb 19 255 102 153'>> fnl_500hpa.gsecho 'set rgb 20 102 0 102'>> fnl_500hpa.gsecho 'set rgb 21 51 0 51'>> fnl_500hpa.gsecho *设置绘图区域>> fnl_500hpa.gsecho 'set lat 10 70'>> fnl_500hpa.gsecho 'set lon 20 150'>> fnl_500hpa.gsecho 'set lev 500'>> fnl_500hpa.gsecho *设置等值线>> fnl_500hpa.gsecho *echo 'setcterp on'>> fnl_500hpa.gsecho *echo 'setcsmooth on'>> fnl_500hpa.gsecho 'set clevs 520 524 528 532 536 540 544 548 552 556 560 564 568 572 576 580 584 588 592'>> fnl_500hpa.gsecho *'set ccolor 18'>> fnl_500hpa.gsecho 'set cmin 4'>> fnl_500hpa.gsecho 'd HGTprs/10'>> fnl_500hpa.gsecho 'set gxout barb'>> fnl_500hpa.gsecho 'd skip(UGRdprs,3);skip(VGRdprs,3)'>> fnl_500hpa.gsecho *设置地图>> fnl_500hpa.gsecho 'set mpdset cnriver' >> fnl_500hpa.gsecho 'set map 4 1 15' >> fnl_500hpa.gsecho 'draw map' >> fnl_500hpa.gsecho 'set mpdset cnworld' >> fnl_500hpa.gsecho 'set map 15 1 5' >> fnl_500hpa.gsecho 'draw map' >> fnl_500hpa.gsecho 'set mpdset xinjiang' >> fnl_500hpa.gsecho 'set map 14 1 7' >> fnl_500hpa.gsecho 'draw map' >> fnl_500hpa.gsecho 'set mproj scaled' >> fnl_500hpa.gsecho 'set mpdset cnworld' >> fnl_500hpa.gsecho *'draw title 2007/07/18 02:00(BJT) 500hPa HGT' >> fnl_500hpa.gs echo 'printim %f%%n%UTC_hgt_wind.png white x1000 y800' >> fnl_500hpa.gs echo ; >> fnl_500hpa.gs rem --------------------------------------------------------------------------- echo 正在绘图，马上完成，请稍等....pushd C:\cd C:\OpenGrADS\Contents\Cygwin\Versions\2.0.a9.oga.1\i686\rem cd C:\OpenGrADS\Contents\Cygwin\Versions\2.0.a9.oga.1\i686\grads -cl H:\test\fnl_500hpa_hgt_wind.gsC:\OpenGrADS\Contents\Cygwin\Versions\2.0.a9.oga.1\i686\grads -cl H:\test\fnl_500hpa.gs将以上复制到记事本，保存为FNL2Graph4GrADS.bat，双击运行即可！。

合集下载

IDL、GRAD、NCL绘图命令对应关系一览表

IDL、GRAD、NCL绘图命令对应关系一览表

IDL、NCL、GRADS、MATLAB绘图命令对应关系一览表1.grads的数据文件与ncl的什么文件对应？grads只支持按照一定顺序存储的二进制数据文件，后缀名以.grd或者.dat或者.grb结束。

这种数据必须以时间为最外层，然后每个变量按照向量形式存储，每个变量由外向内的存储顺序是高度（或等压面）==》纬度==》经度。

这必须注意，否则画图容易出现一堆一堆的乱线条，这就说明你的数据没有按照grads的要求存储。

再看看ncl，可以说ncl支持绝大多数各种数据的读写，包括netcdf，hdf,以及二进制数据甚至ASCII码(如果说是十进制数据或许你会更熟悉)，前两种数据一般都有头文件，不能用C 语言或者Fortran读取，都需要插件才可以读取，ncl可以直接读取，matlab中也可以直接读取NETCDF格式(.nc)的数据。

grads中可以读取按照说明存取的NC数据，这种数据必须又正确的时间说明，也就是说时间必须是真实的，有些模式模拟出来的数据grads的sdfopen 命令是打不开的，因为一般模式都是nonleap run，都是平年，没有设定闰年，造成了时间说明不真实，grads就会报错。

那么grads'如何使用NC数据呢？所以建议使用ncl转换数据，将NC数据，hdf数据或者十进制数据转换成grd数据，供grads使用。

这样说明是在是太空洞了，那么下面我举个例子吧。

eg1) 使用grads将netcdf数据转换成grd数据'reinit'var.1='air';var.2='hgt';var.3='uwnd'var.4='vwnd';var.5='omega';var.6='shum'j=6while (j<=6)'set fwrite/disk3/users/Rao_Jian/ERA-Interim-daily/entropy/'var.j'.daily.1979-2010.grd' 'set gxout fwrite'i=1979while (i<=2010)'sdfopen/disk3/users/lbq/ERA-Interim-daily/pressure/'var.j'.interim.'i'.nc'tt=1if(i=1980|i=1984|i=1988|i=1992|i=1996|i=2000|i=2004|i=2008)while (tt<=366)'set t 'ttzz=1while (zz<=37)'set x 1 240''set y 1 121''set z 'zz'd 'var.j''zz=zz+1endwhilett=tt+1endwhileelsewhile (tt<=365)'set t 'ttzz=1zz=1while (zz<=37)'set x 1 240''set y 1 121''set z 'zz'd 'var.j''zz=zz+1endwhilett=tt+1endwhileendifi=i+1name='/disk3/users/lbq/ERA-Interim-daily/pressure/'var.j'.interim.'i'.nc''close 1'endwhile'disable fwrite'j=j+1endwhileeg2)使用ncl将数据输出成二进制数据，grads可以使用(截取部分)patho = "/disk3/users/Rao_Jian/Hadley/"system("rm "+patho+"ev_ts.grd")system("rm "+patho+"ev_ts2.grd")system("rm "+patho+"ev_ts.txt")system("rm"+patho+"ev_ts2.txt")do nt=0,719fbindirwrite(patho+"ev_ts.grd",ev_ts(:,nt));写成二进制数据，供grads使用end dofbindirwrite(patho+"ev_ts2.grd",ev_ts(time|:,evn|:))asciiwrite(patho+"ev_ts.txt",ev_ts(time|:,evn|:))；写成十进制数据，可以贴到EXCEL中使用asciiwrite(patho+"ev_ts2.txt",ev_ts);print(ev_ts(0,:))此外ncl中还有其他的读写函数，如fbinrecread,fbinrecwrite,fbinwrite,fbinreadncl读取netcdf3/4、hdf4、grib-1/2也是通过文件操作的eg3)ncl支持直接读取格式文件的读法fi = addfile(pathi+"HadISST_sst.nc","r")sst0 =fi->sst(960:1679,:,:) ; load 50 year data duringfrom 1950 to 2009注意：这类文件的后缀名一般为.nc /.hdf/ .grb/.hdfeos/.ccm2.grads中的描述文件与ncl中的什么对应描述文件（.ctl)是一个纯文本文件，我们有数据grads还是不能出图，需要一个描述文件来指定他的存储数据个数，维度（时间长度、层数和经纬度信息）。

(完整word版)GMS-全面+经典教程

(完整word版)GMS-全面+经典教程

地下水模拟软件GMS 教程翻译目录1 GIS模块 (6)1.1简介 (7)1.2未安装ESRI ArcObjects (7)1.3开始 (8)1.4读取Shapefile文件 (8)1.5查看Shapefile文件 (8)1.6查看属性表 (9)1.7文件转换为2D离散点 (9)1.8已安装ESRI ArcObjects (11)1.9启动ArcObjects (11)1.10读出数据 (11)1.11查看shapefile文件 (12)1.12图形文件转化为属性对象 (13)2 属性对象 (17)2.1 目标 (17)2.1简介 (17)2.2开始 (18)2.3属性对象 (18)2.4结论 (24)3 MODFLOW—概念模型法 (25)3.1简介 (26)3.2问题描述 (26)3.3开始 (28)3.4导入模型底图 (28)3.5存储工程文件 (28)3.7定义边界 (29)3.8创建本地源汇图层 (30)3.9定义补给区域 (36)3.11 Locating the Grid Frame定位栅格骨架 (40)3.12 创建栅格 (40)3.13初始化MODFLOW数据 (41)3.14定义活动/不活动区域 (41)3.15插值（篡改）高程值 (41)3.16转换概念模型 (44)3.17检查模拟 (45)3.18运行MODFLOW (45)3.19显示等水位线 (45)3.20侧边浏览水头 (46)3.21查阅水流均衡 (46)3.22结论 (47)4 MODFLOW-Grid Approach (48)4.1 简介 (49)4.2问题描述 (49)4.3开始 (50)4.4单位 (50)4.5创建栅格 (51)4.6创建modflow模拟 (51)4.7单元边界条件直接赋值 (54)4.8 LPF模块 (55)4.9补给（回灌）模块 (57)4.10排水沟模块 (57)4.11井模块 (59)4.12检查模拟 (62)4.13保存模拟 (63)4.14 运行MODFLOW (63)4.15显示结果 (63)4.17结论 (66)5 MODFLOW—层数据插值 (67)5.1简介 (68)5.2开始 (68)5.3层数据插值 (68)5.4案例 (69)5.5例1—完整连续地层 (69)5.6案例2 —裂隙 (71)5.7案例3—出露 (73)5.8案例4—基岩消蚀 (74)5.9 结论 (76)6 MODPATH (78)6.1 简介 (78)6.2 例题描述 (79)6.3开始 (79)6.4导入工程文件 (79)6.5定义孔隙度 (80)6.6定义初始位置 (80)6.7显示设置 (81)6.8粒子示踪属性设置 (81)6.9垃圾填埋场的粒子示踪 (83)6.10 通过区域代码设置颜色 (85)6.11流线长度和时间计算 (85)6.12通过区域代码设定截获区 (86)6.13小结 (86)7 MT3D-概念模型方法 (88)7.1 简介 (88)7.2案例描述 (89)7.4导入工程 (89)7.5确定单位 (89)7.6初始化MT3DMS模拟 (90)7.7 定义含水层特征 (92)7.8 入渗浓度的赋值 (93)7.9概念模型转换 (93)7.10层厚设置 (94)7.11对流包 (94)7.12离散包 (94)7.13查阅结果 (94)7.14源汇混合包 (95)7.15 保存模拟 (95)7.16运行模型 (95)7.17运行MT3DMS (95)7.18查阅结果 (96)7.19物质—时间图 (96)7.20动画显示 (97)7.21吸附和分解的模拟 (97)7.22运行选项 (98)7.23保存模型 (99)7.23.1运行MT3DMS (99)7.24 查看计算结果 (99)7.25随时间变化历史图 (99)7.26小结 (100)1 GIS模块Shapefile文件导入、显示以及转换目标无须使用ArcObjects实现Shapgefile文件的导入和显示。

WRF模式简易操作中文指南

WRF模式简易操作中文指南

WRF 模式操作指南The Institute of Atmospheric Physics， Chinese Academy of Sciences Northeast Institute of Geography and Agroecology, Chinese Academy of Sciences中国科学院大气物理研究所中国科学院东北地理与农业生态研究所二○一七年三月二十日目录1. WRF模式简介 (1)2. WRF模式的安装 (2)2.1 安装环境 (2)2.2 模式源程序 (2)2.3 NetCDF函数库的安装 (2)2.4 标准初始化（SI）的安装 (6)2.5 WRF模式的安装 (9)3. WRF模式与T213模式嵌套 (17)3.1 嵌套方案 (17)3.2 嵌套程序设计 (17)3.3编译嵌套程序 (21)3.4 嵌套的实现 (22)4. WRF模式系统的运行 (29)4.1 理想大气方案 (29)4.2 真实大气方案 (32)5. WRF模式系统作业卡 (47)5.1 源程序 (47)5.2 真实大气方案 (48)6. 模式结果的显示处理 (61)6.1 Vis5D格式 (61)6.2 MICAPS格式 (62)6.2 GrADS格式 (65)附录1. WRF模式参数配置说明 (68)附录2. T213场库参数表 (78)WRF模式系统安装/调试技术报告1. WRF模式简介WRF(Weather Research Forecast)模式系统是由许多美国研究部门及大学的科学家共同参与进行开发研究的新一代中尺度预报模式和同化系统。

WRF模式系统的开发计划是在1997年由NCAR中小尺度气象处、NCEP的环境模拟中心、FSL 的预报研究处和奥克拉荷马大学的风暴分析预报中心四部门联合发起建立的，并由国家自然科学基金和NOAA共同支持。

现在，这项计划，得到了许多其他研究部门及大学的科学家共同参与进行开发研究。

FN1在乳腺癌组织中的表达及其与临床病理的关系

FN1在乳腺癌组织中的表达及其与临床病理的关系

FN1在乳腺癌组织中的表达及其与临床病理的关系陈皓;黄秋霞;韦康来;冯振博【摘要】目的:探讨乳腺癌组织中纤维连接蛋白1(FN1)的表达与临床病理的关系.方法:通过免疫组化SupervisionTM二步法检测79例乳腺癌组织和其中相应的40例癌旁乳腺组织中FN1蛋白的表达.结果:①FN1在乳腺癌组织中的表达明显高于癌旁乳腺组织(P＜0.01).②FN1表达水平随组织学分级、肿瘤大小、淋巴结转移、临床分期的等级升高而增加(均P＜0.05);FN1在HER2阳性中高表达(P＜0.05),在ER/PR阳性中低表达(P＜0.01);FN1的表达在分子亚型中存在差异(P＜0.01),在HER2过表达型中明显高于Luminal型(P＜0.01).③FN1诊断乳腺癌组织的ROC 曲线下面积为0.927,诊断淋巴结转移的ROC曲线下面积为0.631.结论:FN1异常表达参与乳腺癌发生发展过程.【期刊名称】《广西医科大学学报》【年(卷),期】2015(032)002【总页数】4页(P241-244)【关键词】乳腺癌;FN1;分子亚型;转移;免疫组化【作者】陈皓;黄秋霞;韦康来;冯振博【作者单位】广西医科大学第一附属医院病理科南宁530021;广西医科大学第一附属医院病理科南宁530021;广西医科大学第一附属医院病理科南宁530021;广西医科大学第一附属医院病理科南宁530021【正文语种】中文【中图分类】R737.9乳腺癌（breast cancer，BC）已成为目前全球女性中发病率最高的恶性肿瘤，BC 细胞的转移已成为病人复发和死亡的重要原因。

纤维连接蛋白1（Fibronectin 1，FN1）是一种重要的胞外基质（extracellular matrix，ECM）分子，参与了基质重塑，细胞黏附和迁移过程，进而通过调节肌动蛋白聚合而影响细胞的运动［1］。

FN1 在BC 中的表达文献报道并不一致。

多数研究表明FN1 在BC 中高表达［1，2］，但Tiezzi等［3］及Rosman等［4］研究发现FN1在BC细胞系中表达下调。