GrADS描述文件(CTL)解析器 (2)

最新推荐文章于 2024-03-15 17:07:34 发布
最新推荐文章于 2024-03-15 17:07:34 发布
阅读量3.1k 收藏 1
点赞数
分类专栏： Java GrADS 文章标签： email date vector import string
版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。
本文链接：https://blog.csdn.net/BJ0773/article/details/1684883
版权
Java 同时被 2 个专栏收录
27 篇文章 0 订阅
订阅专栏
GrADS
2 篇文章 0 订阅
订阅专栏
gradsCtl.java
01  /**
02    * PACKAGE      : cma.common.grads
03    * DESCRIPTION  :
04    * AUTHOR       : 刘泽军
05    * EMAIL        : BJ0773@gmail.com
06    * Date         : 2004-12-01
07    * Time         : 22:39:59
08    * Relation     : gradsDim.java, gradsCtl.java, gradsVar.java, grads.java
09    * Compile      : javac -d . gradsDim.java gradsCtl.java gradsVar.java grads.java
10    */
11 
12  package  cma.common.grads;
13 
14  import  java.util.Vector;
15 
16  public class  gradsCtl  {
17       public   String  dset;
18       public   String  index;       //暂未支持
19       public   String  dtype;       //暂未支持
20       public   String  options;
21       public   String  title;
22       public   String  undef;
23 
24       public   String      pdef;    //暂未支持
25 
26       public   gradsDim    xdef;
27       public   gradsDim    ydef;
28       public   gradsDim    zdef;
29       public   gradsDim    tdef;
30 
31       public   Vector  vars; //Object of class gradsVar
32       public  gradsCtl () {
33           xdef    =  new  gradsDim () ;
34           ydef    =  new  gradsDim () ;
35           zdef    =  new  gradsDim () ;
36           tdef    =  new  gradsDim () ;
37           vars    =  new  Vector () ;
38       }
39  }

 

gradsVar.java
01  /**
02    * PACKAGE      : cma.common.grads
03    * DESCRIPTION  :
04    * AUTHOR       : 刘泽军
05    * EMAIL        : BJ0773@gmail.com
06    * Date         : 2004-12-01
07    * Time         : 22:40:26
08    * Relation     : gradsDim.java, gradsCtl.java, gradsVar.java, grads.java
09    * Compile      : javac -d . gradsDim.java gradsCtl.java gradsVar.java grads.java
10    */
11 
12  package  cma.common.grads;
13 
14  public class  gradsVar  {
15       public   String  name        =  "" ;
16       public   int      levels      =  0 ;
17       public   String  reserve     =  "" ;
18       public   String  contents    =  "" ;
19       public   String  unicode     =  "" ;
20       public  gradsVar () {
21       }
22       public   String  convertToString () {
23           return (  name +  " "  + String.valueOf ( levels )  +  " "  + reserve +  " "  + contents ) ;
24       }
25  }

 

