[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
[dennou-ruby:001815] Re: GPhys::EP_flux
��Ƿ��Ǥ���
�� �в�����dennou-ruby ML ���äƤޤ����ä���
�����Ϻ�ȤǴ�ĥ�äƤޤ��Τǡ��䤬���˥�뤷�ޤ���
�����ѹ���ݡ��ȡ�
�� SHTLIB (in DCL/MATH2�ˤ������Ĵ�´ؿ�Ÿ���饤�֥�) �ˤĤ���
���ֳ֥ǡ������оݤȤ������Ȥǡ����ֲֳ���������Ƥ���ҴѲ�
�ϥǡ����ν����ˤϤ�äƤ����Ȥ������ȤǻȤ��äǤ�����RubyDCL��
�ޤޤ�Ƥޤ��Τǡ������Ȥ���Τ�̥�ϡ�������RubyDCL��SHTLIB����
�����ե������ˤϤ���äȥХ�������Τ�������ȯ������ľ���Ƥ���
�ޤ��������������ˤϼ����ޤ�ޤ����ʤ⤷�����Ȥ������Ȥ���
��������Ф����˥������ޤ��������ΤȤ��� SHTLIB �ʳ��Υ��å�
�ǡ��ȤϤ���ޤ���Τǡ���
[[�������в�����˼���Ǥ���������롼���� shtint �Ǥϡ�
JM >= (MM+1)/2, IM >= MM+1 �Ǥʤ���Фʤ�ʤ�
�ȽƤ���ޤ�������äƤߤ�� IM, JM �������ꤢ�����پ�����
�Ƥ⤦�ޤ��Ԥ��褦�˸����ޤ�������������������ȼºݤ��Ѵ�������
������äƤ��ޤ������ꤹ��� segmentation fault ���Ǥޤ�����
�ݤΤȤ������ɤ��ʤ�Ǥ��礦���äˡ�MM+1 �� ��+1�פ�ɬ�פǤ�����
�Ȥ���ȡ��㤨��������10ʬ�䤷�����ˡ��ȿ��� 5 �ޤǤǤʤ� 4 ��
�ǤˤʤäƤ��ޤ��ޤ������ʤߤ�¾�Υ롼����Υɥ�����Ȥ˾�����
�㤤��ȯ�����ޤ��������Τ�������������𤬤���Ȼפ��ޤ�]]
(�в�����ؤμ��䤪���)
���ơ�SHTLIB�Ϥɤ���黰���������Ѥʤ褦�Ǥ��������ȿ������̤��
�⾮�����Ȥ�ʤ��Τǡ�EP flux �ǽФƤ���褦���Ӿ�ʿ��ʪ���̡���
�ʤ���ȿ�������ʬ�����˴ؤ�����ʬ���˻Ȥ��ˤ��Ը��Ǥ�������
������������Ʊ����Τԡ����ơʤĤޤ� ary(ny) --> any(nx,ny)
�Τ褦�˹������㤦�˻Ȥ��Ȥ����Τ�ͤ����ΤǤ��������ޤ��̵�̤�
¿��(�� nx=~2*ny)�Τǡ����褦�Ȥ������Ȥˤʤ�ޤ�����
���ȡ�����ʴؿ���Ϳ����ȥ��֥����ݤ��Ǥޤ������ͤޤǡ�������
�ͻ����Ƥޤ����ˤǤ��ð���������٤����֤���Ϥ����ȵ����
���˻פ��ޤ�����äơ�����geometryŪ�˥���ʥ֥�ʴؿ���Ϳ����
���ĵ���geometryŪ�˥���ʥ֥�ʷ��Ƥ�¤������ˤʤ��
���褦�ʵ������ޤ����Ȥ������Ȥ��äϤĤ������ʤ��Ǥ����������
���˶�������ޤ��������ȹͤ��ʤ��ȥȥ�֤�Ȥ����´�����ͭ
����ĺ����Ȼפ��ޤ���
�ƥ��Ȥ��Σ��� g(��,��) = cos�� (�Ӿ�����)
shts2g �� isw=1 �Ȥ���� ��g/�ߦ� �����Ǥ��ޤ������ˤǣ��ˤʤ�
�٤��Ȥ����������ˤʤꡢ���Τ֤��������ӥ��ʥ��֥��γѤ��Ǥޤ�����
������shtg2s �� isw=1 �Ȥ����(cos��)^-1��(g cos��)/�ߦ� ����
�Ǥ��ޤ��������ξ������ꤢ��ޤ��ºݡ�EP flux �����Ӿ�ʿ
���̤��Ф���Ŭ�Ѥ���Τϡ�����Τʤ��ä��ۤ��θ�ԤǤ���
�ƥ��Ȥ��Σ��� g(��,��) = sin�� cos��
���ξ��� ��g/�ߦ� �Ϥ��ޤ��Ԥ���(cos��)^-1��(g cos��)/�ߦ� ��
�㤷�����֥��γѤ��Ǥޤ�����
�� �ǡ���ʬ�Ϥɤ������
���̤Ȥ���geometry���ò������⥸�塼�����ʤ����Ȥˤʤä��Τǡ�
���Ѥ���ʬ�ѥ⥸�塼���Ĥ��롣�⥸�塼��̾�� NumRu::Deriv ��
NumRu::Derivative���оݤ�NArray�ˡ��ޤ�������� GPhys ���Ȥ߹���
����Τ�Ĥ��롣
�Ȥꤢ�����������ʬ��dx_i = x_i+1 - x_i-1 �����סˤ�Ĥ��롣
�������ϼ���ؤ��ġ��ǥե���Ȥ� dx_0 = x_1 - x_0,
dx_N = x_N - x_N-1�������϶������ˤ��碌��������ĥ�����Τ���
x[2..-1] - x[0..-3] ��Ȥ롣�������ˤ��碌����ĥ���åɲ���
�뤳�Ȥǡ���������ѹ����ưפˤ��롣���ʤߤ˺������Ū�˾Ȥ餹�ȡ�
�ˤǤϾ嵭�Υǥե���ȤΤ������ȹͤ�����Τǡ��Ȥꤢ�������Σ�
����Τߤ�Ĥ��롣���ں��֤��¤���Τǡ��������ˤϡ������Ѳ�
�Τ��Ȥ�ͤ����߷פ��Ĥġ�������ɬ����¤��ͤ��ʤࡣ�ס�
> =begin
> = GPhys �� EP-flux ���⥸�塼�����������
>
> ���Υ���, ������¸������Ƿ�⤵�������ǹԤ�줿 EP-flux �⥸�塼������˴ؤ���
> ��Ƿ�⤵��������Ǥ���碌���Ǥ���.
>
> == ����
> 2004/08/03(��) ������ -- ��������.
>
>
> == ��������⥸�塼��
> ̾�Τϲ��Ǥ���.
>
> * ((<module NumRu::Spherical>)) / ((<module NumRu::GPhys::Spherical>))
> * ���̺�ɸ�Ǥ���ʬ�黻��Ԥ��ؿ���. NArray, GPhys ���������Ф��ƺ���.
> * ((<module NumRu::Cartesian>)) / ((<module NumRu::GPhys::Cartesian>))
> * ľ��ľ����ɸ�Ǥ���ʬ�黻��Ԥ��ؿ���. NArray, GPhys ���������Ф��ƺ���.
> * ((<module NumRu::GPhys::EP_flux>))
> * EP_flux �Ф���ؿ�. GPhys ���Ф��ƺ���.
>
> === ��ʬ�黻 �⥸�塼��
>
> * NumRu::Spherical
> * ���̺�ɸ�Ǥ���ʬ�黻�⥸�塼��.
> * xderiv, yderiv, rot, div, xdiv, ydiv, grad, xgrad, ygrad, lapla, ilapla, jacobian, .., ydiv_m0, ygrad_m0, yderiv_m0, etc.
> * ���ֳ֥���åɥǡ���������Ȥ���
> * ����¾�Υ���å�(�ä˥���������)�ˤĤ��ƤϤ��������н�.
> * ISPACK ��Ƴ�����䤹���褦�ռ��������Ȥߤ����
> * ���ڥ��ȥ��Ѵ�, ��ʬ�ˤĤ��Ƥ� DCL �� SHTLIB ������
> * ((<URL|http://www.gfd-dennou.org/arch/ruby/products/ruby-dcl/ruby-dcl-doc/math2/node24.html>))
> * NumRu::Cartesian
> * �ǥ���Ⱥ�ɸ�Ǥ���ʬ�黻�⥸�塼��.
> * �����å����ɤ��Ф���ؿ��ϤȤꤢ�������ʤ�.
> * �Ȥꤢ������ʬˡ
> * cderiv(�����ʬ), fderiv(������ʬ)..etc.
>
> * NumRu::GPhys::Spherical
> * NumRu::Spherical ������
> * �ǡ���������ΰ�Υ����å�(���٤ϲ��٤��鲿�٤ޤǤ�, ���ٷ��٤ʤΤ����٤Τߤ�, �Ȥ�)
> * ��åɲ�����, �������˱����Ƽ���ؤ�����褦�ˤ���
> * ISPACK �� SHTLIB �ǥǡ���������ΰ褬(����ä�)�ۤʤ�
> * ����ñ�̤�����å�(���٤��ɤ���)����, Ƚ��
> * �̤�����˥ǡ��������äƤ�����(NCEP �Ϥ���)��, �ո����ˤ��Ƥ�����
> * ����Ū�ˤ�, (��ǽ�ʸ¤�)�ΰ��Ŭ�礵������褦�ˤ���.
> * ñ��, �ϵ�Ⱦ�¤��θ
> * NumRu::GPhys::Cartesian
> * NumRu::Cartesian ������
> * ��ɸ�ѿ��μ�ư����
> * �� : Ⱦ��������åɤ���������å�
> * ñ�̤��θ
>
> === EP_flux �⥸�塼��(NumRu::GPhys::EP_flux)
> * ���ܤ϶����Ԥ�ʤ�, ������ե륻�åȤ� EP_flux
> * �����Ƴ�����褦�Ȥ�, ��ʬ�黻�⥸�塼��Ϻ��ͤФʤ�ʤ�����, �����Ȥϰ��.
> * ���ץ����Ƕ�������Ѥ��뤫�ɤ����ڤ��ؤ�����褦�ˤ���.
> * �����ѿ�
> * ������, ������, ��ľ��(omega or w), ����(temp or theta)
> * ��ľ��ɸ�ϰ��Ϻ�ɸ(P or logP)
> * �����ѿ�
> * Fy, Fz, uv_bar, ...
> * divF ����åɤ����
> * Fy, Fz ������ˤȤ�
>
> == �������塼��(��ɸ)
>
> (1) NArray ����ʬ�黻�⥸�塼��κ��� -- 8/4(��)
> (2) GPhys ����ʬ�黻�⥸�塼��κ��� -- 8/5(��)
> (3) EP_flux ���Х⥸�塼��κ��� -- 8/6(��), 8/7(��)
> (4) ���塼�ȥꥢ��κ��� -- 8/8(��), 8/8(��)
>
> === ���
> * �����ǥ��ȥƥ��Ȥ�Ʊ���ʹ�
> * �ƥ��Ȥ�����
> * ����Ū��Υǡ����Ƚ��Ϥ���Ӥ���
> * �ƥ��Ȥ����������̤�褦�ˤʤä���, �ɥ�����Ȥ����(���ܸ� and RD)
>
> =end