[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