!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! !! !! !! GNU General Public License !! !! !! !! This file is part of the Flexible Modeling System (FMS). !! !! !! !! FMS is free software; you can redistribute it and/or modify !! !! it and are expected to follow the terms of the GNU General Public !! !! License as published by the Free Software Foundation. !! !! !! !! FMS is distributed in the hope that it will be useful, !! !! but WITHOUT ANY WARRANTY; without even the implied warranty of !! !! MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the !! !! GNU General Public License for more details. !! !! !! !! You should have received a copy of the GNU General Public License !! !! along with FMS; if not, write to: !! !! Free Software Foundation, Inc. !! !! 59 Temple Place, Suite 330 !! !! Boston, MA 02111-1307 USA !! !! or see: !! !! http://www.gnu.org/licenses/gpl.txt !! !! !! !!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!! subroutine set_eta(km, ks, ptop, ak, bk) implicit none integer, intent(in):: km ! vertical dimension integer, intent(out):: ks ! number of pure p layers real, intent(out):: ak(km+1) real, intent(out):: bk(km+1) real, intent(out):: ptop ! model top (Pa) ! local real a24(25),b24(25) ! GFDL AM2L24 real a32(33),b32(33) real a32w(33),b32w(33) real a48(49),b48(49) real a48x(49),b48x(49) real a55(56),b55(56) real a60(61),b60(61) real a64(65),b64(65) real a100(101),b100(101) #ifdef MARS_GCM real:: p0= 6.0E2 real:: pc= 1.0e-3 real a36(37),b36(37) ! For Mars GCM real a28(29),b28(29) ! Mars 28 levels #else real:: p0=1000.E2 real:: pc=200.E2 #endif MARS_GCM real pt, pint, lnpe, dlnp real press(km+1) integer k ! Definition: press(i,j,k) = ak(k) + bk(k) * ps(i,j) #ifdef MARS_GCM !------- Mars vertical levels ----------------------------- data a28 / & 2.0000000e-02,5.7381273e-02,1.9583981e-01,5.9229580e-01, & 1.5660228e+00,2.4454965e+00,2.7683755e+00,2.8851692e+00, & 2.9172228e+00,2.9087039e+00,2.8598939e+00,2.7687652e+00, & 2.6327054e+00,2.4509219e+00,2.2266811e+00,1.9684681e+00, & 1.6894832e+00,1.4055812e+00,1.1324258e+00,8.8289177e-01, & 6.6548467e-01,4.8401020e-01,3.3824119e-01,2.2510704e-01, & 1.3995719e-01,7.7611554e-02,3.3085503e-02,2.0000000e-03, & 0.0000000e+00 / data b28 / & 0.0000000e+00,0.0000000e+00,0.0000000e+00,0.0000000e+00, & 0.0000000e+00,1.9366394e-03,7.4419133e-03,1.6227267e-02, & 2.7075192e-02,4.3641000e-02,6.8106804e-02,1.0280240e-01, & 1.4971954e-01,2.0987133e-01,2.8270233e-01,3.6581610e-01, & 4.5520230e-01,5.4593599e-01,6.3310970e-01,7.1267629e-01, & 7.8196151e-01,8.3977530e-01,8.8620345e-01,9.2223168e-01, & 9.4934533e-01,9.6919618e-01,9.8337259e-01,9.9326941e-01, & 1.0000000e+00 / data a36 / & 2.0000000000e-03, 6.3299931399e-03, 1.2501646444e-02, & 2.4215113043e-02, 4.6000346612e-02, 8.5702012910e-02, & 1.5659441036e-01, 2.8061882660e-01, 4.9318818941e-01, & 8.5008792314e-01, 1.4370449074e+00, 2.0054945771e+00, & 2.3335916338e+00, 2.5221957520e+00, 2.6264602874e+00, & 2.6762480591e+00, 2.6870173757e+00, 2.6657174771e+00, & 2.6140573091e+00, 2.5304052915e+00, 2.4110051828e+00, & 2.2508885698e+00, 2.0446940567e+00, 1.8074358726e+00, & 1.5511877409e+00, 1.2904762231e+00, 1.0396626704e+00, & 8.1055438393e-01, 6.1095195817e-01, 4.4434410651e-01, & 3.1051916267e-01, 2.0665632614e-01, 1.2848513437e-01, & 7.1249514632e-02, 3.0373097709e-02, 2.1040298410e-16, & 0.0000000000e+00 / data b36 / & 0.0000000000e+00, 0.0000000000e+00, 0.0000000000e+00, & 0.0000000000e+00, 0.0000000000e+00, 0.0000000000e+00, & 0.0000000000e+00, 0.0000000000e+00, 0.0000000000e+00, & 0.0000000000e+00, 0.0000000000e+00, 6.2833541661e-04, & 2.5671934940e-03, 6.0923860938e-03, 1.1725248117e-02, & 2.0238695380e-02, 3.2676843179e-02, 5.0373898012e-02, & 7.4959525687e-02, 1.0833815470e-01, 1.5263019875e-01, & 2.1006506655e-01, 2.8281982789e-01, 3.6584598427e-01, & 4.5513763376e-01, 5.4577516959e-01, 6.3285639147e-01, & 7.1233851133e-01, 7.8155013540e-01, 8.3930250332e-01, & 8.8568132228e-01, 9.2167127115e-01, 9.4875609753e-01, & 9.6858585131e-01, 9.8274719475e-01, 9.9326941371e-01, & 1.0000000000e+00 / #endif MARS_GCM !----------------------------------------------- ! GFDL AM2-L24: modified by SJL at the model top !----------------------------------------------- ! data a24 / 100.0000, 1050.0000, 3474.7942, 7505.5556, 12787.2428, & data a24 / 100.0000, 903.4465, 3474.7942, 7505.5556, 12787.2428, & 19111.3683, 21854.9274, 22884.1866, 22776.3058, 21716.1604, & 20073.2963, 18110.5123, 16004.7832, 13877.6253, 11812.5452, & 9865.8840, 8073.9726, 6458.0834, 5027.9899, 3784.6085, & 2722.0086, 1828.9752, 1090.2396, 487.4595, 0.0000 / data b24 / 0.0000000, 0.0000000, 0.0000000, 0.0000000, 0.0000000, & 0.0000000, 0.0435679, 0.1102275, 0.1922249, 0.2817656, & 0.3694997, 0.4532348, 0.5316253, 0.6038733, 0.6695556, & 0.7285176, 0.7808017, 0.8265992, 0.8662148, 0.9000406, & 0.9285364, 0.9522140, 0.9716252, 0.9873523, 1.0000000 / ! High-resolution troposphere setup ! Revised Apr 14, 2004: PINT = 245.027 mb data a32/100.00000, 400.00000, 818.60211, & 1378.88653, 2091.79519, 2983.64084, & 4121.78960, 5579.22148, 7419.79300, & 9704.82578, 12496.33710, 15855.26306, & 19839.62499, 24502.73262, 28177.10152, & 29525.28447, 29016.34358, 27131.32792, & 24406.11225, 21326.04907, 18221.18357, & 15275.14642, 12581.67796, 10181.42843, & 8081.89816, 6270.86956, 4725.35001, & 3417.39199, 2317.75459, 1398.09473, & 632.49506, 0.00000, 0.00000 / data b32/0.00000, 0.00000, 0.00000, & 0.00000, 0.00000, 0.00000, & 0.00000, 0.00000, 0.00000, & 0.00000, 0.00000, 0.00000, & 0.00000, 0.00000, 0.01711, & 0.06479, 0.13730, 0.22693, & 0.32416, 0.42058, 0.51105, & 0.59325, 0.66628, 0.73011, & 0.78516, 0.83217, 0.87197, & 0.90546, 0.93349, 0.95685, & 0.97624, 0.99223, 1.00000 / !--------------------- ! Wilson's 32L settings: !--------------------- ! Top changed to 0.01 mb data a32w/ 1.00, 26.6378, 84.5529, 228.8592, & 539.9597, 1131.7087, 2141.8082, 3712.0454, & 5963.5317, 8974.1873, 12764.5388, 17294.5911, & 20857.7007, 22221.8651, 22892.7202, 22891.1641, & 22286.0724, 21176.0846, 19673.0671, 17889.0989, & 15927.5060, 13877.6239, 11812.5474, 9865.8830, & 8073.9717, 6458.0824, 5027.9893, 3784.6104, & 2722.0093, 1828.9741, 1090.2397, 487.4575, & 0.0000 / data b32w/ 0.0000, 0.0000, 0.0000, 0.0000, & 0.0000, 0.0000, 0.0000, 0.0000, & 0.0000, 0.0000, 0.0000, 0.0000, & 0.0159, 0.0586, 0.1117, 0.1734, & 0.2415, 0.3137, 0.3878, 0.4619, & 0.5344, 0.6039, 0.6696, 0.7285, & 0.7808, 0.8266, 0.8662, 0.9000, & 0.9285, 0.9522, 0.9716, 0.9874, & 1.0000 / ! High trop-resolution (Feb 2004) data a48x/ & 1.00000, 2.69722, 5.17136, & 8.89455, 14.24790, 22.07157, & 33.61283, 50.48096, 74.79993, & 109.40055, 158.00460, 225.44108, & 317.89560, 443.19350, 611.11558, & 833.74392, 1125.83405, 1505.20759, & 1993.15829, 2614.86254, 3399.78420, & 4382.06240, 5600.87014, 7100.73115, & 8931.78242, 11149.97021, 13817.16841, & 17001.20930, 20775.81856, 23967.33875, & 25527.64563, 25671.22552, 24609.29622, & 22640.51220, 20147.13482, 17477.63530, & 14859.86462, 12414.92533, 10201.44191, & 8241.50255, 6534.43202, 5066.17865, & 3815.60705, 2758.60264, 1870.64631, & 1128.33931, 510.47983, 0.00000, & 0.00000 / data b48x/ & 0.00000, 0.00000, 0.00000, & 0.00000, 0.00000, 0.00000, & 0.00000, 0.00000, 0.00000, & 0.00000, 0.00000, 0.00000, & 0.00000, 0.00000, 0.00000, & 0.00000, 0.00000, 0.00000, & 0.00000, 0.00000, 0.00000, & 0.00000, 0.00000, 0.00000, & 0.00000, 0.00000, 0.00000, & 0.00000, 0.00000, 0.01253, & 0.04887, 0.10724, 0.18455, & 0.27461, 0.36914, 0.46103, & 0.54623, 0.62305, 0.69099, & 0.75016, 0.80110, 0.84453, & 0.88127, 0.91217, 0.93803, & 0.95958, 0.97747, 0.99223, & 1.00000 / ! John Wilson 48 stratosphere data a48/ 0.2000, 0.4643, 1.1978, 2.8806, & 6.4919, 13.7758, 27.6477, 52.6975, & 95.7569, 166.4702, 277.7878, 446.2903, & 692.2605, 1039.4396, 1514.4422, 2145.8457, & 2963.0124, 3994.7360, 5267.8238, 6805.7330, & 8627.3676, 10746.1233, 13169.2374, 15897.4710, & 18925.1187, 20785.0282, 21756.8874, 22456.5487, & 22864.2127, 22968.7603, 22767.3528, 22264.7045, & 21472.1263, 20406.4316, 19088.7860, 17543.5673, & 15797.2892, 13877.6239, 11812.5474, 9865.8830, & 8073.9717, 6458.0824, 5027.9893, 3784.6104, & 2722.0093, 1828.9741, 1090.2397, 487.4575, & 0.0000 / data b48/ 0.00000, 0.00000, 0.00000, 0.00000, & 0.00000, 0.00000, 0.00000, 0.00000, & 0.00000, 0.00000, 0.00000, 0.00000, & 0.00000, 0.00000, 0.00000, 0.00000, & 0.00000, 0.00000, 0.00000, 0.00000, & 0.00000, 0.00000, 0.00000, 0.00000, & 0.00000, 0.01436, 0.04016, 0.07108, & 0.10707, 0.14795, 0.19344, 0.24318, & 0.29672, 0.35358, 0.41326, 0.47522, & 0.53893, 0.60387, 0.66956, 0.72852, & 0.78080, 0.82660, 0.86621, 0.90004, & 0.92854, 0.95221, 0.97163, 0.98735, & 1.00000 / data a55/ 1.00000, 2.00000, 3.27000, & 4.75850, 6.60000, 8.93450, & 11.97030, 15.94950, 21.13490, & 27.85260, 36.50410, 47.58060, & 61.67790, 79.51340, 101.94420, & 130.05080, 165.07920, 208.49720, & 262.02120, 327.64330, 407.65670, & 504.68050, 621.68000, 761.98390, & 929.29430, 1127.68880, 1364.33920, & 1645.70720, 1979.15540, 2373.03610, & 2836.78160, 3380.99550, 4017.54170, & 4764.39320, 5638.79380, 6660.33770, & 7851.22980, 9236.56610, 10866.34270, & 12783.70000, 15039.30000, 17693.00000, & 20119.20876, 21686.49129, 22436.28749, & 22388.46844, 21541.75227, 19873.78342, & 17340.31831, 13874.44006, 10167.16551, & 6609.84274, 3546.59643, 1270.49390, & 0.00000, 0.00000 / data b55 /0.00000, 0.00000, 0.00000, & 0.00000, 0.00000, 0.00000, & 0.00000, 0.00000, 0.00000, & 0.00000, 0.00000, 0.00000, & 0.00000, 0.00000, 0.00000, & 0.00000, 0.00000, 0.00000, & 0.00000, 0.00000, 0.00000, & 0.00000, 0.00000, 0.00000, & 0.00000, 0.00000, 0.00000, & 0.00000, 0.00000, 0.00000, & 0.00000, 0.00000, 0.00000, & 0.00000, 0.00000, 0.00000, & 0.00000, 0.00000, 0.00000, & 0.00000, 0.00000, 0.00000, & 0.00696, 0.02801, 0.06372, & 0.11503, 0.18330, 0.27033, & 0.37844, 0.51046, 0.64271, & 0.76492, 0.86783, 0.94329, & 0.98511, 1.00000 / data a60/ & 0.17861, 1.08051, 3.96471, 9.75160, & 19.81658, 36.69595, 62.55057, 99.19962, & 147.92505, 209.47487, 284.22571, 372.41721, & 474.37835, 590.70236, 722.36063, 870.76746, & 1037.81388, 1225.88773, 1437.89246, 1677.27266, & 1948.05064, 2254.87627, 2603.09094, 2998.80592, & 3448.99523, 3961.60289, 4545.66416, 5211.44017, & 5970.56440, 6836.19818, 7823.19060, 8948.23510, & 10230.01066, 11689.28975, 13348.98686, 15234.11106, & 17371.57323, 19789.78458, 22005.56455, 23550.11512, & 24468.58332, 24800.54880, 24582.44507, 23849.99962, & 22640.51974, 20994.73715, 18957.84873, 16579.41323, & 14080.07103, 11753.63092, 9651.69963, 7793.80093, & 6176.90628, 4787.42760, 3605.04975, 2605.98607, & 1766.83282, 1065.61312, 482.26201, 0.00000, & 0.00000 / data b60/ & 0.00000, 0.00000, 0.00000, 0.00000, & 0.00000, 0.00000, 0.00000, 0.00000, & 0.00000, 0.00000, 0.00000, 0.00000, & 0.00000, 0.00000, 0.00000, 0.00000, & 0.00000, 0.00000, 0.00000, 0.00000, & 0.00000, 0.00000, 0.00000, 0.00000, & 0.00000, 0.00000, 0.00000, 0.00000, & 0.00000, 0.00000, 0.00000, 0.00000, & 0.00000, 0.00000, 0.00000, 0.00000, & 0.00000, 0.00000, 0.00506, 0.02008, & 0.04490, 0.07936, 0.12326, 0.17634, & 0.23820, 0.30827, 0.38581, 0.46989, & 0.55393, 0.62958, 0.69642, 0.75458, & 0.80463, 0.84728, 0.88335, 0.91368, & 0.93905, 0.96020, 0.97775, 0.99223, & 1.00000 / data a64/1.00000, 3.90000, 8.70000, & 15.42000, 24.00000, 34.50000, & 47.00000, 61.50000, 78.60000, & 99.13500, 124.12789, 154.63770, & 191.69700, 236.49300, 290.38000, & 354.91000, 431.82303, 523.09300, & 630.92800, 757.79000, 906.45000, & 1079.85000, 1281.00000, 1515.00000, & 1788.00000, 2105.00000, 2470.00000, & 2889.00000, 3362.00000, 3890.00000, & 4475.00000, 5120.00000, 5830.00000, & 6608.00000, 7461.00000, 8395.00000, & 9424.46289, 10574.46880, 11864.80270, & 13312.58890, 14937.03710, 16759.70700, & 18804.78710, 21099.41210, 23674.03710, & 26562.82810, 29804.11720, 32627.31640, & 34245.89840, 34722.28910, 34155.19920, & 32636.50390, 30241.08200, 27101.44920, & 23362.20700, 19317.05270, 15446.17090, & 12197.45210, 9496.39941, 7205.66992, & 5144.64307, 3240.79346, 1518.62134, & 0.00000, 0.00000 / data b64/0.00000, 0.00000, 0.00000, & 0.00000, 0.00000, 0.00000, & 0.00000, 0.00000, 0.00000, & 0.00000, 0.00000, 0.00000, & 0.00000, 0.00000, 0.00000, & 0.00000, 0.00000, 0.00000, & 0.00000, 0.00000, 0.00000, & 0.00000, 0.00000, 0.00000, & 0.00000, 0.00000, 0.00000, & 0.00000, 0.00000, 0.00000, & 0.00000, 0.00000, 0.00000, & 0.00000, 0.00000, 0.00000, & 0.00000, 0.00000, 0.00000, & 0.00000, 0.00000, 0.00000, & 0.00000, 0.00000, 0.00000, & 0.00000, 0.00000, 0.00813, & 0.03224, 0.07128, 0.12445, & 0.19063, 0.26929, 0.35799, & 0.45438, 0.55263, 0.64304, & 0.71703, 0.77754, 0.82827, & 0.87352, 0.91502, 0.95235, & 0.98511, 1.00000 / ! ! Ultra high troposphere resolution data a100/100.00000, 300.00000, 800.00000, & 1762.35235, 3106.43596, 4225.71874, & 4946.40525, 5388.77387, 5708.35540, & 5993.33124, 6277.73673, 6571.49996, & 6877.05339, 7195.14327, 7526.24920, & 7870.82981, 8229.35361, 8602.30193, & 8990.16936, 9393.46399, 9812.70768, & 10248.43625, 10701.19980, 11171.56286, & 11660.10476, 12167.41975, 12694.11735, & 13240.82253, 13808.17600, 14396.83442, & 15007.47066, 15640.77407, 16297.45067, & 16978.22343, 17683.83253, 18415.03554, & 19172.60771, 19957.34218, 20770.05022, & 21559.14829, 22274.03147, 22916.87519, & 23489.70456, 23994.40187, 24432.71365, & 24806.25734, 25116.52754, 25364.90190, & 25552.64670, 25680.92203, 25750.78675, & 25763.20311, 25719.04113, 25619.08274, & 25464.02630, 25254.49482, 24991.06137, & 24674.32737, 24305.11235, 23884.79781, & 23415.77059, 22901.76510, 22347.84738, & 21759.93950, 21144.07284, 20505.73136, & 19849.54271, 19179.31412, 18498.23400, & 17809.06809, 17114.28232, 16416.10343, & 15716.54833, 15017.44246, 14320.43478, & 13627.01116, 12938.50682, 12256.11762, & 11580.91062, 10913.83385, 10255.72526, & 9607.32122, 8969.26427, 8342.11044, & 7726.33606, 7122.34405, 6530.46991, & 5950.98721, 5384.11279, 4830.01153, & 4288.80090, 3760.55514, 3245.30920, & 2743.06250, 2253.78294, 1777.41285, & 1313.88054, 863.12371, 425.13088, & 0.00000, 0.00000 / data b100/0.00000, 0.00000, 0.00000, & 0.00000, 0.00000, 0.00000, & 0.00000, 0.00000, 0.00000, & 0.00000, 0.00000, 0.00000, & 0.00000, 0.00000, 0.00000, & 0.00000, 0.00000, 0.00000, & 0.00000, 0.00000, 0.00000, & 0.00000, 0.00000, 0.00000, & 0.00000, 0.00000, 0.00000, & 0.00000, 0.00000, 0.00000, & 0.00000, 0.00000, 0.00000, & 0.00000, 0.00000, 0.00000, & 0.00000, 0.00000, 0.00000, & 0.00052, 0.00209, 0.00468, & 0.00828, 0.01288, 0.01849, & 0.02508, 0.03266, 0.04121, & 0.05075, 0.06126, 0.07275, & 0.08521, 0.09866, 0.11308, & 0.12850, 0.14490, 0.16230, & 0.18070, 0.20009, 0.22042, & 0.24164, 0.26362, 0.28622, & 0.30926, 0.33258, 0.35605, & 0.37958, 0.40308, 0.42651, & 0.44981, 0.47296, 0.49591, & 0.51862, 0.54109, 0.56327, & 0.58514, 0.60668, 0.62789, & 0.64872, 0.66919, 0.68927, & 0.70895, 0.72822, 0.74709, & 0.76554, 0.78357, 0.80117, & 0.81835, 0.83511, 0.85145, & 0.86736, 0.88286, 0.89794, & 0.91261, 0.92687, 0.94073, & 0.95419, 0.96726, 0.97994, & 0.99223, 1.00000 / select case (km) case (24) ks = 5 do k=1,km+1 ak(k) = a24(k) bk(k) = b24(k) enddo #ifdef MARS_GCM case (28) ks = 4 do k=1,km+1 ak(k) = a28(k) bk(k) = b28(k) enddo case (36) ks = 10 do k=1,km+1 ak(k) = a36(k) bk(k) = b36(k) enddo #endif MARS_GCM case (32) ! ks = 11 ! John Wilson's setup ks = 13 ! high-res trop_32 setup do k=1,km+1 ak(k) = a32(k) bk(k) = b32(k) enddo case (48) ks = 28 do k=1,km+1 ak(k) = a48(k) bk(k) = b48(k) enddo case (55) ks = 41 do k=1,km+1 ak(k) = a55(k) bk(k) = b55(k) enddo case (60) ks = 37 do k=1,km+1 ak(k) = a60(k) bk(k) = b60(k) enddo case (64) ks = 46 do k=1,km+1 ak(k) = a64(k) bk(k) = b64(k) enddo case (100) ks = 38 do k=1,km+1 ak(k) = a100(k) bk(k) = b100(k) enddo case (10) !-------------------------------------------------- ! Pure sigma-coordinate with uniform spacing in "z" !-------------------------------------------------0 ! pt = 2000. ! model top pressure (pascal) press(1) = pt press(km+1) = p0 dlnp = (log(p0) - log(pt)) / real(km) lnpe = log(press(km+1)) do k=km,2,-1 lnpe = lnpe - dlnp press(k) = exp(lnpe) enddo ! Search KS ks = 0 do k=1,km if(press(k) >= pc) then ks = k-1 goto 123 endif enddo 123 continue if(ks /= 0) then do k=1,ks ak(k) = press(k) bk(k) = 0. enddo endif pint = press(ks+1) do k=ks+1,km ak(k) = pint*(press(km)-press(k))/(press(km)-pint) bk(k) = (press(k) - ak(k)) / press(km+1) enddo ak(km+1) = 0. bk(km+1) = 1. do k=2,km ! bk(k) = real(k-1) / real(km) ! ak(k) = pt * ( 1. - bk(k) ) enddo ak(km+1) = 0. bk(km+1) = 1. case default !----------------------------------------------------------------- ! Pure sigma-coordinate with uniform spacing in sigma and ptop = 0 !----------------------------------------------------------------- pt = 0. pint = pt ks = 0 ak(1) = 0. bk(1) = 0. do k=2,km ak(k) = 0. bk(k) = real(k-1) / real(km) enddo ak(km+1) = 0. bk(km+1) = 1. !--- fix the top at 1 mb ak(1) = 100. !--- ! do k=1,km+1 ! write(*,*) k, ak(k), bk(k) ! enddo end select ptop = ak(1) end subroutine set_eta subroutine check_eta(km, ks, ptop, ak, bk) use fv_pack, only: master use fms_mod, only: error_mesg, FATAL implicit none integer, intent(in):: km ! vertical dimension integer, intent(in):: ks ! number of pure p layers real, intent(in):: ak(km+1) real, intent(in):: bk(km+1) real, intent(in):: ptop ! model top (Pa) ! local real dlp, ps_max integer k ps_max = 0. do k=ks+1,km dlp = (ak(k)-ak(k+1))/(bk(k+1)-bk(k)) ! if ( master ) write(*,*) k, dlp*0.01 ps_max = max( ps_max, dlp ) enddo if ( master ) write(*,*) 'Check_eta: allowable min_ps (mb)=', ps_max*0.01 if( ps_max > 450.E2 ) then call error_mesg('check_eta', 'invalid (a,b) in set_eta',FATAL) endif end subroutine check_eta