Remove model drift using control run simulations

Recently, I’ve been trying to remove model drift from my high-resolution CCSM 20th century climate change simulation. The model drift is estimated using the two companion control runs, HRC08 and HRC09. All three runs were spun up from a similar initial condition, with CO2 held fixed in year 2000 level. Following Gupta et al., [2013]1, we tried to identify the model drift by fitting a cubic polynomial to the full record of control runs (nearly 70 years).

Updated Hugo Academic theme to sync the upstream repo

Testing the new features of Hugo-Academic theme…

Using git-latexdiff to track changes in tex-file graphically

It’s been a while since I started using git. The version control system helps to track changes in computer files and coordinating work on those files among multiple people. Composing a manuscript is nothing unlike developing a program. If multiple authors work on the same manuscript, git can be very helpful to avoid conflicts. If you are the only author, git can help you to document what has been changed since the last commit or put your advisor’s comment in another branch.

Fill in the missing data using Python pandas

One of the many advantages of Python is its abundant and often powerful Libraries. For my research, besides plotting maps, I often play with time series. When it comes to manipulating and plotting time series, no other tools can beat python pandas.

Gibbs SeaWater Oceanographic Package for Python

Recently, I’m looking into the Meridional Heat Transport across the 35$^{\circ}$S in our high-resolution CCSM simulation following Dong et al. 2011. To calculate the heat transport, the following formula is used:

\begin{equation} \int \rho C_{p}v\theta dx dz \end{equation}

The tricky part is whether one should use temporal and spatial-varying values of $C_{p}$ (seawater heat capacity) and $\rho$ (seawater density). Since density and heat capacity are not among the standard outputs, I went on looking for some standard libraries to calculate them using the fields I have: temperature, pressure, and salinity. Soon, I found the Python version of Gibbs SeaWater Oceanographic Package.

iPython notebook gallery

Sometimes I got emails requesting codes to plot certain things. Those inquiries remind me that I have an older version of my website on still up and running. That was an answer to the requirement that every student in RSMAS should have their own websites back in 2014. To fill in the blanks, I listed out some of my python works in Ipython notebooks.

How to read the NWS hurricane warning

The University will be closed on both Thursday and Friday in preparation for the approaching Hurricane Matthew. Students and faculties were asked to leave RSMAS by 5 pm on Wednesday since wind shutters would completely seal the MSC building. The threat seems imminent; people are freak out – Publix was like a war zone when I was ordered to pick up some bottled water on my way home.

Automated deployment of Hugo generated site to Github pages

In the last post, I listed the typical workflow of generating a static website with Hugo. Here I would like to share the issues I encountered and the tweaks I’ve made. The source codes for this site can be found here. Basic setting in my config.toml for the Academic Theme Many beautifully designed themes are already available. I decided to use the academic theme, because it’s specifically designed for academic uses, showcasing publications, personal profile, and projects.

Switch from Pelican to Hugo

Being a big fan of all things Python, I built the last two incarnations of this site with Pelican. I even went to the length to set up a Travis CI to automatically build and deploy to GitHub pages. That was the beginning of the summer and only two new posts ever since. Recently, my boss has brought up the topic of building a group website. Also, there is a joint-project website needs to be done.

Workflow for a CCSM4 run

Basic approach, startup run:

  1. Create a new case using the script under ~/CCSM4/scripts, following the template: PATH, compset, resolution and machine file. Compset B is the fully coupled mode.

    ./create_newcase -case ~/ccsm4/$casename \
                           -compset B_2000 \
                           -res 0.9x1.25_gx1v6 \
                           -mach bluefire