gphys 1.1.1a

data/sample/druby_cli2.rb

@@ -0,0 +1,34 @@
+# druby_cli2.rb -- a sample dRuby client
+#
+# USAGE
+#   See ./druby_serv2.rb
+
+require "drb/drb"
+require "numru/ggraph"
+include NumRu
+
+DRb.start_service
+uri = ARGV.shift || raise("Usage: % #{$0} uri")
+gp = DRbObject.new(nil, uri)
+
+p gp.class, gp.name, gp.rank, gp.length
+p gp.grid.axis(0).pos.name
+
+gpmean = gp.mean(0)
+p gpmean.class, gpmean.name, gpmean.rank
+
+class NArray
+  def self._load(o) to_na(*Marshal::load(o)).ntoh end
+end
+
+p ( gpmean_copy = gpmean.copy )
+DCL.gropn(1)
+DCL.sgpset('lfull',true)
+DCL.sgpset('lcntl',false)
+DCL.uzfact(0.7)
+GGraph.set_fig('viewport'=>[0.25,0.7,0.15,0.6], 'itr'=>2 )
+GGraph.contour( gpmean_copy )
+DCL.grcls
+
+print "\nThe following will be a DRbObject rather than GPhys since the size is too large\n(see druby_serv2.rb):\n"
+p ( gpcopy = gp.copy )

data/sample/druby_serv1.rb

