Skip to content

Instantly share code, notes, and snippets.

@bekozi
bekozi / Foo.cpp
Created Oct 25, 2018 — forked from Luthaf/Foo.cpp
Calling C++ from Fortran
View Foo.cpp
#include "Foo.hpp"
#include <iostream>
using namespace std;
Foo::Foo(int _a, int _b): a(_a), b(_b){
cout << "C++ side, constructor" << endl;
}
Foo::~Foo(){
View iofoo.py
import netCDF4 as nc
import ocgis as ESMF
import numpy as np
from ocgis.vmachine.mpi import rank_print
NTIMESTEPS = 3
NELEMPERPET = 4
IO_TAG = 111
View ocgis_regrid_to_camels.py
import ocgis
from ocgis.test.base import create_gridxy_global, create_exact_field
# Path to the input ESRI Shapefile.
# PATH_SHP = '/home/benkoziol/l/data/bekozi-work/i53-jh-camel-watersheds/basin_set_full_res/HCDN_nhru_final_671.shp'
PATH_SHP = '/gcs/camels/basin_dataset_public_v1p2/shapefiles/merge/basinset_gf_nhru.shp'
# Path to example exact data for testing.
DATA_PATH = 'camels_source.nc'
@bekozi
bekozi / logging.bash
Created Mar 1, 2018 — forked from goodmami/logging.bash
Basic logging commands for Linux shell scripts
View logging.bash
#!/bin/bash
##
## Simple logging mechanism for Bash
##
## Author: Michael Wayne Goodman <goodman.m.w@gmail.com>
## Thanks: Jul for the idea to add a datestring. See:
## http://www.goodmami.org/2011/07/simple-logging-in-bash-scripts/#comment-5854
## Thanks: @gffhcks for noting that inf() and debug() should be swapped,
## and that critical() used $2 instead of $1
View spatial_coords_with_time.py
from shapely.geometry import Point
import ocgis
path = '/path/to/data.nc'
lat_variable = 'XLAT'
lon_variable = 'XLONG'
lat_dimension = 'south_north'
lon_dimension = 'west_east'
View xx_calc_modified_v2.py
import numpy
def get_masked_arr(arr, fill_val):
'''
If a masked array is passed, this function does nothing.
If a filled array is passed (fill_value must be passed also), it will be transformed into a masked array.
'''
if isinstance(arr, numpy.ma.MaskedArray): # numpy.ma.MaskedArray