gphys 1.5.0 → 1.5.1

data/sample/ep_flux/demo_NCEP_1.rb

@@ -1,48 +0,0 @@
-=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

@@ -1,57 +0,0 @@
-=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

@@ -1,81 +0,0 @@
-=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

@@ -1,159 +0,0 @@
-## 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

@@ -1,131 +0,0 @@
-## 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