Class: NumRu::Grib

Inherits:
Object
  • Object
show all
Defined in:
lib/numru/grib/grib.rb,
lib/numru/grib/version.rb,
ext/grib.c

Constant Summary collapse

VERSION = 
"0.4.0"
TYPE_UNDEFINED = 
INT2NUM(GRIB_TYPE_UNDEFINED)
TYPE_LONG = 
INT2NUM(GRIB_TYPE_LONG)
TYPE_DOUBLE = 
INT2NUM(GRIB_TYPE_DOUBLE)
TYPE_STRING = 
INT2NUM(GRIB_TYPE_STRING)
TYPE_BYTES = 
INT2NUM(GRIB_TYPE_BYTES)
TYPE_SECTION = 
INT2NUM(GRIB_TYPE_SECTION)
TYPE_LABEL = 
INT2NUM(GRIB_TYPE_LABEL)
TYPE_MISSING = 
INT2NUM(GRIB_TYPE_MISSING)
NEAREST_SAME_GRID = 
ULONG2NUM(GRIB_NEAREST_SAME_GRID)
NEAREST_SAME_DATA = 
ULONG2NUM(GRIB_NEAREST_SAME_DATA)
NEAREST_SAME_POINT = 
ULONG2NUM(GRIB_NEAREST_SAME_POINT)
KEYS_ITERATOR_ALL_KEYS = 
ULONG2NUM(GRIB_KEYS_ITERATOR_ALL_KEYS)
KEYS_ITERATOR_SKIP_READ_ONLY = 
ULONG2NUM(GRIB_KEYS_ITERATOR_SKIP_READ_ONLY)
SUPPORT_BIGMEM = 
Qfalse
@@gribs =

cache due to momery leak in the original grib_api library

Hash.new

Class Method Summary collapse

Instance Method Summary collapse

Constructor Details

#initialize(*args) ⇒ Object

NumRu::Grib#initialize(filename, [, mode])



384
385
386
387
388
389
390
391
392
393
394
395
396
397
398
399
400
401
402
403
404
405
406
407
408
409
410
411
412
413
414
415
416
417
418
419
420
421
422
423
424
425
426
427
428
429
430
431
432
433
434
435
436
437
438
439
440
441
442
443
444
445
446
447
448
449
450
451
452
453
454
455
456
457
458
459
460
461
462
463
464
465
466
467
468
469
470
471
472
 
# File 'ext/grib.c', line 384

static VALUE
rg_file_initialize(int argc, VALUE *argv, VALUE self)
{
  char *fname, *mode;
  VALUE rfname, rmode;
  rb_scan_args(argc, argv, "11", &rfname, &rmode);
  fname = StringValueCStr(rfname);
  if (rmode == Qnil)
    mode = "r";
  else
    mode = StringValueCStr(rmode);

  FILE *file = fopen(fname, mode);
  if (!file)
    rb_raise(rb_eRuntimeError, "unable to open file %s", fname);

  rg_file *gfile;
  Data_Get_Struct(self, rg_file, gfile);
  gfile->file = file;
  gfile->fname = ALLOC_N(char, strlen(fname)+1);
  strcpy(gfile->fname, fname);

  int error;
  grib_handle *handle;
  var_list *var = NULL;
  while ((handle = grib_handle_new_from_file(0, file, &error))) {
    raise_error(error, "grib_handle_new_from_file(0, file, &error)");
    if (handle == NULL)
      rb_raise(rb_eRuntimeError, "unable to create handle in %s", fname);
    push_msg_var(&var, handle, self);
  }
  gfile->var = var;
  var_list *var2;
  const char *vname;
  char buf[MAX_VALUE_LENGTH];
  int first_ltype, first_gtype;
  int i;
  while (var) {
    vname = var->vname;
    var2 = var->next;
    i = 0;
    first_ltype = 1;
    first_gtype = 1;
    while (var2) {
      if (strcmp(var2->vname, vname) == 0) {
	if (var->ltype_id != var2->ltype_id) {
	  if (first_ltype) {
	    if (strcmp(var->ltype,"isobaricInhPa") != 0 && strcmp(var->ltype,"pl") !=0) {
	      if (strcmp(var->ltype, "unknown")==0)
		sprintf(buf, "%s_level%ld", var->vname, var->ltype_id);
	      else
		sprintf(buf, "%s_%s", var->vname, var->ltype);
	      strcpy(var->vname, buf);
	    }
	    first_ltype = 0;
	  }
	  if (strcmp(var2->ltype,"isobaricInhPa") != 0 && strcmp(var2->ltype,"pl") != 0) {
	    if (strcmp(var2->ltype, "unknown")==0)
	      sprintf(buf, "%s_level%ld", var2->vname, var2->ltype_id);
	    else
	      sprintf(buf, "%s_%s", var2->vname, var2->ltype);
	    strcpy(var2->vname, buf);
	  }
	} else if (strcmp(var->gtype, var2->gtype) != 0) {
	  if (first_gtype) {
	    sprintf(buf, "%s_%s", var->vname, var->gtype);
	    strcpy(var->vname, buf);
	    first_gtype = 0;
	  }
	  sprintf(buf, "%s_%s", var2->vname, var2->gtype);
	  strcpy(var2->vname, buf);
	} else {
	  if (i==0) {
	    sprintf(buf, "%s_%d", var->vname, i);
	    strcpy(var->vname, buf);
	    i += 1;
	  }
	  sprintf(buf, "%s_%d", var2->vname, i);
	  strcpy(var2->vname, buf);
	  i += 1;
	}
      }
      var2 = var2->next;
    }
    var = var->next;
  }
  
  return self;
}

