SPMODEL - ISPACK ? gt4f90io ?????????? -

1
SPMODEL- ISPACK ? gt4f90io ?????????? -
?? davis gtool4 ???? ???????

• ????
• SPMODEL ??????
• http//www.gfd-dennou.org/arch/spmodel/

2
?????

• ???????????????
• ???????????????????
• ?????????????????????
• ?????????????
• ??????
• ??????
• ???????????????

3
SPMODEL

• ??????????????????
• GFD ???????????????
• ?????(spml)??????????????
• ????
• ??????????
• ?????
• ??????????????
• ???????????
• ???????????
• ??????? ISPACK, ??? I/O ? gt4f90io

4
SPMODEL ?????(spml)

• ISPACK ? F90 ????????
• F90 ??????????
• ???? spml-0.2.2
• http//www.gfd-dennou.org/arch/spmodel/spml.tar.gz
  
• ?????????
• x86 Linux Fujitsu frt, Intel ifc (6.0, 7.0, 8.0)
• Fujitsu VPP
• NEC SX (?????)

5
?????????

• GFD ????????????????
• 1 ?????
• ???????, KdV ???
• 2 ?????
• ????????????, ??ß??????
• ?????????, ???????
• 3 ?????
• ??????????
• ?? MHD ???
• ?????????(????)

6
??????????????

• ??????
• ????
• ????
• ????????????? gt4f90io ?
• 4 ????????? call ????
• ?????????????????
• netCDF ????????

7
SPMODEL ??

• ??
• ????????????????????
• ?????? g_Zeta
• ???????? s_Zeta
• ????
• (???)_(??)_(???) ??????
• ??????? s_Zetas_g(g_Zeta)
• ?? s_Div_g(g_Zeta)
• ???????????????

8
???????????(1)

• 1 ???????????
• ???? (Euler ??????)
• ???????????????

s_ZetaA s_ZetaB

dt( - Us_Dx_s(s_Zeta)
Ds_Dx_s(s_Dx_s(s_Zeta))
http//www.gfd-dennou.org/arch/spmodel/1d-cyclic-e
/advection-diffusion/sample/f90/advdiff1.f90
9
???????????(2)

• ???????(?????)
• ?????? leapfrog ?????

w_EtaAw_EtaB 2dt
(-
w_Div_xy_xy(xy_w(w_Eta)xy_GradLon_w(w_Chi)/R,

xy_w(w_Eta)xy_GradLat_w(w_Chi)/
R)/R w_Jacobian_w_w(w_Eta, w_Psi)/R2)
http//www.gfd-dennou.org/arch/spmodel/2d-sphere-w
/shallow/shallow-zd/f90/shallow_zd.f90
10
???????????(2)

• ??????
• Leapforg ?? semi-implicit ??????
• ??????????????????

11
???????????(2)

• call get_dHsfc
• w_HsfcA((1 dt2GravH0n(.1))w_HsfcB
  
• 2dtw_dHsfc)/(1 -
  dt2GravH0n(.1))
• End

12
???????????(2)

• soubroutin get_dHsfc
• w_dHsfc- (w_Div_xy_xy(xy_Hsfcxy_GradLon_w(w_Chi)
  /R,
  xy_Hsfcxy_GradLat_w(w_Chi)/R)/R
• - w_Jacobian_w_w(w_Hsfc,w_Psi)/R
  2)
• - H0w_DivB
  
• -H0dt(w_Div_xy_xy(xy_Etaxy_GradLon_w(w_P
  si)/R,
• 
  xy_Etaxy_GradLat_w(w_Psi)/R)/R
• w_Jacobian_w_w(w_Eta,w_Chi)/
  R2
• w_Lapla_w(w_E)/R2)

13
?????????????

• use gt4history
  ! ???????
• call HistoryCreate( ! ????????
• fileoutput_file, titleRun Number,
  sourceshallow institutionGFD Dennou Club
  SPMODEL project, dims(/lon,lat,t
  /),dimsizes(/im,jm,0/), longnames(/longit
  ude,latitude ,time /),
  units(/'deg.','deg.','sec.'/),
  originreal(tinit), intervalreal(intrstdelt
  ) )
• call HistoryAddVariable( varnameeta, !
  ??????
• dims(/lon,lat,t /),
  lognameVorticity, units1/s,xtypedouble)
• do it 1, n ! ?????
• call HistoryPut(eta,xy_Eta) ! ????
• end do
• call HistoryClose ! ?????

14
???? CPU ??

• ??????????
• 15 ????????
• ?t1800 sec, 720 ????(T106 ??t900 sec)
• ??? hhmmss
• ?????????

T21 T42 T63 T106
Xeon 2.4GHz 000003 000033 000109 001352
VPP800 000003 000009 ---------- 000153
15
SPMODEL ????

• ??????????????????
• ?????????
• ??????????????
• ?????????????????
• do ?????????????
• ?????
• ????
• ????????????????
• ????????????????????
• ????? Grid Modeling System ???

16
gt4f90io ????

• ????????????????
• 4 ????????????????
• netCDF ?????????????????

17
SPMODEL ???

• ??????????
• ex). ??????????????
• Grid Modeling System ????
• ?????????????????
• ???????????
• ????????????????
• ?????
• ???????, ???????????
• ???????????????

18
?????

• ?????????
• ??????????????
• 1 ?????????????????????
• ???????
• ?????????
• 3 ???????????
• ??????????????
• ?? MHD ??????

19
?? URL

• SPMODEL
• http//www.gfd-dennou.org/arch/spmodel/
• ISPACK
• http//www.gfd-dennou.org/arch/ispack/
• gt4f90io
• http//www.gfd-dennou.org/arch/gtool4/

20
???

• debian ?????????
• /etc/apt/sources.list ???? 2 ????
• apt ???????
• ISPACK, netCDF ?????????
• ???http//www/gfddennou.org/arch/spmodel

deb ftp//www.gfd-dennou.org/arch/spmodel/debian
woody/ deb-src ftp//www.gfd-dennou.org/arch/spmod
el/debian woody/
apt-get install spml
21
??