gradsDim.java
001  /**
002    * PACKAGE      : cma.common.grads
003    * DESCRIPTION  :
004    * AUTHOR       : 刘泽军
005    * EMAIL        : BJ0773@gmail.com
006    * Date         : 2004-12-01
007    * Time         : 22:40:19
008    * Relation     : gradsDim.java, gradsCtl.java, gradsVar.java, grads.java
009    * Compile      : javac -d . gradsDim.java gradsCtl.java gradsVar.java grads.java
010    */
011 
012  package  cma.common.grads;
013 
014  public class  gradsDim  {
015       public   String  symbol  =  "" ;
016       public   int      count   =  0 ;
017       public   String  type    =  "" ;
018       public   String  start   =  "" ;
019       public   String  increment   =  "" ;
020       public   String  value []      =  { "" } ;
021       public  gradsDim () {
022       }
023       public  boolean  extractValue ( String str []) {
024           if (  str.length <  5  ) {
025               return ( false ) ;
026           }
027           symbol  = str [ 0 ] ;
028           count   = Integer.parseInt ( str [ 1 ]) ;
029           value   =  new  String [ count ] ;
030           for ( int  i= 0 ;i<count;i++ ) {
031               value [ i ]     =  "" ;
032           }
033           type    = str [ 2 ] ;
034           try  {
035               if (  0  == symbol.compareToIgnoreCase ( "tdef" ) ) {
036                   start       = str [ 3 ] ;
037                   increment   = str [ 4 ] ;
038                   int  year=- 1 , month=- 1 , day= 1 , hour= 0 , min= 0 ;
039                   String  monStr []  =  { "Jan" ,  "Feb" ,  "Mar" ,  "Apr" ,  "May" ,  "Jun" ,  "Jul" ,  "Aug" ,  "Sep" ,  "Oct" ,  "Nov" ,  "Dec" } ;
040                   String  mon =  "" ;
041                   String  tmp = str [ 3 ] .substring ( 0 ,str [ 3 ] .length () - 7 ) .toUpperCase () ;
042                   year        = Integer.parseInt ( str [ 3 ] .substring ( str [ 3 ] .length () - 4 )) ;
043                   mon         = str [ 3 ] .substring ( str [ 3 ] .length () - 7 ,str [ 3 ] .length () - 4 ) ;
044                   for ( int  m= 1 ;m<= 12 ;m++ ) {
045                       if (  0  == monStr [ m- 1 ] .compareToIgnoreCase ( mon ) ) {
046                           month   = m;
047                           start   = str [ 3 ] .replaceAll ( mon, monStr [ m- 1 ]) ;
048                       }
049                   }
050                   if (  0  > year ||  9999  <= year ||  1  > month ||  12  < month  ) {
051                       return ( false ) ;
052                   }
053                   try  {
054                       if (  - 1  != tmp.indexOf ( "Z" ) ) {
055                           day = Integer.parseInt ( tmp.substring ( tmp.indexOf ( "Z" ) + 1 )) ;
056                           tmp = tmp.substring ( 0 , tmp.indexOf ( "Z" )) ;
057                       }
058                       if (  - 1  != tmp.indexOf ( ":" ) ) {
059                           hour    = Integer.parseInt ( tmp.substring ( 0 ,tmp.indexOf ( ":" ))) ;
060                           min     = Integer.parseInt ( tmp.substring ( tmp.indexOf ( ":" ) + 1 )) ;
061                       }
062                       else  {
063                           hour    = Integer.parseInt ( tmp ) ;
064                       }
065                   }
066                   catch ( Exception e ) {
067                   }
068                   value [ 0 ]     =
069                       String.valueOf ( year + 10000 ) .substring ( 1 )  +
070                       String.valueOf ( month+ 100   ) .substring ( 1 )  +
071                       String.valueOf ( day  + 100   ) .substring ( 1 )  +
072                       String.valueOf ( hour + 100   ) .substring ( 1 ) ;
073                   int  intIncrement    =  0 ;
074                   String  suffix =  "" ;
075                   if (  str [ 4 ] .toLowerCase () .endsWith ( "mn" ) ) {
076                       intIncrement   = Integer.parseInt ( str [ 4 ] .substring ( 0 ,str [ 4 ] .toLowerCase () .indexOf ( "mn" ))) ;
077                       if (  intIncrement >=  60 ) {
078                           intIncrement    = intIncrement /  60 ;
079                           suffix  =  "hr" ;
080                       }
081                       else  {
082                           suffix  =  "mn" ;
083                       }
084                   }
085                   else if (  str [ 4 ] .toLowerCase () .endsWith ( "hr" ) ) {
086                       intIncrement   = Integer.parseInt ( str [ 4 ] .substring ( 0 ,str [ 4 ] .toLowerCase () .indexOf ( "hr" ))) ;
087                       if (  intIncrement >=  24 ) {
088                           intIncrement    = intIncrement /  24 ;
089                           suffix  =  "dy" ;
090                       }
091                       else  {
092                           suffix  =  "hr" ;
093                       }
094                   }
095                   else if (  str [ 4 ] .toLowerCase () .endsWith ( "dy" ) ) {
096                       intIncrement   = Integer.parseInt ( str [ 4 ] .substring ( 0 ,str [ 4 ] .toLowerCase () .indexOf ( "dy" ))) ;
097                       suffix  =  "dy" ;
098                   }
099                   else if (  str [ 4 ] .toLowerCase () .endsWith ( "mo" ) ) {
100                       intIncrement   = Integer.parseInt ( str [ 4 ] .substring ( 0 ,str [ 4 ] .toLowerCase () .indexOf ( "mo" ))) ;
101                       suffix  =  "mo" ;
102                   }
103                   else if (  str [ 4 ] .toLowerCase () .endsWith ( "yr" ) ) {
104                       intIncrement   = Integer.parseInt ( str [ 4 ] .substring ( 0 ,str [ 4 ] .toLowerCase () .indexOf ( "yr" ))) ;
105                       suffix  =  "yr" ;
106                   }
107                   else  {
108                       try  {
109                           intIncrement   = Integer.parseInt ( str [ 4 ]) ;
110                           suffix  =  "" ;
111                       }
112                       catch ( Exception e ) {
113                       }
114                   }
115                   increment   = String.valueOf ( intIncrement )  + suffix;
116                   for ( int  j= 0 ;j<count;j++ ) {
117                       value [ j ]     = String.valueOf ( 1000000 +j*intIncrement ) .substring ( 1 ) ;
118                       while (  value [ j ] .length ()  >  3  && value [ j ] .startsWith ( "0" ) ) {
119                           value [ j ]     = value [ j ] .substring ( 1 ) ;
120                       }
121                       value [ j ]     = value [ j ]  + suffix;
122                   }
123                   if (  1  == count  ) {
124                       value [ 0 ]     =  "anl" ;
125                   }
126               }
127               else  { //xdef ydef zdef
128                   start   = str [ 3 ] ;
129                   if (  0  == type.compareToIgnoreCase ( "LINEAR" ) ) {
130                       increment   = String.valueOf ( Float.parseFloat ( str [ 4 ])) ;
131                       for ( int  j= 0 ;j<count;j++ ) {
132                           value [ j ]     = start + Float.parseFloat ( increment )  * j;
133                       }
134                   }
135                   else if (  0  == type.compareToIgnoreCase ( "LEVELS" ) ) {
136                       if (  str.length -  3  != count  ) {
137                           return ( false ) ;
138                       }
139                       for ( int  j= 0 ;j<count;j++ ) {
140                           value [ j ]     = String.valueOf ( Float.parseFloat ( str [ 3 +j ])) ;
141                       }
142                   }
143                   else  {
144                       return ( false ) ;
145                   }
146                   if (  symbol.equalsIgnoreCase ( "zdef" ) ) {
147                       //System.out.println(symbol + " ====================== " + String.valueOf(count));
148                       for ( int  i= 0 ;i<count;i++ ) {
149                           //System.out.println("LEVEL = " + value[i]);
150                           if (  - 1  != value [ i ] .indexOf ( "." )  && !value [ i ] .endsWith ( "." )  &&
151                               0  == Integer.parseInt ( value [ i ] .substring ( value [ i ] .indexOf ( "." ) + 1 )) ) {
152                               value [ i ]     = value [ i ] .substring ( 0 ,value [ i ] .indexOf ( "." )) ;
153                               //System.out.println(" ====>>>>  " + value[i]);
154                           }
155                       }
156                   }
157               }
158               return ( true ) ;
159           }
160           catch ( Exception e ) {
161               return ( false ) ;
162           }
163       }
164       public  String convertToString () {
165           String  s = symbol +  " "  + String.valueOf ( count )  +  " "  + type;
166           if (  0  == type.compareToIgnoreCase ( "levels" ) ) {
167               for ( int  i= 0 ;i<count;i++ ) {
168                   s   = s +  "/n"  + value [ i ] ;
169               }
170           }
171           else  {
172               s   = s +  " "  + start +  " "  + increment;
173           }
174  /*
175           if( 0 == symbol.compareToIgnoreCase("tdef") ) {
176               String  s1  = "";
177               for(int i=0;i<count;i++) {
178                   s1  = s1 + "/n" + value[i];
179               }
180               s   = s + s1;
181           }
182  */
183           return ( s ) ;
184       }
185       public static  String newDate ( String yyyymmddhh,  int  incHours ) {
186  /*
187  功能：计算yyyymmddhh加上incHours个小时后对应的日期，仍以yyyymmddhh的形式返回，出错时返回YYYYMMDDHH。
188  作者：刘泽军
189  更新：2004-11-24
190  */
191           if (  10  != yyyymmddhh.length () ) {
192               return ( "YYYYMMDDHH" ) ;
193           }
194           int  yyyy, mm, dd, hh; //开始日期
195           int  yyyy1, mm1, dd1, hh1; //结束日期
196           try  {
197               yyyy    = Integer.parseInt ( yyyymmddhh.substring ( 0 , 4 )) ;
198               mm      = Integer.parseInt ( yyyymmddhh.substring ( 4 , 6 )) ;
199               dd      = Integer.parseInt ( yyyymmddhh.substring ( 6 , 8 )) ;
200               hh      = Integer.parseInt ( yyyymmddhh.substring ( 8 )) ;
201           }
202           catch ( Exception ex ) {
203               return ( "YYYYMMDDHH" ) ;
204           }
205           yyyy1   = yyyy;
206           mm1     = mm;
207           dd1     = dd;
208           hh1     =  ( hh+incHours )  %  24 ;
209           int  days    =  ( hh + incHours - hh1 )  /  24 ;
210           if (  hh1 <  0  ) {
211               days    = days -  1 ;
212               hh1     =  24  + hh1;
213           }
214           //                   12   1   2   3   4   5   6   7   8   9  10  11  12
215           int  daysOfMonth []  =  { 31 ,  31 ,  28 ,  31 ,  30 ,  31 ,  30 ,  31 ,  31 ,  30 ,  31 ,  30 ,  31 } ;
216 
217           if ( ( yyyy %  4 )  ==  0  &&  ( ( yyyy %  100 )  !=  0  ||  ( yyyy %  400 )  ==  0  ) ) { //闰年
218               daysOfMonth [ 2 ]   =  29 ;
219           }
220           if (  yyyy <  0  || yyyy + days <  0  ||
221               yyyy >=  9999  || yyyy + days >=  9999  ||
222               mm <  1  || mm >  12  ||
223               dd <  1  || dd > daysOfMonth [ mm ] ) {
224               return ( "YYYYMMDDHH" ) ;
225           }
226 
227           if (  days >=  0  ) {
228               for ( int  i= 0 ;i<days;i++ ) {
229                   if (  dd1 < daysOfMonth [ mm1 ] ) {
230                       dd1 = dd1 +  1 ;
231                   }
232                   else if  (  mm1 <  12  ) {
233                       mm1 = mm1 +  1 ;
234                       dd1 =  1 ;
235                   }
236                   else  {
237                       yyyy1   = yyyy1 +  1 ;
238                       mm1     =  1 ;
239                       dd1     =  1 ;
240                       if ( ( yyyy1 %  4 )  ==  0  &&  ( ( yyyy1 %  100 )  !=  0  ||  ( yyyy1 %  400 )  ==  0  ) ) { //闰年
241                           daysOfMonth [ 2 ]   =  29 ;
242                       }
243                       else  {
244                           daysOfMonth [ 2 ]   =  28 ;
245                       }
246                   }
247               }
248           }
249           else  {
250               for ( int  i=days;i< 0 ;i++ ) {
251                   if (  dd1 >  1  ) {
252                       dd1 = dd1 -  1 ;
253                   }
254                   else if  (  mm1 >  1  ) {
255                       mm1 = mm1 -  1 ;
256                       dd1 = daysOfMonth [ mm1 ] ;
257                   }
258                   else  {
259                       yyyy1   = yyyy1 -  1 ;
260                       mm1     =  12 ;
261                       dd1     =  31 ;
262                       if ( ( yyyy1 %  4 )  ==  0  &&  ( ( yyyy1 %  100 )  !=  0  ||  ( yyyy1 %  400 )  ==  0  ) ) { //闰年
263                           daysOfMonth [ 2 ]   =  29 ;
264                       }
265                       else  {
266                           daysOfMonth [ 2 ]   =  28 ;
267                       }
268                   }
269               }
270           }
271           return (
272                   String.valueOf ( 10000 +yyyy1 ) .substring ( 1 )  +
273                   String.valueOf ( 100 +mm1 ) .substring ( 1 )  +
274                   String.valueOf ( 100 +dd1 ) .substring ( 1 )  +
275                   String.valueOf ( 100 +hh1 ) .substring ( 1 )
276           ) ;
277       }
278  }
BJ0773
关注
  • 0
    点赞
  • 踩
  • 1
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
专栏目录
Grads使用手册
02-27
中科院大气物理所 大气科学和地球流体力学数值模拟国家重点实验室 必备教材
g2ctl:grads绘图应用 CTL文件转化
01-26
g2ctl:grads绘图应用 CTL文件转化
ctl 的写法_GrADSctl文件编写
weixin_39749501的博客
12-21 1844
DSETdata_filenamebacktotopThisentryspecifiesthefilenameofthedatafilebeingdescribed.Ifthedataandthedescriptorfilearenotinthesamedirectory,thendata_filenamemustincludeafullpa...
Grads保姆级教学——小乌龟考试看这一篇就够了
m0_60394632的博客
06-12 1万+
临近考试根据上课内容和ppt整理的grads笔记，ppt翻来翻去太难受了，就看这篇叭，过考试绝对没问题
GrADS描述文件(CTL)解析器 (1)
气象万千
07-11 1万+
GrADS.java 0001 /**0002  *0003  * 根据863项目《中国气象应用网格CMAG》要求，编写于2004年11-12月，国防科技大学611重点实验室)0004  *0005  * PACKAGE      : cma.common.grads0006  * DESCRIPTION  : 解析 GrADS 的 CTL 描述文件0007  * AUTHOR       : 刘
xgrads：解析和读取GrADS常用的ctl文件
02-04
此python软件包专门用于解析和读取常用的.ctl文件。 现在，它可以解析各种.ctl文件。 然而，只有常用原始二进制数据集4D可以使用读和返回作为其它类型的二进制数据，如dtype是station或者grib ，可能在将来被支持。...
grads处理多个ctl文件和nc文件.pdf
10-04
为了解决这个问题，可以重新写一个 ctl 描述文件，其文件变量都和已知的 ctl 相同，若原来的 n 文件只是时间不同，那么新描述文件的时间维数是所有原文件的时间的和。同样，若其他维数不同时也用同样的方法处理。 ...
GrADS学习资料：实习一 数据文件的转换及数据描述文件的建立.docx
09-21
在GrADS中，数据描述文件（CTL文件）是用于描述数据的文件格式和结构的重要文件。CTL文件中包含了数据的基本信息，例如数据的文件名、数据类型、数据维数、数据范围等信息。在本实习中，我们需要编写数据描述文件uv....
grads绘图时的河流文件
06-24
主要用于grads绘制图形时，将全国主要大江大河的边界给绘制出来，使得对国内的河流分布究竟在何处，以及河流走向、河流长短等一清二楚
气象ARWpost、grads 等使用的ctl和dat格式的grd 二进制文件
最新发布
cwr888的专栏
03-15 728
气象ARWpost、grads 和 Fortran使用的ctl+dat格式的grd 二进制文件，在Python中可以用xgrads包来读取或者转成nc 格式使用。
GrADS地图 GrADS地图
01-13
GrADS地图GrADS地图GrADS地图
grads 最新版2.0.a7.1 for windows
01-04
grads 最新版2.0.a7.1 for windows
java读取grib2文件
09-28
java读取grib2文件，使用netcdfAll，内附jar包与例子，netcdfAll-4.3.jar为jdk6版本，netcdfAll-4.6.jar为jdk7版本
grads 多时次
12-13
grads 多时次 直接可以用 为气象人提供个方便
grads批处理ctl文件和用cmd生成idx文件的步骤
chht的专栏
01-25 3842
 一：grads批处理涉及idx文件的grib资料时，编辑grads批处理所需要的ctl文件：（一般可以把已存在的第一个时次的ctl文件进行更改后另存为一个新的ctl文件即可，比如本次批处理的原来的ctl名字都是200206xxxx.ctl形式，则可以新建一个200206.ctl作为grads将来要调用的ctl文件））grib200206%d2%h2 时间循环 opti
grads 批量处理nc文件
chht的专栏
01-25 7478
若想要提取从1951-2006年56年nc文件中的某些数据，一个一个处理非常麻烦，这里介绍种较为简易的方法。例如想提取6-8月的位势高度资料。reinitt5=1951*作文件名循环while(t5 set gxout fwrite set fwrite D:/sichuan/hgt1/%t5%.dat sdfopen e:/ncep1/hgt/hgt
ctl 里面pdef解说
weixin_33873846的博客
01-08 382
WRF 模式MM5 模式都是目前从网上可以下载的气象软件，因此在国内经常可以见到。但这两种模式的数据特点数据的水平网格都不是标准的经纬度网格。需要在ctl 文件中加入PDEF 定义说明把这种非标准的数据经过GrADS 内部的计算转换成标准的经纬度网格数据使用。 其中 PDEF 说明：pdef 130 120 lcc 31.0 110. 50.0 63.0 30. 60. 110. 30000...
python处理wrf气象数据_气象编程 | Python3之WRF的投影转换
weixin_35903340的博客
02-10 4313
有一种WRF输出的数据采用兰伯特双标准纬线投影，那么除非刚好需要同样的投影，想对这种数据进行处理的话往往要进行投影转换，WRF应该是有一套工具可以进行相关的处理，比如wrf-python，但是作为并不熟悉wrf、仅仅是使用WRF输出数据的小白，去使用WRF系工具的话学习成本就比较高了，如何用熟悉、更通用的工具实现这一投影转换呢？难道不是设置几个投影参数，用常见的投影相关的包就可以实现了吗？对，问题...
Grads 笔记
热门推荐
weixin_40645816的博客
08-27 1万+
注：2012-2016年抄的，已经不用这软件，停止更新 Grads绘图程序几例 例1 （相关及信度检验；中国160站插值绘图） 中国16站插值绘图；打开多个文件；设置画布大小；等值线阴影叠加绘图；改等值线间隔，加粗； 部分数值绘制阴影，自定义颜色；调用长江黄河、中国底图；不打印grads标注；填加标题 'reinit' 'open D:\zother\grid.ctl' 'ope...
GrADS学习资料：第7章 GrADS 描述语言及编程(1) .ppt
05-07
描述语言是GrADS软件包中自带的高级语言，在*.gs文件中可以包括描述语句和指令性语句。描述语句不用引号，需要经过编译，而指令性语句则需要用引号括起来。通过描述语言，如if/else/endif和while/endwhile语句块，...

“相关推荐”对你有帮助么？

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

个

红包个数最小为10个

元

红包金额最低5元

当前余额3.43元 前往充值 >
需支付：10.00元
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付元
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。

余额充值