gphys 1.2.2.1 → 1.4.3
This diff represents the content of publicly available package versions that have been released to one of the supported registries. The information contained in this diff is provided for informational purposes only and reflects changes between package versions as they appear in their respective public registries.
- checksums.yaml +7 -0
- data/.gitignore +8 -17
- data/.rspec +2 -0
- data/.travis.yml +3 -0
- data/ChangeLog +5762 -753
- data/LICENSE.txt +30 -18
- data/Rakefile +1 -0
- data/bin/console +14 -0
- data/bin/gpcat +43 -2
- data/bin/gpcut +16 -0
- data/bin/gpvect +167 -15
- data/bin/gpview +254 -51
- data/bin/setup +7 -0
- data/dim_op.c +1220 -0
- data/doc/attribute.html +19 -0
- data/doc/attributenetcdf.html +15 -0
- data/doc/axis.html +387 -0
- data/doc/coordmapping.html +111 -0
- data/doc/coordtransform.html +36 -0
- data/doc/dclext.html +821 -0
- data/doc/derivative/gphys-derivative.html +100 -0
- data/doc/derivative/index.html +21 -0
- data/doc/derivative/index.rd +14 -0
- data/doc/derivative/math-doc/document.pdf +0 -0
- data/doc/derivative/math-doc/document.tex +158 -0
- data/doc/derivative/math-doc/document/document.css +30 -0
- data/doc/derivative/math-doc/document/document.html +57 -0
- data/doc/derivative/math-doc/document/images.aux +1 -0
- data/doc/derivative/math-doc/document/images.log +385 -0
- data/doc/derivative/math-doc/document/images.pl +186 -0
- data/doc/derivative/math-doc/document/images.tex +364 -0
- data/doc/derivative/math-doc/document/img1.png +0 -0
- data/doc/derivative/math-doc/document/img10.png +0 -0
- data/doc/derivative/math-doc/document/img11.png +0 -0
- data/doc/derivative/math-doc/document/img12.png +0 -0
- data/doc/derivative/math-doc/document/img13.png +0 -0
- data/doc/derivative/math-doc/document/img14.png +0 -0
- data/doc/derivative/math-doc/document/img15.png +0 -0
- data/doc/derivative/math-doc/document/img16.png +0 -0
- data/doc/derivative/math-doc/document/img17.png +0 -0
- data/doc/derivative/math-doc/document/img18.png +0 -0
- data/doc/derivative/math-doc/document/img19.png +0 -0
- data/doc/derivative/math-doc/document/img2.png +0 -0
- data/doc/derivative/math-doc/document/img20.png +0 -0
- data/doc/derivative/math-doc/document/img21.png +0 -0
- data/doc/derivative/math-doc/document/img22.png +0 -0
- data/doc/derivative/math-doc/document/img23.png +0 -0
- data/doc/derivative/math-doc/document/img24.png +0 -0
- data/doc/derivative/math-doc/document/img25.png +0 -0
- data/doc/derivative/math-doc/document/img26.png +0 -0
- data/doc/derivative/math-doc/document/img27.png +0 -0
- data/doc/derivative/math-doc/document/img28.png +0 -0
- data/doc/derivative/math-doc/document/img29.png +0 -0
- data/doc/derivative/math-doc/document/img3.png +0 -0
- data/doc/derivative/math-doc/document/img30.png +0 -0
- data/doc/derivative/math-doc/document/img4.png +0 -0
- data/doc/derivative/math-doc/document/img5.png +0 -0
- data/doc/derivative/math-doc/document/img6.png +0 -0
- data/doc/derivative/math-doc/document/img7.png +0 -0
- data/doc/derivative/math-doc/document/img8.png +0 -0
- data/doc/derivative/math-doc/document/img9.png +0 -0
- data/doc/derivative/math-doc/document/index.html +57 -0
- data/doc/derivative/math-doc/document/labels.pl +13 -0
- data/doc/derivative/math-doc/document/next.png +0 -0
- data/doc/derivative/math-doc/document/next_g.png +0 -0
- data/doc/derivative/math-doc/document/node1.html +238 -0
- data/doc/derivative/math-doc/document/node2.html +75 -0
- data/doc/derivative/math-doc/document/prev.png +0 -0
- data/doc/derivative/math-doc/document/prev_g.png +0 -0
- data/doc/derivative/math-doc/document/up.png +0 -0
- data/doc/derivative/math-doc/document/up_g.png +0 -0
- data/doc/derivative/numru-derivative.html +158 -0
- data/doc/ep_flux/ep_flux.html +469 -0
- data/doc/ep_flux/ggraph_on_merdional_section.html +71 -0
- data/doc/ep_flux/index.html +31 -0
- data/doc/ep_flux/index.rd +24 -0
- data/doc/ep_flux/math-doc/document.pdf +0 -0
- data/doc/ep_flux/math-doc/document.tex +2018 -0
- data/doc/ep_flux/math-doc/document/WARNINGS +1 -0
- data/doc/ep_flux/math-doc/document/contents.png +0 -0
- data/doc/ep_flux/math-doc/document/crossref.png +0 -0
- data/doc/ep_flux/math-doc/document/document.css +30 -0
- data/doc/ep_flux/math-doc/document/document.html +101 -0
- data/doc/ep_flux/math-doc/document/images.aux +1 -0
- data/doc/ep_flux/math-doc/document/images.log +1375 -0
- data/doc/ep_flux/math-doc/document/images.pl +1328 -0
- data/doc/ep_flux/math-doc/document/images.tex +1471 -0
- data/doc/ep_flux/math-doc/document/img1.png +0 -0
- data/doc/ep_flux/math-doc/document/img10.png +0 -0
- data/doc/ep_flux/math-doc/document/img100.png +0 -0
- data/doc/ep_flux/math-doc/document/img101.png +0 -0
- data/doc/ep_flux/math-doc/document/img102.png +0 -0
- data/doc/ep_flux/math-doc/document/img103.png +0 -0
- data/doc/ep_flux/math-doc/document/img104.png +0 -0
- data/doc/ep_flux/math-doc/document/img105.png +0 -0
- data/doc/ep_flux/math-doc/document/img106.png +0 -0
- data/doc/ep_flux/math-doc/document/img107.png +0 -0
- data/doc/ep_flux/math-doc/document/img108.png +0 -0
- data/doc/ep_flux/math-doc/document/img109.png +0 -0
- data/doc/ep_flux/math-doc/document/img11.png +0 -0
- data/doc/ep_flux/math-doc/document/img110.png +0 -0
- data/doc/ep_flux/math-doc/document/img111.png +0 -0
- data/doc/ep_flux/math-doc/document/img112.png +0 -0
- data/doc/ep_flux/math-doc/document/img113.png +0 -0
- data/doc/ep_flux/math-doc/document/img114.png +0 -0
- data/doc/ep_flux/math-doc/document/img115.png +0 -0
- data/doc/ep_flux/math-doc/document/img116.png +0 -0
- data/doc/ep_flux/math-doc/document/img117.png +0 -0
- data/doc/ep_flux/math-doc/document/img118.png +0 -0
- data/doc/ep_flux/math-doc/document/img119.png +0 -0
- data/doc/ep_flux/math-doc/document/img12.png +0 -0
- data/doc/ep_flux/math-doc/document/img120.png +0 -0
- data/doc/ep_flux/math-doc/document/img121.png +0 -0
- data/doc/ep_flux/math-doc/document/img122.png +0 -0
- data/doc/ep_flux/math-doc/document/img123.png +0 -0
- data/doc/ep_flux/math-doc/document/img124.png +0 -0
- data/doc/ep_flux/math-doc/document/img125.png +0 -0
- data/doc/ep_flux/math-doc/document/img126.png +0 -0
- data/doc/ep_flux/math-doc/document/img127.png +0 -0
- data/doc/ep_flux/math-doc/document/img128.png +0 -0
- data/doc/ep_flux/math-doc/document/img129.png +0 -0
- data/doc/ep_flux/math-doc/document/img13.png +0 -0
- data/doc/ep_flux/math-doc/document/img130.png +0 -0
- data/doc/ep_flux/math-doc/document/img131.png +0 -0
- data/doc/ep_flux/math-doc/document/img132.png +0 -0
- data/doc/ep_flux/math-doc/document/img133.png +0 -0
- data/doc/ep_flux/math-doc/document/img134.png +0 -0
- data/doc/ep_flux/math-doc/document/img135.png +0 -0
- data/doc/ep_flux/math-doc/document/img136.png +0 -0
- data/doc/ep_flux/math-doc/document/img137.png +0 -0
- data/doc/ep_flux/math-doc/document/img138.png +0 -0
- data/doc/ep_flux/math-doc/document/img139.png +0 -0
- data/doc/ep_flux/math-doc/document/img14.png +0 -0
- data/doc/ep_flux/math-doc/document/img140.png +0 -0
- data/doc/ep_flux/math-doc/document/img141.png +0 -0
- data/doc/ep_flux/math-doc/document/img142.png +0 -0
- data/doc/ep_flux/math-doc/document/img143.png +0 -0
- data/doc/ep_flux/math-doc/document/img144.png +0 -0
- data/doc/ep_flux/math-doc/document/img145.png +0 -0
- data/doc/ep_flux/math-doc/document/img146.png +0 -0
- data/doc/ep_flux/math-doc/document/img147.png +0 -0
- data/doc/ep_flux/math-doc/document/img148.png +0 -0
- data/doc/ep_flux/math-doc/document/img149.png +0 -0
- data/doc/ep_flux/math-doc/document/img15.png +0 -0
- data/doc/ep_flux/math-doc/document/img150.png +0 -0
- data/doc/ep_flux/math-doc/document/img151.png +0 -0
- data/doc/ep_flux/math-doc/document/img152.png +0 -0
- data/doc/ep_flux/math-doc/document/img153.png +0 -0
- data/doc/ep_flux/math-doc/document/img154.png +0 -0
- data/doc/ep_flux/math-doc/document/img155.png +0 -0
- data/doc/ep_flux/math-doc/document/img156.png +0 -0
- data/doc/ep_flux/math-doc/document/img157.png +0 -0
- data/doc/ep_flux/math-doc/document/img158.png +0 -0
- data/doc/ep_flux/math-doc/document/img159.png +0 -0
- data/doc/ep_flux/math-doc/document/img16.png +0 -0
- data/doc/ep_flux/math-doc/document/img160.png +0 -0
- data/doc/ep_flux/math-doc/document/img161.png +0 -0
- data/doc/ep_flux/math-doc/document/img162.png +0 -0
- data/doc/ep_flux/math-doc/document/img163.png +0 -0
- data/doc/ep_flux/math-doc/document/img164.png +0 -0
- data/doc/ep_flux/math-doc/document/img165.png +0 -0
- data/doc/ep_flux/math-doc/document/img166.png +0 -0
- data/doc/ep_flux/math-doc/document/img167.png +0 -0
- data/doc/ep_flux/math-doc/document/img168.png +0 -0
- data/doc/ep_flux/math-doc/document/img169.png +0 -0
- data/doc/ep_flux/math-doc/document/img17.png +0 -0
- data/doc/ep_flux/math-doc/document/img170.png +0 -0
- data/doc/ep_flux/math-doc/document/img171.png +0 -0
- data/doc/ep_flux/math-doc/document/img172.png +0 -0
- data/doc/ep_flux/math-doc/document/img173.png +0 -0
- data/doc/ep_flux/math-doc/document/img174.png +0 -0
- data/doc/ep_flux/math-doc/document/img175.png +0 -0
- data/doc/ep_flux/math-doc/document/img176.png +0 -0
- data/doc/ep_flux/math-doc/document/img177.png +0 -0
- data/doc/ep_flux/math-doc/document/img178.png +0 -0
- data/doc/ep_flux/math-doc/document/img179.png +0 -0
- data/doc/ep_flux/math-doc/document/img18.png +0 -0
- data/doc/ep_flux/math-doc/document/img180.png +0 -0
- data/doc/ep_flux/math-doc/document/img181.png +0 -0
- data/doc/ep_flux/math-doc/document/img182.png +0 -0
- data/doc/ep_flux/math-doc/document/img183.png +0 -0
- data/doc/ep_flux/math-doc/document/img184.png +0 -0
- data/doc/ep_flux/math-doc/document/img185.png +0 -0
- data/doc/ep_flux/math-doc/document/img186.png +0 -0
- data/doc/ep_flux/math-doc/document/img187.png +0 -0
- data/doc/ep_flux/math-doc/document/img188.png +0 -0
- data/doc/ep_flux/math-doc/document/img189.png +0 -0
- data/doc/ep_flux/math-doc/document/img19.png +0 -0
- data/doc/ep_flux/math-doc/document/img190.png +0 -0
- data/doc/ep_flux/math-doc/document/img191.png +0 -0
- data/doc/ep_flux/math-doc/document/img192.png +0 -0
- data/doc/ep_flux/math-doc/document/img193.png +0 -0
- data/doc/ep_flux/math-doc/document/img194.png +0 -0
- data/doc/ep_flux/math-doc/document/img195.png +0 -0
- data/doc/ep_flux/math-doc/document/img196.png +0 -0
- data/doc/ep_flux/math-doc/document/img197.png +0 -0
- data/doc/ep_flux/math-doc/document/img198.png +0 -0
- data/doc/ep_flux/math-doc/document/img199.png +0 -0
- data/doc/ep_flux/math-doc/document/img2.png +0 -0
- data/doc/ep_flux/math-doc/document/img20.png +0 -0
- data/doc/ep_flux/math-doc/document/img200.png +0 -0
- data/doc/ep_flux/math-doc/document/img21.png +0 -0
- data/doc/ep_flux/math-doc/document/img22.png +0 -0
- data/doc/ep_flux/math-doc/document/img23.png +0 -0
- data/doc/ep_flux/math-doc/document/img24.png +0 -0
- data/doc/ep_flux/math-doc/document/img25.png +0 -0
- data/doc/ep_flux/math-doc/document/img26.png +0 -0
- data/doc/ep_flux/math-doc/document/img27.png +0 -0
- data/doc/ep_flux/math-doc/document/img28.png +0 -0
- data/doc/ep_flux/math-doc/document/img29.png +0 -0
- data/doc/ep_flux/math-doc/document/img3.png +0 -0
- data/doc/ep_flux/math-doc/document/img30.png +0 -0
- data/doc/ep_flux/math-doc/document/img31.png +0 -0
- data/doc/ep_flux/math-doc/document/img32.png +0 -0
- data/doc/ep_flux/math-doc/document/img33.png +0 -0
- data/doc/ep_flux/math-doc/document/img34.png +0 -0
- data/doc/ep_flux/math-doc/document/img35.png +0 -0
- data/doc/ep_flux/math-doc/document/img36.png +0 -0
- data/doc/ep_flux/math-doc/document/img37.png +0 -0
- data/doc/ep_flux/math-doc/document/img38.png +0 -0
- data/doc/ep_flux/math-doc/document/img39.png +0 -0
- data/doc/ep_flux/math-doc/document/img4.png +0 -0
- data/doc/ep_flux/math-doc/document/img40.png +0 -0
- data/doc/ep_flux/math-doc/document/img41.png +0 -0
- data/doc/ep_flux/math-doc/document/img42.png +0 -0
- data/doc/ep_flux/math-doc/document/img43.png +0 -0
- data/doc/ep_flux/math-doc/document/img44.png +0 -0
- data/doc/ep_flux/math-doc/document/img45.png +0 -0
- data/doc/ep_flux/math-doc/document/img46.png +0 -0
- data/doc/ep_flux/math-doc/document/img47.png +0 -0
- data/doc/ep_flux/math-doc/document/img48.png +0 -0
- data/doc/ep_flux/math-doc/document/img49.png +0 -0
- data/doc/ep_flux/math-doc/document/img5.png +0 -0
- data/doc/ep_flux/math-doc/document/img50.png +0 -0
- data/doc/ep_flux/math-doc/document/img51.png +0 -0
- data/doc/ep_flux/math-doc/document/img52.png +0 -0
- data/doc/ep_flux/math-doc/document/img53.png +0 -0
- data/doc/ep_flux/math-doc/document/img54.png +0 -0
- data/doc/ep_flux/math-doc/document/img55.png +0 -0
- data/doc/ep_flux/math-doc/document/img56.png +0 -0
- data/doc/ep_flux/math-doc/document/img57.png +0 -0
- data/doc/ep_flux/math-doc/document/img58.png +0 -0
- data/doc/ep_flux/math-doc/document/img59.png +0 -0
- data/doc/ep_flux/math-doc/document/img6.png +0 -0
- data/doc/ep_flux/math-doc/document/img60.png +0 -0
- data/doc/ep_flux/math-doc/document/img61.png +0 -0
- data/doc/ep_flux/math-doc/document/img62.png +0 -0
- data/doc/ep_flux/math-doc/document/img63.png +0 -0
- data/doc/ep_flux/math-doc/document/img64.png +0 -0
- data/doc/ep_flux/math-doc/document/img65.png +0 -0
- data/doc/ep_flux/math-doc/document/img66.png +0 -0
- data/doc/ep_flux/math-doc/document/img67.png +0 -0
- data/doc/ep_flux/math-doc/document/img68.png +0 -0
- data/doc/ep_flux/math-doc/document/img69.png +0 -0
- data/doc/ep_flux/math-doc/document/img7.png +0 -0
- data/doc/ep_flux/math-doc/document/img70.png +0 -0
- data/doc/ep_flux/math-doc/document/img71.png +0 -0
- data/doc/ep_flux/math-doc/document/img72.png +0 -0
- data/doc/ep_flux/math-doc/document/img73.png +0 -0
- data/doc/ep_flux/math-doc/document/img74.png +0 -0
- data/doc/ep_flux/math-doc/document/img75.png +0 -0
- data/doc/ep_flux/math-doc/document/img76.png +0 -0
- data/doc/ep_flux/math-doc/document/img77.png +0 -0
- data/doc/ep_flux/math-doc/document/img78.png +0 -0
- data/doc/ep_flux/math-doc/document/img79.png +0 -0
- data/doc/ep_flux/math-doc/document/img8.png +0 -0
- data/doc/ep_flux/math-doc/document/img80.png +0 -0
- data/doc/ep_flux/math-doc/document/img81.png +0 -0
- data/doc/ep_flux/math-doc/document/img82.png +0 -0
- data/doc/ep_flux/math-doc/document/img83.png +0 -0
- data/doc/ep_flux/math-doc/document/img84.png +0 -0
- data/doc/ep_flux/math-doc/document/img85.png +0 -0
- data/doc/ep_flux/math-doc/document/img86.png +0 -0
- data/doc/ep_flux/math-doc/document/img87.png +0 -0
- data/doc/ep_flux/math-doc/document/img88.png +0 -0
- data/doc/ep_flux/math-doc/document/img89.png +0 -0
- data/doc/ep_flux/math-doc/document/img9.png +0 -0
- data/doc/ep_flux/math-doc/document/img90.png +0 -0
- data/doc/ep_flux/math-doc/document/img91.png +0 -0
- data/doc/ep_flux/math-doc/document/img92.png +0 -0
- data/doc/ep_flux/math-doc/document/img93.png +0 -0
- data/doc/ep_flux/math-doc/document/img94.png +0 -0
- data/doc/ep_flux/math-doc/document/img95.png +0 -0
- data/doc/ep_flux/math-doc/document/img96.png +0 -0
- data/doc/ep_flux/math-doc/document/img97.png +0 -0
- data/doc/ep_flux/math-doc/document/img98.png +0 -0
- data/doc/ep_flux/math-doc/document/img99.png +0 -0
- data/doc/ep_flux/math-doc/document/index.html +101 -0
- data/doc/ep_flux/math-doc/document/internals.pl +258 -0
- data/doc/ep_flux/math-doc/document/labels.pl +265 -0
- data/doc/ep_flux/math-doc/document/next.png +0 -0
- data/doc/ep_flux/math-doc/document/next_g.png +0 -0
- data/doc/ep_flux/math-doc/document/node1.html +104 -0
- data/doc/ep_flux/math-doc/document/node10.html +164 -0
- data/doc/ep_flux/math-doc/document/node11.html +86 -0
- data/doc/ep_flux/math-doc/document/node12.html +166 -0
- data/doc/ep_flux/math-doc/document/node13.html +897 -0
- data/doc/ep_flux/math-doc/document/node14.html +1065 -0
- data/doc/ep_flux/math-doc/document/node15.html +72 -0
- data/doc/ep_flux/math-doc/document/node16.html +81 -0
- data/doc/ep_flux/math-doc/document/node2.html +82 -0
- data/doc/ep_flux/math-doc/document/node3.html +91 -0
- data/doc/ep_flux/math-doc/document/node4.html +149 -0
- data/doc/ep_flux/math-doc/document/node5.html +330 -0
- data/doc/ep_flux/math-doc/document/node6.html +99 -0
- data/doc/ep_flux/math-doc/document/node7.html +98 -0
- data/doc/ep_flux/math-doc/document/node8.html +83 -0
- data/doc/ep_flux/math-doc/document/node9.html +140 -0
- data/doc/ep_flux/math-doc/document/prev.png +0 -0
- data/doc/ep_flux/math-doc/document/prev_g.png +0 -0
- data/doc/ep_flux/math-doc/document/up.png +0 -0
- data/doc/ep_flux/math-doc/document/up_g.png +0 -0
- data/doc/gdir.html +412 -0
- data/doc/gdir_client.html +16 -0
- data/doc/gdir_connect_ftp-like.html +61 -0
- data/doc/gdir_server.html +45 -0
- data/doc/ggraph.html +1119 -0
- data/doc/gpcat.html +45 -0
- data/doc/gpcut.html +47 -0
- data/doc/gphys.html +624 -0
- data/doc/gphys_fft.html +324 -0
- data/doc/gphys_grads_io.html +69 -0
- data/doc/gphys_grib_io.html +82 -0
- data/doc/gphys_io.html +183 -0
- data/doc/gphys_io_common.html +18 -0
- data/doc/gphys_netcdf_io.html +283 -0
- data/doc/gplist.html +24 -0
- data/doc/gpmath.html +52 -0
- data/doc/gpmaxmin.html +32 -0
- data/doc/gpprint.html +35 -0
- data/doc/gpview.html +349 -0
- data/doc/grads2nc_with_gphys.html +21 -0
- data/doc/grads_gridded.html +307 -0
- data/doc/grib.html +149 -0
- data/doc/grid.html +224 -0
- data/doc/index.html +145 -0
- data/doc/index.rd +138 -0
- data/doc/netcdf_convention.html +136 -0
- data/doc/unumeric.html +176 -0
- data/doc/update +69 -0
- data/doc/update_rdoc +8 -0
- data/doc/varray.html +299 -0
- data/doc/varraycomposite.html +67 -0
- data/ext_init.c +1 -0
- data/extconf.rb +16 -6
- data/gphys.gemspec +33 -26
- data/interpo.c +1 -1
- data/lib/numru/dclext.rb +718 -546
- data/lib/numru/derivative.rb +2 -0
- data/lib/numru/ganalysis.rb +38 -0
- data/lib/numru/ganalysis/beta_plane.rb +103 -0
- data/lib/numru/ganalysis/eof.rb +3 -2
- data/lib/numru/ganalysis/fitting.rb +559 -0
- data/lib/numru/ganalysis/histogram.rb +36 -19
- data/lib/numru/ganalysis/log_p.rb +130 -0
- data/lib/numru/ganalysis/met.rb +396 -2
- data/lib/numru/ganalysis/met_z.rb +300 -0
- data/lib/numru/ganalysis/planet.rb +17 -7
- data/lib/numru/ganalysis/qg.rb +685 -0
- data/lib/numru/ganalysis/sigma_coord.rb +90 -0
- data/lib/numru/gdir.rb +2 -1
- data/lib/numru/ggraph.rb +204 -60
- data/lib/numru/ggraph_on_merdional_section.rb +1 -1
- data/lib/numru/gphys.rb +6 -0
- data/lib/numru/gphys/assoccoords.rb +18 -3
- data/lib/numru/gphys/axis.rb +209 -8
- data/lib/numru/gphys/derivative.rb +11 -0
- data/lib/numru/gphys/gphys.rb +539 -48
- data/lib/numru/gphys/gphys_dim_op.rb +331 -0
- data/lib/numru/gphys/gphys_fft.rb +48 -2
- data/lib/numru/gphys/gphys_io.rb +241 -13
- data/lib/numru/gphys/gphys_netcdf_io.rb +77 -39
- data/lib/numru/gphys/gphys_nusdas_io.rb +3 -0
- data/lib/numru/gphys/grib.rb +133 -54
- data/lib/numru/gphys/grib_params.rb +26 -3
- data/lib/numru/gphys/grid.rb +75 -34
- data/lib/numru/gphys/interpolate.rb +24 -10
- data/lib/numru/gphys/mdstorage.rb +160 -0
- data/lib/numru/gphys/netcdf_convention.rb +4 -2
- data/lib/numru/gphys/subsetmapping.rb +0 -1
- data/lib/numru/gphys/unumeric.rb +50 -5
- data/lib/numru/gphys/varray.rb +15 -30
- data/lib/numru/gphys/varraycomposite.rb +107 -24
- data/lib/numru/gphys/varraynetcdf.rb +9 -3
- data/lib/numru/gphys/version.rb +5 -0
- data/sample/druby_cli1.rb +2 -0
- data/sample/druby_cli2.rb +0 -6
- data/sample/druby_serv2.rb +0 -13
- data/spec/gphys_spec.rb +11 -0
- data/spec/spec_helper.rb +2 -0
- data/test/test_assoccoords.rb +102 -0
- data/test/test_axis.rb +61 -0
- data/test/test_fitting.rb +116 -0
- data/test/test_gphys.rb +20 -0
- data/test/test_met_z.rb +96 -0
- data/test/test_sigma_coord.rb +50 -0
- data/{test → test_old}/eof_slp.rb +0 -0
- data/{test → test_old}/mltbit.dat +0 -0
- data/{test → test_old}/test_ep_flux.rb +0 -0
- data/{test → test_old}/test_multibitIO.rb +0 -0
- metadata +530 -191
- data/README.md +0 -29
- data/lib/gphys.rb +0 -2
- data/lib/numru/dclext_datetime_ax.rb +0 -220
- data/lib/version.rb +0 -3
|
Binary file
|
|
Binary file
|
|
Binary file
|
|
Binary file
|
|
@@ -0,0 +1,158 @@
|
|
|
1
|
+
<?xml version="1.0" ?>
|
|
2
|
+
<!DOCTYPE html
|
|
3
|
+
PUBLIC "-//W3C//DTD XHTML 1.0 Transitional//EN"
|
|
4
|
+
"http://www.w3.org/TR/xhtml1/DTD/xhtml1-transitional.dtd">
|
|
5
|
+
<html xmlns="http://www.w3.org/1999/xhtml">
|
|
6
|
+
<head>
|
|
7
|
+
<title>../../lib/numru/derivative.rb</title>
|
|
8
|
+
</head>
|
|
9
|
+
<body>
|
|
10
|
+
<h1><a name="label-0" id="label-0">module NumRu::Derivative in derivative.rb</a></h1><!-- RDLabel: "module NumRu::Derivative in derivative.rb" -->
|
|
11
|
+
<h2><a name="label-1" id="label-1">todo</a></h2><!-- RDLabel: "todo" -->
|
|
12
|
+
<ul>
|
|
13
|
+
<li>decide argument of b_expand_linear_ext is Symbol or Numeric.
|
|
14
|
+
<ul>
|
|
15
|
+
<li>now is Numeric.</li>
|
|
16
|
+
<li>it denpends the treatment of dRuby.</li>
|
|
17
|
+
</ul></li>
|
|
18
|
+
<li>support other boundary conditions.</li>
|
|
19
|
+
</ul>
|
|
20
|
+
<h2><a name="label-2" id="label-2">Index</a></h2><!-- RDLabel: "Index" -->
|
|
21
|
+
<ul>
|
|
22
|
+
<li><a href="#label-3">module NumRu::Derivative</a>
|
|
23
|
+
<ul>
|
|
24
|
+
<li><a href="#label-4">threepoint_O2nd_deriv</a>
|
|
25
|
+
<ul>
|
|
26
|
+
<li>First derivative (2nd Order difference use three point.)</li>
|
|
27
|
+
</ul></li>
|
|
28
|
+
<li><a href="#label-5">cderiv</a>
|
|
29
|
+
<ul>
|
|
30
|
+
<li>First derivative (center difference use two point.)</li>
|
|
31
|
+
</ul></li>
|
|
32
|
+
<li><!-- Reference, RDLabel "deriv2d" doesn't exist --><em class="label-not-found">deriv2d</em><!-- Reference end -->
|
|
33
|
+
<ul>
|
|
34
|
+
<li>2nd derivative</li>
|
|
35
|
+
</ul></li>
|
|
36
|
+
<li><a href="#label-7">b_expand_linear_ext</a>
|
|
37
|
+
<ul>
|
|
38
|
+
<li>return array extended boundaries with linear extention.</li>
|
|
39
|
+
</ul></li>
|
|
40
|
+
<li><a href="#label-8">cdiff</a>
|
|
41
|
+
<ul>
|
|
42
|
+
<li>return difference. (center difference)</li>
|
|
43
|
+
</ul></li>
|
|
44
|
+
</ul></li>
|
|
45
|
+
</ul>
|
|
46
|
+
<h1><a name="label-3" id="label-3">module NumRu::Derivative</a></h1><!-- RDLabel: "module NumRu::Derivative" -->
|
|
47
|
+
<p>Module functions of Derivative Operater for NArray.</p>
|
|
48
|
+
<dl>
|
|
49
|
+
<dt><h4><a name="label-4" id="label-4"><code>threepoint_O2nd_deriv(<var>z</var>, <var>x</var>, <var>dim</var>, <var>bc</var>=<var>LINEAR_EXT</var>)</code></a></h4></dt><!-- RDLabel: "threepoint_O2nd_deriv" -->
|
|
50
|
+
<dd>
|
|
51
|
+
<p>Derivate of <code>z</code> with respect to <code>dim</code> th dim using a 2nd
|
|
52
|
+
order 3-point differentiation valid for non-uniform grid: </p>
|
|
53
|
+
<pre>(s**2*z_{i+1} + (t**2 - s**2)*f_{i} - t**2*f_{i-1}) / (s*t*(s + t))</pre>
|
|
54
|
+
<p>Here, s represents (x_{i} - x_{i-1}) ,t represents (x_{i+1} - x_{i})
|
|
55
|
+
and _{i} represents the suffix of {i} th element in the <!-- Reference, RDLabel "dim" doesn't exist --><em class="label-not-found">dim</em><!-- Reference end --> th
|
|
56
|
+
dimension of the array. ).</p>
|
|
57
|
+
<p>ARGUMENTS</p>
|
|
58
|
+
<ul>
|
|
59
|
+
<li>z (NArray): a NArray which you want to derivative.</li>
|
|
60
|
+
<li>x (NArray): a NArray represents the dimension which derivative respect to.
|
|
61
|
+
z.rank must be 1.</li>
|
|
62
|
+
<li>dim (Numeric): a Numeric represents the dimention which derivative respect to.
|
|
63
|
+
you can give number count backward (<!-- Reference, RDLabel "dim" doesn't exist --><em class="label-not-found">dim</em><!-- Reference end --><0), but <!-- Reference, RDLabel "z.rank ��dim" doesn't exist --><em class="label-not-found">z.rank ��dim</em><!-- Reference end --> must be > 0. </li>
|
|
64
|
+
<li>bc (Numeric) : a Numeric to represent boundary condition.
|
|
65
|
+
Supported conditions are expressed in terms of boundary extension
|
|
66
|
+
applied before differentiation:
|
|
67
|
+
<ul>
|
|
68
|
+
<li>LINEAR_EXT(=1): default value. linear extension using two
|
|
69
|
+
grid points at the boundaries.</li>
|
|
70
|
+
<li>CYCLIC_EXT(=2): cyclic extension as k,0,1,..,k,0 for indices 0,1,..,k</li>
|
|
71
|
+
<li>MIRROR_A(=3): for indices 0,1,..,k, extends as 0,0,1,..,k,k
|
|
72
|
+
(good for cell interior values with the Neumann condition)</li>
|
|
73
|
+
<li>MIRROR_B(=4): for indices 0,1,..,k, extends as 1,0,1,..,k,k-1
|
|
74
|
+
(good for cell boundary values with the Neumann condition)</li>
|
|
75
|
+
</ul></li>
|
|
76
|
+
</ul>
|
|
77
|
+
<p>RETURN VALUE</p>
|
|
78
|
+
<ul>
|
|
79
|
+
<li>O2nd_deriv_data (NArray): (s**2*z_{i+1} + (t**2 - s**2)*f_{i} - t**2*f_{i-1}) / (s*t*(s + t))</li>
|
|
80
|
+
</ul></dd>
|
|
81
|
+
<dt><h4><a name="label-5" id="label-5"><code>cderiv(<var>z</var>, <var>x</var>, <var>dim</var>, <var>bc</var>=<var>LINEAR_EXT</var>)</code></a></h4></dt><!-- RDLabel: "cderiv" -->
|
|
82
|
+
<dd>
|
|
83
|
+
<p>Derivate of <code>z</code> with respect to <code>dim</code> th dim using centeral
|
|
84
|
+
differenciation: (z_{i+1} - z_{i-1}) / (x_{i+1} - x_{i-1})</p>
|
|
85
|
+
<p>ARGUMENTS</p>
|
|
86
|
+
<ul>
|
|
87
|
+
<li>z (NArray): a NArray which you want to derivative.</li>
|
|
88
|
+
<li>x (NArray): a NArray represents the dimension which derivative respect
|
|
89
|
+
to. z.rank must be 1.</li>
|
|
90
|
+
<li>dim (Numeric): a Numeric represents the dimention which derivative
|
|
91
|
+
respect to. you can give number count backward (<!-- Reference, RDLabel "dim" doesn't exist --><em class="label-not-found">dim</em><!-- Reference end --><0), but
|
|
92
|
+
<!-- Reference, RDLabel "z.rank ��dim" doesn't exist --><em class="label-not-found">z.rank ��dim</em><!-- Reference end --> must be > 0. </li>
|
|
93
|
+
<li>bc (Numeric) : a Numeric to represent boundary condition.
|
|
94
|
+
See <a href="#label-4">threepoint_O2nd_deriv</a> for supported conditions.</li>
|
|
95
|
+
</ul>
|
|
96
|
+
<p>RETURN VALUE</p>
|
|
97
|
+
<ul>
|
|
98
|
+
<li>cderiv_data (NArray): (z_{i+1} - z_{i-1}) / (x_{i+1} - x_{i-1})</li>
|
|
99
|
+
</ul></dd>
|
|
100
|
+
<dt><h4><a name="label-6" id="label-6"><code>deriv2nd(<var>z</var>, <var>x</var>, <var>dim</var>, <var>bc</var>=<var>LINEAR_EXT</var>)</code></a></h4></dt><!-- RDLabel: "deriv2nd" -->
|
|
101
|
+
<dd>
|
|
102
|
+
<p>2nd Derivate of <code>z</code> with respect to <code>dim</code>-th dim
|
|
103
|
+
covering non-uniform grids. Based on:</p>
|
|
104
|
+
<pre>( (z_{i+1}-z_{i})/(x_{i+1}-x_{i}) - (z_{i}-z_{i-1})/(x_{i}-x_{i-1}) )
|
|
105
|
+
/ ((x_{i+1}-x_{i-1})/2)</pre>
|
|
106
|
+
<p>ARGUMENTS</p>
|
|
107
|
+
<ul>
|
|
108
|
+
<li>z (NArray): a NArray which you want to derivative.</li>
|
|
109
|
+
<li>x (NArray): a NArray represents the dimension which derivative respect
|
|
110
|
+
to. z.rank must be 1.</li>
|
|
111
|
+
<li>dim (Numeric): a Numeric represents the dimention which derivative
|
|
112
|
+
respect to. you can give number count backward (<!-- Reference, RDLabel "dim" doesn't exist --><em class="label-not-found">dim</em><!-- Reference end --><0), but
|
|
113
|
+
<!-- Reference, RDLabel "z.rank ��dim" doesn't exist --><em class="label-not-found">z.rank ��dim</em><!-- Reference end --> must be > 0. </li>
|
|
114
|
+
<li>bc (Numeric) : a Numeric to represent boundary condition.
|
|
115
|
+
See <a href="#label-4">threepoint_O2nd_deriv</a> for supported conditions.</li>
|
|
116
|
+
</ul>
|
|
117
|
+
<p>RETURN VALUE</p>
|
|
118
|
+
<ul>
|
|
119
|
+
<li>cderiv_data (NArray): (z_{i+1} - z_{i-1}) / (x_{i+1} - x_{i-1})</li>
|
|
120
|
+
</ul></dd>
|
|
121
|
+
<dt><h4><a name="label-7" id="label-7"><code>b_expand_linear_ext(<var>z</var>, <var>dim</var>)</code></a></h4></dt><!-- RDLabel: "b_expand_linear_ext" -->
|
|
122
|
+
<dd>
|
|
123
|
+
<p>expand boundary with linear value. extend array with 1 grid at each
|
|
124
|
+
boundary with <!-- Reference, RDLabel "dim" doesn't exist --><em class="label-not-found">dim</em><!-- Reference end --> th dimension, and assign th value which diffrential
|
|
125
|
+
value between a grid short of boundary and boundary grid in original array.
|
|
126
|
+
(on other wards, 2*z_{0}-z_{1} or 2*z_{n-1}-z_{n-2}: now _{i} represents the </p>
|
|
127
|
+
<pre>suffix of {i} th element in the ((<dim>)) th dimension of array. ).</pre>
|
|
128
|
+
<p>ARGUMENTS</p>
|
|
129
|
+
<ul>
|
|
130
|
+
<li>z (NArray): a NArray which you want to expand boundary.</li>
|
|
131
|
+
<li>dim (Numeric): a Numeric represents the dimention which derivative
|
|
132
|
+
respect to. you can give number count backward (<!-- Reference, RDLabel "dim" doesn't exist --><em class="label-not-found">dim</em><!-- Reference end --><0), but
|
|
133
|
+
<!-- Reference, RDLabel "z.rank ��dim" doesn't exist --><em class="label-not-found">z.rank ��dim</em><!-- Reference end --> must be > 0. </li>
|
|
134
|
+
</ul>
|
|
135
|
+
<p>RETURN VALUE</p>
|
|
136
|
+
<ul>
|
|
137
|
+
<li>expand_data (NArray): </li>
|
|
138
|
+
</ul></dd>
|
|
139
|
+
<dt><h4><a name="label-8" id="label-8"><code>cdiff(<var>x</var>, <var>dim</var>)</code></a></h4></dt><!-- RDLabel: "cdiff" -->
|
|
140
|
+
<dd>
|
|
141
|
+
<p>Diffrence operater. return an NArray which a difference <!-- Reference, RDLabel "x" doesn't exist --><em class="label-not-found">x</em><!-- Reference end -->
|
|
142
|
+
( in other wards, (x_{i+1} - x_{i-1}): now _{i} represents the suffix of
|
|
143
|
+
{i} th element in the <!-- Reference, RDLabel "dim" doesn't exist --><em class="label-not-found">dim</em><!-- Reference end --> th dimension of array. ).</p>
|
|
144
|
+
<p>ARGUMENTS</p>
|
|
145
|
+
<ul>
|
|
146
|
+
<li>x (NArray): a NArray which you want to get difference.</li>
|
|
147
|
+
<li>dim (Numeric): a Numeric representing the dimention which derivative
|
|
148
|
+
respect to. you can give number count backward (<!-- Reference, RDLabel "dim" doesn't exist --><em class="label-not-found">dim</em><!-- Reference end --><0), but
|
|
149
|
+
<!-- Reference, RDLabel "z.rank ��dim" doesn't exist --><em class="label-not-found">z.rank ��dim</em><!-- Reference end --> must be > 0. </li>
|
|
150
|
+
</ul>
|
|
151
|
+
<p>RETURN VALUE</p>
|
|
152
|
+
<ul>
|
|
153
|
+
<li>cdiff_data (NArray): (x_{i+1} - x_{i-1})</li>
|
|
154
|
+
</ul></dd>
|
|
155
|
+
</dl>
|
|
156
|
+
|
|
157
|
+
</body>
|
|
158
|
+
</html>
|
|
@@ -0,0 +1,469 @@
|
|
|
1
|
+
<?xml version="1.0" ?>
|
|
2
|
+
<!DOCTYPE html
|
|
3
|
+
PUBLIC "-//W3C//DTD XHTML 1.0 Transitional//EN"
|
|
4
|
+
"http://www.w3.org/TR/xhtml1/DTD/xhtml1-transitional.dtd">
|
|
5
|
+
<html xmlns="http://www.w3.org/1999/xhtml">
|
|
6
|
+
<head>
|
|
7
|
+
<title>../../lib/numru/gphys/ep_flux.rb</title>
|
|
8
|
+
</head>
|
|
9
|
+
<body>
|
|
10
|
+
<h1><a name="label:0" id="label:0">module NumRu::GPhys::EP_Flux in ep_flux.</a></h1><!-- RDLabel: "module NumRu::GPhys::EP_Flux in ep_flux." -->
|
|
11
|
+
<h2><a name="label:1" id="label:1">Testprogram</a></h2><!-- RDLabel: "Testprogram" -->
|
|
12
|
+
<p>Test script path is 'test/test_ep_flux.rb' in expand dir of gphys.tar.gz.</p>
|
|
13
|
+
<h2><a name="label:2" id="label:2">Index</a></h2><!-- RDLabel: "Index" -->
|
|
14
|
+
<ul>
|
|
15
|
+
<li><a href="#label:3">module NumRu::GPhys::EP_Flux</a>
|
|
16
|
+
<ul>
|
|
17
|
+
<li><a href="#label:4">Functions:</a>
|
|
18
|
+
<ul>
|
|
19
|
+
<li><a href="#label:5">ep_full_sphere</a>
|
|
20
|
+
<ul>
|
|
21
|
+
<li>calculate EP Flux with full set equations on spherical coordinate.</li>
|
|
22
|
+
</ul></li>
|
|
23
|
+
<li><a href="#label:6">div_sphere</a>
|
|
24
|
+
<ul>
|
|
25
|
+
<li>calculate divergence on spherical coordinate, Not only for EP Flux!</li>
|
|
26
|
+
</ul></li>
|
|
27
|
+
<li><a href="#label:34">strm_rmean</a>
|
|
28
|
+
<ul>
|
|
29
|
+
<li>calculate residual mass stream function in spherical coordinate
|
|
30
|
+
(it might be good precision, but enough test has not done yet.)</li>
|
|
31
|
+
</ul></li>
|
|
32
|
+
<li><a href="#label:9">scale_height</a>
|
|
33
|
+
<ul>
|
|
34
|
+
<li>get the scale height.</li>
|
|
35
|
+
</ul></li>
|
|
36
|
+
<li><a href="#label:10">scale_height=</a>
|
|
37
|
+
<ul>
|
|
38
|
+
<li>set the scale height.</li>
|
|
39
|
+
</ul></li>
|
|
40
|
+
<li><a href="#label:11">radius</a>
|
|
41
|
+
<ul>
|
|
42
|
+
<li>get the radus of planet.</li>
|
|
43
|
+
</ul></li>
|
|
44
|
+
<li><a href="#label:12">radius=</a>
|
|
45
|
+
<ul>
|
|
46
|
+
<li>set the radus of planet.</li>
|
|
47
|
+
</ul></li>
|
|
48
|
+
<li><a href="#label:13">rot_period</a>
|
|
49
|
+
<ul>
|
|
50
|
+
<li>get the rotation period of planet.</li>
|
|
51
|
+
</ul></li>
|
|
52
|
+
<li><a href="#label:14">rot_period=</a>
|
|
53
|
+
<ul>
|
|
54
|
+
<li>set the rotation period of planet.</li>
|
|
55
|
+
</ul></li>
|
|
56
|
+
<li><a href="#label:15">g_forces</a>
|
|
57
|
+
<ul>
|
|
58
|
+
<li>get the gravitational acceleration in surface.</li>
|
|
59
|
+
</ul></li>
|
|
60
|
+
<li><a href="#label:16">g_forces=</a>
|
|
61
|
+
<ul>
|
|
62
|
+
<li>set the gravitational acceleration in surface.</li>
|
|
63
|
+
</ul></li>
|
|
64
|
+
<li><a href="#label:17">p00</a>
|
|
65
|
+
<ul>
|
|
66
|
+
<li>get the reference surface pressure.</li>
|
|
67
|
+
</ul></li>
|
|
68
|
+
<li><a href="#label:18">p00=</a>
|
|
69
|
+
<ul>
|
|
70
|
+
<li>set the reference surface pressure.</li>
|
|
71
|
+
</ul></li>
|
|
72
|
+
<li><a href="#label:19">cp</a>
|
|
73
|
+
<ul>
|
|
74
|
+
<li>get the the specific heat at constant pressure of the atmosphere.</li>
|
|
75
|
+
</ul></li>
|
|
76
|
+
<li><a href="#label:20">cp=</a>
|
|
77
|
+
<ul>
|
|
78
|
+
<li>set the the specific heat at constant pressure of the atmosphere.</li>
|
|
79
|
+
</ul></li>
|
|
80
|
+
<li><a href="#label:21">gas_const</a>
|
|
81
|
+
<ul>
|
|
82
|
+
<li>get the specific heat at constant pressure of the atmosphere.</li>
|
|
83
|
+
</ul></li>
|
|
84
|
+
<li><a href="#label:22">gas_const=</a>
|
|
85
|
+
<ul>
|
|
86
|
+
<li>set the specific heat at constant pressure of the atmosphere.</li>
|
|
87
|
+
</ul></li>
|
|
88
|
+
<li><a href="#label:23">get_constants</a>
|
|
89
|
+
<ul>
|
|
90
|
+
<li>get the module variables</li>
|
|
91
|
+
</ul></li>
|
|
92
|
+
<li><a href="#label:24">set_constants</a>
|
|
93
|
+
<ul>
|
|
94
|
+
<li>set the module variables</li>
|
|
95
|
+
</ul></li>
|
|
96
|
+
<li><a href="#label:25">make_gphys</a>
|
|
97
|
+
<ul>
|
|
98
|
+
<li>convert <!-- Reference, RDLabel "Axis" doesn't exist --><em class="label-not-found">Axis</em><!-- Reference end --> to <!-- Reference, RDLabel "GPhys" doesn't exist --><em class="label-not-found">GPhys</em><!-- Reference end -->.</li>
|
|
99
|
+
</ul></li>
|
|
100
|
+
<li><a href="#label:26">to_w_if_omega</a>
|
|
101
|
+
<ul>
|
|
102
|
+
<li>convert to velocity if <!-- Reference, RDLabel "gp" doesn't exist --><em class="label-not-found">gp</em><!-- Reference end --> is pressure velocity.</li>
|
|
103
|
+
</ul></li>
|
|
104
|
+
<li><a href="#label:27">to_z_if_pressure</a>
|
|
105
|
+
<ul>
|
|
106
|
+
<li>convert to altitude if <!-- Reference, RDLabel "gp" doesn't exist --><em class="label-not-found">gp</em><!-- Reference end --> is pressure.</li>
|
|
107
|
+
</ul></li>
|
|
108
|
+
<li><a href="#label:28">to_p_if_altitude</a>
|
|
109
|
+
<ul>
|
|
110
|
+
<li>convert to pressure if <!-- Reference, RDLabel "gp" doesn't exist --><em class="label-not-found">gp</em><!-- Reference end --> is altitude.</li>
|
|
111
|
+
</ul></li>
|
|
112
|
+
<li><a href="#label:29">to_theta_if_temperature</a>
|
|
113
|
+
<ul>
|
|
114
|
+
<li>convert to potential temperature of temperature if <!-- Reference, RDLabel "flag" doesn't exist --><em class="label-not-found">flag</em><!-- Reference end --> is true.</li>
|
|
115
|
+
</ul></li>
|
|
116
|
+
<li><a href="#label:30">to_rad_if_deg</a>
|
|
117
|
+
<ul>
|
|
118
|
+
<li>convert to radian if <!-- Reference, RDLabel "gp" doesn't exist --><em class="label-not-found">gp</em><!-- Reference end --> is degrees.</li>
|
|
119
|
+
</ul></li>
|
|
120
|
+
<li><a href="#label:31">eddy_products</a>
|
|
121
|
+
<ul>
|
|
122
|
+
<li>calculate eddy flux respect to <!-- Reference, RDLabel "dimname" doesn't exist --><em class="label-not-found">dimname</em><!-- Reference end -->.</li>
|
|
123
|
+
</ul></li>
|
|
124
|
+
<li><a href="#label:32">remove_0_at_poles</a>
|
|
125
|
+
<ul>
|
|
126
|
+
<li>set value if the cos(phi) is 0 at poles (phi is latitude).</li>
|
|
127
|
+
</ul></li>
|
|
128
|
+
<li><a href="#label:33">preparate_for_vector_on_merdional_section</a>
|
|
129
|
+
<ul>
|
|
130
|
+
<li>preparate for <!-- Reference, RDLabel "GGraph::vector_on_merdional_section" doesn't exist --><em class="label-not-found">GGraph::vector_on_merdional_section</em><!-- Reference end -->
|
|
131
|
+
in vector_on_merdional_section.rb</li>
|
|
132
|
+
</ul></li>
|
|
133
|
+
</ul></li>
|
|
134
|
+
<li><a href="#label:35">Constants:</a>
|
|
135
|
+
<ul>
|
|
136
|
+
<li><a href="#label:36">Deriv_methods</a>
|
|
137
|
+
<ul>
|
|
138
|
+
<li>derivative method names.</li>
|
|
139
|
+
</ul></li>
|
|
140
|
+
</ul></li>
|
|
141
|
+
</ul></li>
|
|
142
|
+
</ul>
|
|
143
|
+
<h1><a name="label:3" id="label:3">module NumRu::GPhys::EP_Flux</a></h1><!-- RDLabel: "module NumRu::GPhys::EP_Flux" -->
|
|
144
|
+
<p>Module functions of EP_Flux Operater for GPhys.</p>
|
|
145
|
+
<h2><a name="label:4" id="label:4">Functions:</a></h2><!-- RDLabel: "Functions:" -->
|
|
146
|
+
<dl>
|
|
147
|
+
<dt><a name="label:5" id="label:5"><code>ep_full_sphere(<var>gp_u</var>, <var>gp_v</var>, <var>gp_w_or_omega</var>, <var>gp_temp_or_theta</var>, <var>flag_temp_or_theta</var>, <var>xyzdims</var>=[<var>0</var>,<var>1</var>,<var>2</var>])</code></a></dt><!-- RDLabel: "ep_full_sphere" -->
|
|
148
|
+
<dd>
|
|
149
|
+
<p>Calculate Eliassen-Palm Flux(EP-Flux) from full set equations on the
|
|
150
|
+
spherical coordinate. this method calculates EP-Flux from 4 GPhys objects,
|
|
151
|
+
zonal-wind velocity(U), merdional-wind velocity(V), vertical-wind velocity(W)
|
|
152
|
+
or pressure velocity(Omega), and temperature(T) or potential(Theta)
|
|
153
|
+
temperature. check the equations on documents. </p>
|
|
154
|
+
<p>Furthermore, Residual mean merdional circulation (0, v*, w*) can be calculated.</p>
|
|
155
|
+
<p>ARGUMENTS</p>
|
|
156
|
+
<ul>
|
|
157
|
+
<li>gp_u (GPhys): a GPhys which data is U.</li>
|
|
158
|
+
<li>gp_v (GPhys): a GPhys which data is V.</li>
|
|
159
|
+
<li>gp_w_or_omega (GPhys): a GPhys which data is W or Omega. if you give
|
|
160
|
+
gp_omega, convert to W in this method and calculate EP-Flux.</li>
|
|
161
|
+
<li>gp_temp_or_theta (GPhys): a GPhys which data is T or Theta.</li>
|
|
162
|
+
<li><p>xyzdims (Array): an Array which represents location of upper gphyses's </p>
|
|
163
|
+
<pre>coordinate. if coordinate configuration is
|
|
164
|
+
(longitude, latitude, z), then xyzdims = [0, 1, 2].
|
|
165
|
+
else if coordinate configuration is
|
|
166
|
+
(z, latitude, longitude), then xyzdims = [2, 1, 0].</pre></li>
|
|
167
|
+
</ul>
|
|
168
|
+
<p>RETURN VALUE</p>
|
|
169
|
+
<ul>
|
|
170
|
+
<li>epflx_y (GPhys): EP-Flux y-component. it is on the merdional hoge.</li>
|
|
171
|
+
<li>epflx_z (GPhys): EP-Flux z-component.</li>
|
|
172
|
+
<li>v_rmean (GPhys): residual zonal mean V.</li>
|
|
173
|
+
<li>w_rmean (GPhys): residual zonal mean W.</li>
|
|
174
|
+
<li>gp_lat (GPhys): latitude (its units is radian)</li>
|
|
175
|
+
<li>gp_z (GPhys): from vertical axis (z)</li>
|
|
176
|
+
<li>u_mean (GPhys): zonal-mean U.</li>
|
|
177
|
+
<li>theta_mean (GPhys): zonal-mean Theta.</li>
|
|
178
|
+
<li>uv_dash (GPhys): zonal mean of zonal-eddy products U and V.</li>
|
|
179
|
+
<li>vt_dash (GPhys): zonal mean of zonal-eddy products V and Theta.</li>
|
|
180
|
+
<li>uw_dash (GPhys): zonal mean of zonal-eddy products U and W.</li>
|
|
181
|
+
<li>dtheta_dz (GPhys): zonal mean Theta derivate with z.</li>
|
|
182
|
+
</ul></dd>
|
|
183
|
+
<dt><a name="label:6" id="label:6"><code>div_sphere(<var>gp_y</var>, <var>gp_z</var>)</code></a></dt><!-- RDLabel: "div_sphere" -->
|
|
184
|
+
<dd>
|
|
185
|
+
<p>Calculate divergence on the spherical coordinate. it is exclusive to in
|
|
186
|
+
merdional cross section.</p>
|
|
187
|
+
<p>ARGUMENTS</p>
|
|
188
|
+
<ul>
|
|
189
|
+
<li><p>gp_y (GPhys): a GPhys which is merdional component you want to calculate</p>
|
|
190
|
+
<pre>divergence.</pre></li>
|
|
191
|
+
<li><p>gp_z (GPhys): a GPhys which is vertical component you want to calculate </p>
|
|
192
|
+
<pre>divergence.</pre></li>
|
|
193
|
+
</ul>
|
|
194
|
+
<p>RETURN VALUE</p>
|
|
195
|
+
<ul>
|
|
196
|
+
<li>gp_div (GPhys): a GPhys which is divergence on the spherical coordinate.</li>
|
|
197
|
+
</ul></dd>
|
|
198
|
+
<dt><a name="label:7" id="label:7"><code>set_deriv_method( <var>method_name</var> )</code></a></dt><!-- RDLabel: "set_deriv_method" -->
|
|
199
|
+
<dd>
|
|
200
|
+
<p>Set derivative method. methods are defined in <!-- Reference, RDLabel "GPhys::derivative" doesn't exist --><em class="label-not-found">GPhys::derivative</em><!-- Reference end -->.
|
|
201
|
+
Now <!-- Reference, RDLabel "cderiv" doesn't exist --><em class="label-not-found">cderiv</em><!-- Reference end --> and <!-- Reference, RDLabel "threepoint_O2nd_deriv" doesn't exist --><em class="label-not-found">threepoint_O2nd_deriv</em><!-- Reference end --> supported</p>
|
|
202
|
+
<p>ARGUMENTS</p>
|
|
203
|
+
<ul>
|
|
204
|
+
<li>method_name (String): derivative method name.</li>
|
|
205
|
+
</ul>
|
|
206
|
+
<p>RETURN VALUE</p>
|
|
207
|
+
<ul>
|
|
208
|
+
<li>nil</li>
|
|
209
|
+
</ul></dd>
|
|
210
|
+
<dt><a name="label:8" id="label:8"><code>deriv( *<var>args</var> )</code></a></dt><!-- RDLabel: "deriv" -->
|
|
211
|
+
<dd>
|
|
212
|
+
<p>Call derivative method defined in <!-- Reference, RDLabel "GPhys::Derivative" doesn't exist --><em class="label-not-found">GPhys::Derivative</em><!-- Reference end --> refer to
|
|
213
|
+
<!-- Reference, RDLabel "@@deriv_method" doesn't exist --><em class="label-not-found">@@deriv_method</em><!-- Reference end -->. </p>
|
|
214
|
+
<p>ARGUMENTS</p>
|
|
215
|
+
<ul>
|
|
216
|
+
<li>args : Option for derivative method. Pleaase see <!-- Reference, RDLabel "GPhys::Derivative" doesn't exist --><em class="label-not-found">GPhys::Derivative</em><!-- Reference end -->.</li>
|
|
217
|
+
</ul>
|
|
218
|
+
<p>RETURN VALUE</p>
|
|
219
|
+
<ul>
|
|
220
|
+
<li>nil</li>
|
|
221
|
+
</ul></dd>
|
|
222
|
+
<dt><a name="label:9" id="label:9"><code>scale_height</code></a></dt><!-- RDLabel: "scale_height" -->
|
|
223
|
+
<dd>
|
|
224
|
+
<p>return a value of the scale height on the log-pressure coordinate.
|
|
225
|
+
default value is "7000 m".</p>
|
|
226
|
+
<pre>RETURN VALUE
|
|
227
|
+
* scale height (UNumeric)</pre></dd>
|
|
228
|
+
<dt><a name="label:10" id="label:10"><code>scale_height=(<var>h</var>)</code></a></dt><!-- RDLabel: "scale_height=" -->
|
|
229
|
+
<dd>
|
|
230
|
+
<p>set a value of the scale height on the log-pressure coordinate.</p>
|
|
231
|
+
<pre>RETURN VALUE
|
|
232
|
+
* nil</pre></dd>
|
|
233
|
+
<dt><a name="label:11" id="label:11"><code>radius</code></a></dt><!-- RDLabel: "radius" -->
|
|
234
|
+
<dd>
|
|
235
|
+
<p>return a value of the radius of the planet. default value is "6.37E6 m".</p>
|
|
236
|
+
<pre>RETURN VALUE
|
|
237
|
+
* radius (UNumeric)</pre></dd>
|
|
238
|
+
<dt><a name="label:12" id="label:12"><code>radius=(<var>a</var>)</code></a></dt><!-- RDLabel: "radius=" -->
|
|
239
|
+
<dd>
|
|
240
|
+
<p>set a value of the radius of the planet.</p>
|
|
241
|
+
<pre>RETURN VALUE
|
|
242
|
+
* nil</pre></dd>
|
|
243
|
+
<dt><a name="label:13" id="label:13"><code>rot_period</code></a></dt><!-- RDLabel: "rot_period" -->
|
|
244
|
+
<dd>
|
|
245
|
+
<p>return a value of the rotation period of the planet.
|
|
246
|
+
default value is "8.64E4 s".</p>
|
|
247
|
+
<pre>RETURN VALUE
|
|
248
|
+
* rotation period (UNumeric)</pre></dd>
|
|
249
|
+
<dt><a name="label:14" id="label:14"><code>rot_period=(<var>rp</var>)</code></a></dt><!-- RDLabel: "rot_period=" -->
|
|
250
|
+
<dd>
|
|
251
|
+
<p>set a value of the rotation period of the planet.</p>
|
|
252
|
+
<pre>RETURN VALUE
|
|
253
|
+
* nil</pre></dd>
|
|
254
|
+
<dt><a name="label:15" id="label:15"><code>g_forces</code></a></dt><!-- RDLabel: "g_forces" -->
|
|
255
|
+
<dd>
|
|
256
|
+
<p>return a value of the gravitational acceleration on the surface.
|
|
257
|
+
default value is "9.81 m/s2".</p>
|
|
258
|
+
<pre>RETURN VALUE
|
|
259
|
+
* rotation period (UNumeric)</pre></dd>
|
|
260
|
+
<dt><a name="label:16" id="label:16"><code>g_forces=(<var>g</var>)</code></a></dt><!-- RDLabel: "g_forces=" -->
|
|
261
|
+
<dd>
|
|
262
|
+
<p>set a value of the gravitational acceleration on the surface.</p>
|
|
263
|
+
<pre>RETURN VALUE
|
|
264
|
+
* nil</pre></dd>
|
|
265
|
+
<dt><a name="label:17" id="label:17"><code>p00</code></a></dt><!-- RDLabel: "p00" -->
|
|
266
|
+
<dd>
|
|
267
|
+
<p>return a value of the reference surface pressure.
|
|
268
|
+
default value is "1.0E5 Pa".</p>
|
|
269
|
+
<pre>RETURN VALUE
|
|
270
|
+
* reference surface pressure (UNumeric)</pre></dd>
|
|
271
|
+
<dt><a name="label:18" id="label:18"><code>p00=(<var>p00</var>)</code></a></dt><!-- RDLabel: "p00=" -->
|
|
272
|
+
<dd>
|
|
273
|
+
<p>set a value of the reference surface pressure.</p>
|
|
274
|
+
<pre>RETURN VALUE
|
|
275
|
+
* nil</pre></dd>
|
|
276
|
+
<dt><a name="label:19" id="label:19"><code>cp</code></a></dt><!-- RDLabel: "cp" -->
|
|
277
|
+
<dd>
|
|
278
|
+
<p>return a value of the specific heat at constant pressure of the atmosphere.
|
|
279
|
+
default value is "1004.0[J.K-1.kg-1]"</p>
|
|
280
|
+
<pre>RETURN VALUE
|
|
281
|
+
* reference surface pressure (UNumeric)</pre></dd>
|
|
282
|
+
<dt><a name="label:20" id="label:20"><code>cp=(<var>cp</var>)</code></a></dt><!-- RDLabel: "cp=" -->
|
|
283
|
+
<dd>
|
|
284
|
+
<p>set a value of the specific heat at constant pressure of the atmosphere.</p>
|
|
285
|
+
<pre>RETURN VALUE
|
|
286
|
+
* nil</pre></dd>
|
|
287
|
+
<dt><a name="label:21" id="label:21"><code>gas_const</code></a></dt><!-- RDLabel: "gas_const" -->
|
|
288
|
+
<dd>
|
|
289
|
+
<p>return a value of the gas constant divided by molecular mass.
|
|
290
|
+
default value is "287.0[J.K-1.kg-1]".</p>
|
|
291
|
+
<pre>RETURN VALUE
|
|
292
|
+
* reference surface pressure (UNumeric)</pre></dd>
|
|
293
|
+
<dt><a name="label:22" id="label:22"><code>gas_const=(<var>r</var>)</code></a></dt><!-- RDLabel: "gas_const=" -->
|
|
294
|
+
<dd>
|
|
295
|
+
<p>set a value of the gas constant divided by molecular mass.</p>
|
|
296
|
+
<pre>RETURN VALUE
|
|
297
|
+
* nil</pre></dd>
|
|
298
|
+
<dt><a name="label:23" id="label:23"><code>get_constants</code></a></dt><!-- RDLabel: "get_constants" -->
|
|
299
|
+
<dd>
|
|
300
|
+
<p>return values of the scale height, radius, rotation period,
|
|
301
|
+
gravitational acceleration, reference surface pressure, specific heat,
|
|
302
|
+
gas constant.</p>
|
|
303
|
+
<pre>RETURN VALUE
|
|
304
|
+
* scale height (UNumeric)
|
|
305
|
+
* radius (UNumeric)
|
|
306
|
+
* rotation period (UNumeric)
|
|
307
|
+
* gravitational acceleration (UNumeric)
|
|
308
|
+
* reference surface pressure (UNumeric)
|
|
309
|
+
* specific heat at constant pressure of the atmosphere (UNumeric)
|
|
310
|
+
* gas constant divided by molecular mass (UNumeric)</pre></dd>
|
|
311
|
+
<dt><a name="label:24" id="label:24"><code>set_constants(<var>scale_height</var>, <var>radius</var>, <var>rot_period</var>, <var>g_forces</var>, <var>p00</var>, <var>cp</var>, <var>gas_const</var>)</code></a></dt><!-- RDLabel: "set_constants" -->
|
|
312
|
+
<dd>
|
|
313
|
+
<p>set values of the scale height, radius, rotation period, and gravitational
|
|
314
|
+
acceleration.</p>
|
|
315
|
+
<pre>ARGUMENTS
|
|
316
|
+
* scale height (UNumeric)
|
|
317
|
+
* radius (UNumeric)
|
|
318
|
+
* rotation period (UNumeric)
|
|
319
|
+
* gravitational acceleration (UNumeric)
|
|
320
|
+
* reference surface pressure (UNumeric)
|
|
321
|
+
* specific heat at constant pressure of the atmosphere (UNumeric)
|
|
322
|
+
* gas constant divided by molecular mass (UNumeric)
|
|
323
|
+
|
|
324
|
+
RETURN VALUE
|
|
325
|
+
* nil</pre></dd>
|
|
326
|
+
<dt><a name="label:25" id="label:25"><code>make_gphys(*<var>ax_ary</var>)</code></a></dt><!-- RDLabel: "make_gphys" -->
|
|
327
|
+
<dd>
|
|
328
|
+
<p>make GPhys objects from Axis or VArray. data components is VArray of
|
|
329
|
+
<!-- Reference, RDLabel "Axis.pos" doesn't exist --><em class="label-not-found">Axis.pos</em><!-- Reference end -->.</p>
|
|
330
|
+
<p>ARGUMENTS</p>
|
|
331
|
+
<ul>
|
|
332
|
+
<li>ax_ary (Array): an Array each objects are <!-- Reference, RDLabel "Axis" doesn't exist --><em class="label-not-found">Axis</em><!-- Reference end --> or <!-- Reference, RDLabel "VArray" doesn't exist --><em class="label-not-found">VArray</em><!-- Reference end -->.</li>
|
|
333
|
+
</ul>
|
|
334
|
+
<p>RETURN VALUE</p>
|
|
335
|
+
<ul>
|
|
336
|
+
<li>gp_ary (Array): an Array each objects are <!-- Reference, RDLabel "GPhys" doesn't exist --><em class="label-not-found">GPhys</em><!-- Reference end -->.</li>
|
|
337
|
+
</ul></dd>
|
|
338
|
+
<dt><a name="label:26" id="label:26"><code>to_w_if_omega(<var>gp</var>, <var>z</var>)</code></a></dt><!-- RDLabel: "to_w_if_omega" -->
|
|
339
|
+
<dd>
|
|
340
|
+
<p>convert to velocity(W) if <!-- Reference, RDLabel "gp" doesn't exist --><em class="label-not-found">gp</em><!-- Reference end --> is pressure velocity(Omega).
|
|
341
|
+
decide from units <!-- Reference, RDLabel "gp.data.units" doesn't exist --><em class="label-not-found">gp.data.units</em><!-- Reference end -->. if it compatible with "m/s" then
|
|
342
|
+
deem it <!-- Reference, RDLabel "W" doesn't exist --><em class="label-not-found">W</em><!-- Reference end -->, else if "Pa/s" then deem it <!-- Reference, RDLabel "Omega" doesn't exist --><em class="label-not-found">Omega</em><!-- Reference end -->.</p>
|
|
343
|
+
<p>ARGUMENTS</p>
|
|
344
|
+
<ul>
|
|
345
|
+
<li>gp(GPhys): a GPhys which data represents velocity or pressure velocity.</li>
|
|
346
|
+
<li>z(GPhys): a GPhys which data represents z-coordinate.</li>
|
|
347
|
+
</ul>
|
|
348
|
+
<p>RETURN VALUE</p>
|
|
349
|
+
<ul>
|
|
350
|
+
<li>gp_w(GPhys): a GPhys which data represents velocity or pressure velocity</li>
|
|
351
|
+
</ul></dd>
|
|
352
|
+
<dt><a name="label:27" id="label:27"><code>to_z_if_pressure(<var>gp</var>)</code></a></dt><!-- RDLabel: "to_z_if_pressure" -->
|
|
353
|
+
<dd>
|
|
354
|
+
<p>convert to altitude(z) if <!-- Reference, RDLabel "gp" doesn't exist --><em class="label-not-found">gp</em><!-- Reference end --> is pressure coordinate (p).
|
|
355
|
+
decide from units <!-- Reference, RDLabel "gp.data.units" doesn't exist --><em class="label-not-found">gp.data.units</em><!-- Reference end -->. if it compatible with "Pa" then
|
|
356
|
+
deem it (p).</p>
|
|
357
|
+
<p>ARGUMENTS</p>
|
|
358
|
+
<ul>
|
|
359
|
+
<li>gp(GPhys): a GPhys which data represents z or pressure coordinate.</li>
|
|
360
|
+
</ul>
|
|
361
|
+
<p>RETURN VALUE</p>
|
|
362
|
+
<ul>
|
|
363
|
+
<li>gp_z(GPhys): a GPhys which data represents z-coordinate.</li>
|
|
364
|
+
</ul></dd>
|
|
365
|
+
<dt><a name="label:28" id="label:28"><code>to_p_if_altitude(<var>gp</var>)</code></a></dt><!-- RDLabel: "to_p_if_altitude" -->
|
|
366
|
+
<dd>
|
|
367
|
+
<p>convert to pressure(p) if <!-- Reference, RDLabel "gp" doesn't exist --><em class="label-not-found">gp</em><!-- Reference end --> is altitude(z).
|
|
368
|
+
decide from units <!-- Reference, RDLabel "gp.data.units" doesn't exist --><em class="label-not-found">gp.data.units</em><!-- Reference end -->. if it compatible with "m" then
|
|
369
|
+
deem it (z).</p>
|
|
370
|
+
<p>ARGUMENTS</p>
|
|
371
|
+
<ul>
|
|
372
|
+
<li>gp(GPhys): a GPhys which data represents z or pressure coordinate.</li>
|
|
373
|
+
</ul>
|
|
374
|
+
<p>RETURN VALUE</p>
|
|
375
|
+
<ul>
|
|
376
|
+
<li>gp_p(GPhys): a GPhys which data represents p-coordinate.</li>
|
|
377
|
+
</ul></dd>
|
|
378
|
+
<dt><a name="label:29" id="label:29"><code>to_theta_if_temperature(<var>gp_t</var>, <var>z</var>, <var>flag_temp_or_theta</var>=<var>true</var>)</code></a></dt><!-- RDLabel: "to_theta_if_temperature" -->
|
|
379
|
+
<dd>
|
|
380
|
+
<p>convert <!-- Reference, RDLabel "gp" doesn't exist --><em class="label-not-found">gp</em><!-- Reference end --> to potential temperature(\theta) if
|
|
381
|
+
<!-- Reference, RDLabel "flag_temp_or_theta" doesn't exist --><em class="label-not-found">flag_temp_or_theta</em><!-- Reference end --> is true. </p>
|
|
382
|
+
<p>ARGUMENTS</p>
|
|
383
|
+
<ul>
|
|
384
|
+
<li><p>gp_t(GPhys): a GPhys which data represents potential temperature or </p>
|
|
385
|
+
<pre>temperature.</pre></li>
|
|
386
|
+
<li>z(GPhys) : a GPhys which data represents z-coordinate.</li>
|
|
387
|
+
<li>flag_temp_or_theta(True or False): a flagment if <!-- Reference, RDLabel "gp_t" doesn't exist --><em class="label-not-found">gp_t</em><!-- Reference end --> convert to.</li>
|
|
388
|
+
</ul>
|
|
389
|
+
<p>RETURN VALUE</p>
|
|
390
|
+
<ul>
|
|
391
|
+
<li>gp_theta(GPhys): a GPhys which data represents potential temperature.</li>
|
|
392
|
+
</ul></dd>
|
|
393
|
+
<dt><a name="label:30" id="label:30"><code>to_rad_if_deg(<var>gp</var>)</code></a></dt><!-- RDLabel: "to_rad_if_deg" -->
|
|
394
|
+
<dd>
|
|
395
|
+
<p>convert to radian if <!-- Reference, RDLabel "gp.data.units" doesn't exist --><em class="label-not-found">gp.data.units</em><!-- Reference end --> is degrees. </p>
|
|
396
|
+
<p>ARGUMENTS</p>
|
|
397
|
+
<ul>
|
|
398
|
+
<li>gp(GPhys): a GPhys which represents angle (radian or degree).</li>
|
|
399
|
+
</ul>
|
|
400
|
+
<p>RETURN VALUE</p>
|
|
401
|
+
<ul>
|
|
402
|
+
<li>gp_rad(GPhys): a GPhys which units is radian.</li>
|
|
403
|
+
</ul></dd>
|
|
404
|
+
<dt><a name="label:31" id="label:31"><code>eddy_products(<var>gp_u</var>, <var>gp_v</var>, <var>gp_w</var>, <var>gp_t</var>, <var>dimname</var>)</code></a></dt><!-- RDLabel: "eddy_products" -->
|
|
405
|
+
<dd>
|
|
406
|
+
<p>calculate eddy products along "dimname" dimension. now in this documents,
|
|
407
|
+
' means eddy from zonal mean, and () means zonal mean.</p>
|
|
408
|
+
<p>ARGUMENTS</p>
|
|
409
|
+
<ul>
|
|
410
|
+
<li>gp_u(GPhys): a GPhys which data represents zonal-wind velocity(m/s).</li>
|
|
411
|
+
<li>gp_v(GPhys): a GPhys which data represents merdional-wind velocity(m/s).</li>
|
|
412
|
+
<li>gp_w(GPhys): a GPhys which data represents vertical-wind velocity(m/s).</li>
|
|
413
|
+
<li>gp_t(GPhys): a GPhys which data represents temperature(K).</li>
|
|
414
|
+
</ul>
|
|
415
|
+
<p>RETURN VALUE</p>
|
|
416
|
+
<ul>
|
|
417
|
+
<li>uv_dash(GPhys): a GPhys which represents (gp_u'*gp_v').</li>
|
|
418
|
+
<li>vt_dash(GPhys): a GPhys which represents (gp_v'*gp_t').</li>
|
|
419
|
+
<li>uw_dash(GPhys): a GPhys which represents (gp_u'*gp_w').</li>
|
|
420
|
+
</ul></dd>
|
|
421
|
+
<dt><a name="label:32" id="label:32"><code>remove_0_at_poles(<var>cos_gp</var>)</code></a></dt><!-- RDLabel: "remove_0_at_poles" -->
|
|
422
|
+
<dd>
|
|
423
|
+
<p>set value if the cos(latitude) is nearly equal to 0 (|x|< 1e-6) at poles.
|
|
424
|
+
at North pole, new value is <a href="(a_cos_lat.val[0] + a_cos_lat.val[1])">2</a>
|
|
425
|
+
and at South pole <a href="(a_cos_lat.val[-1] + a_cos_lat.val[-2])">2</a></p>
|
|
426
|
+
<pre>ARGUMENTS
|
|
427
|
+
* cos_gp(GPhys): a GPhys which represents latitude.
|
|
428
|
+
|
|
429
|
+
RETURN VALUE
|
|
430
|
+
* new_cos_gp(GPhys): a GPhys which value at poles displaceed.</pre></dd>
|
|
431
|
+
<dt><a name="label:33" id="label:33"><code>preparate_for_vector_on_merdional_section(<var>xax</var>, <var>yax</var>)</code></a></dt><!-- RDLabel: "preparate_for_vector_on_merdional_section" -->
|
|
432
|
+
<dd>
|
|
433
|
+
<p>preparate for <!-- Reference, RDLabel "GGraph::vector_on_merdional_section" doesn't exist --><em class="label-not-found">GGraph::vector_on_merdional_section</em><!-- Reference end --> in
|
|
434
|
+
vector_on_merdional_section. </p>
|
|
435
|
+
<ol>
|
|
436
|
+
<li>check <!-- Reference, RDLabel "yax" doesn't exist --><em class="label-not-found">yax</em><!-- Reference end --> if it is proportional to p </li>
|
|
437
|
+
<li>get axis ( a*phi, z ) </li>
|
|
438
|
+
</ol>
|
|
439
|
+
<pre>ARGUMENTS
|
|
440
|
+
* xax(VArray): a VArray which represents x axis.
|
|
441
|
+
* yax(VArray): a VArray which represents y axis.
|
|
442
|
+
|
|
443
|
+
RETURN VALUE
|
|
444
|
+
* va_aphi(VArray): a VArray which represents x-coordinate(radius * phi).
|
|
445
|
+
* va_z(VArray): a VArray which represents z-coordinate.
|
|
446
|
+
* was_proportional_to_p(True or False): flag original axis proportional to
|
|
447
|
+
pressure or z.</pre></dd>
|
|
448
|
+
<dt><a name="label:34" id="label:34"><code>strm_rmean(<var>v_rmean</var>, <var>yzdims</var>=[<var>0</var>,<var>1</var>])</code></a></dt><!-- RDLabel: "strm_rmean" -->
|
|
449
|
+
<dd>
|
|
450
|
+
<p>Calculate mass stream function for residual zonal mean circulation. </p>
|
|
451
|
+
<p>ARGUMENTS</p>
|
|
452
|
+
<ul>
|
|
453
|
+
<li>v_rmean (GPhys): a GPhys which is residual zonal mean V.</li>
|
|
454
|
+
<li>yzdims (Array): an Array which represents axis.</li>
|
|
455
|
+
</ul>
|
|
456
|
+
<p>RETURN VALUE</p>
|
|
457
|
+
<ul>
|
|
458
|
+
<li>gp_strm (GPhys): a GPhys which is mass stream function on merdional section.</li>
|
|
459
|
+
</ul></dd>
|
|
460
|
+
</dl>
|
|
461
|
+
<h2><a name="label:35" id="label:35">Constants:</a></h2><!-- RDLabel: "Constants:" -->
|
|
462
|
+
<dl>
|
|
463
|
+
<dt><a name="label:36" id="label:36"><code>Deriv_methods</code></a></dt><!-- RDLabel: "Deriv_methods" -->
|
|
464
|
+
<dd>
|
|
465
|
+
<p>derivative method name [ 'cderiv', 'threepoint_O2nd_deriv' ]</p></dd>
|
|
466
|
+
</dl>
|
|
467
|
+
|
|
468
|
+
</body>
|
|
469
|
+
</html>
|