$Id: modificationhistory.txt 238 2012-11-12 04:11:13Z yamasaki $

--------------------------------------------------------------------------------
5th Nov. 2004

 *=*=*=*=*  m_FFT.F90  *=*=*=*=*
ASLFFTのkimg=1の場合への適用。
NECの小久保さんから貰ったソースを、山本さんがテスト計算しつつ修正。山崎がさらにほんの少し修正。

 *=*=*=*=*  m_XC_Potential.f90  *=*=*=*=*
ASLFFTのkimg=1の場合への適用に伴って、必要な修正を行った。
山本さんがテスト計算しつつ修正。山崎がさらにほんの少し修正。

 *=*=*=*=*  makefile  *=*=*=*=*
SXでASLFFTを使うように変更（山本）。

--------------------------------------------------------------------------------
2nd Nov. 2004

 *=*=*=*=*  Convergence_Check.F90  *=*=*=*=*
1. nfstopや標準入力のControl{max_iteration = n}で指定した回数で終るように修正（宇佐見）。
これまでは、nfstopや標準入力のmax_iterationで指定した回数より１回多く計算していた。
2. Empricalな計算が出来るように(make -f makefile.ep intelが実行できるように）修正。

 *=*=*=*=*  Forces.f90  *=*=*=*=*
1. call m_IS_evaluate_v_verlet(forc_l)を挿入。m_IS_evaluate_v_verletは新規作成（宇佐見）。
2. コメントアウトされていた行を削除。

 *=*=*=*=*  Forces_ep.f90  *=*=*=*=*
1. call m_IS_evaluate_v_verlet(forc_l)を挿入。m_IS_evaluate_v_verletは新規作成（宇佐見）。

 *=*=*=*=*  InputData_Analysis.F90  *=*=*=*=*
1. Empricalな計算が出来るように(make -f makefile.ep intelが実行できるように）修正。

 *=*=*=*=*  Preparation_ep.f90  *=*=*=*=*
1. Intelコンパイラがinternal errorを引き起こすのを抑止するためにmoduleのuseの順を変更。

 *=*=*=*=*  WriteDownData_onto_Files.F90  *=*=*=*=*
1. Empricalな計算が出来るように(make -f makefile.ep intelが実行できるように）修正。

 *=*=*=*=*  m_Control_Parameters.F90  *=*=*=*=*
1. Empricalな計算が出来るように(make -f makefile.ep intelが実行できるように）修正。

 *=*=*=*=*  m_Force_ep.f90  *=*=*=*=*
1. Empricalな計算が出来るように(make -f makefile.ep intelが実行できるように）修正。
 subroutine m_Force_wd_force_and_cps を加えた。

 *=*=*=*=*  m_Ionic_System.F90  *=*=*=*=*
1. Empricalな計算が出来るように(make -f makefile.ep intelが実行できるように）修正。
2. EtotalとEkinetcの和が一定値からずれる問題を解決（宇佐見）。
3. m_IS_evaluate_v_verletを作成（宇佐見）。subroutine Forcesから呼ばれる。


--------------------------------------------------------------------------------
1st Nov. 2004
 
 *=*=*=*=*  ChargeDensity_Mixing.f90  *=*=*=*=*
1. iprixxxの値をmypeの値によって設定するように（mypeが０の場合は指定した値を、そうでない
場合には、０を設定）変更したことが原因で、もとのままだとdeadlockが発生することがあった。
deadlockを防ぐための変更。

>   if( ipritotalcharge >= 2 .or. ( ipritotalcharge >= 1 .and. &
>        &   ( nspin == 2 .and. af == 0))) call m_CD_check(nfout)
を
<   call m_CD_check(nfout)
に変更した。

 *=*=*=*=*  Convergence_Check.F90  *=*=*=*=*
1. 回数制限で止るときに、標準入力ファイルnfpos.data起因なのか、nfstop.dataファイル起因なのかが分るようにした。
2. 標準出力への書出しに、if(printable)を付けた。

 *=*=*=*=*  Initialization.F90  *=*=*=*=*
1. 標準出力ファイル名を変更した(output00x、あるいはfile_names.dataで指定する場合は任意。
mype==0以外のPEの場合は、output00x_00yあるいは、somedesignatedfile_00y）ことに伴う変更。
標準出力ファイルを、call m_Files_rd_file_names_data()より後でopen。

2. 標準出力ファイルへの書き出しに if(printable) thenあるいは、if(printable)を追加。

 *=*=*=*=*  InputData_Analysis.F90  *=*=*=*=*
1. 標準出力ファイルへの書き出しに if(printable) thenあるいは、if(printable)を追加。
2. 継続ファイルのあるなしを調べた結果を保持している変数 file_existance_contfileを使用。

 *=*=*=*=*  IterationNumbers_Setting.f90  *=*=*=*=*
1. 計算状況報告用ファイルへの書き出し。出力形式の制御も行う。

 *=*=*=*=*  Postprocessing.F90  *=*=*=*=*
1. 標準出力ファイルへの書き出しに if(ipri>=n) thenあるいは、if(ipri>=n)を追加。

 *=*=*=*=*  Preparation.F90  *=*=*=*=*
1. 標準出力ファイルへの書き出しに if(ipri>=n) thenあるいは、if(ipri>=n)を追加。

 *=*=*=*=*  Preparation_for_mpi.F90  *=*=*=*=*
1. 標準出力ファイルへの書き出しの制御に、m_Control_Parametersのprintableを使用。
m_Parallel_xxxxへ引数でprintableを渡す。

 *=*=*=*=*  PseudoPotential_Construction.F90  *=*=*=*=*
1. 標準出力ファイルへの書き出しに if(ipripp>=1|2) thenあるいは、if(ipripp>=1|2)を追加。

 *=*=*=*=*  Real_space_integ.f90  *=*=*=*=*
1. $id $をくっつけた。

 *=*=*=*=*  Renewal_of_WaveFunctions.F90  *=*=*=*=*
1. 標準出力ファイルへの書き出しに if(ipri>=n) thenあるいは、if(ipri>=n)を追加。

 *=*=*=*=*  WriteDownData_onto_Files.F90  *=*=*=*=*
1. 計算状況報告用ファイルへのジョブ終了時の書き出し。出力形式の制御も行う。

 *=*=*=*=*  b_Kpoints.f90  *=*=*=*=*
1. logical :: openを logical :: open_checkに変更。

 *=*=*=*=*  b_Ldos_f77.F  *=*=*=*=*
1. $id $を付加。

 *=*=*=*=*  b_PseudoPotential.F90  *=*=*=*=*
1. 標準出力ファイルへの書き出しに if(ipri>=n) thenあるいは、if(ipri>=n)を追加。
2. 標準出力ファイルへの書き出し内容の一部を修正。

 *=*=*=*=*  b_PseudoPotential_f77.f  *=*=*=*=*
1. 字下げを少々変更。

 *=*=*=*=*  crotylm.f90  *=*=*=*=*
1. $id $を付加。

 *=*=*=*=*  heap_sort.F90  *=*=*=*=*
1. SR用の自動SMP並列抑制のコンパイラ指示行を追加。
#ifdef HIUX
*OPTION MP(P(0))
#endif
を付加。これに伴い、ファイル名もheap_sort.f90からheap_sort.F90に変更。
2. $id $を追加。

 *=*=*=*=*  input_interface.f90  *=*=*=*=*
1. 標準出力ファイルへの書き出しに if(printable) thenあるいは、if(printable)を追加。

 *=*=*=*=*  jrcat_fft.F  *=*=*=*=*
1. SR用の自動SMP並列抑制のコンパイラ指示行を追加。
#ifdef HIUX
*OPTION MP(P(0))
#endif
を付加。
2. call eqivvlをコメントアウト

1,3d0
< #ifdef HIUX
< *OPTION MP(P(0))
< #endif

 *=*=*=*=*  m_Charge_Density.F90  *=*=*=*=*
1. m_CD_checkの大幅な変更。
2. 標準出力ファイルへの書き出しに if(iprixxxx >= n) thenあるいは、if(iprixxxx >= n)を追加。
3. SR用のSMP並列化を促進するための次のコンパイラ指示行を追加。
#ifdef HIUX
*POPTION PARALLEL
#endif

 *=*=*=*=*  m_Const_Parameters.f90  *=*=*=*=*
1. 定数、AUTOMATIC = -1を追加。
2. 定数、ITERATIVE = 1, TAG_FORMAT = 0, TAG_LINE = 1, TABLE = 2を追加。

 *=*=*=*=*  m_Control_Parameters.F90  *=*=*=*=*
1. icond = AUTOMATICの場合の制御を追加。
2. 標準出力ファイルへの出力制御変数printableを追加。
3. 標準出力ファイルへの書き出しに if(iprixxxx >= n) thenあるいは、if(iprixxxx >= n)を追加。
4. 計算状況報告用ファイル(jobstatus00x)への書き出し制御変数の管理。


 *=*=*=*=*  m_Crystal_Structure.f90  *=*=*=*=*
1. 標準出力ファイルへの書き出しに if(printable) thenあるいは、if(printable)を追加。

 *=*=*=*=*  m_ELF.f90  *=*=*=*=*
1. $id $を追加。

 *=*=*=*=*  m_ES_IO.f90  *=*=*=*=*
1. 標準出力ファイルへの書き出しに if(printable) thenあるいは、if(printable)や
 if(ipri >= n) then あるいは if(ipri >= n)を付加。

 *=*=*=*=*  m_ES_WF_by_RMM.F90  *=*=*=*=*
1. 標準出力ファイルへの書き出しに if(printable) thenあるいは、if(printable)を付加。
2. 配列の割り当てのsubroutine構成を変更。
3. RMM実行時のメモリ削減。波動関数の残差を保持する配列を３分の２に削減。
4. 残差用ファイルをallocateするのに失敗したときにメッセージを出力するように修正。

 *=*=*=*=*  m_ES_WF_by_SDorCG.F90  *=*=*=*=*
1. 配列の割り当てのsubroutine構成を変更。
2. 直前の波動関数用配列zaj_oldを単精度から倍精度へ。

 *=*=*=*=*  m_ES_WF_by_submat.F90  *=*=*=*=*
1. 標準出力ファイルへの書き出しに if(iprisubmat >= 1 ) thenあるいは、if(iprisubmat >= 1)を付加。

 *=*=*=*=*  m_ES_dos.f90  *=*=*=*=*
1. sw_pdos == ONのときの処理(pdos = 0.d0; sumpdos = 0.d0)の間違いを修正。

356,358c356
<     if(sw_pdos == ON) then
<        pdos=0.d0; sumpdos=0.d0
<     end if
---
>     if(sw_pdos == ON) pdos=0.d0; sumpdos=0.d0

 *=*=*=*=*  m_ES_occup.F90  *=*=*=*=*
1. 標準出力ファイルへの書き出しに if(iprioccup >= 1 )を付加。

 *=*=*=*=*  m_Electronic_Structure.F90  *=*=*=*=*
1. 配列の割り当てのsubroutine構成を変更。mgs_allocからm_ES_mgs_allocへ。mgs_deallocからm_ES_mgs_deallocへ。
2. 標準出力ファイルへの書き出しに if(printable) thenあるいは、if(printable)を付加。

 *=*=*=*=*  m_FFT.F90  *=*=*=*=*
1. 標準出力ファイルへの書き出しに if(ipri >= n) thenあるいは、if(ipri >= n)を付加。
2. #ifndef GOEDECKER_FFT ... #endif を挿入。
3. init_fft_coefficients_arraysを init_fft_coefficients_arrays_WFとinit_fft_coefficients_arrays_CDに分離。特にSRで実行時のエラーへの対処。

 *=*=*=*=*  m_Files.f90  *=*=*=*=*
1. 計算状況報告用ファイル(jobstatus00x)名の管理。当該ファイルのopenとclose。
2. 継続用ファイルの存在確認用変数の管理。
3. 標準出力用ファイル名の指定。そのopen。

 *=*=*=*=*  m_Force.F90  *=*=*=*=*
1. 標準出力ファイルへの書き出しに if(ipriforce>= 1) を付加。

618c618,619
<     if(ipriforce >= 1) write(6,'(" !D forcmx = ",d20.12)') forcmx
---
> !!$    if(mype==0) print '(" !D forcmx = ",d20.12)', forcmx
>     write(6,'(" !D forcmx = ",d20.12)') forcmx

 *=*=*=*=*  m_Ionic_System.F90  *=*=*=*=*
1. 標準出力ファイルへの書き出しに if(printable) thenあるいは、if(printable)を付加。

 *=*=*=*=*  m_Kpoints.F90  *=*=*=*=*
1. 標準出力ファイルへの書き出しに if(printable) thenあるいは、if(printable)、
あるいは if(ipri_kp >= n) を付加。

 *=*=*=*=*  m_Ldos.F90  *=*=*=*=*
1. call check_naldos_range -> m_CtrlP_check_naldos_range
2. 標準出力ファイルへの書き出しに if(printable) thenあるいは、if(printable)を付加。

 *=*=*=*=*  m_NonLocal_Potential.f90  *=*=*=*=*
1. 標準出力ファイルへの書き出しに if(iprisnl>=1) thenを付加。

 *=*=*=*=*  m_Parallelization.F90  *=*=*=*=*
1. SRの場合の自動SMP並列の抑止。
2. m_Parallel_init_comm_worldの新規作成。
3. m_Parallel_init_mpi_firstの削除。
4. m_Parallel_init_mpiの新規作成。
5. 標準出力ファイル名の定義を削除。
6. 標準出力ファイルへの書き出しに if(printable) thenあるいは、if(printable)を付加。

1,3d0
< #ifdef HIUX
< *option MP(P(0))
< #endif
< 

 *=*=*=*=*  m_Phonon.f90  *=*=*=*=*
1. 標準出力ファイルへの書き出しに if(printable) thenあるいは、if(printable)を付加。
あるいは、if(ipriphonon >=1) thenまたはif(ipriphonon >= 1)を付加。
2. $id $を追加。

 *=*=*=*=*  m_PlaneWaveBasisSet.F90  *=*=*=*=*
1. 標準出力ファイルへの書き出しに if(printable) thenあるいは、if(printable)を付加。
あるいは、if(ipri>=1) thenまたはif(ipri>=1)を付加。

 *=*=*=*=*  m_PseudoPotential.F90  *=*=*=*=*
1. 標準出力ファイルへの書き出しに if(ipripp>=n) thenまたはif(ipripp>=n)を付加。
if(ipri >= n)はif(ipripp >= n)に変更。
2. このモジュールからの標準出力ファイルへの書き出し時には、" !PP"を行頭に付ける。

 *=*=*=*=*  m_Stress.f90  *=*=*=*=*
1. 標準出力ファイルへの書き出しに if(ipri >= n) thenまたはif(ipri >= n)を付加。

 *=*=*=*=*  m_Timing.F90  *=*=*=*=*
1. SR(HIUX)で実行時に時間測定subroutine xclockを用いる。これに伴いm_Timing.f90をm_Timing.F90に変更。
2. 標準出力ファイルへの書き出しに if(printable) thenあるいは、if(printable)を付加。
3. 現在時刻を書き出すsubroutine m_Timing_wd_timenow()を新規作成。
4. 時刻測定用変数の初期化用subroutine m_Timing_init_timer()を新規作成。
5. 現在のiterationまでのCPU時間、総elapse時間を書き出すsubroutine m_Timing_wd_statusを作成。


 *=*=*=*=*  m_Total_Energy.F90  *=*=*=*=*
2. 標準出力ファイルへの書き出しに if(printable) thenあるいは、if(printable)を付加。

 *=*=*=*=*  m_XC_Potential.f90  *=*=*=*=*
1. 不必用な定義nd3pのコメントアウト。

 *=*=*=*=*  makefile  *=*=*=*=*
1. SR11000用の定義を加えた。


 *=*=*=*=*  makefile.HI-UX  *=*=*=*=*
1. SMP並列用コンパイルオプションを加えた。

 *=*=*=*=*  makefile.UNIX_System_V  *=*=*=*=*
1. VPP300用のコンパイラオプションをコメント中に示した。

 *=*=*=*=*  miscellaneous.F90  *=*=*=*=*
1. 時間測定subroutine system_clockの呼び出し後の変数の書き出し部分をコメントアウト。


 *=*=*=*=*  rmmsubs.f  *=*=*=*=*
1. 字下げ

 *=*=*=*=*  spg+tetra.f  *=*=*=*=*
1. 一部のsubroutine中のigen,jgen（これらは引数で渡されている）の定義を
igen(3), jgen(2,3,3)からigen(ngen),jgen(2,3,ngen)に修正。
2. 標準出力ファイルへの書き出しに if(ipri_kp >= n) thenまたはif(ipri_kp >= n)を付加。


--------------------------------------------------------------------------------
10th June 2004

 *=*=*=*=*  Convergence_Check.F90  *=*=*=*=*
 1. function <EigenValues_are_Converged> において use m_Const_Parameters にGRIDを追加
 2. function <EigenValues_are_Converged> において
    >      call m_ES_cp_eko_l_to_eko_ek()
を
    <      if(ekmode == ON) call m_ES_cp_eko_l_to_eko_ek()
と変更
 3. function <AllKpoints_are_Calculated>において
    >   use m_Control_Parameters, only : nspin, m_CtrlP_set_iconvergence
の順序を変更。Intel compilerのinternal errorを回避するため

 *=*=*=*=*  Ewald_and_Structure_Factor.f90  *=*=*=*=*
 1. 不必要なuse文を消去
>   use m_Charge_Density,    only : m_CD_initial_CD_by_Gauss_func

 *=*=*=*=*  Finalization_of_mpi.F90  *=*=*=*=*
1. use文のカラム位置を調整
2. #ifdef _POSITRON_ 
   #endif
  を挿入
3. call m_pWF_deallocate_pzaj_etc()の追加

 *=*=*=*=*  Initial_Electronic_Structure.f90  *=*=*=*=*
subroutine Initial_WaveFunctions_ek において
  write(6,'("!! icond, nk_in_the_process, first_kpoint_in_this_job = ",3i8)') &
       icond, nk_in_the_process, first_kpoint_in_this_job
とフォーマットをi5からi8へ変更

 *=*=*=*=*  Initialization.F90  *=*=*=*=*
1. subroutine Initialization_set_ekmode_GRID()を追加作成

 *=*=*=*=*  InputData_Analysis.F90  *=*=*=*=*
1. 次のcall文を追加
<      call m_CtrlP_check_matm(nfout,natm) ! rmm_precal_phase_matm
2. ekmode がGRIDになる場合があることに対する対処

 *=*=*=*=*  IterationNumbers_Setting.f90  *=*=*=*=*
次のサブルーチンの追加作成。
1. subroutine IterationNumber_Setting_g()
2. subroutine IterationNumber_reset()

 *=*=*=*=*  Postprocessing.F90  *=*=*=*=*
1. ファイルの拡張子をf90からF90へ変更
2. Local (atomicおよびlayer) DOSを計算するための修正
3. subroutine Postprocessing_kの追加作成
4. 31文字を越えるsubroutine名の変更

 *=*=*=*=*  Preparation.F90  *=*=*=*=*
1. ファイルの拡張子をf90からF90へ変更
2. Local (atomicおよびlayer) DOSを計算するための修正
3. ekmode がGRIDになる場合があることに対する対処
4. subroutine Preparation_gridの追加作成
5. subroutine Preparation_grid2の追加作成

 *=*=*=*=*  Renewal_of_WaveFunctions.F90  *=*=*=*=*
1. 次の関数の追加作成。
   function EigenValue(ilevel)

 *=*=*=*=*  Renewal_of_pWaveFunctions.F90  *=*=*=*=*
1. 陽電子波動関数のソルバー。かなり変更。

x *=*=*=*=*  WriteDownData_onto_Files.F90  *=*=*=*=*
x1. 余計な行を削除

 *=*=*=*=*  WriteDownData_onto_Files_ek.f90  *=*=*=*=*
次の２つのサブルーチンを追加。
1. subroutine WriteDown_totalcpu()
2. subroutine write_energies(mnk,nk,kxyz,ek,mnt,nt,tk,ie)

 *=*=*=*=*  b_Kpoints.f90  *=*=*=*=*
1.subroutine readk0中のipriをipri_kpに変更
2. 同じく allocate(trmat(3,3)); allocate(trbp(3,2)); ..., を
    allocate(trmat(3,3),trbp(3,3),...) と変更
3. 同じく deallocateに関しても、一行に。
4. call get_trmatのsubroutine get_trmatをcontainしているものから、その外側に出し、
  他のsubroutineからも呼べるように変更した。

 *=*=*=*=*  b_Ldos_f77.F  *=*=*=*=*
1. LDOS用ファイル。新規追加

 *=*=*=*=*  bottom_Subroutines.F90  *=*=*=*=*
1. subroutine k_plus_G_vectors()において、単位BZの外にあるk点を単位BZ内に戻して
  運動エネルギーを計算。

 *=*=*=*=*  decfft_ent.F  *=*=*=*=*
positron用のentryを追加

 *=*=*=*=*  ekmain.f90  *=*=*=*=*
 LDOS用の
  call Postprocessing_k()を追加

 *=*=*=*=*  input_interface.f90  *=*=*=*=*
1. function <setUnit>の引数 usysの属性を intent(in)からintent(inout)に変更
2. function <f_readUnitCell>の引数 avecの属性を intent(in)からintent(inout)に変更

 *=*=*=*=*  m_BP_Properties.f90  *=*=*=*=*
1. 31文字を越すsubroutine名の変更
2. write フォーマットの修正（カンマの挿入）

 *=*=*=*=*  m_BerryPhase.F90  *=*=*=*=*
1. subroutine <add_deficit_term>中の use m_PseudoPotential, only : 
  の中にindex_lmt1_lmt2が二回出現するうちの一方を削除
2. subroutine <constract_of_ftq>中の use m_PseudoPotential, only : 
  の中にindex_lmt1_lmt2が二回出現するうちの一方を削除

 *=*=*=*=*  m_Charge_Density.F90  *=*=*=*=*
LDOS用の追加サブルーチン多数

 *=*=*=*=*  m_Const_Parameters.f90  *=*=*=*=*
1. LDOS用定数追加
    REGULAR_INTERVALS = 1, EAZY = 1, BY_ATOMIC_POSITIONS = 2, BY_ATOMS = 2
    SOFTPART = 1, HARDPART = 2
 :: ALDOS = 1, LAYERDOS = 2
2. 定数 GRID追加
3. 陽電子計算関係の定数追加
    ELECTRON = 0, POSITRON = 1
    TOTAL = 0

 *=*=*=*=*  m_Control_Parameters.F90  *=*=*=*=*
1. LDOS関連の変更
2. 陽電子関係のsubroutineの追加

 *=*=*=*=*  m_ELF.f90  *=*=*=*=*
1. use文中の ngabc, igfp_lが重複。一方を消去

 *=*=*=*=*  m_ES_Intgr_VlhxcQlm.f90  *=*=*=*=*
1. subroutine m_pwBS_sphrp2の引数の仕様が修正されたことに伴う修正

<        call m_pwBS_sphrp2(ilm3,rltv,ista_kngp,iend_kngp,ylm_t)
---
>        call m_pwBS_sphrp2(ilm3,rltv,ylm_t)

 *=*=*=*=*  m_ES_WF_by_RMM.F90  *=*=*=*=*
1.Renew_WF_RMM_alloc の修正
2. rmm_printout -> iprirmm

 *=*=*=*=*  m_ES_WF_by_SDorCG.F90  *=*=*=*=*
1. iprisolverの追加。
  isolver_coreを書き出すときに使う

 *=*=*=*=*  m_ES_dos.f90  *=*=*=*=*
1. LDOS計算用にいろいろと変更

 *=*=*=*=*  m_ES_occup.F90  *=*=*=*=*
1. check_if_metalic_flagの追加
2. subroutine check_if_metalic_ek をcheck_if_metalicに変更
   subroutine m_ESoc_check_if_metalicの追加
  状況(ekmode)に応じて処理を変えるようにした(Phaseからもekcalからも使う。
      m_ESoc_check_if_metalicを通じて）

 *=*=*=*=*  m_Electronic_Structure.F90  *=*=*=*=*
1. eko_lとeko1_lをallocateしたあとで初期化
2. subroutine <Vnonlocal_W_part_sum_over_lmt1>中でsnlをsnl_tに複写していたのを
  止めた。処理の高速化
3. subroutine <m_FFT_W_Vlocal_W>の引数が修正されたことに伴う修正
その他いろいろ

 *=*=*=*=*  m_FFT.F90  *=*=*=*=*
1. positron用の変更

 *=*=*=*=*  m_Files.f90  *=*=*=*=*
1. ldos用ファイルの追加

 *=*=*=*=*  m_Force.F90  *=*=*=*=*
1. subroutine m_pwBS_sphrp2の引数の変更

 *=*=*=*=*  m_Force_ep.f90  *=*=*=*=*
1. 著者(author)の変更

 *=*=*=*=*  m_Ionic_System.F90  *=*=*=*=*
1. LDOS関係の変更

 *=*=*=*=*  m_IterationNumbers.F90  *=*=*=*=*
1. subroutine m_Iter_positrion_set の作成

 *=*=*=*=*  m_Kpoints.F90  *=*=*=*=*
1. wd_kpointsの引数を加えた
2. subroutine m_Kp_cp_kxyz_to_vkxyzの作成
3. subroutine m_Kp_set_kv3の作成

 *=*=*=*=*  m_Ldos.F90  *=*=*=*=*
1. 新規作成モジュール。LDOS用。

 *=*=*=*=*  m_NonLocal_Potential.f90  *=*=*=*=*
1. 一部のwrite文を、iprisnl>=2 という条件下で実行できるように変更

 *=*=*=*=*  m_Parallelization.F90  *=*=*=*=*
1. map_eの書き出し個所を整理
2. map_ekの書き出し個所を修正

 *=*=*=*=*  m_Phonon.f90  *=*=*=*=*
1. 31文字を越える語を短く修正
2. writeフォーマットを修正
など

 *=*=*=*=*  m_PlaneWaveBasisSet.F90  *=*=*=*=*
1. subroutine m_pwBS_sphrp2の引数に配列ylmの大きさを加えた
    ista_kngp -> ista_ylm
    iend_kngp -> iend_ylm
2. 陽電子用のfft-boxへの写像配列などの追加
3. subroutine m_pwBS_increase_kg1の作成追加
など

 *=*=*=*=*  m_Positron_Wave_Functions.F90  *=*=*=*=*
大変更

 *=*=*=*=*  m_PseudoPotential.F90  *=*=*=*=*
1. stop 文の内容変更
<              stop ' istress /= 0 or 1'
---
>              stop ' istress \= 0 or 1'
2. 31文字を越えるファイル名の修正
3. subroutine read_ps_xctypeにおいて語の比較を大文字小文字の区別をしないように変更

 *=*=*=*=*  m_Stress.f90  *=*=*=*=*
1. subroutine m_pwBS_sphrp2の引数変更

 *=*=*=*=*  m_XC_Potential.f90  *=*=*=*=*
1. subroutine m_pwBS_sphrp2の引数変更

 *=*=*=*=*  makefile  *=*=*=*=*
1. LDOS計算用ファイルの追加
2. 陽電子計算用ファイルの追加
など

 *=*=*=*=*  makefile.HI-UX  *=*=*=*=*
1. コメント行追加
  #FFLAG = -c -64 -Oss -procnum=8 -i,P -msg=w

 *=*=*=*=*  makefile.OSF1  *=*=*=*=*
1. コンパイルオプションの変更
< FFLAG = -O3 -c -fast

 *=*=*=*=*  makefile.ek  *=*=*=*=*
1. LDOS計算用ファイルの追加
など

 *=*=*=*=*  mdmain.f90  *=*=*=*=*
1. 陽電子計算用の処理の実行

 *=*=*=*=*  spg+tetra.f  *=*=*=*=*
1. subroutine <wtetra>中のwrite文の内容を少し変更

