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

最新推荐文章于 2024-03-15 17:07:34 发布

BJ0773

最新推荐文章于 2024-03-15 17:07:34 发布

阅读量1w

点赞数

分类专栏： GrADS Java 文章标签： null string file border 脚本 vector

本文链接：https://blog.csdn.net/BJ0773/article/details/1684881

版权

Java 同时被 2 个专栏收录

27 篇文章 0 订阅

订阅专栏

GrADS

2 篇文章 0 订阅

订阅专栏

GrADS.java

0001 /**
0002 *
0003 * 根据863项目《中国气象应用网格CMAG》要求，编写于2004年11-12月，国防科技大学611重点实验室)
0004 *
0005 * PACKAGE      : cma.common.grads
0006 * DESCRIPTION  : 解析 GrADS 的 CTL 描述文件
0007 * AUTHOR       : 刘泽军
0008 * EMAIL        : BJ0773@gmail.com
0009 * Date         : 2004-12-01
0010 * Time         : 21:43:32
0011 * Relation     : gradsDim.java, gradsCtl.java, gradsVar.java, grads.java
0012 * Compile      : javac -d . gradsDim.java gradsCtl.java gradsVar.java grads.java
0013 */
0014
0015 package cma.common.grads;
0016
0017 import java.io.*;
0018 import java.lang.*;
0019 import java.util.*;
0020 import java.text.NumberFormat;
0021 import java.text.DecimalFormat;
0022
0023 import cma.common.dataio.*;
0024
0025 public class GrADS {
0026
0027 //测试
0028      public static void main ( String argv []) throws Exception {
0029
0030          if ( 1 == argv.length ) {
0031              GrADS   g = new GrADS () ;
0032              boolean enabled = g.loadFromFile ( argv [ 0 ]) ;
0033              if ( g.enabled ) {
0034                  enabled = g.writeToFile ( "tmp_GrADS.ctl" ) ;
0035                  String  env [] = { "PATH=$PATH:/usr/local/bin" } ;
0036                  try {
0037                      g.procScript ( env, "/home/grapes/grapes_grads/" , "/home/grapes/Result" , "h" , "500" , "1" ) ;
0038                      g.overlapping   = true ;
0039                      g.procScript ( env, "/home/grapes/grapes_grads/" , "/home/grapes/Result" , "t" , "850" , "1" ) ;
0040                  }
0041                  catch ( Exception e ) {
0042                      System.out.println ( e.toString ()) ;
0043                      e.printStackTrace () ;
0044                  }
0045              }
0046          }
0047      }
0048
0049 //
0050      public   int          width       = 1024 ;      //缺省的图像宽度
0051      public   int          height      = 768 ;       //缺省的图像高度
0052      public   int          border      = 0 ;         //裁掉的边框格点数
0053      public   String      background  = "white" ;   //缺省的图像背景色
0054      public   String      model       = "" ;        //模式名，主要用于图像文件名的前缀
0055
0056      private double       LEFT        = 104.0 ;     //广西区域暂时未用
0057      private double       TOP         = 26.5 ;
0058      private double       RIGHT       = 112.5 ;
0059      private double       BOTTOM      = 20.0 ;
0060
0061      public   double       left        = 104.0 ;
0062      public   double       top         = 26.5 ;
0063      public   double       right       = 112.5 ;
0064      public   double       bottom      = 20.0 ;
0065
0066      public   gradsCtl    ctl         = new gradsCtl () ;
0067      public   boolean      overlapping = false ;     //是否迭加
0068      public   boolean      shading     = true ;      //是否填充阴影
0069
0070      public   String      title       = "" ;
0071
0072      private String      pathname    = "." ;       //当前解释的 CTL 文件所在的目录名
0073      private String      filename    = "" ;        //当前解释的 CTL 文件名
0074      private String      element     = "" ;        //要素
0075      private String      level       = "" ;        //层次
0076      private String      valid       = "" ;        //时次
0077      private String      message     = "" ;        //执行 GrADS 命令时输出的内容
0078      private String      command     = "/public/local/grads/gradsc" ; //
0079      private String      scriptfile  = "gimdr_grads.gs" ;
0080      private boolean      enabled     = false ;
0081
0082      private static   int size        = 4096 ;      //一个 GrADS 的描述文件最大不能超过的字节数
0083 //
0084      public   boolean isEnabled () {
0085          return ( enabled ) ;
0086      }
0087 //获得经向格点数
0088      public int getLongitudeNumber () {
0089          if ( enabled ) {
0090              return ( ctl.xdef.count ) ;
0091          }
0092          return 0 ;
0093      }
0094 //获得纬向格点数
0095      public int getLatitudeNumber () {
0096          if ( enabled ) {
0097              return ( ctl.ydef.count ) ;
0098          }
0099          return 0 ;
0100      }
0101 //获得层次数
0102      public int getLevelNumber () {
0103          if ( enabled ) {
0104              return ( ctl.zdef.count ) ;
0105          }
0106          return 0 ;
0107      }
0108 //获得要素数
0109      public int getElementNumber () {
0110          if ( enabled ) {
0111              return ( ctl.vars.size () ) ;
0112          }
0113          return 0 ;
0114      }
0115 //获得要素的中文说明
0116      public String getElementUnicode () {
0117          return ( getElementUnicode ( element )) ;
0118      }
0119
0120 //获得要素的中文说明
0121      public String getElementUnicode ( String _name ) {
0122          if ( null == _name || "" .equals ( _name.trim ()) ) {
0123              return ( "" ) ;
0124          }
0125          String  eGrapesUnicode [][]   = { //GRAPES
0126              { "u" ,        "风U分量" } ,         { "v" ,        "风V分量" } ,                 { "t" ,        "温度" } ,               { "h" ,            "高度" } ,           { "qv" ,       "比湿" } ,
0127              { "qc" ,       "云水" } ,       { "qr" ,       "雨水" } ,               { "qv" ,       "比湿" } ,               { "q1" ,           "比湿" } ,           { "q2" ,       "比湿" } ,
0128              { "q3" ,       "比湿" } ,       { "q4" ,       "比湿" } ,               { "w" ,        "垂直速度" } ,             { "ps" ,           "地面气压" } ,         { "psl" ,      "海平面气压" } ,
0129              { "rainc" ,    "对流降水" } ,     { "rainnc" ,   "非对流降水" } ,            { "ts" ,       "地面气温" } ,             { "glw" ,          "地面长波辐射" } ,   { "gsw" ,      "地面短波辐射" } ,
0130              { "hfx" ,      "地面热通量" } ,    { "qfx" ,      "地面水汽通量" } ,       { "q2m" ,      "2米比湿" } ,                 { "t2m" ,          "2米温度" } ,             { "u10m" ,     "10米风U分量" } ,
0131              { "v10m" ,     "10米风V分量" } , { "lu" ,       "陆面" } ,               { "zs" ,       "地形" } ,
0132              { "u;v" ,      "风场" } ,       { "uv" ,       "流场" } ,               { "rain" ,     "降水(对流+非对流)" } ,   { "u10m;v10m" ,    "10米风场" }
0133          } ;
0134
0135          String  eMM5Unicode [][] = { //MM5
0136              { "rc" ,       "对流降水" } ,     { "rn" ,       "非对流降水" } ,            { "ter" ,      "地形高度" } ,             { "xlat" ,         "纬度" } ,           { "xlon" ,     "经度" } ,
0137              { "lu" ,       "陆面" } ,       { "t2m" ,      "2米温度" } ,                 { "q2m" ,      "2米比湿" } ,                 { "u10" ,          "10米风U分量" } ,      { "v10" ,      "10米风V分量" } ,
0138              { "pslv" ,     "海平面气压" } ,    { "cref" ,     "" } ,                     { "iclw" ,     "整层云水" } ,             { "irnw" ,         "整层雨水" } ,         { "pwat" ,     "可降水量" } ,
0139              { "clfrlo" ,   "低云量" } ,      { "clfrmi" ,   "中云量" } ,              { "clfrhi" ,   "高云量" } ,              { "u" ,            "风U分量" } ,             { "v" ,        "风V分量" } ,
0140              { "w" ,        "垂直速度" } ,     { "pp" ,       "" } ,                     { "t" ,        "温度" } ,               { "q" ,            "比湿" } ,           { "clw" ,      "云水" } ,
0141              { "rnw" ,      "雨水" } ,       { "rtnd" ,     "辐射趋势" } ,             { "z" ,        "地表高度" } ,             { "h" ,            "高度" } ,           { "td" ,       "露点温度" } ,
0142              { "rh" ,       "相对湿度" } ,     { "vor" ,      "涡度" } ,               { "div" ,      "散度" } ,
0143              { "u;v" ,      "风场" } ,       { "uv" ,       "流场" } ,               { "rain" ,     "降水(对流+非对流)" } ,   { "u10;v10" ,      "10米风场" }
0144          } ;
0145
0146          String  rainName []   = { "rainc" , "rainnc" , "rc" , "rn" , "rain" } ;
0147          String  prefix      = ctl.tdef.increment;
0148          String  dateNames [] = { "yr" , "mo" , "dy" , "hr" , "mn" , "we" } ;
0149          String  dateValues [] = { "年" , "月" , "日" , "小时" , "分钟" , "星期" } ;
0150          String  _name_tmp   = _name.trim () ;
0151          if ( _name_tmp.startsWith ( "_" ) ) { //如果要素名以 "_" 开始，则为对应要素的当前预报时段与上一预报时段的差值
0152              _name_tmp   = _name.trim () .substring ( 1 ) ;
0153              for ( int i= 0 ;i<dateNames.length;i++ ) {
0154                  prefix  = prefix.replaceAll ( dateNames [ i ] , dateValues [ i ]) ;
0155              }
0156          }
0157          else {
0158              _name_tmp   = _name.trim () ;
0159              prefix      = "" ;
0160              for ( int i= 0 ;i<rainName.length;i++ ) {
0161                  if ( _name_tmp.equalsIgnoreCase ( rainName [ i ]) ) {
0162                      prefix      = "累积" ;
0163                  }
0164              }
0165          }
0166
0167          if ( "grapes" .equalsIgnoreCase ( model ) ) { //grapes 模式的 GrADS 描述文件
0168              for ( int i= 0 ;i<eGrapesUnicode.length;i++ ) {
0169                  if ( 0 == eGrapesUnicode [ i ][ 0 ] .compareToIgnoreCase ( _name_tmp ) ) {
0170                      return ( prefix + eGrapesUnicode [ i ][ 1 ]) ;
0171                  }
0172              }
0173          }
0174          else if ( "mm5" .equalsIgnoreCase ( model ) ) { //mm5 模式的 GrADS 描述文件
0175              for ( int i= 0 ;i<eMM5Unicode.length;i++ ) {
0176                  if ( 0 == eMM5Unicode [ i ][ 0 ] .compareToIgnoreCase ( _name_tmp ) ) {
0177                      return ( prefix+ eMM5Unicode [ i ][ 1 ]) ;
0178                  }
0179              }
0180          }
0181          else { //其它模式
0182              for ( int i= 0 ;i<eGrapesUnicode.length;i++ ) {
0183                  if ( 0 == eGrapesUnicode [ i ][ 0 ] .compareToIgnoreCase ( _name_tmp ) ) {
0184                      return ( prefix + eGrapesUnicode [ i ][ 1 ] ) ;
0185                  }
0186              }
0187              for ( int i= 0 ;i<eMM5Unicode.length;i++ ) {
0188                  if ( 0 == eMM5Unicode [ i ][ 0 ] .compareToIgnoreCase ( _name_tmp ) ) {
0189                      return ( prefix + eMM5Unicode [ i ][ 1 ] ) ;
0190                  }
0191              }
0192          }
0193          return ( getElementDescription ( _name_tmp )) ;
0194      }
0195
0196 //获得时效数
0197      public int getValidNumber () {
0198          if ( enabled ) {
0199              return ( ctl.tdef.count ) ;
0200          }
0201          return 0 ;
0202      }
0203 //设置 GrADS 的命令行
0204      public   boolean setCommand ( String cmd ) {
0205          if ( null != cmd &&
0206              ( cmd.trim () .startsWith ( "grads" ) || (( new File ( cmd )) .getName () .startsWith ( "grads" )) ) ) {
0207              command = cmd.trim () ;
0208              return ( true ) ;
0209          }
0210          return ( false ) ;
0211      }
0212 //获得 GrADS 命令的输出
0213      public   String  getMessage () {
0214          return ( message ) ;
0215      }
0216 //根据缺省参数获得生成的图像文件的说明，以作为网页中<img>的title
0217      public String getTitle ( boolean unicode ) {
0218          //System.out.println("getTitle: " + element + " " + level + " " + valid);
0219          String  tmpTitle =
0220              ( unicode? "初始场：" : "Start    : " ) + ctl.tdef.start + "/n" +
0221              ( unicode? "预报时效：" : "Validate : " ) + getValid () + "/n" +
0222              ( unicode? "要素名：" : "Element  : " ) + getElement () + ( unicode? "(" +getElementUnicode ( element ) + ")" :getElementDescription ()) + "/n" +
0223              ( unicode? "层　　次：" : "Level    : " ) + ( checkElement ( element ) > 1 ?getLevel () : "" ) ;
0224          return ( tmpTitle ) ;
0225      }
0226 //根据指定的参数获得生成的图像文件的说明，以作为网页中<img>的title
0227      public String getTitle ( String e, String l, String t ) {
0228          if ( !enabled || !setElement ( e ) || !setLevel ( l ) || !setValid ( t ) ) {
0229              return ( "" ) ;
0230          }
0231          return ( getTitle ( false )) ;
0232      }
0233 //根据缺省的参数获得在 gs 中生成的图像文件名
0234      public   String  getImageFilename () {
0235          String  imageFile   = getImageFilename ( model+ctl.xdef.increment, ".jpg" ) ;
0236          return ( imageFile ) ;
0237      }
0238 //根据指定的参数获得在 gs 中生成的图像文件名
0239      public   String  getImageFilename ( String prefix, String ext ) {
0240          String  imageFile       = "" ;
0241          if ( !enabled ) {
0242              return ( imageFile ) ;
0243          }
0244          String  prefixString    = null ==prefix|| "" .equals ( prefix.trim ()) ?
0245              model+ctl.xdef.increment : //"grapes" :
0246              prefix.trim () ;
0247          String  extString       = null ==ext   || "" .equals ( ext.trim ())    ||!ext.trim () .startsWith ( "." ) ?
0248              ".jpg" :
0249              ext.trim () ;
0250          imageFile               =
0251              prefixString    + "_" + //文件名前缀
0252              ctl.tdef.start  + "_" + //开始时间
0253              getElement () .replaceAll ( ";" , "" ) + "_" + //要素
0254              ( checkElement ( element ) > 1 ?getLevel () + "_" : "" ) +    //层次
0255              getValid ()       + "_" + //预报时效
0256              width + "x" + height + "_" +
0257              ( background.equalsIgnoreCase ( "white" ) ? "W" : "B" ) +
0258              ( overlapping? "O" : "" ) +
0259              ( shading? "S" : "" ) +
0260              extString;               //扩展名
0261          return ( imageFile ) ;
0262      }
0263 //与要素有关的方法
0264      public   String  getElement () { //获得当前指定的要素
0265          return ( element ) ;
0266      }
0267      public   int      checkElement ( String e ) { //检测要素名是否存在，返回值：0不存在、1单层要素、>1多层要素
0268          if ( !enabled || null == e || "" .equals ( e.trim ()) ) {
0269              return ( 0 ) ;
0270          }
0271          if ( "uv" .equals ( e ) ) {
0272              return ( 2 ) ;
0273          }
0274
0275          //原 CTL 中列出的要素
0276          for ( int i= 0 ;i<ctl.vars.size () ;i++ ) {
0277              if ( 0 == (( gradsVar ) ctl.vars.get ( i )) .name.compareToIgnoreCase ( e ) ) {
0278                  return ((( gradsVar ) ctl.vars.get ( i )) .levels== 0 ? 1 : (( gradsVar ) ctl.vars.get ( i )) .levels ) ;
0279              }
0280          }
0281
0282          //扩展要素，在原 CTL 中列出的要素前加 _ 符号，用于显示该要素的当前时段与上一时段的差值
0283          for ( int i= 0 ;i<ctl.vars.size () ;i++ ) {
0284              if ( 0 == ( "_" + (( gradsVar ) ctl.vars.get ( i )) .name ) .compareToIgnoreCase ( e ) ) {
0285                  return ((( gradsVar ) ctl.vars.get ( i )) .levels== 0 ? 1 : (( gradsVar ) ctl.vars.get ( i )) .levels ) ;
0286              }
0287          }
0288
0289          //扩展要素，用于显示 rainc+rainnc 之和
0290          if ( "rain" .equalsIgnoreCase ( e ) && "grapes" .equalsIgnoreCase ( model ) &&
0291              checkElement ( "rainc" ) > 0   && checkElement ( "rainnc" ) > 0 ) { //GRAPES
0292              return ( checkElement ( "rainc" )) ;
0293          }
0294
0295          //扩展要素，用于显示 rainc+rainnc 之和与上一时段的相应值之差
0296          if ( "_rain" .equalsIgnoreCase ( e ) && "grapes" .equalsIgnoreCase ( model ) &&
0297              checkElement ( "rainc" ) > 0    && checkElement ( "rainnc" ) > 0 ) { //GRAPES
0298              return ( checkElement ( "rainc" )) ;
0299          }
0300
0301          //扩展要素，用于显示 rc+rn 之和
0302          if ( "rain" .equalsIgnoreCase ( e ) && "mm5" .equalsIgnoreCase ( model ) &&
0303              checkElement ( "rc" ) > 0      && checkElement ( "rn" ) > 0 ) { //MM5
0304              return ( checkElement ( "rc" )) ;
0305          }
0306
0307          //扩展要素，用于显示 rc+rn 之和与上一时段的相应值之差
0308          if ( "_rain" .equalsIgnoreCase ( e ) && "mm5" .equalsIgnoreCase ( model ) &&
0309              checkElement ( "rc" ) > 0       && checkElement ( "rn" ) > 0 ) { //MM5
0310              return ( checkElement ( "rc" )) ;
0311          }
0312
0313          //组合要素
0314          String  combineVar []     = { ";" , "+" , "-" , "*" , "/" } ;
0315          for ( int i= 0 ;i<combineVar.length;i++ ) {
0316              if ( 0 < e.indexOf ( combineVar [ i ]) )  {
0317                  String  e1 = "" , e2 = "" ;
0318                  boolean b1 = false, b2 = false ;
0319                  e1  = e.substring ( 0 ,e.indexOf ( combineVar [ i ])) ;
0320                  e2  = e.substring ( e.indexOf ( combineVar [ i ]) + 1 ) ;
0321                  for ( int j= 0 ;j<ctl.vars.size () ;j++ ) {
0322                      b1  = (( gradsVar ) ctl.vars.get ( j )) .name.equals ( e1 ) ? true :b1;
0323                      b2  = (( gradsVar ) ctl.vars.get ( j )) .name.equals ( e2 ) ? true :b1;
0324                  }
0325                  if ( b1 && b2 ) {
0326                      element = e;
0327                      return ( checkElement ( e1 )) ;
0328                  }
0329              }
0330          }
0331          if ( e.startsWith ( "_" ) ) {
0332              return ( checkElement ( e.substring ( 1 ))) ;
0333          }
0334
0335          return ( 0 ) ;
0336      }
0337      public   String  getElementDescription () { //获得要素的说明
0338          for ( int i= 0 ;i<ctl.vars.size () ;i++ ) {
0339              if ( 0 == (( gradsVar ) ctl.vars.get ( i )) .name.compareToIgnoreCase ( element ) ) {
0340                  return (
0341                      "" .equals ((( gradsVar ) ctl.vars.get ( i )) .contents ) ?
0342                      "" :
0343                      "(" + (( gradsVar ) ctl.vars.get ( i )) .contents + ")" ) ;
0344              }
0345          }
0346          return ( "" ) ;
0347      }
0348      public   String  getElementDescription ( String e ) { //获得要素的说明
0349          for ( int i= 0 ;i<ctl.vars.size () ;i++ ) {
0350              if ( 0 == (( gradsVar ) ctl.vars.get ( i )) .name.compareToIgnoreCase ( e ) ) {
0351                  return (
0352                      "" .equals ((( gradsVar ) ctl.vars.get ( i )) .contents ) ?
0353                      "" :
0354                      "(" + (( gradsVar ) ctl.vars.get ( i )) .contents + ")" ) ;
0355              }
0356          }
0357          return ( "" ) ;
0358      }
0359      public   boolean setElement ( String e ) { //设置要素名
0360          if ( checkElement ( e ) >= 0 ) {
0361              element = e;
0362              return ( true ) ;
0363          }
0364          System.out.println ( "ERROR: setElement( " + e + " )" ) ;
0365          return ( false ) ;
0366      }
0367      public   boolean setElement ( int e ) { //索引
0368          if ( !enabled ) {
0369              return ( false ) ;
0370          }
0371          else if ( enabled && e >= 0 && e < ctl.vars.size () ) {
0372              element = (( gradsVar ) ctl.vars.get ( e )) .name;
0373              return ( true ) ;
0374          }
0375          else {
0376              return ( false ) ;
0377          }
0378      }
0379      public   int createElementScript ( Vector vector ) { //根据当前给定的要素生成 gs 脚本
0380          int len = vector.size () ;
0381          if ( "u;v" .equals ( element ) ) { //风矢
0382              vector.add ( "'set gxout barb'" ) ;
0383              vector.add ( "'d skip(u,8);skip(v,8)'" ) ;
0384          }
0385          else if ( "uv" .equals ( element ) ) { //流场
0386              vector.add ( "'set gxout stream'" ) ;
0387              vector.add ( "'d skip(u,8);skip(v,8)'" ) ;
0388          }
0389          else if ( "h" .equals ( element ) ) {
0390              vector.add ( "'" + "set cint 4'" ) ;
0391              vector.add ( "'" + "d " + element + "/10'" ) ;
0392          }
0393          else if ( element.startsWith ( "t" ) || element.startsWith ( "td" ) ) {
0394              vector.add ( "'" + "set cint 4'" ) ;
0395              vector.add ( "'" + "d " + element + "-273.15'" ) ;
0396          }
0397          else if ( element.startsWith ( "p" ) ) {
0398              vector.add ( "'" + "set cint 2.5'" ) ;
0399              vector.add ( "'" + "d " + element + "'" ) ;
0400          }
0401          else if ( //累积降水量
0402              "rainc" .equals ( element ) || "rainnc" .equals ( element ) ||
0403              "rc" .equals ( element ) || "rn" .equals ( element )    ) {
0404              vector.add ( "'set cmin 0.01'" ) ;
0405              vector.add ( "'set clevs 0 10 25 50 100 250'" ) ;
0406              vector.add ( "'" + "d " + element + ( "mm5" .equals ( model ) ? "*10" : "" ) + "'" ) ;
0407          }
0408          else if ( "rain" .equalsIgnoreCase ( element ) ) { //总累积降水量rain=rainc+rainnc或rain=rc+rn
0409              String  e1 = "" , e2 = "" ;
0410              if ( "grapes" .equalsIgnoreCase ( model ) && checkElement ( "rainc" ) > 0 && checkElement ( "rainnc" ) > 0 ) {
0411                  e1 = "rainc" ;
0412                  e2 = "rainnc" ;
0413              }
0414              else if ( "mm5" .equalsIgnoreCase ( model ) && checkElement ( "rc" ) > 0 && checkElement ( "rn" ) > 0 ) {
0415                  e1 = "rc*10" ;
0416                  e2 = "rn*10" ;
0417              }
0418              else {
0419                  return ( vector.size () - len ) ;
0420              }
0421              vector.add ( "'set clevs 0 10 25 50 100 250'" ) ;
0422              vector.add ( "'" + "d " + e1 + "+" + e2 + "'" ) ;
0423          }
0424          else if ( element.startsWith ( "_rain" ) ) { //总降水量(当前预报时段与上一预报时段的差值)
0425              int tIndex = Integer.parseInt ( valid ) ;
0426              int tCount = "_rain" .equalsIgnoreCase ( element ) ? 1 : Integer.parseInt ( element.substring ( "_rain" .length ())) ;
0427              vector.add ( "'set clevs 0 10 25 50 100 250'" ) ;
0428 //          vector.add("'set ccols 2 3 4 5 6 7 8'");
0429 //          vector.add("'set cmin 0.01'");
0430              String  e1 = "" , e2 = "" ;
0431              if ( "grapes" .equalsIgnoreCase ( model )     &&
0432                      checkElement ( "rainc" )    > 0      &&
0433                      checkElement ( "rainnc" )   > 0      ) {
0434                  e1 = "rainc" ;
0435                  e2 = "rainnc" ;
0436              }
0437              else if ( "mm5" .equalsIgnoreCase ( model )   &&
0438                      checkElement ( "rc" )   > 0          &&
0439                      checkElement ( "rn" )   > 0          ) {
0440                  e1 = "10*rc" ;
0441                  e2 = "10*rn" ;
0442              }
0443              else {
0444                  return ( vector.size () - len ) ;
0445              }
0446              if ( tIndex > tCount ) {
0447                  vector.add ( "'" + "d " + element + "'" ) ;
0448                  vector.add (
0449                      "'" + "d " + e1 + "(t=" + ( tIndex ) + ")-" + e1 + "(t=" + ( tIndex-tCount ) + ")+" + e2 + "(t=" + ( tIndex ) + ")-" + e2 + "(t=" + ( tIndex-tCount ) + ")'" ) ;
0450              }
0451              else {
0452                  vector.add ( "'" + "d " + e1 + "+" + e2 + "'" ) ;
0453              }
0454          }
0455          else if ( element.startsWith ( "_006" ) || element.startsWith ( "_012" ) || element.startsWith ( "_024" ) || element.startsWith ( "_048" ) ) { //当前预报时段与指定的前某一预报时段的差值，一般用于统计时段内累积降水量
0456              int tIndex = Integer.parseInt ( valid ) ;
0457              int tCount = Integer.parseInt ( element.substring ( 1 , 4 )) ;
0458              vector.add ( "'set clevs 0 10 25 50 100 250'" ) ;
0459              if ( tIndex > tCount ) {
0460                  vector.add ( "'" + "d " + element.substring ( 4 ) + "(t=" + ( tIndex ) + ")-" + element.substring ( 4 ) + "(t=" + ( tIndex-tCount ) + ")'" ) ;
0461              }
0462              else {
0463                  vector.add ( "'" + "d " + element.substring ( 4 ) + "'" ) ;
0464              }
0465          }
0466          else if ( element.startsWith ( "_" ) ) { //当前预报时段与上一预报时段的差值，一般用于降水量
0467              int tIndex = Integer.parseInt ( valid ) ;
0468              vector.add ( "'set clevs 0 10 25 50 100 250'" ) ;
0469              if ( tIndex > 1 ) {
0470                  vector.add ( "'" + "d " + element.substring ( 1 ) + "(t=" + ( tIndex ) + ")-" + element.substring ( 1 ) + "(t=" + ( tIndex- 1 ) + ")'" ) ;
0471              }
0472              else {
0473                  vector.add ( "'" + "d " + element.substring ( 1 ) + "'" ) ;
0474              }
0475          }
0476          else { //其它
0477              vector.add ( "'" + "d " + element + "'" ) ;
0478          }
0479          return ( vector.size () - len ) ;
0480      }
0481 //与层次有关的方法
0482      public   String  getLevel () { //获得当前指定的层次
0483          return ( level ) ;
0484      }
0485      public   boolean checkLevel ( String l ) { //检测层次
0486          if ( !enabled || null == l || "" .equals ( l ) ) {
0487              return ( false ) ;
0488          }
0489          for ( int i= 0 ;i<ctl.zdef.count;i++ ) {
0490              if ( ctl.zdef.value [ i ] .equals ( l ) ) {
0491                  return ( true ) ;
0492              }
0493          }
0494          return ( false ) ;
0495      }
0496      public   boolean setLevel ( String l ) { //设置层次(500、700、850、…………)
0497          if ( !enabled ) {
0498              return ( false ) ;
0499          }
0500          if ( null == l || "" .equals ( l.trim ()) ) {
0501              level   = ctl.zdef.value [ 0 ] ;
0502              return ( false ) ;
0503          }
0504          for ( int i= 0 ;i<ctl.zdef.count;i++ ) {
0505              if ( ctl.zdef.value [ i ] .equals ( l ) ) {
0506                  level   = checkElement ( element ) > 1 ?l:ctl.zdef.value [ 0 ] ;
0507                  return ( true ) ;
0508              }
0509          }
0510          return ( false ) ;
0511      }
0512      public   boolean setLevel ( int l ) { //设置层次(第几层)
0513          if ( !enabled ) {
0514              return ( false ) ;
0515          }
0516          else if ( 0 <= l && l < ctl.zdef.count ) {
0517              level   = checkElement ( element ) > 1 ?ctl.zdef.value [ l ] :ctl.zdef.value [ 0 ] ;
0518              return ( true ) ;
0519          }
0520          else {
0521              level   = ctl.zdef.value [ 0 ] ;
0522              return ( false ) ;
0523          }
0524      }
0525 //与时效有关的方法
0526      public   int getValidIndex () { //获得当前指定的时效
0527          try {
0528              int index = Integer.parseInt ( valid ) - 1 ;
0529              if ( 0 <= index && index < ctl.tdef.count ) {
0530                  return ( index ) ;
0531              }
0532          }
0533          catch ( Exception e ) {
0534          }
0535          return ( 0 ) ;
0536      }
0537      public   String  getValid () { //获得当前指定的时效
0538          try {
0539              int index = Integer.parseInt ( valid ) - 1 ;
0540              if ( 0 <= index && index < ctl.tdef.count ) {
0541                  return ( ctl.tdef.value [ index ]) ;
0542              }
0543          }
0544          catch ( Exception e ) {
0545          }
0546          return ( valid ) ;
0547      }
0548      public   boolean checkValid ( String t ) { //检测时效
0549          if ( !enabled || null == t || "" .equals ( t ) ) {
0550              return ( false ) ;
0551          }
0552          for ( int i= 0 ;i<ctl.tdef.count;i++ ) {
0553              if ( ctl.zdef.value [ i ] .equals ( t ) ) {
0554                  return ( true ) ;
0555              }
0556          }
0557          for ( int i= 0 ;i<ctl.tdef.count;i++ ) {
0558              if ( String.valueOf ( i+ 1 ) .equals ( t ) ) {
0559                  return ( true ) ;
0560              }
0561          }
0562          return ( false ) ;
0563      }
0564      public   boolean setValid ( String t ) { //设置时效
0565          if ( !enabled ) {
0566              return ( false ) ;
0567          }
0568          else if ( null == t || "" .equals ( t.trim ()) ) {
0569              valid   = "1" ;
0570              return ( false ) ;
0571          }
0572          for ( int i= 0 ;i<ctl.tdef.count;i++ ) { //传入的是 000hr、006hr、012hr ...
0573              if ( ctl.tdef.value [ i ] .equals ( t.trim ()) ) {
0574                  valid   = String.valueOf ( i+ 1 ) ;
0575                  return ( true ) ;
0576              }
0577          }
0578          for ( int i= 0 ;i<ctl.tdef.count;i++ ) { //传入的是 1、2、3 ...
0579              if ( String.valueOf ( i+ 1 ) .equals ( t.trim ()) ) {
0580                  valid   = String.valueOf ( i+ 1 ) ;
0581                  return ( true ) ;
0582              }
0583          }
0584          return ( false ) ;
0585      }
0586      public   boolean setValid ( int t ) { //设置时效(第几时效)
0587          if ( !enabled ) {
0588              return ( false ) ;
0589          }
0590          else if ( 0 <= t && t < ctl.tdef.count ) {
0591              valid   = String.valueOf ( t + 1 ) ;
0592              return ( true ) ;
0593          }
0594          return ( false ) ;
0595      }
0596 //构造函数
0597      public GrADS () {
0598          ctl     = new gradsCtl () ;
0599      }
0600 //将给定的 line 字符串按空格、回车、换行、TAB键分隔开
0601      public String [] separateLine ( String line ) {
0602          if ( null == line || "" .equals ( line.trim ()) ) {
0603              return ( null ) ;
0604          }
0605          String  s = line.trim () ;
0606          Vector  v = new Vector () ;
0607 //注：可以考虑使用 StringTokenizer，见 Diamond04.java 的 loadFromFile() 方法
0608          s   = s.replaceAll ( "/r" , " " ) ; //回车
0609          s   = s.replaceAll ( "/n" , " " ) ; //换行
0610          s   = s.replaceAll ( "/t" , " " ) ; //TAB键
0611          while ( - 1 != s.indexOf ( "  " ) ) { //双空格换为单空格
0612              s   = s.replaceAll ( "  " , " " ) .trim () ;
0613          }
0614          int      index = - 1 ;
0615          while ( - 1 != ( index = s.indexOf ( " " )) ) {
0616              v.add ( s.substring ( 0 , index )) ;
0617              s   = s.substring ( index+ 1 ) .trim () ;
0618          }
0619          if ( ! "" .equals ( s ) ) {
0620              v.add ( s ) ;
0621          }
0622          String  str []    = new String [ v.size ()] ;
0623          for ( int i= 0 ;i<v.size () ;i++ ) {
0624              str [ i ]   = ( String ) v.get ( i ) ;
0625          }
0626          return ( str ) ;
0627      }
0628 //判断该行是否以 GrADS 的关键字开头
0629      public boolean isSymbol ( String line ) {
0630          String symbol [] = { "dset" , "options" , "title" , "undef" , "pdef" , "xdef" , "ydef" , "zdef" , "tdef" , "vars" , "endvars" } ;
0631          for ( int i= 0 ;i<symbol.length;i++ ) {
0632              if ( line.toLowerCase () .startsWith ( symbol [ i ]) ) {
0633                  return ( true ) ;
0634              }
0635          }
0636          return ( false ) ;
0637      }
0638 //判断 file 是否是一个有效的 GrADS 的 CTL 文件
0639      public static boolean isValid ( String file ) {
0640
0641          File    f = new File ( file ) ;
0642          if ( f.exists ()       && f.isFile ()            && //确保文件存在
0643              f.length () > 32 && f.length () <= size ) { //文件大小符合要求
0644              try {
0645                  FileReader      fr  = new FileReader ( file ) ;
0646                  BufferedReader  br  = new BufferedReader ( fr ) ;
0647                  String  stringLine  = null ;
0648                  Vector  v           = new Vector () ;
0649                  boolean match       = false ;
0650
0651                  while ( null != ( stringLine = br.readLine ()) ) {
0652                      if ( ! "" .equals ( stringLine.trim ())    &&
0653                          stringLine.toLowerCase () .startsWith ( "dset" ) ) { //数据文件所在行
0654                          br.close () ;
0655                          fr.close () ;
0656                          String  fData   = stringLine.substring ( 4 ) .trim () ;
0657                          if ( fData.startsWith ( "^" ) ) {
0658                              fData   = f.getParent () + "/" + fData.substring ( 1 ) ;
0659                          }
0660                          //System.out.println(fData);
0661                          File    d   = new File ( fData ) ;
0662                          return ( d.exists () && d.isFile ()) ;
0663                      }
0664                  }
0665                  br.close () ;
0666                  fr.close () ;
0667              }
0668              catch ( Exception e ) {
0669              }
0670          }
0671          return ( false ) ;
0672      }
0673 //在指定的目录中列出以 prefix 开头且以 suffix 结尾的有效的 GrADS 描述文件名(不含目录)
0674      public static String [] listFiles ( String path, String prefix, String suffix ) {
0675          File    f   = new File ( path ) ;
0676          Vector  v   = new Vector () ;
0677          if ( f.isDirectory () && f.exists () ) {
0678              //System.out.println(path + " , " + prefix + " , " + suffix);
0679              File    files [] = f.listFiles () ;
0680              String  fname = "" ;
0681              for ( int i= 0 ;i<files.length;i++ ) {
0682                  //System.out.println(files[i]);
0683                  fname   = files [ i ] .getName () ;
0684                  if (
0685                      ( null == prefix || "" .equals ( prefix ) || fname.startsWith ( prefix ))    &&
0686                      ( null == suffix || "" .equals ( suffix ) || fname.endsWith ( suffix )  )    &&
0687                      isValid ( files [ i ] .getAbsolutePath ()) ) {
0688                      //System.out.println("true");
0689                      v.add ( fname ) ;
0690                  }
0691              }
0692          }
0693          if ( v.size () == 0 ) {
0694              return ( null ) ;
0695          }
0696          else {
0697              String  fnames [] = new String [ v.size ()] ;
0698              for ( int i= 0 ;i<v.size () ;i++ ) {
0699                  fnames [ i ]    = ( String ) v.get ( i ) ;
0700              }
0701              java.util.Arrays.sort ( fnames ) ;
0702              return ( fnames ) ;
0703          }
0704      }
0705 //读取 GrADS 的 CTL 说明文件
0706      public boolean loadFromFile ( String file ) throws IOException {
0707          enabled = false ;
0708          ctl.vars.clear () ;
0709          pathname    = "." ;
0710          filename    = "" ;
0711          File    f = new File ( file ) ;
0712          if ( !f.exists () ) {
0713              return ( false ) ;
0714          }
0715          pathname    = f.getPath () ;
0716          if ( null == pathname || "" .equals ( pathname ) ) {
0717              pathname    = "." ;
0718          }
0719          filename    = f.getName () ;
0720          try {
0721          FileReader      fr          = new FileReader ( file ) ;
0722          BufferedReader  br          = new BufferedReader ( fr ) ;
0723          String          stringLine  = null ;
0724          Vector          lines       = new Vector () ;
0725          boolean          vars        = false ;
0726          while ( null != ( stringLine = br.readLine ()) ) {
0727              stringLine  = stringLine.trim () ;
0728              if ( ! "" .equals ( stringLine ) ) {
0729                  if ( stringLine.startsWith ( "vars" ) ) {
0730                      vars    = true ;
0731                  }
0732                  if ( isSymbol ( stringLine ) || vars ) {
0733                      lines.add ( stringLine ) ;
0734                  }
0735                  else if ( lines.size () == 1 ) { //某些 CTL 文件的 options 行(第2行)并不带有 options 标识
0736                      lines.add ( "options " + stringLine ) ;
0737                  }
0738                  else {
0739                      if ( 0 == lines.size () ) {
0740                          lines.add ( stringLine ) ;
0741                      }
0742                      else {
0743                          lines.set (
0744                              lines.size () - 1 ,
0745                              ( String ) lines.get ( lines.size () - 1 ) + " " + stringLine ) ;
0746                      }
0747                  }
0748              }
0749          }
0750          br.close () ;
0751          fr.close () ;
0752          enabled = extractStrings ( lines ) ;
0753          return ( enabled ) ;
0754          }
0755          catch ( Exception e ) {
0756              return ( false ) ;
0757          }
0758      }
0759 //解析 GrADS 的 CTL 文件内容
0760      private boolean extractStrings ( Vector strings ) {
0761          String  line = "" ;
0762          String  str [] = null ;
0763          try {
0764              for ( int i= 0 ;i<strings.size () ;i++ ) {
0765                  line    = ( String ) strings.get ( i ) ;
0766                  str = separateLine ( line ) ;
0767                  if ( line.startsWith ( "dset" ) ) {
0768                      ctl.dset    = str [ 1 ] ;
0769                  }
0770                  else if ( line.startsWith ( "options" ) ) {
0771                      ctl.options = line.substring ( "options" .length ()) .trim () ;
0772                  }
0773                  else if ( line.startsWith ( "title" ) ) {
0774                      ctl.title   = line.substring ( "title" .length ()) .trim () ;
0775                  }
0776                  else if ( line.startsWith ( "undef" ) ) {
0777                      ctl.undef   = line.substring ( "undef" .length ()) .trim () ;
0778                  }
0779                  else {
0780                      if ( line.startsWith ( "xdef" ) ) {
0781                          ctl.xdef.extractValue ( str ) ;
0782                      }
0783                      else if ( line.startsWith ( "ydef" ) ) {
0784                          ctl.ydef.extractValue ( str ) ;
0785                      }
0786                      else if ( line.startsWith ( "zdef" ) ) {
0787                          ctl.zdef.extractValue ( str ) ;
0788                      }
0789                      else if ( line.startsWith ( "tdef" ) ) {
0790                          ctl.tdef.extractValue ( str ) ;
0791                      }
0792                      else if ( line.startsWith ( "vars" ) ) {
0793                          int varCount    = Integer.parseInt ( line.substring ( "vars" .length ()) .trim ()) ;
0794                          for ( int j= 0 ;j<varCount;j++ ) {
0795                              gradsVar    var = new gradsVar () ;
0796                              line    = ( String ) strings.get ( i+j+ 1 ) ;
0797                              str = separateLine ( line ) ;
0798                              var.name    = str [ 0 ] ;
0799                              var.levels  = Integer.parseInt ( str [ 1 ]) ;
0800                              var.reserve = str [ 2 ] ;
0801                              for ( int k= 3 ;k<str.length;k++ ) {
0802                                  var.contents    = var.contents + str [ k ] + " " ;
0803                              }
0804                              var.contents    = var.contents.trim () ;
0805                              var.unicode     = getElementUnicode ( var.name ) ;
0806                              if ( "" .equals ( var.unicode ) ) {
0807                                  var.unicode = var.contents;
0808                              }
0809                              ctl.vars.add ( var ) ;
0810                          }
0811                      }
0812                  }
0813              }
0814          }
0815          catch ( Exception e ) {
0816              return ( false ) ;
0817          }
0818          enabled = ( 0 < ctl.vars.size () && 0 < ctl.tdef.count && 0 < ctl.zdef.count ) ;
0819          setElement ( "" ) ;
0820          setLevel   ( "" ) ;
0821          setValid   ( "" ) ;
0822          return ( enabled ) ;
0823      }
0824 //生成 GrADS 的 CTL 说明文件
0825      public boolean writeToFile ( String file ) {
0826          if ( !enabled ) {
0827              return ( false ) ;
0828          }
0829          try {
0830              FileWriter      fw  = new FileWriter ( file ) ;
0831              BufferedWriter  bw  = new BufferedWriter ( fw ) ;
0832              bw.write ( "dset " + ctl.dset + "/n" ) ;
0833              bw.write ( "options " + ctl.options + "/n" ) ;
0834              bw.write ( "title " + ctl.title + "/n" ) ;
0835              bw.write ( "undef " + ctl.undef + "/n" ) ;
0836              bw.write ( ctl.xdef.convertToString () + "/n" ) ;
0837              bw.write ( ctl.ydef.convertToString () + "/n" ) ;
0838              bw.write ( ctl.zdef.convertToString () + "/n" ) ;
0839              bw.write ( ctl.tdef.convertToString () + "/n" ) ;
0840              bw.write ( "vars " + String.valueOf ( ctl.vars.size ()) + "/n" ) ;
0841              for ( int i= 0 ;i<ctl.vars.size () ;i++ ) {
0842                  bw.write ((( gradsVar ) ctl.vars.get ( i )) .convertToString ()) ;
0843                  bw.newLine () ;
0844              }
0845              bw.write ( "endvars" + "/n" ) ;
0846              bw.close () ;
0847              fw.close () ;
0848              return ( true ) ;
0849          }
0850          catch ( Exception e ) {
0851              System.out.println ( e.getMessage ()) ;
0852              e.printStackTrace () ;
0853              return ( false ) ;
0854          }
0855      }
0856 /** 根据给定的参数生成 GrADS 的 gs 脚本
0857 参数:
0858      path    - 存放 gs 脚本文件及图像文件的路径
0859      e, l, t - 变量
0860 返回:
0861      脚本串
0862 */
0863      public String getScript ( String path, String e, String l, String t ) {
0864          if ( !setElement ( e ) || !setLevel ( l ) || !setValid ( t ) ) {
0865              return ( "" ) ;
0866          }
0867 /*
0868          System.out.println("element = " + element);
0869          System.out.println("level   = " + level);
0870          System.out.println("valid   = " + valid);
0871 */
0872          String  gsString    = "" ;
0873          String  tmpTitle        = getElement () + "_" + ( checkElement ( element ) > 1 ?getLevel () + "_" : "" ) + ctl.tdef.start + "+" + getValid () ; //getElementDescription();
0874
0875          Vector  gsVector        = new Vector () ;
0876          int      originLength    = 0 ;
0877          File    f = new File ( path + ( path.endsWith ( "/" ) ? "" : "/" ) + scriptfile ) ;
0878          String  strLine = null ;
0879          if ( overlapping && //迭加
0880              f.exists () && f.isFile () ) { //文件存在
0881              try {
0882                  FileReader      fr  = new FileReader ( f ) ;
0883                  BufferedReader  br  = new BufferedReader ( fr ) ;
0884                  while ( null != ( strLine = br.readLine ()) ) {
0885                      if ( ! "" .equals ( strLine.trim ()) ) {
0886                          gsVector.add ( strLine.trim ()) ;
0887                      }
0888                  }
0889                  br.close () ;
0890                  fr.close () ;
0891                  for ( int i=gsVector.size () - 1 ;i> 0 ;i-- ) {
0892                      strLine = ( String ) gsVector.get ( i ) ;
0893                      if ( strLine.startsWith ( "'" + "draw title" ) ) {
0894                          tmpTitle   = strLine.substring (( "'" + "draw title" ) .length () ,strLine.length () - 1 ) .trim () + "/" + tmpTitle;
0895                          gsVector.set ( i, "" ) ;
0896                      }
0897                      if ( strLine.startsWith ( "'" + "draw xlab" ) ) {
0898                          gsVector.set ( i, "" ) ;
0899                      }
0900                      if ( strLine.startsWith ( "'" + "set string" ) ) {
0901                          gsVector.set ( i, "" ) ;
0902                      }
0903                      if ( strLine.startsWith ( "'" + "draw string" ) ) {
0904                          gsVector.set ( i, "" ) ;
0905                      }
0906                      if ( strLine.startsWith ( "'" + "set lon" ) ) {
0907                          gsVector.set ( i, "" ) ;
0908                      }
0909                      if ( strLine.startsWith ( "'" + "set lat" ) ) {
0910                          gsVector.set ( i, "" ) ;
0911                      }
0912                  }
0913                  for ( int i=gsVector.size () - 1 ;i> 0 ;i-- ) {
0914                      strLine = ( String ) gsVector.get ( i ) ;
0915                      if ( ( "'" + "quit'" ) .equalsIgnoreCase ( strLine ) ) {
0916                          gsVector.set ( i, "" ) ;
0917                      }
0918                  }
0919                  for ( int i= 0 ;i<gsVector.size () ;i++ ) {
0920                      System.out.println (( String ) gsVector.get ( i )) ;
0921                  }
0922 //              gsVector.add("'set lon " + String.valueOf(Math.min(left, right)) + " " + String.valueOf(Math.max(left, right)) + "'");
0923 //              gsVector.add("'set lat " + String.valueOf(Math.min(top, bottom)) + " " + String.valueOf(Math.max(top, bottom)) + "'");
0924              }
0925              catch ( Exception ex ) {
0926              }
0927          }
0928          String  img     = path + ( path.endsWith ( "/" ) ? "" : "/" ) + getImageFilename () ;
0929          originLength    = gsVector.size () ;
0930          if ( - 1 == gsVector.indexOf ( "'" + "reinit" + "'" ) ) {
0931              gsVector.add ( "'" + "reinit" + "'" ) ;
0932              gsVector.add ( "'set grads off'" ) ;
0933          }
0934          gsVector.add ( "'set mpdset cnworld'" ) ; //画中国的省界，需要将 cnworld 文件复制到 /usr/local/lib/grads 路径
0935          gsVector.add ( "'set rgb 17 160 160 160'" ) ;
0936          gsVector.add ( "'set map 17 0 1'" ) ; //棕色缺省线型线宽为1
0937          if ( null != pathname && "" .equals ( pathname ) ) {
0938              gsVector.add ( "'" + "! cd " + pathname + "'" ) ;
0939          }
0940          gsVector.add ( "'" + "open " + filename + "'" ) ;
0941
0942          if ( !overlapping && shading ) { //非迭加时显示填充图
0943              gsVector.add ( "'" + "set gxout shaded" + "'" ) ;
0944              gsVector.add ( "'" + "set lev " + level + "'" ) ;
0945              gsVector.add ( "'" + "set t " + valid + "'" ) ;
0946              gsVector.add ( "'" + "set x " + ( border> 0 ?border: 1 ) + " " + ( border> 0 ?ctl.xdef.count - border:ctl.xdef.count ) + "'" ) ;
0947              gsVector.add ( "'" + "set y " + ( border> 0 ?border: 1 ) + " " + ( border> 0 ?ctl.ydef.count - border:ctl.ydef.count ) + "'" ) ;
0948
0949 //          gsVector.add("'set lon " + String.valueOf(Math.min(left, right)) + " " + String.valueOf(Math.max(left, right)) + "'");
0950 //          gsVector.add("'set lat " + String.valueOf(Math.min(top, bottom)) + " " + String.valueOf(Math.max(top, bottom)) + "'");
0951
0952              createElementScript ( gsVector ) ;
0953
0954          }
0955          gsVector.add ( "'" + "set gxout contour" + "'" ) ;
0956          gsVector.add ( "'" + "set lev " + level + "'" ) ;
0957          gsVector.add ( "'" + "set t " + valid + "'" ) ;
0958              gsVector.add ( "'" + "set x " + ( border> 0 ?border: 1 ) + " " + ( border> 0 ?ctl.xdef.count - border:ctl.xdef.count ) + "'" ) ;
0959              gsVector.add ( "'" + "set y " + ( border> 0 ?border: 1 ) + " " + ( border> 0 ?ctl.ydef.count - border:ctl.ydef.count ) + "'" ) ;
0960
0961 //      gsVector.add("'set lon " + String.valueOf(Math.min(left, right)) + " " + String.valueOf(Math.max(left, right)) + "'");
0962 //      gsVector.add("'set lat " + String.valueOf(Math.min(top, bottom)) + " " + String.valueOf(Math.max(top, bottom)) + "'");
0963
0964          createElementScript ( gsVector ) ;
0965
0966          gsVector.add ( "'" + "set font 0'" ) ;
0967 //      gsVector.add("'" + "set string 19'");
0968          gsVector.add ( "'" + "draw title " + ( "" .equals ( title ) ?tmpTitle:title ) + "'" ) ;
0969 //      String  author = "广西气象减灾研究所";
0970          String  author = "Guangxi Institude of Meteorology and Disaster-reducing Research" ;
0971 //      gsVector.add("'" + "set xlopts 17 1 1'");
0972 //      gsVector.add("'" + "draw xlab " + author + "'");
0973
0974          gsVector.add ( "'" + "set font 2'" ) ;
0975          gsVector.add ( "'" + "set string 17 bl'" ) ;
0976          gsVector.add ( "'" + "draw string 0 0.1 " + author + "'" ) ;
0977          gsVector.add ( "'draw map'" ) ; //画省界
0978          gsVector.add ( "'" + "printim " + img + " x" + width + " y" + height + " " + background + "'" ) ;
0979          gsVector.add ( "'" + "close 1" + "'" ) ;
0980          gsVector.add ( "'" + "quit" + "'" ) ;
0981
0982          for ( int i= 0 ;i<gsVector.size () ;i++ ) {
0983              strLine = (( String ) gsVector.get ( i )) .trim () ;
0984              if ( 0 < strLine.length () ) {
0985                  gsString    =
0986                      gsString +
0987                      (( String ) gsVector.get ( i )) .trim () +
0988                      "/n" ;
0989              }
0990          }
0991          gsString    = gsString + "/n/n" ;
0992
0993          return ( gsString ) ;
0994      }
0995 /**
0996 功能：合并两个 gs 文件
0997 */
0998      public   String  combineScript ( String gs1, String gs2 ) {
0999          return ( "" ) ;
1000      }
1001 /**
1002 功能：调用 GrADS 脚本
1003 参数：
1004      env[]       - 运行 GrADS 所需要的环境，如：PATH=$PATH:/usr/local/bin，一般地，需要将 grads 所在的目录加入到 PATH 中
1005      resultPath  - resultPath 所调用的 CTL 文件所在的路径
1006      imagePath   - 所生成的图像的存放路径
1007      e, l, t     - 变量
1008
1009 返回值：grads 的 gs 脚本
1010 */
1011      public String procScript ( String env [] , String resultPath, String imagePath, String e, String l, String t ) {
1012          //System.out.println("GrADS #1025: procScript(env, resultPath, imagePath, " + e + ", " + l + ", " + t + ")");
1013          message = "" ;
1014          File    rPath   = new File ( resultPath ) ;
1015          File    iPath   = new File ( imagePath ) ;
1016          if ( !enabled        ||
1017              !rPath.exists () || !rPath.isDirectory () ||
1018              !iPath.exists () || !iPath.isDirectory () ||
1019              !setElement ( e )   || !setLevel ( l )          || !setValid ( t ) ) {
1020              System.out.println ( "ERROR: Grapes -> procScript : parameter(s)" ) ;
1021              return ( "ERROR: Grapes -> procScript : parameter(s)" ) ;
1022          }
1023          try {
1024              String  imageFile   =
1025                          imagePath + ( imagePath.endsWith ( "/" ) ? "" : "/" ) + getImageFilename () ;
1026              File    f = new File ( imageFile ) ;
1027              if ( f.exists () ) {
1028                  f.delete () ;
1029              }
1030              String          gs          = imagePath  + ( imagePath.endsWith ( "/" ) ? "" : "/" ) + scriptfile;
1031              String          gsString;
1032              //生成 GrADS 的 gs 脚本
1033              gsString    = getScript ( imagePath,element,level,valid ) ;
1034              FileWriter      fw          = new FileWriter ( gs ) ;
1035              BufferedWriter  bw          = new BufferedWriter ( fw ) ;
1036              bw.write ( gsString ) ;
1037              bw.close () ;
1038              fw.close () ;
1039
1040              message = procScript ( env, resultPath, imagePath, gs ) ;
1041          }
1042          catch ( IOException ex ) {
1043              System.out.println ( ex.getMessage ()) ;
1044              ex.printStackTrace () ;
1045          }
1046          return ( message ) ;
1047      } //调用 GrADS 脚本
1048 /**
1049 * 功能 : 调用 GrADS 脚本
1050 * 参数 :
1051 *      env[]   - 运行环境
1052 *      path    - 开始路径，一般地，设为 GrADS 数据文件所在的路径
1053 *      gs      - GrADS 的 gs 脚本文件名
1054 * 返回值      : GrADS 运行结果
1055 */
1056      public String procScript ( String env [] , String pathCtl, String pathShell, String gs ) {
1057          //System.out.println("GrADS #1070: procScript(env, " + pathCtl + ", " + pathShell + ", " + gs + ")");
1058          message = "" ;
1059          String  cmdString   = command + " -blc /"run " + gs + "/"" ;
1060          String  cmdStrings []     = { command, "-blc" , "/"run " + gs + "/"" } ;
1061          //System.out.println(cmdStrings[0] + " " + cmdStrings[1] + " " + cmdStrings[2]);
1062
1063          try {
1064              File    f = new File ( pathShell ) ;
1065              if ( f.isDirectory () && f.exists () ) {
1066                  //System.out.println("GrADS #1079: " + pathShell);
1067                  FileWriter      fw = new FileWriter ( pathShell + "/gradsc.sh" ) ;
1068                  BufferedWriter  bw = new BufferedWriter ( fw ) ;
1069                  bw.write ( "date" ) ;
1070                  bw.newLine () ;
1071                  bw.write ( "source /etc/profile" ) ; //确保 .sh 的执行环境与人工执行时一样
1072                  bw.newLine () ;
1073                  bw.write ( "source $HOME/.bash_profile" ) ; //确保 .sh 的执行环境与人工执行时一样
1074                  bw.newLine () ;
1075                  bw.write ( "cd " + pathCtl ) ;
1076                  bw.newLine () ;
1077                  bw.write ( cmdString ) ;
1078                  bw.newLine () ;
1079                  bw.close () ;
1080                  fw.close () ;
1081                  String stringLine;
1082                  Process cmd = Runtime.getRuntime () .exec ( "/bin/chmod +x " + pathShell + "/gradsc.sh" , env, f ) ;
1083                  cmd.waitFor () ;
1084 /*
1085                  cmd = Runtime.getRuntime().exec("/bin/ls -ltr " + pathShell, env, f);//测试
1086                  BufferedReader  bufferedReader = new BufferedReader( new InputStreamReader(cmd.getInputStream()));
1087                  while(null != (stringLine = bufferedReader.readLine()) ) {
1088                      System.out.println(stringLine);
1089                  }
1090 System.out.println("/bin/bash " + pathShell + "/gradsc.sh");
1091 */
1092                  cmd = Runtime.getRuntime () .exec ( "/bin/bash " + pathShell + "/gradsc.sh" , env, f ) ;
1093                  BufferedReader  bufferedReader = new BufferedReader ( new InputStreamReader ( cmd.getInputStream ())) ;
1094                  while ( null != ( stringLine = bufferedReader.readLine ()) ) {
1095                      //stringLine  = stringLine.replaceAll(" ", " ");
1096                      //System.out.println(stringLine);
1097                      message     = message + stringLine + "/n" ;
1098                  }
1099 //              message = message.replaceAll("/n", "<br>");
1100                  bufferedReader.close () ;
1101 //              cmd = Runtime.getRuntime().exec("/bin/rm -f " + pathShell + "/grads.sh", env, f);
1102              }
1103          }
1104          catch ( Exception ex ) {
1105              System.out.println ( "GrADS Exception: " + ex.getMessage ()) ;
1106              ex.printStackTrace () ;
1107          }
1108          return ( message ) ;
1109      } //调用 GrADS 脚本
1110 //删除图像文件，在 JSP 中当获得的要素、时次、层次等三个参数的值均为空时执行本函数，以避免保留过多的垃圾文件。
1111      public static int   deleteImageFile ( String path, String ext ) { //suffix=.ext
1112          if ( null == ext || "" .equals ( ext.trim ()) ) {
1113              return ( 0 ) ;
1114          }
1115          boolean isImageFile = false ;
1116          String  exts [] = { ".jpg" , ".jpeg" , ".gif" , ".bmp" , ".dib" , ".png" , ".ico" } ;
1117          for ( int i= 0 ;i<exts.length;i++ ) { //确保删除的是图像文件
1118              if ( 0 == ext.compareToIgnoreCase ( exts [ i ]) ) {
1119                  i   = exts.length;
1120                  isImageFile = true ;
1121              }
1122          }
1123          if ( false == isImageFile ) {
1124              return ( 0 ) ;
1125          }
1126          File    f   = new File ( path ) ;
1127          int      iCount  = 0 ;
1128          if ( f.exists () && f.isDirectory () ) {
1129              File    listFiles [] = f.listFiles () ;
1130              for ( int i= 0 ;i<listFiles.length;i++ ) {
1131                  if ( !listFiles [ i ] .isDirectory () && listFiles [ i ] .getName () .endsWith ( ext ) ) {
1132                      iCount  = iCount + ( listFiles [ i ] .delete () ? 1 : 0 ) ;
1133                  }
1134              }
1135          }
1136          return ( iCount ) ;
1137      }
1138 /**
1139 * 功能：列出所有文件
1140 * 参数：
1141 *      result      - 存放搜索结果
1142 *      path        - 目录
1143 *      prefix      - 文件名前缀
1144 *      suffix      - 文件后缀
1145 *      recrise     - 是否包括子目录
1146 * 返回值：
1147 *      找到的文件总数
1148 */
1149      public static int listGrADS ( Vector result, String path, String prefix, String suffix, boolean recuise ) {
1150          File    f   = new File ( path ) ;
1151          int      resultCount = result.size () ;
1152          if ( f.exists () && f.isDirectory () ) {
1153              File    listFiles [] = f.listFiles () ;
1154              for ( int i= 0 ;i<listFiles.length;i++ ) {
1155                  //System.out.println(listFiles[i].getAbsolutePath());
1156                  if ( recuise && listFiles [ i ] .isDirectory () ) {
1157                      listGrADS ( result, listFiles [ i ] .getAbsolutePath () , prefix, suffix, recuise ) ;
1158                  }
1159                  if ( !listFiles [ i ] .isDirectory () &&
1160                          listFiles [ i ] .getName () .startsWith ( prefix ) &&
1161                          listFiles [ i ] .getName () .endsWith ( suffix ) ) {
1162                      result.add ( listFiles [ i ] .getAbsolutePath ()) ;
1163                  }
1164              }
1165          }
1166          return ( result.size () - resultCount ) ;
1167      }
1168 //获得数据 - 未完成
1169      public int getData ( String file, float [][] data, int lon, int lat ) {
1170          int      count   = 0 ;
1171          byte []   b       = new byte [ 4 ] ;
1172          try {
1173              RandomAccessFile    raf = new RandomAccessFile ( file, "r" ) ;
1174              DecimalFormat       df  = new DecimalFormat ( "0.000000" ) ;
1175
1176              String  s   = "                " ;
1177              String  s1  = "" ;
1178              for ( int j= 0 ;j<lat;j++ ) {
1179                  for ( int i= 0 ;i<lon;i++ ) {
1180                      raf.read ( b ) ;
1181                      data [ i ][ j ]   = DataConverter.getFloat ( b ) ;
1182                      s1          = s + df.format ( data [ i ][ j ]) ;
1183                      System.out.print (( i% 10 == 0 ? "/r/n" : "" ) + s1.substring ( s1.length () - 16 )) ;
1184                      count++;
1185                  }
1186                  System.out.println ( "" ) ;
1187              }
1188              raf.close () ;
1189          }
1190          catch ( Exception e ) {
1191              System.out.println ( e.getMessage ()) ;
1192              e.printStackTrace () ;
1193          }
1194          return ( count ) ;
1195      }
1196 }