module zahyo
  implicit none
  integer,private :: i
  integer,parameter :: nlon=128,nlat=64,nlev=124
  real,parameter ::lon(nlon)=(/((i-1)*360./128.,i=1,128)/)
  real,parameter ::lat(nlat)=(/&
       & 87.8638  , 85.096527, 82.312912, 79.525604, 76.7369  ,&
       & 73.947517, 71.157753,  68.36776, 65.577606, 62.787354,&
       & 59.997021, 57.206631, 54.416199, 51.625732, 48.835239,&
       & 46.044727, 43.254196,  40.46365, 37.673088, 34.882523,&
       & 32.091946, 29.301359, 26.510769, 23.720175, 20.929575,&
       & 18.138971, 15.348365, 12.557756, 9.7671452, 6.9765334,&
       & 4.1859207, 1.3953069,-1.3953069,-4.1859207,-6.9765334,&
       &-9.7671452,-12.557756,-15.348365,-18.138971,-20.929575,&
       &-23.720175,-26.510769,-29.301359,-32.091946,-34.882523,&
       &-37.673088, -40.46365,-43.254196,-46.044727,-48.835239,&

       &-51.625732,-54.416199,-57.206631,-59.997021,-62.787354,&
       &-65.577606, -68.36776,-71.157753,-73.947517,  -76.7369,&
       &-79.525604,-82.312912,-85.096527,-87.8638/)
  real,parameter :: pi=acos(-1.)
  real,parameter :: tal(nlat)=lat(nlat:1:-1),phi(nlat)=(/(lat(i)*pi/180,i=1,nlat)/)
  real,parameter :: lev(nlev)=(/&
       &997.49927    ,991.49854    ,982.99695    ,971.99561    ,958.49304    ,&
       &941.98975    ,922.48578    ,900.48254    ,875.97699    ,848.97247    ,&
       &819.96735    ,786.95099    ,746.92285    ,698.88483    ,642.83325    ,&
       &581.80328    ,519.27246    ,460.30438    ,412.37866    ,374.90823    ,&
       &343.42712    ,305.1062     ,261.19775    ,223.60966    ,191.4277     ,&
       &163.8775     ,140.29442    ,120.10387    ,102.81879    ,88.021881    ,&
       &75.354477    ,64.510246    ,55.225895    ,47.278111    ,40.47448     ,&
       &34.649654    ,29.662722    ,25.393799    ,21.739408    ,18.610655    ,&
       &15.93258     ,13.639717    ,11.676627    ,9.9961815    ,8.5575638    ,&
       &7.3259969    ,6.2717037    ,5.3691502    ,4.5964541    ,3.9349279    ,&

       &3.3686352    ,2.8838825    ,2.4688289    ,2.1134844    ,1.8093534    ,&
       &1.5489892    ,1.326049     ,1.1351857    ,0.97181851   ,0.83197844   ,&
       &0.71224427   ,0.60974282   ,0.52199066   ,0.44686964   ,0.38256067   ,&
       &0.3275052    ,0.28037396   ,0.24002287   ,0.20547763   ,0.17590386   ,&
       &0.15059212   ,0.12892254   ,0.11036584   ,0.09448228   ,0.080885492  ,&
       &0.069244854  ,0.05927949   ,0.050748114  ,0.043444846  ,0.037192326  ,&
       &0.031839669  ,0.027257456  ,0.02333477   ,0.019976677  ,0.017101737  ,&
       &0.014640503  ,0.012533514  ,0.010729859  ,0.009185764  ,0.0078637609 ,&
       &0.0067320229 ,0.0057632183 ,0.0049338117 ,0.0042237719 ,0.0036159067 ,&
       &0.0030955311 ,0.0026500493 ,0.0022686739 ,0.0019422075 ,0.0016626571 ,&

       &0.0014233802 ,0.0012185762 ,0.0010431918 ,0.00089304463,0.00076452183,&
       &0.00065449567,0.00056030322,0.00047967111,0.00041064052,0.00035154296,&
       &0.00030094912,0.00025763526,0.00022056192,0.00018609798,0.00015180476,&
       &0.00011893146,8.8784742e-05,6.2559826e-05,4.1141917e-05,2.4917852e-05,&
       &1.3679726e-05,6.6807565e-06,2.8386589e-06,1.0222096e-06/)
  !100hPa->29,10hPa->44,1hPa->59,0.1hPa->74,0.01hPa->88,0.001hPa->104!
  real,parameter :: vel(nlev)=lev(nlev:1:-1)
  real,parameter :: g=9.81,H=7000.,Ps=1000.,z(nlev)=-H*log(lev/Ps),rho(nlev)=100*lev/(g*H)
end module zahyo