CM1

このページは、CM1(Bryan's Cloud Model)に関するメモです。

Last modified: 2013/08/14

本家HOME PAGE：

http://www.mmm.ucar.edu/people/bryan/cm1/

Reference:

Bryan, G. H., and J. M. Fritsch, 2002: A Benchmark Simulation for Moist Nonhydrostatic Numerical Models. Mon. Wea. Rev., 130, 2917-2928. LINK

---

# 以下のメモは ver. 1.16 で確認したものである。

初期擾乱をカスタマイズする場合：

src/init3d.F を編集する。

・iinit - 3D initialization option:

1 = warm bubble

2 = cold pool

3 = line of warm bubbles

4 = initialization for moist benchmark

5 = cold blob

7 = Rotunno-Emanuel tropical cyclone

8 = line thermal with random perturbations

9 = forced convergence (Loftus et al 2008)

気温・湿度及び、風速の鉛直分布をカスタマイズする場合：

src/base.F を編集する。

・isnd - Base-state sounding: 1 = Dry adiabatic

2 = Dry isothermal

3 = Dry, constant dT/dz

4 = Saturated neutrally stable (BF02 sounding)

5 = Weisman-Klemp analytic sounding

7 = External file (named 'input_sounding')

(see isnd=7 section of base.F for info)

(some soundings are available at http://www.mmm.ucar.edu/people/bryan/cm1)

8 = Dry, constant d(theta)/dz

9 = Dry, constant Brunt-Vaisala frequency

10 = Saturated, constant Brunt-Vaisala frequency

11 = Saturated, constant equiv. pot. temp.

・iwnd - Base-state wind profile: (ignored if isnd=7)

0 = zero winds

1 = RKW-type profile

2 = Weisman-Klemp supercell

3 = multicell

4 = Weisman-Klemp multicell

5 = Dornbrack etal analytic profile

6 = constant wind

地形をカスタマイズする場合：

src/init_terrain.F を編集する。

・itern - Initial topography specifications:

1 = bell-shaped hill

2 = Schaer test case

3 = (case from T. Lane and J. Doyle)

4 = specified in external GrADS file

CM1で使用できる 成形済みnamelist.input, input_sounding：

http://www.mmm.ucar.edu/people/bryan/cm1/namelists/

http://www.mmm.ucar.edu/people/bryan/cm1/soundings/

注意点：

・terrain_flag を使用する場合、出力時に地形のマスクをかけるには output_interp = 1 としておく。

このとき、stat_out の出力が0だと、12998 のエラーコードで計算が止まる。

そこで、stat_out の出力を1以上に設定しておく必要がある。

Z方向のストレッチ座標にしても地形のマスクはかからない。

12998 のエラーコードは、バイナリ出力時に stat_out の出力が一つもない場合にも生じる。

・output_interp = 1 にした場合、out_s の出力が無駄なので、以下のように変更すると便利。

○ writeout.F ! 346行目

!IF(s_out.ge.1)THEN

IF(s_out.ge.1 .and. output_interp.eq.0)THEN

○ cm1.F ! 857行目