Class Method Details

._new ⇒ Object 



13
 
# File 'lib/numru/grib/grib.rb', line 13

alias :_new :new

.is_a_Grib?(fname) ⇒ Boolean

Returns:

  • (Boolean) 


6
7
8
9
10
11
12
 
# File 'lib/numru/grib/grib.rb', line 6

def is_a_Grib?(fname)
  is = false
  File.open(fname) do |file|
    is = (file.read(4) == "GRIB")
  end
  return is
end

.multi=(flag) ⇒ Object

NumRu::Grib.multi=(flag) -> True/False



369
370
371
372
373
374
375
376
377
378
379
 
# File 'ext/grib.c', line 369

static VALUE
rg_multi(VALUE self, VALUE flag)
{
  if (flag == Qtrue)
    grib_multi_support_on(0);
  else if (flag == Qfalse)
    grib_multi_support_off(0);
  else
    rb_raise(rb_eArgError, "flag must be true or false");
  return flag;
}

.new(fname, mode = "r") ⇒ Object Also known as: open 



14
15
16
17
18
19
20
21
22
23
24
 
# File 'lib/numru/grib/grib.rb', line 14

def new(fname, mode="r")
  key = [fname,mode]
  if (obj = @@gribs[key])
    obj[1] += 1
    return obj[0]
  else
    grib = _new(fname, mode)
    @@gribs[key] = [grib, 1]
    return grib
  end
end

Instance Method Details

#_close ⇒ Object

class << self



28
 
# File 'lib/numru/grib/grib.rb', line 28

alias :_close :close

#close ⇒ Object

NumRu::Grib#close() -> nil



476
477
478
479
480
481
482
483
484
485
486
487
488
 
# File 'ext/grib.c', line 476

def close
  @@gribs.each do |k, val|
    if self == val[0]
      val[1] -= 1
      if val[1] < 1
        _close
        @@gribs.delete(k)
        return true
      end
    end
  end
  return false
end

#inspect ⇒ Object 



43
44
45
 
# File 'lib/numru/grib/grib.rb', line 43

def inspect
  "Grib: #{path}"
end

#path ⇒ Object

NumRu::Grib#path() -> String



488
489
490
491
492
493
494
 
# File 'ext/grib.c', line 488

static VALUE
rg_file_path(VALUE self)
{
  rg_file *gfile;
  Data_Get_Struct(self, rg_file, gfile);
  return rb_str_new2(gfile->fname);
}

#var(rvname) ⇒ Object

NumRu::Grib#var(String) -> NumRu::GribVar



517
518
519
520
521
522
523
524
525
526
527
528
529
530
531
532
533
534
535
536
537
 
# File 'ext/grib.c', line 517

static VALUE
rg_file_var(VALUE self, VALUE rvname)
{
  char *vname = StringValueCStr(rvname);
  rg_file *gfile;
  Data_Get_Struct(self, rg_file, gfile);
  var_list *var_list = gfile->var;
  while (var_list) {
    if (strcmp(var_list->vname, vname) == 0) {
      VALUE rvar = var_alloc(cVar);
      rg_var *var;
      Data_Get_Struct(rvar, rg_var, var);
      var->file = self;
      var->var_list = var_list;
      rb_funcall(rvar, id_init, 0);
      return rvar;
    }
    var_list = var_list->next;
  }
  rb_raise(rb_eArgError, "cannot find variable: %s", vname);
}

#var_names ⇒ Object

NumRu::Grib#var_names() -> Array



499
500
501
502
503
504
505
506
507
508
509
510
511
 
# File 'ext/grib.c', line 499

static VALUE
rg_file_var_names(VALUE self)
{
  rg_file *gfile;
  Data_Get_Struct(self, rg_file, gfile);
  VALUE ary = rb_ary_new();
  var_list *var = gfile->var;
  while (var) {
    rb_ary_push(ary, rb_str_new2(var->vname));
    var = var->next;
  }
  return ary;
}