#!/usr/bin/env ruby
#
# 表題  Mitchell and Vallis (2010) 追試 : Figure 2 (U/mass stream function)
#
# 履歴  2026/10/13  竹広 真一
#
require "numru/ggraph"
include NumRu

TimeCut=1000

gtURL = sprintf('U_rank*.nc@U,time=%f',TimeCut)
u = GPhys::IO.open_gturl( gtURL ).mean('lon')

gtURL = sprintf('Msf_rank*.nc@Msf,time=%f',TimeCut)
msf = GPhys::IO.open_gturl( gtURL ).mean('lon')

IntU=10.0

DCL.gropn(1)
GGraph.tone msf, true
GGraph.contour u, false, 'int'=>IntU
DCL.grcls