@@ -0,0 +1,30 @@
+=begin
+=druby_serv1.rb
+
+a sample dRuby server
+
+==USAGE
+
+To start the server:
+
+     % ruby druby_serv1.rb
+
+This will print the address of the server such as
+druby://[hostname]:[port] (For exaple druby://horihost:39391
+
+Then run the client as (if the address is druby://horihost:39391):
+
+     % ruby druby_cli1.rb druby://horihost:39391
+=end
+
+require "drb/drb"
+require "numru/gphys"
+include NumRu
+
+file = NetCDF.open("../testdata/T.jan.nc")
+gp = GPhys::NetCDF_IO.open(file,"T")
+
+DRb.start_service(nil, gp)
+puts 'URI: '+DRb.uri
+puts '[return] to exit'
+gets

data/sample/druby_serv2.rb

@@ -0,0 +1,64 @@
+=begin
+=druby_serv2.rb
+
+a sample dRuby server. The difference from druby_serv1.rb is that
+a couple of methods are added to NArray, which enables us to serialize
+NArray with Marshal.dump. Therefore, a GPhys object can be transferred by 
+value (not by reference) if all the VArray in it is based on NArray.
+This can be done with the copy method of GPhys. See druby_cli2.rb
+for an example.
+
+==USAGE
+
+To start the server:
+
+     % ruby druby_serv2.rb [port]
+
+This will print the address of the server such as
+druby://[hostname]:[port] (For exaple druby://horihost:39391
+
+Then run the client as (if the address is druby://horihost:39391):
+
+     % ruby druby_cli2.rb druby://horihost:39391
+=end
+
+require "drb/drb"
+require "numru/gphys"
+include NumRu
+
+def usage
+  print <<-EOS
+
+  USAGE:
+    % ruby #{$0} [port]
+
+    Here, port (optional) is the port number to assign.
+  EOS
+  raise RuntimeError
+end
+
+port = ARGV.shift
+usage if port && port.to_i.to_s != port
+usage if ARGV.length > 0
+
+class NArray
+  DUMP_SIZE_LIMIT = 5000
+  def self._load(o) to_na(*Marshal::load(o)).ntoh end
+  def _dump(limit) 
+    if size <= DUMP_SIZE_LIMIT
+      Marshal::dump([hton.to_s, typecode, *shape]) 
+    else
+      raise "size of the NArray (#{size}) is too large to dump "+
+            "(limit: #{DUMP_SIZE_LIMIT})"
+    end
+  end
+end
+
+file = NetCDF.open("../testdata/T.jan.nc")
+gp = GPhys::NetCDF_IO.open(file,"T")
+
+uri_seed = ( port ? 'druby://:'+port : nil )   # 'druby://:'+port, or nil
+DRb.start_service(uri_seed, gp)
+puts 'URI: '+DRb.uri
+puts '[return] to exit'
+gets

data/sample/ep_flux/demo_NCEP_1.rb

@@ -0,0 +1,48 @@
+=begin
+=demo_NCEP_1.rb 
+
+  a demo program for GPhys::EP_Flux. show sample code to for calculating
+  EP-Flux and other physical values, and save as a NetCDF file.
+
+USAGE
+  % ruby demo_NCEP_1.rb
+
+Note
+  This program requires testdata. The URL is 
+
+    http://dennou-k.gfd-dennou.org/arch/ruby/products/ep_flux/testdata.tar.gz
+
+=end
+
+require 'numru/gphys'
+require 'numru/gphys/ep_flux'
+include NumRu
+
+datadir = "./testdata/"
+gp_u =     GPhys::IO.open( datadir+"UWND_NCEP.nc",  "uwnd")
+gp_v =     GPhys::IO.open( datadir+"VWND_NCEP.nc",  "vwnd")
+gp_omega = GPhys::IO.open( datadir+"OMEGA_NCEP.nc", "omega")
+gp_t =     GPhys::IO.open( datadir+"TEMP_NCEP.nc",  "temp")
+
+epflx_y, epflx_z, v_rmean, w_rmean, gp_lat, = ary =
+          GPhys::EP_Flux::ep_full_sphere(gp_u, gp_v, gp_omega, gp_t, true)
+gp_lat.rename('phi')
+
+ary.each{|gp|                                  #  This part will not
+  gp.data.att_names.each{|nm|                  #  be needed in future.
+    gp.data.del_att(nm) if /^valid_/ =~ nm     #  (Even now, it is not
+  }                                            #  needed if the valid
+}                                              #  range is wide enough)
+
+epflx_div = GPhys::EP_Flux::div_sphere(epflx_y, epflx_z)
+strm_rmean = GPhys::EP_Flux::strm_rmean(v_rmean)
+
+ofile = NetCDF.create("epflx_NCEP.nc")
+
+ary.each{|gp|
+  GPhys::IO.write(ofile, gp)
+}
+GPhys::IO.write(ofile, epflx_div)
+GPhys::IO.write(ofile, strm_rmean)
+
+ofile.close

data/sample/ep_flux/demo_NCEP_2.rb

@@ -0,0 +1,57 @@
+=begin
+=demo_NCEP_2.rb 
+
+  a demo program with GPhys::NetCDF_IO::each_along_dims_write, which 
+  is Iterator to process GPhys objects too big to read on memory at once.
+  this program show same performance with demo_NCEP_1.rb.
+
+USAGE
+  % ruby demo_NCEP_2.rb
+
+Note
+  This program requires testdata. The URL is written in demo_NCEP_1.rb
+
+=end
+
+require 'numru/gphys'
+require 'numru/gphys/ep_flux'
+include NumRu
+
+datadir = "./testdata/"
+gp_u =     GPhys::IO.open( datadir+"UWND_NCEP.nc",  "uwnd")
+gp_v =     GPhys::IO.open( datadir+"VWND_NCEP.nc",  "vwnd")
+gp_omega = GPhys::IO.open( datadir+"OMEGA_NCEP.nc", "omega")
+gp_t =     GPhys::IO.open( datadir+"TEMP_NCEP.nc",  "temp")
+
+ofile = NetCDF.create("epflx_NCEP.nc")
+
+nt = gp_u.shape[-1]
+i = 0
+GPhys::IO.each_along_dims_write([gp_u, gp_v, gp_omega, gp_t], ofile, -1){
+  |u, v, omega, t|
+  i += 1
+  print "processing #{i} / #{nt} ..\n" if (i % (nt/20+1))==1
+
+  epflx_y, epflx_z, v_rmean, w_rmean, gp_lat, gp_z, u_mean, theta_mean,
+        uv_dash, vt_dash, uw_dash, dtheta_dz = ary =
+                GPhys::EP_Flux::ep_full_sphere(u, v, omega, t, true)
+
+  ary.each{|gp|                                  #  This part will not
+    gp.data.att_names.each{|nm|                  #  be needed in future.
+      gp.data.del_att(nm) if /^valid_/ =~ nm     #  (Even now, it is not
+    }                                            #  needed if the valid
+  }                                              #  range is wide enough)
+
+  epflx_div = GPhys::EP_Flux::div_sphere(epflx_y, epflx_z)
+  strm_rmean = GPhys::EP_Flux::strm_rmean(v_rmean)
+
+  if i==1    # time independent => write only once
+    gp_lat.rename('phi')
+    GPhys::IO.write(ofile, gp_lat)
+    GPhys::IO.write(ofile, gp_z)
+  end
+  [ epflx_y, epflx_z, v_rmean, w_rmean, strm_rmean, epflx_div, u_mean, theta_mean,
+    uv_dash, vt_dash, uw_dash, dtheta_dz ]
+}
+
+ofile.close

data/sample/ep_flux/demo_NCEP_3.rb

@@ -0,0 +1,81 @@
+=begin
+=demo_NCEP_3.rb 
+
+  a demo program for drawing EP-Flux vector and Residual mean circulation.
+
+USAGE
+  % ruby demo_NCEP_3.rb
+
+Note
+  This program requires testdata. The URL is written in demo_NCEP_1.rb
+
+=end
+
+require 'numru/gphys'
+require 'numru/gphys/ep_flux'
+require 'numru/ggraph_on_merdional_section'
+include NumRu
+
+epflx_fnm = './epflx_NCEP.nc'
+
+if File::exist?(epflx_fnm)
+
+  epflx_y =   GPhys::IO.open(epflx_fnm,  'epflx_y')
+  epflx_z =   GPhys::IO.open(epflx_fnm,  'epflx_z')
+  epflx_div = GPhys::IO.open(epflx_fnm,  'epflx_div')
+  v_rmean =   GPhys::IO.open(epflx_fnm,  'v_rmean')
+  w_rmean =   GPhys::IO.open(epflx_fnm,  'w_rmean')
+  strm_rmean =   GPhys::IO.open(epflx_fnm,  'strm_rmean')
+
+else
+
+  datadir = "./testdata/"
+  gp_u =     GPhys::IO.open( datadir+"UWND_NCEP.nc",  "uwnd")
+  gp_v =     GPhys::IO.open( datadir+"VWND_NCEP.nc",  "vwnd")
+  gp_omega = GPhys::IO.open( datadir+"OMEGA_NCEP.nc", "omega")
+  gp_t =     GPhys::IO.open( datadir+"TEMP_NCEP.nc",  "temp")
+  
+  epflx_y, epflx_z, v_rmean, w_rmean= ary =
+    GPhys::EP_Flux::ep_full_sphere(gp_u, gp_v, gp_omega, gp_t, true)
+
+  ary.each{|gp|                                  #  This part will not
+    gp.data.att_names.each{|nm|                  #  be needed in future.
+      gp.data.del_att(nm) if /^valid_/ =~ nm     #  (Even now, it is not
+    }                                            #  needed if the valid
+  }                                              #  range is wide enough)
+
+  epflx_div = GPhys::EP_Flux::div_sphere(epflx_y, epflx_z)
+  strm_rmean = GPhys::EP_Flux::strm_rmean(v_rmean)
+  
+end
+  
+
+DCL.gropn(1)
+DCL::sglset('LFULL', true)                       # use full area
+DCL::slrat(1.0, 0.85)                            # set aspect ratio of drwable area
+DCL.sgpset('lcntl', false)                       # don't rede control character.
+DCL.sgpset('lfprop',true)                        # use prportional font
+DCL.uzfact(0.6)                                  # set character size
+
+## show vector (Fy, Fz)
+fy = epflx_y.mean('time')                        
+fz = epflx_z.mean('time')                        
+epdiv = epflx_div.mean('time')                   
+
+GGraph::set_fig('view'=>[0.15, 0.85, 0.25, 0.55])
+GGraph::set_fig('itr'=>2)                        # contour && tone only
+GGraph::tone(epdiv)                              # tone
+GGraph::contour(epdiv,false)                           # contour
+xfact, yfact = GGraph::vector_on_merdional_section(fy, fz, false,
+   'fact'=>3.0,'xintv'=>1,'unit'=>true, 'annot'=>false
+				    )
+
+## show residual mean merdional circulation (cut nearly surface and Antarctica)
+vrm =  v_rmean.cut('lat'=>-70..90, 'level'=>850..100).mean('time')
+wrm =  w_rmean.cut('lat'=>-70..90, 'level'=>850..100).mean('time')
+strm = strm_rmean.cut('lat'=>-70..90, 'level'=>850..100).mean('time')
+GGraph::contour(strm, true, 'nlev'=>25)
+GGraph::vector_on_merdional_section(vrm, wrm, false,
+	  	        'fact'=>2.0,'xintv'=>1,'unit'=>true, 'annot'=>false
+				    )
+DCL.grcls

data/sample/ggraph_latlon_labelling_dr002690.rb

@@ -0,0 +1,159 @@
+## from [dennou-ruby:002690] by koshiro
+## Sample of axis-drwaing options: 'xmaplabel', 'ymaplabel' in GGraph,
+## which labels axes bettern than in the default drwaing when itr==10
+## uses 90E, 90W, etc instead of 90, 270, etc.
+
+require 'numru/ggraph'
+include NumRu
+
+wsn = ( ARGV[0] ? ARGV[0].to_i : 1 )
+
+path = '../testdata/T.jan.nc'
+var = 'T'
+gp = GPhys::IO.open(path,var)
+
+DCL.swpset('ldump',true) if wsn == 4
+DCL.gropn(wsn)
+DCL.sldiv('y',3,2)
+#DCL.uzfact(0.8)
+DCL.sgpset('lcntl',false)
+
+GGraph.set_map('coast_world'=>true)
+
+gpjpn = gp.cut(110..160,10..70,false)
+
+itr = 10
+
+# contour
+GGraph.next_fig('itr'=>itr)
+GGraph.contour(gp, true,'map_axes'=>true,'title'=>'default axes')
+
+GGraph.next_fig('itr'=>itr)
+GGraph.next_axes('xmaplabel'=>'lon','ymaplabel'=>'lat')
+GGraph.contour(gp, true,'map_axes'=>true,'title'=>'lon/lat axes')
+
+GGraph.next_fig('itr'=>itr)
+GGraph.next_axes('xmaplabel'=>'lon','ymaplabel'=>'lat','xtickint'=>30, 'ytickint'=>10, 'xlabelint'=>90, 'ylabelint'=>30)
+GGraph.contour(gp, true,'map_axes'=>true,'title'=>'lon/lat axes (set interval)')
+
+# tone
+GGraph.next_fig('itr'=>itr)
+GGraph.tone(gp, true,'map_axes'=>true,'title'=>'default axes')
+
+GGraph.next_fig('itr'=>itr)
+GGraph.next_axes('xmaplabel'=>'lon','ymaplabel'=>'lat')
+GGraph.tone(gp, true,'map_axes'=>true,'title'=>'lon/lat axes')
+
+GGraph.next_fig('itr'=>itr)
+GGraph.next_axes('xmaplabel'=>'lon','ymaplabel'=>'lat','xtickint'=>30, 'ytickint'=>10, 'xlabelint'=>90, 'ylabelint'=>30)
+GGraph.tone(gp, true,'map_axes'=>true,'title'=>'lon/lat axes (set interval)')
+
+
+# tone & contour
+GGraph.next_fig('itr'=>itr)
+GGraph.tone(gpjpn, true,'map_axes'=>true,'title'=>'default axes')
+GGraph.contour( gpjpn, false )
+
+GGraph.next_fig('itr'=>itr)
+GGraph.next_axes('xmaplabel'=>'lon','ymaplabel'=>'lat')
+GGraph.tone(gpjpn, true,'map_axes'=>true,'title'=>'lon/lat axes')
+GGraph.contour( gpjpn, false )
+
+GGraph.next_fig('itr'=>itr)
+GGraph.next_axes('xmaplabel'=>'lon','ymaplabel'=>'lat','xtickint'=>10, 'ytickint'=>10, 'xlabelint'=>30, 'ylabelint'=>30)
+GGraph.tone(gpjpn, true,'map_axes'=>true,'title'=>'lon/lat axes (set interval)')
+GGraph.contour( gpjpn, false )
+
+
+# vector ('flow_vect'=>false)
+GGraph.next_fig('itr'=>itr)
+GGraph.vector(gpjpn,gpjpn, true,'map_axes'=>true,'title'=>'default axes','flow_vect'=>false,'unit_vect'=>true)
+
+GGraph.next_fig('itr'=>itr)
+GGraph.next_axes('xmaplabel'=>'lon','ymaplabel'=>'lat')
+GGraph.vector(gpjpn,gpjpn, true,'map_axes'=>true,'title'=>'lon/lat axes','flow_vect'=>false,'unit_vect'=>true)
+
+GGraph.next_fig('itr'=>itr)
+GGraph.next_axes('xmaplabel'=>'lon','ymaplabel'=>'lat','xtickint'=>10, 'ytickint'=>10, 'xlabelint'=>30, 'ylabelint'=>30)
+GGraph.vector(gpjpn,gpjpn, true,'map_axes'=>true,'title'=>'lon/lat axes (set interval)','flow_vect'=>false,'unit_vect'=>true)
+
+# vector ('flow_vect'=>true)
+GGraph.next_fig('itr'=>itr)
+GGraph.vector(gpjpn,gpjpn, true,'map_axes'=>true,'title'=>'default axes','unit_vect'=>true)
+
+GGraph.next_fig('itr'=>itr)
+GGraph.next_axes('xmaplabel'=>'lon','ymaplabel'=>'lat')
+GGraph.vector(gpjpn,gpjpn, true,'map_axes'=>true,'title'=>'lon/lat axes','unit_vect'=>true)
+
+GGraph.next_fig('itr'=>itr)
+GGraph.next_axes('xmaplabel'=>'lon','ymaplabel'=>'lat','xtickint'=>10, 'ytickint'=>10, 'xlabelint'=>30, 'ylabelint'=>30)
+GGraph.vector(gpjpn,gpjpn, true,'map_axes'=>true,'title'=>'lon/lat axes (set interval)','unit_vect'=>true)
+
+DCL.grfrm
+DCL.grfrm
+DCL.grfrm
+
+### DCLExt::lon_ax, DCLExt::lat_ax
+
+DCL.grfrm
+DCL.grswnd(0.0, 360.0, -90.0, 90.0)
+DCL.grsvpt(0.15, 0.95, 0.3, 0.7)
+DCL.grstrn(1)
+DCL.grstrf
+DCLExt.lon_ax
+DCLExt.lat_ax
+DCLExt.lon_ax('cside'=>'t')
+DCLExt.lat_ax('cside'=>'r')
+
+DCL.grfrm
+DCL.grswnd(0.0, 360.0, -90.0, 90.0)
+DCL.grsvpt(0.15, 0.95, 0.3, 0.7)
+DCL.grstrn(1)
+DCL.grstrf
+DCLExt.lon_ax('dtick1'=>30, 'dtick2'=>90)
+DCLExt.lat_ax('dtick1'=>10, 'dtick2'=>30)
+DCLExt.lon_ax('cside'=>'t', 'dtick1'=>30, 'dtick2'=>90)
+DCLExt.lat_ax('cside'=>'r', 'dtick1'=>10, 'dtick2'=>30)
+
+DCL.grfrm
+DCL.grswnd(-90.0, 90.0, 0.0, 360.0)
+DCL.grsvpt(0.15, 0.95, 0.3, 0.7)
+DCL.grstrn(1)
+DCL.grstrf
+DCLExt.lon_ax('yax'=>true,'dtick1'=>30, 'dtick2'=>90)
+DCLExt.lat_ax('xax'=>true,'dtick1'=>10, 'dtick2'=>30)
+DCLExt.lon_ax('yax'=>true,'cside'=>'r', 'dtick1'=>30, 'dtick2'=>90)
+DCLExt.lat_ax('xax'=>true,'cside'=>'t', 'dtick1'=>10, 'dtick2'=>30)
+
+DCL.grfrm
+DCL.grswnd(90.0, -90.0, 360.0, 0.0)
+DCL.grsvpt(0.15, 0.95, 0.3, 0.7)
+DCL.grstrn(1)
+DCL.grstrf
+DCLExt.lon_ax('yax'=>true,'dtick1'=>30, 'dtick2'=>90)
+DCLExt.lat_ax('xax'=>true,'dtick1'=>10, 'dtick2'=>30)
+DCLExt.lon_ax('yax'=>true,'cside'=>'r', 'dtick1'=>30, 'dtick2'=>90)
+DCLExt.lat_ax('xax'=>true,'cside'=>'t', 'dtick1'=>10, 'dtick2'=>30)
+
+DCL.grfrm
+DCL.grswnd(-90.0, 90.0, 1000.0, 1.0)
+DCL.grsvpt(0.15, 0.95, 0.3, 0.7)
+DCL.grstrn(2)
+DCL.grstrf
+DCL.ulpset('iytype',3)
+DCL.ulylog('l', 3, 9)
+DCL.ulylog('r', 3, 9)
+DCLExt.lat_ax('xax'=>true,'dtick1'=>10, 'dtick2'=>30)
+DCLExt.lat_ax('xax'=>true,'cside'=>'t', 'dtick1'=>10, 'dtick2'=>30)
+
+DCL.grfrm
+DCL.grswnd(-180.0, 180.0, -90.0, 90.0)
+DCL.grsvpt(0.15, 0.95, 0.3, 0.7)
+DCL.grstrn(1)
+DCL.grstrf
+DCLExt.lon_ax('dtick1'=>30, 'dtick2'=>90)
+DCLExt.lat_ax('dtick1'=>10, 'dtick2'=>30)
+DCLExt.lon_ax('cside'=>'t', 'dtick1'=>30, 'dtick2'=>90)
+DCLExt.lat_ax('cside'=>'r', 'dtick1'=>10, 'dtick2'=>30)
+
+DCL.grcls

data/sample/ggraph_mapfit-axes_dr002687.rb

@@ -0,0 +1,131 @@
+## from [dennou-ruby:002687] by koshiro
+## simultaneous map-fitting and axes-drawing when itr=10
+
+require 'numru/ggraph'
+include NumRu
+
+wsn = ( ARGV[0] ? ARGV[0].to_i : 1 )
+
+path = '../testdata/T.jan.nc'
+var = 'T'
+gp = GPhys::IO.open(path,var)
+
+DCL.swpset('ldump',true) if wsn == 4
+DCL.gropn(wsn)
+DCL.sldiv('t',3,2)
+DCL.uzfact(0.8)
+DCL.sgpset('lcntl',false)
+
+GGraph.set_map('coast_world'=>true)
+
+gpjpn = gp.cut(110..160,10..70,false)
+
+itr = 10
+
+# contour
+GGraph.next_fig('itr'=>itr)
+GGraph.contour(gp,true,'title'=>'map_fit=T, map_axes=F (default)')
+
+GGraph.next_fig('itr'=>itr)
+GGraph.contour(gp, true,'map_axes'=>true,'title'=>'map_fit=T map_axes=T')
+
+GGraph.next_fig('itr'=>itr,'map_fit'=>false)
+GGraph.contour(gp,true,'map_axes'=>false,'title'=>'map_fit=F map_axes=F' )
+
+GGraph.next_fig('itr'=>itr,'map_fit'=>false)
+GGraph.contour(gp,true,'map_axes'=>true,'title'=>'map_fit=F map_axes=T')
+
+GGraph.next_fig('itr'=>itr,'map_fit'=>false,'map_axis'=>[0,0,0],'map_window'=>[-120,120,-60,60])
+GGraph.contour(gp,true,'map_axes'=>false,'title'=>'map_fit=F map_axes=F map_axis & map_window' )
+
+GGraph.next_fig('itr'=>itr,'map_fit'=>false,'map_axis'=>[0,0,0],'map_window'=>[-120,120,-60,60])
+GGraph.contour(gp,true,'map_axes'=>true,'title'=>'map_fit=F map_axes=T map_axis & map_window')
+
+
+# tone
+GGraph.next_fig('itr'=>itr)
+GGraph.tone(gp,true,'title'=>'map_fit=T, map_axes=F (default)')
+
+GGraph.next_fig('itr'=>itr)
+GGraph.tone(gp, true,'map_axes'=>true,'title'=>'map_fit=T map_axes=T')
+
+GGraph.next_fig('itr'=>itr,'map_fit'=>false)
+GGraph.tone(gp,true,'map_axes'=>false,'title'=>'map_fit=F map_axes=F' )
+
+GGraph.next_fig('itr'=>itr,'map_fit'=>false)
+GGraph.tone(gp,true,'map_axes'=>true,'title'=>'map_fit=F map_axes=T')
+
+GGraph.next_fig('itr'=>itr,'map_fit'=>false,'map_axis'=>[0,0,0],'map_window'=>[-120,120,-60,60])
+GGraph.tone(gp,true,'map_axes'=>false,'title'=>'map_fit=F map_axes=F map_axis & map_window' )
+
+GGraph.next_fig('itr'=>itr,'map_fit'=>false,'map_axis'=>[0,0,0],'map_window'=>[-120,120,-60,60])
+GGraph.tone(gp,true,'map_axes'=>true,'title'=>'map_fit=F map_axes=T map_axis & map_window')
+
+
+# tone & contour
+GGraph.next_fig('itr'=>itr)
+GGraph.tone(gpjpn,true,'title'=>'map_fit=T, map_axes=F (default)')
+GGraph.contour( gpjpn, false )
+
+GGraph.next_fig('itr'=>itr)
+GGraph.tone(gpjpn, true,'map_axes'=>true,'title'=>'map_fit=T map_axes=T')
+GGraph.contour( gpjpn, false )
+
+GGraph.next_fig('itr'=>itr,'map_fit'=>false)
+GGraph.tone(gpjpn,true,'map_axes'=>false,'title'=>'map_fit=F map_axes=F' )
+GGraph.contour( gpjpn, false )
+
+GGraph.next_fig('itr'=>itr,'map_fit'=>false)
+GGraph.tone(gpjpn,true,'map_axes'=>true,'title'=>'map_fit=F map_axes=T')
+GGraph.contour( gpjpn, false )
+
+GGraph.next_fig('itr'=>itr,'map_fit'=>false,'map_axis'=>[0,0,0],'map_window'=>[90,180,0,90])
+GGraph.tone(gpjpn,true,'map_axes'=>false,'title'=>'map_fit=F map_axes=F map_axis & map_window' )
+GGraph.contour( gpjpn, false )
+
+GGraph.next_fig('itr'=>itr,'map_fit'=>false,'map_axis'=>[0,0,0],'map_window'=>[90,180,0,90])
+GGraph.tone(gpjpn,true,'map_axes'=>true,'title'=>'map_fit=F map_axes=T map_axis & map_window')
+GGraph.contour( gpjpn, false )
+
+
+# vector ('flow_vect'=>false)
+GGraph.next_fig('itr'=>itr)
+GGraph.vector(gpjpn,gpjpn,true,'title'=>'map_fit=T, map_axes=F (default)','flow_vect'=>false,'unit_vect'=>true)
+
+GGraph.next_fig('itr'=>itr)
+GGraph.vector(gpjpn,gpjpn, true,'map_axes'=>true,'title'=>'map_fit=T map_axes=T','flow_vect'=>false,'unit_vect'=>true)
+
+GGraph.next_fig('itr'=>itr,'map_fit'=>false)
+GGraph.vector(gpjpn,gpjpn,true,'map_axes'=>false,'title'=>'map_fit=F map_axes=F' ,'flow_vect'=>false,'unit_vect'=>true)
+
+GGraph.next_fig('itr'=>itr,'map_fit'=>false)
+GGraph.vector(gpjpn,gpjpn,true,'map_axes'=>true,'title'=>'map_fit=F map_axes=T','unit_vect'=>true,'flow_vect'=>false,'unit_vect'=>true)
+
+GGraph.next_fig('itr'=>itr,'map_fit'=>false,'map_axis'=>[0,0,0],'map_window'=>[90,180,0,90])
+GGraph.vector(gpjpn,gpjpn,true,'map_axes'=>false,'title'=>'map_fit=F map_axes=F map_axis & map_window','flow_vect'=>false,'unit_vect'=>true)
+
+GGraph.next_fig('itr'=>itr,'map_fit'=>false,'map_axis'=>[0,0,0],'map_window'=>[90,180,0,90])
+GGraph.vector(gpjpn,gpjpn,true,'map_axes'=>true,'title'=>'map_fit=F map_axes=T map_axis & map_window','flow_vect'=>false,'unit_vect'=>true)
+
+
+# vector ('flow_vect'=>true)
+GGraph.next_fig('itr'=>itr)
+GGraph.vector(gpjpn,gpjpn,true,'title'=>'map_fit=T, map_axes=F (default)','unit_vect'=>true)
+
+GGraph.next_fig('itr'=>itr)
+GGraph.vector(gpjpn,gpjpn, true,'map_axes'=>true,'title'=>'map_fit=T map_axes=T','unit_vect'=>true)
+
+GGraph.next_fig('itr'=>itr,'map_fit'=>false)
+GGraph.vector(gpjpn,gpjpn,true,'map_axes'=>false,'title'=>'map_fit=F map_axes=F','unit_vect'=>true)
+
+GGraph.next_fig('itr'=>itr,'map_fit'=>false)
+GGraph.vector(gpjpn,gpjpn,true,'map_axes'=>true,'title'=>'map_fit=F map_axes=T','unit_vect'=>true,'unit_vect'=>true)
+
+GGraph.next_fig('itr'=>itr,'map_fit'=>false,'map_axis'=>[0,0,0],'map_window'=>[90,180,0,90])
+GGraph.vector(gpjpn,gpjpn,true,'map_axes'=>false,'title'=>'map_fit=F map_axes=F map_axis & map_window','unit_vect'=>true)
+
+GGraph.next_fig('itr'=>itr,'map_fit'=>false,'map_axis'=>[0,0,0],'map_window'=>[90,180,0,90])
+GGraph.vector(gpjpn,gpjpn,true,'map_axes'=>true,'title'=>'map_fit=F map_axes=T map_axis & map_window','unit_vect'=>true)
+
+
+DCL.grcls