Skip to content

Instantly share code, notes, and snippets.

View raphaeldussin's full-sized avatar

raphael dussin raphaeldussin

98 followers · 37 following
View GitHub Profile
@raphaeldussin
raphaeldussin / func2_roms2zarr.py
Last active December 3, 2020 02:20
function for the roms2zarr tutorial
def select_interior(ds):
"""
discard "exterior" u,v,rho-points to build a symetric grid
Parameters:
ds (xarray.Dataset): ROMS dataset
"""
ds = ds.isel(xi_rho=slice(1,-1), eta_rho=slice(1,-1))
if 'xi_v' in ds.dims:
ds = ds.isel(xi_v=slice(1,-1))
@raphaeldussin
raphaeldussin / func1_roms2zarr.py
Created December 3, 2020 02:09
function for the roms2zarr tutorial
def remove_all_small_variables(ds):
""" remove all the variables that have less than 3 dimensions
Parameters:
ds (xarray.Dataset): ROMS dataset
"""
for v in ds.variables:
if v not in ['Cs_r', 'Cs_w', 'hc', 'Vtransform', 'ocean_time']:
ds = ds.drop_vars(v) if len(ds[v].dims) < 3 else ds
return ds
@raphaeldussin
raphaeldussin / install_FRE-NCtools_workstation.bash
Created September 1, 2020 20:26
#!/bin/bash
# for PPAN, replace gfdl-ws by gfdl
git clone https://github.com/NOAA-GFDL/FRE-NCtools
pushd FRE-NCtools
here=$(pwd)
export CONFIG_SITE=${here}/site-configs/gfdl-ws/config.site
source ${here}/site-configs/gfdl-ws/env.sh
@raphaeldussin
raphaeldussin / mom6_slurm_job.bash
Created July 6, 2020 20:51
run script for MOM6 using SLURM batch scheduler
#!/bin/bash
#SBATCH -n 240
#SBATCH -J MOM6
#SBATCH --error=MOM6.err
#SBATCH --output=MOM6.out
#SBATCH --exclusive
njobs=40
#--------------------------------- system settings ----- ---------------------------
@raphaeldussin
raphaeldussin / roms_to_zarr.py
Last active June 8, 2020 13:50
convert ROMS output to zarr
#!/usr/bin/env python
# useful ref:
# http://xarray.pydata.org/en/stable/examples/ROMS_ocean_model.html
import xarray as xr
import glob
# list of files to convert
@raphaeldussin
raphaeldussin / gist:e4cda7dbe54cf6a1033cd0662b3e74df
Created December 11, 2019 19:58
using gsw with xarray
# this is a wrong & stupid computation, just here for example
pt = xr.apply_ufunc(gsw.pt0_from_t, ds.so, ds.thetao, ds.z_l, dask='parallelized', output_dtypes=[ds.so.dtype])
@raphaeldussin
raphaeldussin / gist:49f00a3bd9214d93acb2b8e87411dca5
Created October 18, 2019 14:17
set time axis with cftime
pic_start = cftime.num2date(0, 'days since 0001-01-01')
pic_end = cftime.num2date(1200*365.25, 'days since 0001-01-01')
plt.figure(figsize=[8,6])
ax1 = plt.subplot(211)
plt.plot(ds_nadw_CM4_piC['time'], ds_nadw_CM4_piC['vol_nadw_atl'], 'k')
ax1.set_xlim([pic_start, pic_end])
ax2 = plt.subplot(212)
plt.plot(ds_nadw_CESM2_piC['time'], ds_nadw_CESM2_piC['vol_nadw_atl'], 'r')
ax2.set_xlim([pic_start, pic_end])
@raphaeldussin
raphaeldussin / apple_keyboard_linux.txt
Last active April 26, 2019 17:30
french apple keyboard on macbook pro running debian
# this works with the laptop's keyboard fine:
>>> more /etc/default/keyboard
# KEYBOARD CONFIGURATION FILE
# Consult the keyboard(5) manual page.
XKBMODEL="macintosh"
XKBLAYOUT="fr"
@raphaeldussin
raphaeldussin / keybase.md
Created April 23, 2019 21:47

Keybase proof

I hereby claim:

  • I am raphaeldussin on github.
  • I am raphaeldussin (https://keybase.io/raphaeldussin) on keybase.
  • I have a public key ASDelxZLvjXQhqoc0gM86v23mA7H7Z3QzcvuWKWvo51O0Ao

To claim this, I am signing this object:

@raphaeldussin
raphaeldussin / convert_float2double_mitgcm.py
Created December 16, 2018 23:20
#!/usr/bin/env python
import numpy as np
import struct
def write_new_data(data,fileout):
fid = open(fileout, "wb")
flatdata = data.flatten()
for kk in np.arange(len(flatdata)):
tmp = struct.pack('>d',flatdata[kk])