これだけでは、まだ１つのプログラムとしては稼働はしませんが、近日中に
更に公開できる部分を増やすべく、（他の）共同研究者とも相談する予定です。
目標としては、一つのプログラムとして稼働できる状態としての公開を目指
します。
そして、これを元にして、新たにバンド計算プログラムを作ってみようと思 います。構想としては、既存の公開されたルーチンを可能な限り再利用する。 言語仕様はFOTRAN77（できればFORTRAN90以上）とする。ＶＰＰでの並列化 （ｋ点並列）が可能である。可能な限り、汎用性、可読性、効率性（特に、速 度やメモリーに関して）を考慮する（汎用性、可読性と効率性は互いに相反す る面があることに注意）。などを考慮しつつ作成します（それ以外にもいろい ろ考えていますが、まだ完全に煮詰まっていません）。
なるべく話しをどんどん進めていった方が良いと 判断しましたが、事実上頓座しています。
昨日（４／４、１９９７）、日大の石田先生から、プログラム公開に関して、 問題無しの回答を頂きました（既に融合研の森川氏からは許可を貰っています）。 両名には快諾していただき深く感謝する次第です。
これによってこれまで公開していなかったサブルーチンの大部分が公開可能 となりました。実質的に、既に公開しているものも含めると、公開する（した） サブルーチンだけから１本のバンド計算プログラムの作成がほほ可能であるこ とがわかりました（厳密に言うと、まだ完全ではありません）。
従って、取り敢えず、新たに公開可能になったサブルーチンの説明から行な いたいと思います。各サブルーチンの説明は、むしろ各論編で行なった方が良 いと思います。[ 各論編]へ行く。
４／１６、１９９８、朝、各論編での公開、説明一段落着きました。いよい よ本題に入りたいと思います。
このサブルーチンリストで公開されていないもので、（１）本バンド計算プ ログラムに必須だが公開していない、できないもの。（２）必須でなく公開し ていないもの。があります。
（１）としては、 各論編でも述べた、対角化（固有値問題）や、その他にもＦＦＴ部分（Ｍ ＦＦＴはインターネットからソースを 入手可能）、特殊関数（球べッセル関数、誤差関数、Γ関数など）などを 使っている部分は公開できません。
これらは、各自自力で入手してもらうしかないと言え ます（手段は選んで入手して下さい）。
またこのようなもの以外に、重要なサブルーチンが３つ欠けています。それ
はHPSORTとFERMI（とWIDTH2）です。
HPSORTは作った逆格子ベクトルを小さい値（つまりＧ＝０）から、大きい方
へ並べ替える（ヒープ法を使っています）サブルーチンなのですが、作成者が
外国人で、現在日本にいないため、公開の許可がとれていません。
FERMI（とWIDTH2）は、エネルギー固有値（バンド）からフェルミエネルギー
（フェルミ面）を求めるサブルーチンです。現在、公開できるように作業を進
めています（既に公開、５／１７、１９９９）。
並べ替えに関しては、大変効率は悪いですが、バブル法によっても可能で、
この方法自体は、サブルーチン
ＭＳＤのバンドの並べ替えのところで使っています。多分、これを参考に
すれば逆格子ベクトルの並べ替え用のルーチンも作れると思います（但し、バ
ブルソートはソートすべき数Ｎの二乗で計算量が増えるので、大規模な計算に
なると計算量が莫大なものになってしまいます。既に、これも公開、５／１７、１９９９）。
フェルミエネルギーを求めるサブルーチンは取り敢えず後回しにして話しを
進めていきたいと思います。
「やさしい（？）バンド計算プログラムの
作り方」で既に述べられているように、（主に）筆者の作ったプログラム
は汎用性も可読性もない、公開には元々不向きなものです。
そこで次のような方針を立ててみました。
従って、メインプログラムの構成を変えれば、ある程度自由な、プログラム
実行に関しての仕様変更が可能です（例えば、電子状態計算と分子動力学部分
との計算ステップ数の割り振り〔１回毎に原子を動かすか、１０回毎に原子を
動かすか〕など）。
更に、メインプログラム内に新しい機能を持たせたサブルーチンを呼び出す
ようにするだけで、プログラムの拡張ができます（そういうようなプログラム
にする）。
そこで筆者としては、まずメインプログラムに最重点を置き、その改良、拡 張を目指していくことにします。
むしろ問題なのは、インクルード用のプログラムPACVPP
で定義されているCOMMON文（とその〔配列〕変数）です。作った本人にとって
は、COMMON文による配列定義の多用は比較的便利（時間が経つと便利でなくな
る）なのですが、汎用性、可読性を考えると、現行のPACVPP内での定義の仕方
は、大変好ましくないものと言えます。
従って、PACVPP内のCOMMON文による定義部分にも大幅に手を加える予定でい
ます（いまのところCOMMON文による定義の仕方そのものは変更しない方針です）。
一番最初はリスト（項目２）にもあるように、現行のま まの、とにかく公開されたものを繋げただけのものを 公開したいと思います。
今回、公開したものは色分け等は施していません。最低限必要なＨＴＭＬ用
のタグ以外は何も付けていないので、テキスト形式でそのままダウンロードし
てＨＴＭＬタグの付いた最初の数行と、最後の数行を取り除けば、即FORTRAN
プログラムとして利用可能です。
全体で４９６１行もあり３００ｋｂ以上ありますので、ダウンロードの際は
気を付けて下さい。
但し、公開できない、していない部分があるので、これ単独ではバンド計算 プログラムとして稼働しません。その他に、必要なものとしてインクルード用 ファイルPACVPP とREINIT、 入力データファイルI.DAT、 プログラム制御用入力ファイルCONTL.DAT などがあります。
また忘れていけないものに、使うべき擬ポテンシャルデータが必要です（Ｎ ＣＰＳ９５〔９７〕の利用が一番簡単）。
このプログラムは公開を全く想定していません。ある時点で実際に使用して いたものを提示しています。従って、余計な注釈行や、デバッグ用の制御命令 が不用意に挿入されている可能性があります（筆者にも把握しきれていない ^^;）。
公開できない、していない部分が何とかなったとして、このプログラムの簡 単な動かし方について説明したいと思います。
（補足２、６／２６、１９９８）
これらＳＧＩ、ＳＵＮ、ＩＢＭで動かしたプログラムrevpe_d.fは、本ペー
ジで公開しているものと若干バージョンが異なります。
まずＳＧＩでは組み込み関数であるIMAGを
DIMAGに変更しています。ＳＵＮではDABS、DEXPをABS、EXPに変更して
います。またＳＵＮでは倍精度化がコンパイラオプションとして存在しないの
で、netlibにある、s2dなるツールを使って倍精度化しています。ＩＢＭでは
特に留意事項はありません（但し、オプションには扱う計算の規模に依存した
ものもあり、計算によって値が異なるものもあります。またオプションも非常
に多いですが、中には実行には必ずしも必要のないものもあるかと思われます）。
またいわゆるPC-UNIX環境では、フォートランのコ ンパイラが無い可能性があります。本ソースコードはf2c+ccでも正常に動作す るはずです。
f2c -r8 revpe_d.f
cc -I /usr/local/include/f2c -lF77 -lI77 -lm revpe_d.c -lf2c
又は、
gcc -lm revpe_d.c -lf2c
上記で、-r8は倍精度化のオプションです。その他のオプションは環境によっ
て全く異なる（必要ない）可能性があります。またf2cでrevpe_d.cが生成され
ます。-I /usr/local/include/f2cや-lI77、-lF77は必要ない場合もあります。
またcc、gcc、いずれでもコンパイル可能だと思われます。
いま公開しているrevpe_d.fではないですが、ほぼ同じバージョンのソース
コードをf2cでＣ言語に変換し、それをＣコンパイラでコンパイルし、正常に
動作することを確認しています。
動作確認をしたPC-UNIX環境としては、Linux
MLDII、Turbo Linux(ver 2.0)、PANIX(ver4.0)です。
最近、GNUがg77なるフォートランコンパイラを世に出しましたが、現時点で
は動きません（倍精度化ができない、前述の倍精度化を行なうソフトで倍精度
化した後では動く可能性がありますが、動作未確認です）。
（１／２５、２０００）現在、富士通から富士通Fortran&C@Linux
Packageが購入可能です。いくつかのPC-UNIXに対応しており、TurboLinux
（ver 4.0や4.2）も対応しています。revpe_d.fもプログラムfrt -Ad
revpe_d.fでコンパイル可能です。
まずは、メインプログラム部分に注目し、プログラム全体としての見通しの
悪いこの部分を思い切ってすっきりさせて行きたいと思います。
これ以降は、比較的規模の大きなソースコードに対する改良を行なっていく
ので、作業の進行速度が著しく遅くなります。気長に待って下さい。
現在、メインプログラムの解析中
■■■■■■■■■■■■＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿＿
0% 10% 20% 30% 40% 50% 60% 70% 80% 90% 100%（３／１１、１９９８、事実上中断状態です） 。多分、この作業は続行が困難と考えられます。作業開始の目処は立っ ていません。
現在（５／２０、１９９７）、メインプログ ラムを以下の方針で改良中です（上に同じく中断 中）。
誤差関数に関しては、「FORTRAN-77による科学技術計算サブルーチンライブ ラリ」、黒瀬、松島、松尾著、啓学出版に、プログラムリストが載っています。
また、書籍「Numerical Recipes」のFortran版にも球ベッセル関数、誤差関 数、エルミート行列の固有値問題など、多数の有用なルーチンが収録されてい ます。[関連サイト]
現在、これらのサブルーチン、関数をrevpe_d.fにくっつけて動作のテスト を行なっています（全てをチェックしている訳ではない）。既にエルミート行 列の対角化（固有値問題）のサブルーチンは正常に動作することを確認しまし た。現在も他の関数のルーチンの動作テストを行なうべく準備を進めています。
尚、これらのサブルーチン、関数の入手は、前述の文献「Fortran77による
数値計算ソフトウェア」の購入を強く勧めます。
価格は筆者が購入した当時（数年前）で７００４円（消費税３％を含む。当時
はまだ５％でないことに注意）でした。
高いと思われるかもしれませんが、これは収録されているサブルーチン、関
数も豊富で、内容も非常に良い本です。また、このような形で科学技術計算関
連のサブルーチン、関数群が入手できる例は他にあまりなく、この意味でも著
者の先生方（の御尽力）に敬意を払う（購入する）必要があると思います。
当然のことですが、購入以外の手段での入手は、著作権法が認めるごく内輪 （この場合研究室内というのはごく内輪には該当しな い可能性があります、丸善に問い合わせる必要があると思います）の 場合、著者（及び丸善）の許諾を得た場合を除けば、明らかな犯罪ですので、 くれぐれも注意してください。
現在（１０／２７、１９９７）、更に前述の文献「Fortran77による数値計 算ソフトウェア」を利用して、誤差関数DERF、球ベッセル関数DSJBESがテスト 計算において正しく動作することを確認しました。加えて、逆格子ベクトルを 小さい順に並べ替えるサブルーチン（バブルソートによる方法を採用、このた め大変遅いです）も用意しました。並べ替えに関しては本当はヒープ法による ものがあるのですが、作成者に公開の可否について確認していないので、今回 は公開できません（御容赦下さい）。
現在、高速フーリエ変換サブルーチン[ＭＦＦＴ]の最 新版が旧バージョンと同じように、revpep_d.fで正常に動作するかどうか確認 できれば、（まだ対称性には対応できませんが）一応単独（計算機やメーカー 依存の科学技術計算ライブラリに頼らないで動くこと）で動くバンド計算プロ グラムが提供できます（但し、文献「Fortran77による数値計算ソフトウェア」 の購入またはそのサブルーチンソースの合法的 な入手が前提）。
これらのサブルーチン、関数使用に関しては今のところ、ごく限定的なテス
ト計算しか行なっていません。既に対角化（固有値問題）サブルーチンに関し
て、文献「Fortran77による数値計算ソフトウェア」によるものは、縮退のあ
る系ではうまく計算でいきない可能性があるという指摘を頂いております（情
報提供に感謝）。ただ、これがこのサブルーチン固有の問題なのか、それとも
組み込まれたプログラムとの相性の問題なのかは今のところはっきりしません
（現在、調査準備中）。
このように、今後本格的な計算を行なうようになると、システム（ＯＳや機
種）や扱う系に依存するような問題（大きな系では動かなくなる、コンパイラー
のバージョンや機種によって動かない、ある特定の系の時に正常に動かなくな
るなど）に遭遇する可能性があります。
現在（１１／６、１９９７）、ベッセル関数、誤差関数に関して、文献 「Fortran77による数値計算ソフトウェア」のフロッピーに収録されているも のが利用可能（正常に稼働すること）が確認できています。更に、最新版の [ＭＦＦＴ] をくっつけても正常に稼働するかどうかをテストしましたが、これは動く場合 と全く動かない場合（システム）とに分れました。
ＤＥＣ上では、オプションの追加で対応可能であることがわかりました。プ ログラムrevpe_d.fに対して、古いＭＦＦＴルーチンから新しいＭＦＦＴルー チンへ付け替えて、コンパイルを、
f77 -r8 -O5 revpe_d.f
から、
f77 -r8 -i8 -O5 revpe_d.f
とするだけで、テスト計算では問題なく動作することを確認しました。古い ＭＦＦＴルーチンは、開発途中で倍精度用の改良がなされており、オリジナル なものとは内容が変わっています。従って、逆に新しいＭＦＦＴをそのまま付 け替えて、全く従来と同じコンパイル方法を用いただけではプログラムは動き ません。
動かすためには、新たに-i8というオプションが必要になります。実は古い ＭＦＦＴも改良する前はこのオプションを付けていました。
そして問題なのは、ＤＥＣ上では比較的すんなり（本当は全く別の理由で 問題はありますが）動きましたが、それ以外の場合には大分問題がありそ うです。
まずＰＣ（ＤＯＳ／Ｖ互換機、本当はＩＢＭ ＰＣ‐ＡＴ互換）上で動く富 士通のフォートラン９０コンパイラでは正常に動作しません。現時点であらゆ る操作を試みていますが、正しい計算結果を得られません。この様子からＶＰ Ｐ等の富士通、ファコム系スパコンでも動作しない可能性があります（まだス パコン上ではテスト計算していない）。
いろいろ試行錯誤した（１１／１１、１９９７）結果、ようやくＰＣ上の富 士通FORTRAN90コンパイラで動くようになりました（まだテスト計算のでのみ 確認）。ここまで来るためには、いろいろなことを試しました。
まずrevpe_d.fに新しいＭＦＦＴルーチンをくっつけ（旧ＭＦＦＴは削除）、 従来のコンパイル（frt -Ad revpe_d.f）ではまるで動かないので、オプショ ンをいろいろ変えて試してみました。試したのは、-AAd、/Cp、/Cp4、/Cp8な どです（それぞれのオプションの意味は各自で調べて下さい。どれも精度拡張 に関係したものです）。しかしながら結果としてはどの場合も計算途中で止まっ たり、全エネルギーが発散したりして全くまるで動きませんでした。
次に、revpe_d.f内のsin,cos,exp,cmplx,conjgなどの組み込み関数で、倍精 度表示になっているものを単精度にしたり、ＭＦＦＴ用の作業配列のサイズを 変えたりと、いろいろな場合を試してみましたが、どれも全く状況の改善につ ながりませんでした。
そうこうしている内に、ＤＥＣ側ではオプションを-r8から-r8 -i8と、-i8 という整数変数を８バイト（普通の既定値は４バイト）にするオプションを付 けると正しく動くこと、富士通FORTRANコンパイラのマニュアルを良く見てみ ると、このコンパイラにはＤＥＣの-i8に相当するオプションはないが、 INTEGER*8（８バイト整数）の宣言が可能であることがわかりました。
ひょっとしたらと思い、revpe_d.f（新ＭＦＦＴを含む）のINTEGER宣言を全 てINTEGER*8にしてみました。これは正しく動きませんでしたが、新しいＭＦ ＦＴ部分内のみのINTEGER宣言をINTEGER*8にし、frt /Cp4 revpe_d.fの条件で コンパイルすると正しく動くことが確認できました。
オプション/Cp4は、単精度→倍精度（精度拡張を行なう）、倍精度はそのま まで、記憶合わせを行なうものです。但し、INTEGER*8という宣言は、コンパ イラの仕様によって宣言できない場合もあります。
また最近、PC-UNIX系（筆者はLinuxを使用、普通の UNIX、例えばDigital UNIXでも可能）で、f2cを使ってまずrevpe_d.fをＣ言語 版に変換し（f2c -r8 revpe_d.f）、それをgcc（gnuのＣコンパイラ）で翻訳 すると、従来のＭＦＦＴを使っている場合は問題なく動きますが、新しいＭＦ ＦＴではどうしても動かないことを確認しました（実行途中で止まってしまい ます）。先の方法で回避できるかと思い、ＭＦＦＴ部分のINTEGER宣言を INTEGER*8に書き変えてみましたが、これもうまく動作しませんでした。今の ところ、新ＭＦＦＴでは、f2cによる方法での実行はできていません（１１／ ２１、１９９７）。
gnuのg77（フォートランコンパイラ）は筆者の使用しているバージョンでは 全くコンパイルできませんでした。
ヒューレットパッカード（ＨＰ）上では、旧ＭＦＦＴを新ＭＦＦＴに挿し替 え、問題なくコンパイル、実行できましたが、実行結果に問題がありました。 まず、従来のＭＦＦＴでの実行結果（全エネルギー）が、-48.7833775a.u.で した。一方、新ＭＦＦＴの方は、-48.7833740a.u.でした。エネルギー差は、 0.0000035a.u.（0.0952meV）です。一見、この位なら問題ないだろうと思いそ うですが、話はそう簡単ではありません。このＨＰ上の新ＭＦＦＴ付きのコー ドをそのままＤＥＣ側にコピーしてコンパイル、実行させると全エネルギーは ちゃんと-48.7833775a.u.と出てきます（-48.7833740a.u.にならない）。従って、新ＭＦＦＴ に関してＨＰ上では何らかの問題があることが判明しましたが、全エネルギー の違いは0.1meV以下と非常に小さく、それ以外に異常終了や、明らかにおかし な計算結果が出ている訳ではないので、現状で問題点を突き止め、修正を行な うことは非常に難しいと考えられます（５／１７、１９９９現在、ＨＰではＣ ＰＵが６４ビット化したことに伴い、フォートランコンパイラもバージョンアッ プしているため、これ以上の原因調査い続行は無理と判断しました）。
因みにＨＰでは、INTEGER*8はコンパイルエラーによりコンパイルできませ んでした（但し、コンパイラのバージョンは古いです。最新版でどうなってい るかは不明です。〔５／１７、１９９９現在いまだ不明〕）。
以上をまとめてみると、新ＭＦＦＴが動作する環境は（[目次]に戻る）、
f77 -r8 -i8 revpe_d.f
で実行可能。
frt /Cp4 revpe_d.f
で実行可能。
f77 +autodblpad +O3 +FPD revpe_d.f
で実行は可能であるが、精度に問題あり。
f2c -r8 revpe_d.f
gcc -lm -lF77 -lI77 revpe_d.f -lf2c または gcc -lm revpe_d.f -lf2c
で、旧ＭＦＴＴでは稼働するが、新ＭＦＦＴでは全く動かない。
全くコンパイルできない（コンパイルエラーで実行ファイルが作られない）。
frt -Ad revpe_d.f
で、コンパイル、実行可能。
次に、現在（１１／２２、１９９７）公開可能な追加サブルーチン（二つ） を示したいと思います。
まず、バブルソートを行なうサブルーチンBUBBLE です。そしてもう一つは、ベッセル関数呼び出しのためのサブルーチンDNV です。
バブルソートのルーチンは、各自力のある人はヒープ法や、クィックソート
にルーチンを改良変更することを強く勧めます。
バブルソートのままでは系が大きくなると大変遅くなります。
ベッセル関数呼び出しルーチンでは、文献「Fortran77による数値計算ソフ
トウェア」にある球ベッセル関数DSJBESを呼び出しています。このルーチンは
先の文献付属のフロッピー内に収録されています
（この文献を持ってない人は購入しましょう）。
（３／１１、１９９８）それでは、これらDVN,BUBBLEを付けた形での、公開 可能な全体 プログラムを示します。
この段階で、バンド計算にとって基本的に必要な部分は、全て揃ったと言え ます。後は、先の文献「Fortran77による数値計算ソフトウェア」とＭＦＦＴ ルーチンさえあれば、対称性を考慮しない場合での電子状態計算が可能です。 但し、このままではコンパイル時に対称性関係のサブルーチンが存在しないと いうエラーが出てくる可能性があります。その場合は、サブルーチンSYMM 内のサブルーチンの呼び出し（CALL文）を注釈にすれば良いと思います。
対称性を考慮しない場合の指定は、CONTL.DAT での、上から１０行目の、0 0 2 0 4、という記述を 0 0 2 0 1、とすれば対称性なしのバルクの計算ができます（表 面系では1を0に すればよい）。
サブルーチンKPMWBZ で、バルクでの対称性を考慮しない場合のｋ点の設定を行なっています（KPMSF は表面系用です）。この中で、Γ点のみの計算から等間隔でのｋ点（少しずら して設定するなど）の設定などが可能です。ただこれはプログラム内で設定す るようになっているため、変更後再コンパイルが必要です（まだ汎用性がない 好例）。
この全体プログラムは以前示した ものと、細かな内容が異なる可能性がありますが、基本的な内容、性能は 全く変わらないと思って良いと思います。
但し、これに関しての改良の目処は今のところ立っていません。対称性計算 のために必要な変数データを見つけ出し、それを出力できるようにプログラム を変更し（出力用の専用プログラムを作る）、更にデータ入力側としてのプロ グラムを新たに作成、公開する必要が生じます。つまり入力されたデータを処 理するルーチンを加えたプログラムを作らなければなりません。
システムや、ブラウザの種類、バージョンによって
はダウンロードできない可能性がありますが、残念ながらこちらでは対応（の
しようがない）できません。自力で解決されることを強く期待します。
特に、最新のバージョンのブラウザ環境では、その
ままではダウンロードできない可能性が高いです。
圧縮形式は、まずファイル群をtarでアーカイブし、それをgzipで圧縮して います（ＯＳはUNIX上であることを想定しています）。多分、先のリンク部分 をクリックすればファイルをダウンロードするためのウィンドウが開くと思い ます。アーカイブ、圧縮したのはＤＥＣのマシン上ですが、他のUNIXマシンで もtarとgzipが使えれば、問題なく解凍できると思います。
解凍の仕方は、ダウンロードしたディレクトリ上で、
gzip -d Al.tar.gz
tar -xvf Al.tar
で必要なファイルが解凍、展開されるはずです。
ＤＯＳ／Ｖ（ＩＢＭ ＰＣ‐ＡＴ互換）マシン上でも、ＤＯＳ用のtarやgzip がフリーソフトとして入手が可能です。多分、ダウンロード時にバイナリとし て正しく自分のマシン上に転送されていれば、それらＤＯＳ用のソフトでの解 凍、展開は可能かと思います。
（注意）但し、このパッケージに含まれるプ ログラムには、前述の通り、EIGCH、DSJBES、DERF、FTT、対称性関連等のサブ ルーチンは付いていません。従って、解凍 したソースをコンパイルしてもおそらくリンク段階でエラーが出て実行可能ファ イルは作られません。
実行させるためには、「とにかく[単]独で動くように する」部分を熟読して下さい。必要なファイル（サブルーチン）の入手方法、 コンパイル、実行の仕方などが書かれています。
（注意）上記ソフトなどは自力で入手して下さ
い。不正な手段での入手は厳禁です。当然、こち
らから提供することもできません。
（免責）
また、このプログラムパッケージをダウンロード、解凍、展開、コン
パイル、実行、運用等を行なった過程で生じた、如何なる障害、損害に対して、
筆者及び、無機材研は一切責任を負わないものとします。
（禁止）本パッケージの再配布は禁止します（いわゆる常識の範囲内での内輪での
使用は可）。直接、間接に関わらず悪用すること（それを目的とした使用、利
用全てを含む）も禁止します。
（補足）上記３事項に関連するより詳しい
説明があります（ご一読下さい）。
どうしても会員になりたいという奇特な方は、E-mail（のみ）で、入会希望 のメイルを筆者（メイルアドレス削除）に送って下さい（形式は、個々の創意 工夫に任せます）。筆者の独断と偏見で入 会許可を行なう可能性もあります。
最近（１０／２、１９９８）、新たに超イオン伝導体（Superionic Conductor）がSICと略されることが分かった。
（世界のＮＣＰＳ）
（８／１、１９９８）探してみると、世界中にＮＣＰＳ
が存在することがわかった（AltaVistaで検索）。その代表的なものを
以下に挙げる。
North Country Photograph Society (NCPS)
The North Carolina Poety Society (NCPS)
The Northern California Psychiatric Society (NCPS)
Novell Cross-Platform Services (NCPS)
The National Crime Prevention Strategy (NCPS)
National Commission on Public Service (NCPS)
など多数存在することが判明した。
また、Network Control Protocolsのことを略 して、NCPsと書くようです。