Path: | anvargetnum.f90 |
Last Update: | Fri Sep 14 00:29:33 JST 2007 |
Authors: | Yasuhiro MORIKAWA, Eizi TOYODA |
Version: | $Id: anvargetnum.f90,v 1.7 2007/09/13 15:29:33 morikawa Exp $ |
Tag Name: | $Name: gt4f90io-20080211 $ |
Copyright: | Copyright (C) GFD Dennou Club, 2000-2005. All rights reserved. |
License: | See COPYRIGHT |
以下のサブルーチン、関数は an_generic から gtdata_generic#Get として提供されます。
Subroutine : | |||
var : | type(AN_VARIABLE), intent(in) | ||
start(:) : | integer, intent(in) | ||
cnt(:) : | integer, intent(in) | ||
stride(:) : | integer, intent(in) | ||
imap(:) : | integer, intent(in)
| ||
siz : | integer, intent(in) | ||
value(siz) : | real(DP), intent(out) | ||
iostat : | integer, intent(out) |
subroutine ANVarGetDouble(var, start, cnt, stride, imap, siz, value, iostat) use an_types, only: AN_VARIABLE use an_vartable, only: AN_VARIABLE_ENTRY, vtable_lookup use netcdf_f77, only: NF_NOERR, NF_EINVAL, NF_EINDEFINE, nf_get_varm_Double, nf_get_var1_Double, nf_redef, nf_enddef use dc_types, only: DP use dc_trace, only: BeginSub, EndSub, DbgMessage implicit none type(AN_VARIABLE), intent(in):: var integer, intent(in):: start(:) integer, intent(in):: cnt(:) integer, intent(in):: stride(:) integer, intent(in):: imap(:) ! NetCDF変数と内部データ配列のメモリ内構 ! 造間のマッピングを指定する整数ベクトル. ! 詳しくは NetCDF マニュアル ! (NF_PUT_VARM_type 等 を参照のこと) integer, intent(in):: siz real(DP), intent(out):: value(siz) integer, intent(out):: iostat integer:: nd ! var が保持する変数が依存する次元変数の数 integer:: ipos, i integer:: defstat type(AN_VARIABLE_ENTRY):: ent integer, allocatable:: istart(:), istride(:), iimap(:) character(*), parameter:: subname = 'ANVarGetDouble' continue call BeginSub(subname, fmt='varmap=%d, start=%*d, cnt=%*d, stride=%*d, imap=%*d siz=%d', i=(/var%id, start(:), cnt(:), stride(:), imap(:), siz/), n=(/size(start), size(cnt), size(stride), size(imap)/)) iostat = vtable_lookup(var, ent) if (iostat /= NF_NOERR) goto 999 ! --- nd check --- nd = 0 if (associated(ent%dimids)) nd = size(ent%dimids) if (min(size(start), size(cnt), size(stride), size(imap)) < nd) then iostat = NF_EINVAL goto 999 endif if (nd == 0) then iostat = nf_get_var1_Double(ent%fileid, ent%varid, start, value(1)) goto 999 endif ! --- stride ovarwrite buffer --- allocate(istart(nd), istride(nd), iimap(nd)) istart(1:nd) = start(1:nd) istride(1:nd) = stride(1:nd) iimap(1:nd) = imap(1:nd) ipos = 1 ! --- do read --- if (ent%varid <= 0 .or. count(cnt(1:nd) == 0) > 0) then call BeginSub('fake_map_get') call fake_map_get call EndSub('fake_map_get', 'iostat=%d', i=(/iostat/)) elseif (count(stride(1:nd) < 0) <= 0) then iostat = nf_get_varm_Double(ent%fileid, ent%varid, istart, cnt, istride, iimap, value) if (iostat == NF_EINDEFINE) then defstat = nf_enddef(ent%fileid) iostat = nf_get_varm_Double(ent%fileid, ent%varid, istart, cnt, istride, iimap, value) defstat = nf_redef(ent%fileid) end if else ! negative stride is not allowed for netcdf do, i = 1, nd if (stride(i) > 0) cycle ipos = ipos + (cnt(i) - 1) * imap(i) istart(i) = start(i) + (cnt(i) - 1) * stride(i) istride(i) = -stride(i) iimap(i) = -imap(i) call DbgMessage('dim %d negate: stride->%d start->%d map->%d', i=(/i, istride(i), istart(i), iimap(i)/)) enddo iostat = nf_get_varm_Double(ent%fileid, ent%varid, istart, cnt, istride, iimap, value(ipos)) endif deallocate(istart, istride, iimap) 999 continue call EndSub(subname, 'iostat=%d', i=(/iostat/)) contains subroutine fake_map_get ! ! ent%varid <= 0 の場合, または整数型配列 cnt に 0 が含まれる ! 場合に呼ばれるサブルーチン. ! ! 本来は nf_get_varm_Double で行う 1 次元データの入力を, ! nf_get_var1_Double で 1 要素ごとに入力し, value に格納する. ! ! iostat にはデータ入出力時のステータスが返る. ! use netcdf_f77, only: NF_EINDEFINE, nf_enddef, nf_redef, nf_get_var1_Double implicit none integer:: ofs(nd), here(nd) integer:: j, defstat continue iostat = NF_NOERR ofs(1:nd) = 0 do j = ipos + dot_product(ofs(1:nd), imap(1:nd)) here(1:nd) = istart(1:nd) + ofs(1:nd) * istride(1:nd) if (j < 1 .or. j > siz) then iostat = NF_EINVAL call DbgMessage('nf_get_var1_Double(ncid=%d, varid=%d, indx=[%*d], out-ofs=%d)', i=(/ent%fileid, ent%varid, here(1:nd), j/), n=(/nd/)) return endif if (ent%varid == 0) then value(j) = j iostat = NF_NOERR else iostat = nf_get_var1_Double(ent%fileid, ent%varid, here(1), value(j)) if (iostat == NF_EINDEFINE) then defstat = nf_enddef(ent%fileid) iostat = nf_get_var1_Double(ent%fileid, ent%varid, here(1), value(j)) defstat = nf_redef(ent%fileid) end if endif if (iostat /= NF_NOERR) return ofs(1) = ofs(1) + 1 do, j = 1, nd - 1 if (ofs(j) < cnt(j)) exit ofs(j) = 0 ofs(j + 1) = ofs(j + 1) + 1 enddo if (ofs(nd) >= cnt(nd)) exit enddo end subroutine fake_map_get end subroutine ANVarGetDouble
Subroutine : | |||
var : | type(AN_VARIABLE), intent(in) | ||
start(:) : | integer, intent(in) | ||
cnt(:) : | integer, intent(in) | ||
stride(:) : | integer, intent(in) | ||
imap(:) : | integer, intent(in)
| ||
siz : | integer, intent(in) | ||
value(siz) : | integer, intent(out) | ||
iostat : | integer, intent(out) |
subroutine ANVarGetInt(var, start, cnt, stride, imap, siz, value, iostat) use an_types, only: AN_VARIABLE use an_vartable, only: AN_VARIABLE_ENTRY, vtable_lookup use netcdf_f77, only: NF_NOERR, NF_EINVAL, NF_EINDEFINE, nf_get_varm_Int, nf_get_var1_Int, nf_redef, nf_enddef use dc_types, only: DP use dc_trace, only: BeginSub, EndSub, DbgMessage implicit none type(AN_VARIABLE), intent(in):: var integer, intent(in):: start(:) integer, intent(in):: cnt(:) integer, intent(in):: stride(:) integer, intent(in):: imap(:) ! NetCDF変数と内部データ配列のメモリ内構 ! 造間のマッピングを指定する整数ベクトル. ! 詳しくは NetCDF マニュアル ! (NF_PUT_VARM_type 等 を参照のこと) integer, intent(in):: siz integer, intent(out):: value(siz) integer, intent(out):: iostat integer:: nd ! var が保持する変数が依存する次元変数の数 integer:: ipos, i integer:: defstat type(AN_VARIABLE_ENTRY):: ent integer, allocatable:: istart(:), istride(:), iimap(:) character(*), parameter:: subname = 'ANVarGetInt' continue call BeginSub(subname, fmt='varmap=%d, start=%*d, cnt=%*d, stride=%*d, imap=%*d siz=%d', i=(/var%id, start(:), cnt(:), stride(:), imap(:), siz/), n=(/size(start), size(cnt), size(stride), size(imap)/)) iostat = vtable_lookup(var, ent) if (iostat /= NF_NOERR) goto 999 ! --- nd check --- nd = 0 if (associated(ent%dimids)) nd = size(ent%dimids) if (min(size(start), size(cnt), size(stride), size(imap)) < nd) then iostat = NF_EINVAL goto 999 endif if (nd == 0) then iostat = nf_get_var1_Int(ent%fileid, ent%varid, start, value(1)) goto 999 endif ! --- stride ovarwrite buffer --- allocate(istart(nd), istride(nd), iimap(nd)) istart(1:nd) = start(1:nd) istride(1:nd) = stride(1:nd) iimap(1:nd) = imap(1:nd) ipos = 1 ! --- do read --- if (ent%varid <= 0 .or. count(cnt(1:nd) == 0) > 0) then call BeginSub('fake_map_get') call fake_map_get call EndSub('fake_map_get', 'iostat=%d', i=(/iostat/)) elseif (count(stride(1:nd) < 0) <= 0) then iostat = nf_get_varm_Int(ent%fileid, ent%varid, istart, cnt, istride, iimap, value) if (iostat == NF_EINDEFINE) then defstat = nf_enddef(ent%fileid) iostat = nf_get_varm_Int(ent%fileid, ent%varid, istart, cnt, istride, iimap, value) defstat = nf_redef(ent%fileid) end if else ! negative stride is not allowed for netcdf do, i = 1, nd if (stride(i) > 0) cycle ipos = ipos + (cnt(i) - 1) * imap(i) istart(i) = start(i) + (cnt(i) - 1) * stride(i) istride(i) = -stride(i) iimap(i) = -imap(i) call DbgMessage('dim %d negate: stride->%d start->%d map->%d', i=(/i, istride(i), istart(i), iimap(i)/)) enddo iostat = nf_get_varm_Int(ent%fileid, ent%varid, istart, cnt, istride, iimap, value(ipos)) endif deallocate(istart, istride, iimap) 999 continue call EndSub(subname, 'iostat=%d', i=(/iostat/)) contains subroutine fake_map_get ! ! ent%varid <= 0 の場合, または整数型配列 cnt に 0 が含まれる ! 場合に呼ばれるサブルーチン. ! ! 本来は nf_get_varm_Int で行う 1 次元データの入力を, ! nf_get_var1_Int で 1 要素ごとに入力し, value に格納する. ! ! iostat にはデータ入出力時のステータスが返る. ! use netcdf_f77, only: NF_EINDEFINE, nf_enddef, nf_redef, nf_get_var1_Int implicit none integer:: ofs(nd), here(nd) integer:: j, defstat continue iostat = NF_NOERR ofs(1:nd) = 0 do j = ipos + dot_product(ofs(1:nd), imap(1:nd)) here(1:nd) = istart(1:nd) + ofs(1:nd) * istride(1:nd) if (j < 1 .or. j > siz) then iostat = NF_EINVAL call DbgMessage('nf_get_var1_Int(ncid=%d, varid=%d, indx=[%*d], out-ofs=%d)', i=(/ent%fileid, ent%varid, here(1:nd), j/), n=(/nd/)) return endif if (ent%varid == 0) then value(j) = j iostat = NF_NOERR else iostat = nf_get_var1_Int(ent%fileid, ent%varid, here(1), value(j)) if (iostat == NF_EINDEFINE) then defstat = nf_enddef(ent%fileid) iostat = nf_get_var1_Int(ent%fileid, ent%varid, here(1), value(j)) defstat = nf_redef(ent%fileid) end if endif if (iostat /= NF_NOERR) return ofs(1) = ofs(1) + 1 do, j = 1, nd - 1 if (ofs(j) < cnt(j)) exit ofs(j) = 0 ofs(j + 1) = ofs(j + 1) + 1 enddo if (ofs(nd) >= cnt(nd)) exit enddo end subroutine fake_map_get end subroutine ANVarGetInt
Subroutine : | |||
var : | type(AN_VARIABLE), intent(in) | ||
start(:) : | integer, intent(in) | ||
cnt(:) : | integer, intent(in) | ||
stride(:) : | integer, intent(in) | ||
imap(:) : | integer, intent(in)
| ||
siz : | integer, intent(in) | ||
value(siz) : | real, intent(out) | ||
iostat : | integer, intent(out) |
subroutine ANVarGetReal(var, start, cnt, stride, imap, siz, value, iostat) use an_types, only: AN_VARIABLE use an_vartable, only: AN_VARIABLE_ENTRY, vtable_lookup use netcdf_f77, only: NF_NOERR, NF_EINVAL, NF_EINDEFINE, nf_get_varm_Real, nf_get_var1_Real, nf_redef, nf_enddef use dc_types, only: DP use dc_trace, only: BeginSub, EndSub, DbgMessage implicit none type(AN_VARIABLE), intent(in):: var integer, intent(in):: start(:) integer, intent(in):: cnt(:) integer, intent(in):: stride(:) integer, intent(in):: imap(:) ! NetCDF変数と内部データ配列のメモリ内構 ! 造間のマッピングを指定する整数ベクトル. ! 詳しくは NetCDF マニュアル ! (NF_PUT_VARM_type 等 を参照のこと) integer, intent(in):: siz real, intent(out):: value(siz) integer, intent(out):: iostat integer:: nd ! var が保持する変数が依存する次元変数の数 integer:: ipos, i integer:: defstat type(AN_VARIABLE_ENTRY):: ent integer, allocatable:: istart(:), istride(:), iimap(:) character(*), parameter:: subname = 'ANVarGetReal' continue call BeginSub(subname, fmt='varmap=%d, start=%*d, cnt=%*d, stride=%*d, imap=%*d siz=%d', i=(/var%id, start(:), cnt(:), stride(:), imap(:), siz/), n=(/size(start), size(cnt), size(stride), size(imap)/)) iostat = vtable_lookup(var, ent) if (iostat /= NF_NOERR) goto 999 ! --- nd check --- nd = 0 if (associated(ent%dimids)) nd = size(ent%dimids) if (min(size(start), size(cnt), size(stride), size(imap)) < nd) then iostat = NF_EINVAL goto 999 endif if (nd == 0) then iostat = nf_get_var1_Real(ent%fileid, ent%varid, start, value(1)) goto 999 endif ! --- stride ovarwrite buffer --- allocate(istart(nd), istride(nd), iimap(nd)) istart(1:nd) = start(1:nd) istride(1:nd) = stride(1:nd) iimap(1:nd) = imap(1:nd) ipos = 1 ! --- do read --- if (ent%varid <= 0 .or. count(cnt(1:nd) == 0) > 0) then call BeginSub('fake_map_get') call fake_map_get call EndSub('fake_map_get', 'iostat=%d', i=(/iostat/)) elseif (count(stride(1:nd) < 0) <= 0) then iostat = nf_get_varm_Real(ent%fileid, ent%varid, istart, cnt, istride, iimap, value) if (iostat == NF_EINDEFINE) then defstat = nf_enddef(ent%fileid) iostat = nf_get_varm_Real(ent%fileid, ent%varid, istart, cnt, istride, iimap, value) defstat = nf_redef(ent%fileid) end if else ! negative stride is not allowed for netcdf do, i = 1, nd if (stride(i) > 0) cycle ipos = ipos + (cnt(i) - 1) * imap(i) istart(i) = start(i) + (cnt(i) - 1) * stride(i) istride(i) = -stride(i) iimap(i) = -imap(i) call DbgMessage('dim %d negate: stride->%d start->%d map->%d', i=(/i, istride(i), istart(i), iimap(i)/)) enddo iostat = nf_get_varm_Real(ent%fileid, ent%varid, istart, cnt, istride, iimap, value(ipos)) endif deallocate(istart, istride, iimap) 999 continue call EndSub(subname, 'iostat=%d', i=(/iostat/)) contains subroutine fake_map_get ! ! ent%varid <= 0 の場合, または整数型配列 cnt に 0 が含まれる ! 場合に呼ばれるサブルーチン. ! ! 本来は nf_get_varm_Real で行う 1 次元データの入力を, ! nf_get_var1_Real で 1 要素ごとに入力し, value に格納する. ! ! iostat にはデータ入出力時のステータスが返る. ! use netcdf_f77, only: NF_EINDEFINE, nf_enddef, nf_redef, nf_get_var1_Real implicit none integer:: ofs(nd), here(nd) integer:: j, defstat continue iostat = NF_NOERR ofs(1:nd) = 0 do j = ipos + dot_product(ofs(1:nd), imap(1:nd)) here(1:nd) = istart(1:nd) + ofs(1:nd) * istride(1:nd) if (j < 1 .or. j > siz) then iostat = NF_EINVAL call DbgMessage('nf_get_var1_Real(ncid=%d, varid=%d, indx=[%*d], out-ofs=%d)', i=(/ent%fileid, ent%varid, here(1:nd), j/), n=(/nd/)) return endif if (ent%varid == 0) then value(j) = j iostat = NF_NOERR else iostat = nf_get_var1_Real(ent%fileid, ent%varid, here(1), value(j)) if (iostat == NF_EINDEFINE) then defstat = nf_enddef(ent%fileid) iostat = nf_get_var1_Real(ent%fileid, ent%varid, here(1), value(j)) defstat = nf_redef(ent%fileid) end if endif if (iostat /= NF_NOERR) return ofs(1) = ofs(1) + 1 do, j = 1, nd - 1 if (ofs(j) < cnt(j)) exit ofs(j) = 0 ofs(j + 1) = ofs(j + 1) + 1 enddo if (ofs(nd) >= cnt(nd)) exit enddo end subroutine fake_map_get end subroutine ANVarGetReal
Subroutine : |
ent%varid <= 0 の場合, または整数型配列 cnt に 0 が含まれる 場合に呼ばれるサブルーチン.
本来は nf_get_varm_Double で行う 1 次元データの入力を, nf_get_var1_Double で 1 要素ごとに入力し, value に格納する.
iostat にはデータ入出力時のステータスが返る.
subroutine fake_map_get ! ! ent%varid <= 0 の場合, または整数型配列 cnt に 0 が含まれる ! 場合に呼ばれるサブルーチン. ! ! 本来は nf_get_varm_Double で行う 1 次元データの入力を, ! nf_get_var1_Double で 1 要素ごとに入力し, value に格納する. ! ! iostat にはデータ入出力時のステータスが返る. ! use netcdf_f77, only: NF_EINDEFINE, nf_enddef, nf_redef, nf_get_var1_Double implicit none integer:: ofs(nd), here(nd) integer:: j, defstat continue iostat = NF_NOERR ofs(1:nd) = 0 do j = ipos + dot_product(ofs(1:nd), imap(1:nd)) here(1:nd) = istart(1:nd) + ofs(1:nd) * istride(1:nd) if (j < 1 .or. j > siz) then iostat = NF_EINVAL call DbgMessage('nf_get_var1_Double(ncid=%d, varid=%d, indx=[%*d], out-ofs=%d)', i=(/ent%fileid, ent%varid, here(1:nd), j/), n=(/nd/)) return endif if (ent%varid == 0) then value(j) = j iostat = NF_NOERR else iostat = nf_get_var1_Double(ent%fileid, ent%varid, here(1), value(j)) if (iostat == NF_EINDEFINE) then defstat = nf_enddef(ent%fileid) iostat = nf_get_var1_Double(ent%fileid, ent%varid, here(1), value(j)) defstat = nf_redef(ent%fileid) end if endif if (iostat /= NF_NOERR) return ofs(1) = ofs(1) + 1 do, j = 1, nd - 1 if (ofs(j) < cnt(j)) exit ofs(j) = 0 ofs(j + 1) = ofs(j + 1) + 1 enddo if (ofs(nd) >= cnt(nd)) exit enddo end subroutine fake_map_get
Subroutine : |
ent%varid <= 0 の場合, または整数型配列 cnt に 0 が含まれる 場合に呼ばれるサブルーチン.
本来は nf_get_varm_Int で行う 1 次元データの入力を, nf_get_var1_Int で 1 要素ごとに入力し, value に格納する.
iostat にはデータ入出力時のステータスが返る.
subroutine fake_map_get ! ! ent%varid <= 0 の場合, または整数型配列 cnt に 0 が含まれる ! 場合に呼ばれるサブルーチン. ! ! 本来は nf_get_varm_Int で行う 1 次元データの入力を, ! nf_get_var1_Int で 1 要素ごとに入力し, value に格納する. ! ! iostat にはデータ入出力時のステータスが返る. ! use netcdf_f77, only: NF_EINDEFINE, nf_enddef, nf_redef, nf_get_var1_Int implicit none integer:: ofs(nd), here(nd) integer:: j, defstat continue iostat = NF_NOERR ofs(1:nd) = 0 do j = ipos + dot_product(ofs(1:nd), imap(1:nd)) here(1:nd) = istart(1:nd) + ofs(1:nd) * istride(1:nd) if (j < 1 .or. j > siz) then iostat = NF_EINVAL call DbgMessage('nf_get_var1_Int(ncid=%d, varid=%d, indx=[%*d], out-ofs=%d)', i=(/ent%fileid, ent%varid, here(1:nd), j/), n=(/nd/)) return endif if (ent%varid == 0) then value(j) = j iostat = NF_NOERR else iostat = nf_get_var1_Int(ent%fileid, ent%varid, here(1), value(j)) if (iostat == NF_EINDEFINE) then defstat = nf_enddef(ent%fileid) iostat = nf_get_var1_Int(ent%fileid, ent%varid, here(1), value(j)) defstat = nf_redef(ent%fileid) end if endif if (iostat /= NF_NOERR) return ofs(1) = ofs(1) + 1 do, j = 1, nd - 1 if (ofs(j) < cnt(j)) exit ofs(j) = 0 ofs(j + 1) = ofs(j + 1) + 1 enddo if (ofs(nd) >= cnt(nd)) exit enddo end subroutine fake_map_get
Subroutine : |
ent%varid <= 0 の場合, または整数型配列 cnt に 0 が含まれる 場合に呼ばれるサブルーチン.
本来は nf_get_varm_Real で行う 1 次元データの入力を, nf_get_var1_Real で 1 要素ごとに入力し, value に格納する.
iostat にはデータ入出力時のステータスが返る.
subroutine fake_map_get ! ! ent%varid <= 0 の場合, または整数型配列 cnt に 0 が含まれる ! 場合に呼ばれるサブルーチン. ! ! 本来は nf_get_varm_Real で行う 1 次元データの入力を, ! nf_get_var1_Real で 1 要素ごとに入力し, value に格納する. ! ! iostat にはデータ入出力時のステータスが返る. ! use netcdf_f77, only: NF_EINDEFINE, nf_enddef, nf_redef, nf_get_var1_Real implicit none integer:: ofs(nd), here(nd) integer:: j, defstat continue iostat = NF_NOERR ofs(1:nd) = 0 do j = ipos + dot_product(ofs(1:nd), imap(1:nd)) here(1:nd) = istart(1:nd) + ofs(1:nd) * istride(1:nd) if (j < 1 .or. j > siz) then iostat = NF_EINVAL call DbgMessage('nf_get_var1_Real(ncid=%d, varid=%d, indx=[%*d], out-ofs=%d)', i=(/ent%fileid, ent%varid, here(1:nd), j/), n=(/nd/)) return endif if (ent%varid == 0) then value(j) = j iostat = NF_NOERR else iostat = nf_get_var1_Real(ent%fileid, ent%varid, here(1), value(j)) if (iostat == NF_EINDEFINE) then defstat = nf_enddef(ent%fileid) iostat = nf_get_var1_Real(ent%fileid, ent%varid, here(1), value(j)) defstat = nf_redef(ent%fileid) end if endif if (iostat /= NF_NOERR) return ofs(1) = ofs(1) + 1 do, j = 1, nd - 1 if (ofs(j) < cnt(j)) exit ofs(j) = 0 ofs(j + 1) = ofs(j + 1) + 1 enddo if (ofs(nd) >= cnt(nd)) exit enddo end subroutine fake_map_get