#!/usr/bin/ruby
# -*- coding: utf-8 -*-
# 重力を変えたときに表面温度がどのくらい変化するか, の絵
# 具体的には
# Ts_{g=1.46g*} - Ts_{g=1.0g*} の絵を書く

require "numru/ggraph"
require 'fileutils'
include NumRu
WSN=2
gp1 = GPhys::IO.open('E_T42L26_sync_R1.46/SurfTemp.nc','SurfTemp')
gp2 = GPhys::IO.open('E_T42L26_sync_1.46Rg/SurfTemp.nc','SurfTemp')
gp1 = gp1.cut('time'=>725..1095).mean('time')
gp2 = gp2.cut('time'=>300..670).mean('time')
gp= gp2 - gp1

DCL.gropn(WSN)
DCL.sgpset('lcntl',false)
GGraph.set_fig('itr'=>1, 'viewport'=>[0.2,0.8,0.3,0.6] )
GGraph.tone_and_contour(gp,true,'max'=>5,'min'=>-5,'annotate'=>false,'color_bar'=>true)
DCL.grcls