bennyistanto

xkcd-style Pettitt's illustration on Homogeneity testing

The script below try to simulate how Pettitt test detects inhomogeneous or homogeneous on precipitation data. Here's the updated approach:

1. Daily Precipitation Data:

   • We have daily precipitation.

2. Generate Monthly Maximum Precipitation Values:

• Create synthetic monthly maxima with distinct differences to ensure inhomogeneity.

bennyistanto

import xarray as xr
import numpy as np
import os
import glob
import pandas as pd
from tqdm import tqdm
from scipy.stats import genextreme

# Define directories and filenames for threshold calculations
threshold_calculations = [

bennyistanto

import xarray as xr
import numpy as np
import os
import glob
import pandas as pd
from tqdm import tqdm

# Define the input and output directories
max_calculations = [
    {"days": 1, "input_dir": "/mnt/d/temp/imerg/data/fdd_1day", "output_dir": "/mnt/d/temp/imerg/data/fddmaxd_1day"},

bennyistanto

import xarray as xr
import os
import glob
from tqdm import tqdm
import pandas as pd

# Define directories and filenames for max value calculations
max_calculations = [
    {"days": 1, "input_dir": "/mnt/d/temp/imerg/data/fdd_1day", "output_dir": "/mnt/d/temp/imerg/data/fddmaxm_1day"},
    {"days": 2, "input_dir": "/mnt/d/temp/imerg/data/fdd_2day", "output_dir": "/mnt/d/temp/imerg/data/fddmaxm_2day"},

bennyistanto

import xarray as xr
import numpy as np
import os
import glob
from tqdm import tqdm
import pandas as pd

# Define the rolling day output directories
rolling_days = [2, 3, 4, 5]
input_dir = "/mnt/d/temp/imerg/data/fdd_1day"

bennyistanto

import xarray as xr
import os
import glob
from tqdm import tqdm

# Define the input and output directories
input_dir = "/mnt/d/temp/imerg/data/nc4_final"
output_dir = "/mnt/d/temp/imerg/data/fdd_1day"

# Check if input folder exists

bennyistanto

Custom xkcd-style map for Indonesia

Label setting for capital, country name and the ocean, and the red-dot is manually set, the value is calculate from top-left main map area.

You can paste the below code into an online Python compiler like https://python-fiddle.com/ and grab the result instantly.

bennyistanto

xkcd style for country map

This is example of Python script that could generate a xkcd-style for Country Map.

Example:

bennyistanto

xkcd style for LSEQM illustration

This is example of Python script that could generate a xkcd graph-style to illustrate concept of Linear Scaling and Empirical Quantile Mapping (LSEQM) method for daily precipitation bias correction.

You can paste the below code into an online Python compiler like https://python-fiddle.com/ and grab the result instantly.

version a

bennyistanto

<?xml version="1.0" encoding="utf-8"?>
<!-- Modified by: Benny Istanto/KLI/G2501222008, bennyistanto@apps.ipb.ac.id -->
<!-- English version available on: https://gist.github.com/bennyistanto/ef6b1a2bb1ed303a7832f789ac7e7a70 -->
<style xmlns="http://purl.org/net/xbiblio/csl" class="in-text" version="1.0" demote-non-dropping-particle="sort-only" default-locale="en-US">
  <!-- changed default-locale to en-US -->
  <info>
    <title>Institut Pertanian Bogor: Pedoman Penulisan Karya Ilmiah Edisi ke-4 (English)</title>
    <!-- updated title -->
    <title-short>IPB PPKI 4</title-short>
    <id>http://www.zotero.org/styles/institut-pertanian-bogor</id>