擾乱場の初期値を netCDF ファイルから取得する.
real(8), allocatable :: fs_VelXDistrb(:,:) !速度 u real(8), allocatable :: sf_VelZDistrb(:,:) !速度 w real(8), allocatable :: ss_ExnerDistrb(:,:) !無次元圧力 real(8), allocatable :: ss_PotTempDistrb(:,:) !温位
擾乱場の初期値を netCDF ファイルから取得する. netCDF から値を取得するために gt4f90io を利用する.
subroutine distrbset_init(distrbfile)
use dc_trace, only: BeginSub, EndSub use gt4_history use bcset, only: ss_BC, sf_BC, fs_BC use gridset, only: DimXMin, DimXMax, DimZMin, DimZMax
character(*), intent(in) :: distrbfile !擾乱場のファイル
name = "VelXDistrb" call HistoryGet( distrbfile, name, fs_VelXDistrb ) name = "VelZDistrb" call HistoryGet( distrbfile, name, sf_VelZDistrb ) name = "ExnerDistrb" call HistoryGet( distrbfile, name, ss_ExnerDistrb ) name = "PotTempDistrb" call HistoryGet( distrbfile, name, ss_PotTempDistrb )