module zahyo_high implicit none integer,private :: i integer,parameter :: nlon=1920,nlat=960,nlev=340 real,parameter :: lon(nlon)=(/(i*360./nlon,i=0,nlon-1)/) real,parameter :: lat(nlat)=(/& & 89.90625, 89.71875, 89.53125, 89.34375, 89.15625, 88.96875, 88.78125, 88.59375, 88.40625, 88.21875,& & 88.03125, 87.84375, 87.65625, 87.46875, 87.28125, 87.09375, 86.90625, 86.71875, 86.53125, 86.34375,& & 86.15625, 85.96875, 85.78125, 85.59375, 85.40625, 85.21875, 85.03125, 84.84375, 84.65625, 84.46875,& & 84.28125, 84.09375, 83.90625, 83.71875, 83.53125, 83.34375, 83.15625, 82.96875, 82.78125, 82.59375,& & 82.40625, 82.21875, 82.03125, 81.84375, 81.65625, 81.46875, 81.28125, 81.09375, 80.90625, 80.71875,& & 80.53125, 80.34375, 80.15625, 79.96875, 79.78125, 79.59375, 79.40625, 79.21875, 79.03125, 78.84375,& & 78.65625, 78.46875, 78.28125, 78.09375, 77.90625, 77.71875, 77.53125, 77.34375, 77.15625, 76.96875,& & 76.78125, 76.59375, 76.40625, 76.21875, 76.03125, 75.84375, 75.65625, 75.46875, 75.28125, 75.09375,& & 74.90625, 74.71875, 74.53125, 74.34375, 74.15625, 73.96875, 73.78125, 73.59375, 73.40625, 73.21875,& & 73.03125, 72.84375, 72.65625, 72.46875, 72.28125, 72.09375, 71.90625, 71.71875, 71.53125, 71.34375,& & 71.15625, 70.96875, 70.78125, 70.59375, 70.40625, 70.21875, 70.03125, 69.84375, 69.65625, 69.46875,& & 69.28125, 69.09375, 68.90625, 68.71875, 68.53125, 68.34375, 68.15625, 67.96875, 67.78125, 67.59375,& & 67.40625, 67.21875, 67.03125, 66.84375, 66.65625, 66.46875, 66.28125, 66.09375, 65.90625, 65.71875,& & 65.53125, 65.34375, 65.15625, 64.96875, 64.78125, 64.59375, 64.40625, 64.21875, 64.03125, 63.84375,& & 63.65625, 63.46875, 63.28125, 63.09375, 62.90625, 62.71875, 62.53125, 62.34375, 62.15625, 61.96875,& & 61.78125, 61.59375, 61.40625, 61.21875, 61.03125, 60.84375, 60.65625, 60.46875, 60.28125, 60.09375,& & 59.90625, 59.71875, 59.53125, 59.34375, 59.15625, 58.96875, 58.78125, 58.59375, 58.40625, 58.21875,& & 58.03125, 57.84375, 57.65625, 57.46875, 57.28125, 57.09375, 56.90625, 56.71875, 56.53125, 56.34375,& & 56.15625, 55.96875, 55.78125, 55.59375, 55.40625, 55.21875, 55.03125, 54.84375, 54.65625, 54.46875,& & 54.28125, 54.09375, 53.90625, 53.71875, 53.53125, 53.34375, 53.15625, 52.96875, 52.78125, 52.59375,& & 52.40625, 52.21875, 52.03125, 51.84375, 51.65625, 51.46875, 51.28125, 51.09375, 50.90625, 50.71875,& & 50.53125, 50.34375, 50.15625, 49.96875, 49.78125, 49.59375, 49.40625, 49.21875, 49.03125, 48.84375,& & 48.65625, 48.46875, 48.28125, 48.09375, 47.90625, 47.71875, 47.53125, 47.34375, 47.15625, 46.96875,& & 46.78125, 46.59375, 46.40625, 46.21875, 46.03125, 45.84375, 45.65625, 45.46875, 45.28125, 45.09375,& & 44.90625, 44.71875, 44.53125, 44.34375, 44.15625, 43.96875, 43.78125, 43.59375, 43.40625, 43.21875,& & 43.03125, 42.84375, 42.65625, 42.46875, 42.28125, 42.09375, 41.90625, 41.71875, 41.53125, 41.34375,& & 41.15625, 40.96875, 40.78125, 40.59375, 40.40625, 40.21875, 40.03125, 39.84375, 39.65625, 39.46875,& & 39.28125, 39.09375, 38.90625, 38.71875, 38.53125, 38.34375, 38.15625, 37.96875, 37.78125, 37.59375,& & 37.40625, 37.21875, 37.03125, 36.84375, 36.65625, 36.46875, 36.28125, 36.09375, 35.90625, 35.71875,& & 35.53125, 35.34375, 35.15625, 34.96875, 34.78125, 34.59375, 34.40625, 34.21875, 34.03125, 33.84375,& & 33.65625, 33.46875, 33.28125, 33.09375, 32.90625, 32.71875, 32.53125, 32.34375, 32.15625, 31.96875,& & 31.78125, 31.59375, 31.40625, 31.21875, 31.03125, 30.84375, 30.65625, 30.46875, 30.28125, 30.09375,& & 29.90625, 29.71875, 29.53125, 29.34375, 29.15625, 28.96875, 28.78125, 28.59375, 28.40625, 28.21875,& & 28.03125, 27.84375, 27.65625, 27.46875, 27.28125, 27.09375, 26.90625, 26.71875, 26.53125, 26.34375,& & 26.15625, 25.96875, 25.78125, 25.59375, 25.40625, 25.21875, 25.03125, 24.84375, 24.65625, 24.46875,& & 24.28125, 24.09375, 23.90625, 23.71875, 23.53125, 23.34375, 23.15625, 22.96875, 22.78125, 22.59375,& & 22.40625, 22.21875, 22.03125, 21.84375, 21.65625, 21.46875, 21.28125, 21.09375, 20.90625, 20.71875,& & 20.53125, 20.34375, 20.15625, 19.96875, 19.78125, 19.59375, 19.40625, 19.21875, 19.03125, 18.84375,& & 18.65625, 18.46875, 18.28125, 18.09375, 17.90625, 17.71875, 17.53125, 17.34375, 17.15625, 16.96875,& & 16.78125, 16.59375, 16.40625, 16.21875, 16.03125, 15.84375, 15.65625, 15.46875, 15.28125, 15.09375,& & 14.90625, 14.71875, 14.53125, 14.34375, 14.15625, 13.96875, 13.78125, 13.59375, 13.40625, 13.21875,& & 13.03125, 12.84375, 12.65625, 12.46875, 12.28125, 12.09375, 11.90625, 11.71875, 11.53125, 11.34375,& & 11.15625, 10.96875, 10.78125, 10.59375, 10.40625, 10.21875, 10.03125, 9.84375, 9.65625, 9.46875,& & 9.28125, 9.09375, 8.90625, 8.71875, 8.53125, 8.34375, 8.15625, 7.96875, 7.78125, 7.59375,& & 7.40625, 7.21875, 7.03125, 6.84375, 6.65625, 6.46875, 6.28125, 6.09375, 5.90625, 5.71875,& & 5.53125, 5.34375, 5.15625, 4.96875, 4.78125, 4.59375, 4.40625, 4.21875, 4.03125, 3.84375,& & 3.65625, 3.46875, 3.28125, 3.09375, 2.90625, 2.71875, 2.53125, 2.34375, 2.15625, 1.96875,& & 1.78125, 1.59375, 1.40625, 1.21875, 1.03125, 0.84375, 0.65625, 0.46875, 0.28125, 0.09375,& &- 0.09375,- 0.28125,- 0.46875,- 0.65625,- 0.84375,- 1.03125,- 1.21875,- 1.40625,- 1.59375,- 1.78125,& &- 1.96875,- 2.15625,- 2.34375,- 2.53125,- 2.71875,- 2.90625,- 3.09375,- 3.28125,- 3.46875,- 3.65625,& &- 3.84375,- 4.03125,- 4.21875,- 4.40625,- 4.59375,- 4.78125,- 4.96875,- 5.15625,- 5.34375,- 5.53125,& &- 5.71875,- 5.90625,- 6.09375,- 6.28125,- 6.46875,- 6.65625,- 6.84375,- 7.03125,- 7.21875,- 7.40625,& &- 7.59375,- 7.78125,- 7.96875,- 8.15625,- 8.34375,- 8.53125,- 8.71875,- 8.90625,- 9.09375,- 9.28125,& &- 9.46875,- 9.65625,- 9.84375,-10.03125,-10.21875,-10.40625,-10.59375,-10.78125,-10.96875,-11.15625,& &-11.34375,-11.53125,-11.71875,-11.90625,-12.09375,-12.28125,-12.46875,-12.65625,-12.84375,-13.03125,& &-13.21875,-13.40625,-13.59375,-13.78125,-13.96875,-14.15625,-14.34375,-14.53125,-14.71875,-14.90625,& &-15.09375,-15.28125,-15.46875,-15.65625,-15.84375,-16.03125,-16.21875,-16.40625,-16.59375,-16.78125,& &-16.96875,-17.15625,-17.34375,-17.53125,-17.71875,-17.90625,-18.09375,-18.28125,-18.46875,-18.65625,& &-18.84375,-19.03125,-19.21875,-19.40625,-19.59375,-19.78125,-19.96875,-20.15625,-20.34375,-20.53125,& &-20.71875,-20.90625,-21.09375,-21.28125,-21.46875,-21.65625,-21.84375,-22.03125,-22.21875,-22.40625,& &-22.59375,-22.78125,-22.96875,-23.15625,-23.34375,-23.53125,-23.71875,-23.90625,-24.09375,-24.28125,& &-24.46875,-24.65625,-24.84375,-25.03125,-25.21875,-25.40625,-25.59375,-25.78125,-25.96875,-26.15625,& &-26.34375,-26.53125,-26.71875,-26.90625,-27.09375,-27.28125,-27.46875,-27.65625,-27.84375,-28.03125,& &-28.21875,-28.40625,-28.59375,-28.78125,-28.96875,-29.15625,-29.34375,-29.53125,-29.71875,-29.90625,& &-30.09375,-30.28125,-30.46875,-30.65625,-30.84375,-31.03125,-31.21875,-31.40625,-31.59375,-31.78125,& &-31.96875,-32.15625,-32.34375,-32.53125,-32.71875,-32.90625,-33.09375,-33.28125,-33.46875,-33.65625,& &-33.84375,-34.03125,-34.21875,-34.40625,-34.59375,-34.78125,-34.96875,-35.15625,-35.34375,-35.53125,& &-35.71875,-35.90625,-36.09375,-36.28125,-36.46875,-36.65625,-36.84375,-37.03125,-37.21875,-37.40625,& &-37.59375,-37.78125,-37.96875,-38.15625,-38.34375,-38.53125,-38.71875,-38.90625,-39.09375,-39.28125,& &-39.46875,-39.65625,-39.84375,-40.03125,-40.21875,-40.40625,-40.59375,-40.78125,-40.96875,-41.15625,& &-41.34375,-41.53125,-41.71875,-41.90625,-42.09375,-42.28125,-42.46875,-42.65625,-42.84375,-43.03125,& &-43.21875,-43.40625,-43.59375,-43.78125,-43.96875,-44.15625,-44.34375,-44.53125,-44.71875,-44.90625,& &-45.09375,-45.28125,-45.46875,-45.65625,-45.84375,-46.03125,-46.21875,-46.40625,-46.59375,-46.78125,& &-46.96875,-47.15625,-47.34375,-47.53125,-47.71875,-47.90625,-48.09375,-48.28125,-48.46875,-48.65625,& &-48.84375,-49.03125,-49.21875,-49.40625,-49.59375,-49.78125,-49.96875,-50.15625,-50.34375,-50.53125,& &-50.71875,-50.90625,-51.09375,-51.28125,-51.46875,-51.65625,-51.84375,-52.03125,-52.21875,-52.40625,& &-52.59375,-52.78125,-52.96875,-53.15625,-53.34375,-53.53125,-53.71875,-53.90625,-54.09375,-54.28125,& &-54.46875,-54.65625,-54.84375,-55.03125,-55.21875,-55.40625,-55.59375,-55.78125,-55.96875,-56.15625,& &-56.34375,-56.53125,-56.71875,-56.90625,-57.09375,-57.28125,-57.46875,-57.65625,-57.84375,-58.03125,& &-58.21875,-58.40625,-58.59375,-58.78125,-58.96875,-59.15625,-59.34375,-59.53125,-59.71875,-59.90625,& &-60.09375,-60.28125,-60.46875,-60.65625,-60.84375,-61.03125,-61.21875,-61.40625,-61.59375,-61.78125,& &-61.96875,-62.15625,-62.34375,-62.53125,-62.71875,-62.90625,-63.09375,-63.28125,-63.46875,-63.65625,& &-63.84375,-64.03125,-64.21875,-64.40625,-64.59375,-64.78125,-64.96875,-65.15625,-65.34375,-65.53125,& &-65.71875,-65.90625,-66.09375,-66.28125,-66.46875,-66.65625,-66.84375,-67.03125,-67.21875,-67.40625,& &-67.59375,-67.78125,-67.96875,-68.15625,-68.34375,-68.53125,-68.71875,-68.90625,-69.09375,-69.28125,& &-69.46875,-69.65625,-69.84375,-70.03125,-70.21875,-70.40625,-70.59375,-70.78125,-70.96875,-71.15625,& &-71.34375,-71.53125,-71.71875,-71.90625,-72.09375,-72.28125,-72.46875,-72.65625,-72.84375,-73.03125,& &-73.21875,-73.40625,-73.59375,-73.78125,-73.96875,-74.15625,-74.34375,-74.53125,-74.71875,-74.90625,& &-75.09375,-75.28125,-75.46875,-75.65625,-75.84375,-76.03125,-76.21875,-76.40625,-76.59375,-76.78125,& &-76.96875,-77.15625,-77.34375,-77.53125,-77.71875,-77.90625,-78.09375,-78.28125,-78.46875,-78.65625,& &-78.84375,-79.03125,-79.21875,-79.40625,-79.59375,-79.78125,-79.96875,-80.15625,-80.34375,-80.53125,& &-80.71875,-80.90625,-81.09375,-81.28125,-81.46875,-81.65625,-81.84375,-82.03125,-82.21875,-82.40625,& &-82.59375,-82.78125,-82.96875,-83.15625,-83.34375,-83.53125,-83.71875,-83.90625,-84.09375,-84.28125,& &-84.46875,-84.65625,-84.84375,-85.03125,-85.21875,-85.40625,-85.59375,-85.78125,-85.96875,-86.15625,& &-86.34375,-86.53125,-86.71875,-86.90625,-87.09375,-87.28125,-87.46875,-87.65625,-87.84375,-88.03125,& &-88.21875,-88.40625,-88.59375,-88.78125,-88.96875,-89.15625,-89.34375,-89.53125,-89.71875,-89.90625/) real,parameter,private :: pi=acos(-1.) real,parameter :: tal(nlat)=lat(nlat:1:-1),phi(nlat)=lat*pi/180 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 ,316.94592 ,294.69244 ,275.60632 ,259.14319 ,244.86353 ,232.42278 ,221.40031 ,211.18596 ,201.44164 ,& &192.15227 ,183.28781 ,174.8334 ,166.76891 ,159.07445 ,151.73491 ,144.73541 ,138.06084 ,131.69125 ,125.61667 ,& &119.82204 ,114.2924 ,109.01775 ,103.9881 ,99.191414 ,94.615715 ,90.251007 ,86.08728 ,82.115547 ,78.327301 ,& &74.714035 ,71.267265 ,67.979492 ,64.843697 ,61.852394 ,58.999084 ,56.277264 ,53.680935 ,51.204601 ,48.842258 ,& &46.588909 ,44.439548 ,42.389187 ,40.433819 ,38.568439 ,36.789059 ,35.092171 ,33.473278 ,31.928881 ,30.455978 ,& &29.051071 ,27.710659 ,26.432245 ,25.212826 ,24.049906 ,22.940475 ,21.882048 ,20.872614 ,19.909679 ,18.991238 ,& &18.115299 ,17.279854 ,16.482405 ,15.721957 ,14.996502 ,14.304552 ,13.645096 ,13.015635 ,12.415176 ,11.842212 ,& &11.295752 ,10.774784 ,10.277819 ,9.8036985 ,9.3512793 ,8.9199085 ,8.5084352 ,8.1159105 ,7.7414861 ,7.3843598 ,& &7.0436821 ,6.7187037 ,6.4087749 ,6.1131439 ,5.8311129 ,5.5620804 ,5.3054981 ,5.0607638 ,4.8272796 ,4.6045947 ,& &4.3922086 ,4.1895719 ,3.9962847 ,3.8119469 ,3.6360583 ,3.4683201 ,3.3083303 ,3.1556904 ,3.0101001 ,2.8712094 ,& &2.7387683 ,2.6124265 ,2.4918845 ,2.3769424 ,2.2672994 ,2.1627066 ,2.062963 ,1.9677691 ,1.8769755 ,1.7903811 ,& &1.7077868 ,1.6290419 ,1.5538968 ,1.4822016 ,1.4138061 ,1.3485605 ,1.2863648 ,1.2270185 ,1.1704226 ,1.1164256 ,& &1.0648793 ,1.0157727 ,0.96893567 ,0.92423368 ,0.88159651 ,0.84092414 ,0.80212677 ,0.7651242 ,0.72982657 ,0.69615883 ,& &0.66404593 ,0.63341296 ,0.60418987 ,0.57631677 ,0.54972845 ,0.52436519 ,0.50017679 ,0.47710332 ,0.45509481 ,0.43410119 ,& &0.41407248 ,0.39496878 ,0.37674999 ,0.35937116 ,0.34279224 ,0.3269783 ,0.3118943 ,0.29750526 ,0.28378117 ,0.27068701 ,& &0.25819787 ,0.24628866 ,0.23492439 ,0.22408514 ,0.21375082 ,0.20389147 ,0.19448209 ,0.18550768 ,0.17695327 ,0.1687938 ,& &0.16100429 ,0.15357481 ,0.14649028 ,0.13973072 ,0.13328616 ,0.12713654 ,0.12127198 ,0.1156773 ,0.11033768 ,0.10524801 ,& &0.10039435 ,0.095764138 ,0.091345936 ,0.087131716 ,0.083111979 ,0.079278238 ,0.075620979 ,0.072132215 ,0.068804935 ,0.065630637 ,& &0.062602848 ,0.059715036 ,0.056960218 ,0.054332394 ,0.051825557 ,0.049434718 ,0.047154367 ,0.044979014 ,0.042904153 ,0.040924784 ,& &0.039036911 ,0.037236027 ,0.035518143 ,0.033879753 ,0.032316856 ,0.030825952 ,0.029403547 ,0.028047139 ,0.026753223 ,0.025518807 ,& &0.024341885 ,0.023218958 ,0.022147529 ,0.021126099 ,0.020151662 ,0.019221723 ,0.018334782 ,0.017488839 ,0.016681891 ,0.015912445 ,& &0.015178492 ,0.014478038 ,0.013810084 ,0.013173126 ,0.012565667 ,0.011986204 ,0.011433241 ,0.010905776 ,0.010402308 ,0.009922293 ,& &0.0094647715 ,0.0090281516 ,0.0086116781 ,0.0082144048 ,0.0078354301 ,0.0074739545 ,0.007129177 ,0.0068002986 ,0.0064865691 ,0.006187289 ,& &0.0059018582 ,0.0056296266 ,0.0053699436 ,0.0051222099 ,0.0048858756 ,0.0046604406 ,0.004445455 ,0.0042403685 ,0.0040447316 ,0.0038581439 ,& &0.0036801558 ,0.0035104172 ,0.0033484777 ,0.003193988 ,0.0030466476 ,0.0029061073 ,0.002772066 ,0.0026441743 ,0.0025221827 ,0.0024058404 ,& &0.0022948477 ,0.0021889547 ,0.0020879614 ,0.001991668 ,0.0018997737 ,0.00181213 ,0.001728535 ,0.0016487908 ,0.0015727455 ,0.0015001503 ,& &0.0014309552 ,0.0013649593 ,0.0013019636 ,0.0012419177 ,0.0011846212 ,0.0011299752 ,0.0010778785 ,0.0010281317 ,0.000980685 ,0.00093544793,& &0.00089229079,0.00085112854,0.00081186608,0.00077441358,0.00073868595,0.00070460828,0.00067210541,0.00064109743,0.00061151938,0.00058331131,& &0.00055640307,0.00053073477,0.00050625135,0.00048289291,0.00046061442,0.00043936583,0.0004190972 ,0.00039976346,0.0003813197 ,0.00036373088,& &0.00034695194,0.00033094303,0.00031567406,0.00030111001,0.00028722096,0.00027397182,0.00026133267,0.00024927847,0.00023777921,0.00022680995,& &0.00021634564,0.00020636129,0.00019593128,0.00018411466,0.00017087616,0.00015625046,0.00014034733,0.00012339713,0.00010575509,8.789701e-05 ,& &7.0421178e-05,5.3995624e-05,3.9283099e-05,2.6840047e-05,1.7013586e-05,9.8620822e-06,5.1385573e-06,2.358194e-06 ,9.3059265e-07,2.8699617e-07/) real,parameter :: vel(nlev)=lev(nlev:1:-1),z(nlev)=7000*log(1000/lev) end module zahyo_high