=begin
= Pendergrass and Willoughby 2009 SEQ ǥ
ܥץ, Pendergrass and Willoughby 2009 (MWR) ˤ 2 оΥ䡼ꥢåǥ (SEQ) Ǥ. 

ܥǥ꼰ʤɤξܺ٤ϸʸ򻲾 (Ȥǥޥ˥奢벽ޤ).  () αʬۤȴĶβٹ¤, ɱꤷǮǮʬۤͿ, ϳ, Х󥹤 2 ۴ (ե + 徺ή) ʬۤ׻ (붯).

ե빽ϰʲΤ褦ˤʤäƤ.

 sound_make.f90 -> Ķαľ 1 ¤ǡƥȷǺ.
                   sounding.dat Ȥ̾ΤΥե뤬.
 initial_make.f90 -> ®ζʬۤ, sounding.dat β٥ǡ
                     ϳ, Х󥹤褦˲٤ζʬۤ׻.
                     ޤ, ǮǮʬۤ,
                      initial.nc Ȥ NetCDF ΥեˤޤȤ.
 SEQ.f90 -> Хץ.
            initial.nc ǡ, ι¤, ˱ 2 ۴Ĺ¤
            ׻, result_initial.nc Ȥ NetCDF ΥǡˤޤȤ.

== ץμ¹

=== ¸ץ

* ɬ

 STPK, NetCDF, gtool5

* ѤΥץ

 DCLAF90, DCL

NetCDF ǽϤΤ, 褹ġ뤬, Ǥ OK.


=== Makefile Խ
* Makefile 򼫿ȤδĶ˹碌ưʲιܤԽƲ.

 FC      = ifort  # STPK 򥳥ѥ뤷ݤΥѥ
 LIBNCDIR= /home/satoki/usr/netcdf-4.0/lib  # NetCDF  lib Υ󥹥ȡ
 LIBGT5D = /home/satoki/usr/gtool5/lib  # gtool5 (lib) Υ󥹥ȡ
 LIBGT5IN= /home/satoki/usr/gtool5/include  # gtool5 (include) Υ󥹥ȡ


=== Mkinclude Խ
* ܥǥ쥯ȥ 2 ľˤ Mkinclude Խޤ.

 ISTDIR	= ${HOME}/usr/local  # STPK Υ󥹥ȡ


=== ץΥѥ

 $ make

ǥץץå, ᥤ󥽥Фѥ뤵,

 sound_make, initial_make, SEQ

ѥ뤵.

ޤ, ץȤ, ѥץ (DCL, DCLAF90 ɬ) 򥳥ѥ뤹,

 $ make draw

 draw ץब.

=== ץμ¹
==== ץץåμ¹

 sound_make, initial_make

ϽǡΤɬפʤǤ, ȼ˥ǡѰդϤ¹ԤɬפϤʤ.

㤨, ǥ󥰤Ѱդ, initial_make Ϥ뤳ȤǽǤ. Ǥ, 2 ΥǡѰդƤ, ľ SEQ μ¹Ԥ˰ܤäƤ褤.

Ǥ, ʤ֤Ϥ뤳ȤꤷƤ.

==== namelist Խ

Ʊե namelist Ȥե뤬äƤȻפ. ,  3 ץǶ̤Ѥ͡ꥹȥեǤ.

 nr = 751  # ʿʻҿ
 nz = 21  # ľʻҿ
 dr = 2000.0  # ʿʻҲ [m]
 dz = 1000.0  # ľʻҲ [m]
 bc = '1112'  #  sign (Poisson_Gausei ڡ)
 fname = 'initial.nc'  # 2 ͥǡΥե̾
 sound_name = 'sounding.dat'  # 1 ǥ󥰥ǡΥե̾

ǥեȤǤ, Pendergrass and Willoughby (2009) ͤꤵƤ.

==== ºݤμ¹

ʲν֤Ǽ¹Ԥ.

 $ ./sound_make < namelist
 $ ./initial_make < namelist
 $ ./SEQ < namelist

ʾǷ̤ result_initial.nc ˺Ƥ.

ǡξܺ٤ NetCDF Υإå󻲾.

=== 襪ץ

ץȤƱƤץ draw μ¹ԤˤĤƤ⵭ҤƤ.

==== namelist Խ

draw_?? Ȥ̾ΤΥե draw ѤΥ͡ꥹȥեǤ.

 nr = 751  # ʿʻҿ
 nz = 21  # ľʻҿ
 rmin = 0.0  # ΰ躸ü [m]
 rmax = 150.0e3  # ΰ豦ü [m]
 zmin = 0.0  # ΰ貼ü [m]
 zmax = 20.0e3  # ΰü [m]
 fname = 'result_initial.nc'  # ɤ߽Фǡե̾
 value_name1 = 'u'  # 󥿡ɽѿ
 value_name2 = 'w'  # 顼ɤɽѿ
 cmin = -4.0  # 󥿡κǾ
 cmax = 4.0  # 󥿡κ
 smin = -2.0  # ɤκǾ
 smax = 2.0  # ɤκ

ꤷ,

 $ ./draw < draw_??

¹Ԥ, 㤨, ʲΤ褦ʥդ褵.

((<¹|URL:./sample1.png>))

((:<hr>:))

((<|URL:../index.html>))

=end
