! +-======-+ ! Copyright (c) 2003-2007 United States Government as represented by ! the Admistrator of the National Aeronautics and Space Administration. ! All Rights Reserved. ! ! THIS OPEN SOURCE AGREEMENT ("AGREEMENT") DEFINES THE RIGHTS OF USE, ! REPRODUCTION, DISTRIBUTION, MODIFICATION AND REDISTRIBUTION OF CERTAIN ! COMPUTER SOFTWARE ORIGINALLY RELEASED BY THE UNITED STATES GOVERNMENT AS ! REPRESENTED BY THE GOVERNMENT AGENCY LISTED BELOW ("GOVERNMENT AGENCY"). ! THE UNITED STATES GOVERNMENT, AS REPRESENTED BY GOVERNMENT AGENCY, IS AN ! INTENDED THIRD-PARTY BENEFICIARY OF ALL SUBSEQUENT DISTRIBUTIONS OR ! REDISTRIBUTIONS OF THE SUBJECT SOFTWARE. ANYONE WHO USES, REPRODUCES, ! DISTRIBUTES, MODIFIES OR REDISTRIBUTES THE SUBJECT SOFTWARE, AS DEFINED ! HEREIN, OR ANY PART THEREOF, IS, BY THAT ACTION, ACCEPTING IN FULL THE ! RESPONSIBILITIES AND OBLIGATIONS CONTAINED IN THIS AGREEMENT. ! ! Government Agency: National Aeronautics and Space Administration ! Government Agency Original Software Designation: GSC-15354-1 ! Government Agency Original Software Title: GEOS-5 GCM Modeling Software ! User Registration Requested. Please Visit http://opensource.gsfc.nasa.gov ! Government Agency Point of Contact for Original Software: ! Dale Hithon, SRA Assistant, (301) 286-2691 ! ! +-======-+ subroutine set_eta ( km,ptop,ak,bk ) implicit none integer k, km real ak(km+1),bk(km+1) real ptop ! model top (Pa) real a20(21), b20(21) real a32(33), b32(33) real a55(56), b55(56) real a72(73), b72(73) real a91(92), b91(92) data a20 / & 737.44866, 320.17959, 186.17679, & 135.92544, 96.25746, 87.67601, & 61.47969, 66.97722, 42.58292, & 56.11863, 30.33517, 49.84043, & 21.48820, 46.05084, 14.61322, & 43.75791, 8.98483, 42.43210, & 4.19583, 41.76551, 0.00000 / data b20 / & 0.00000, 0.05000, 0.10000, & 0.15000, 0.20000, 0.25000, & 0.30000, 0.35000, 0.40000, & 0.45000, 0.50000, 0.55000, & 0.60000, 0.65000, 0.70000, & 0.75000, 0.80000, 0.85000, & 0.90000, 0.95000, 1.00000 / data a32/40.00000, 100.00000, 200.00000, & 370.00000, 630.00000, 1000.00000, & 1510.00000, 2160.00000, 2900.00000, & 3680.00000, 4535.00000, 5505.00000, & 6607.26750, 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 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.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 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 / ! 72-Level (prior to 5/11/2005) ! ----------------------------- #if 0 data a72 / & 1.00000, 2.00000, 3.27000, 4.75850, 6.60000, 8.93450, & 11.9703, 15.9495, 21.1349, 27.8526, 36.5041, 47.5806, & 61.6779, 79.5134, 101.944, 130.051, 165.079, 208.497, & 262.021, 327.643, 407.657, 504.680, 621.680, 761.984, & 929.294, 1127.69, 1364.34, 1645.71, 1979.16, 2373.04, & 2836.78, 3381.00, 4017.54, 4764.39, 5638.79, 6660.34, & 7851.23, 9236.57, 10866.3, 12783.7, 15039.3, 17693.0, & 20119.2, 21686.5, 22436.3, 22389.8, 21877.6, 21215.0, & 20325.9, 19309.7, 18161.9, 16960.9, 15626.0, 14291.0, & 12869.6, 11409.0, 9936.52, 8909.99, 7883.42, 6856.90, & 5805.32, 5169.61, 4533.90, 3898.20, 3257.08, 2609.20, & 1961.31, 1313.48, 659.375, 332.086, 4.80469, 0.00000, & 0.00000 / data b72 / & 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.00696000, 0.0280100, 0.0637200, 0.113602, 0.156224, 0.200350, & 0.246741, 0.294403, 0.343381, 0.392891, 0.443740, 0.494590, & 0.546304, 0.598410, 0.650635, 0.685900, 0.721166, 0.756431, & 0.791947, 0.813304, 0.834661, 0.856018, 0.877429, 0.898908, & 0.920387, 0.941865, 0.963406, 0.974179, 0.984952, 0.992500, & 1.00000 / #endif ! 72-Level (as of 5/11/2005) ! -------------------------- #if 1 data a72 / & 1.0000000, 2.0000002, 3.2700005, 4.7585009, 6.6000011, & 8.9345014, 11.970302, 15.949503, 21.134903, 27.852606, & 36.504108, 47.580610, 61.677911, 79.513413, 101.94402, & 130.05102, 165.07903, 208.49704, 262.02105, 327.64307, & 407.65710, 504.68010, 621.68012, 761.98417, 929.29420, & 1127.6902, 1364.3402, 1645.7103, 1979.1604, 2373.0405, & 2836.7806, 3381.0007, 4017.5409, 4764.3911, 5638.7912, & 6660.3412, 7851.2316, 9236.5722, 10866.302, 12783.703, & 15039.303, 17693.003, 20119.201, 21686.501, 22436.301, & 22389.800, 21877.598, 21214.998, 20325.898, 19309.696, & 18161.897, 16960.896, 15625.996, 14290.995, 12869.594, & 11895.862, 10918.171, 9936.5219, 8909.9925, 7883.4220, & 7062.1982, 6436.2637, 5805.3211, 5169.6110, 4533.9010, & 3898.2009, 3257.0809, 2609.2006, 1961.3106, 1313.4804, & 659.37527, 4.8048257, 0.0000000 / data b72 / & 0.0000000, 0.0000000, 0.0000000, 0.0000000, 0.0000000, & 0.0000000, 0.0000000, 0.0000000, 0.0000000, 0.0000000, & 0.0000000, 0.0000000, 0.0000000, 0.0000000, 0.0000000, & 0.0000000, 0.0000000, 0.0000000, 0.0000000, 0.0000000, & 0.0000000, 0.0000000, 0.0000000, 0.0000000, 0.0000000, & 0.0000000, 0.0000000, 0.0000000, 0.0000000, 0.0000000, & 0.0000000, 0.0000000, 0.0000000, 0.0000000, 0.0000000, & 0.0000000, 0.0000000, 0.0000000, 0.0000000, 0.0000000, & 0.0000000, 8.1754130e-09, 0.0069600246, 0.028010041, 0.063720063, & 0.11360208, 0.15622409, 0.20035011, 0.24674112, 0.29440312, & 0.34338113, 0.39289115, 0.44374018, 0.49459020, 0.54630418, & 0.58104151, 0.61581843, 0.65063492, 0.68589990, 0.72116594, & 0.74937819, 0.77063753, 0.79194696, 0.81330397, 0.83466097, & 0.85601798, 0.87742898, 0.89890800, 0.92038701, 0.94186501, & 0.96340602, 0.98495195, 1.0000000 / #endif ! 72-Level (as of 3/28/2007) Thin Layer Values ! -------------------------------------------- #if 0 data a72 / & 1.0000000, 2.0000002, 3.2700005, 4.7585009, 6.6000011, & 8.9345014, 11.970302, 15.949503, 21.134903, 27.852606, & 36.504108, 47.580610, 61.677911, 79.513413, 101.94402, & 130.05102, 165.07903, 208.49704, 262.02105, 327.64307, & 407.65710, 504.68010, 621.68012, 761.98417, 929.29420, & 1127.6902, 1364.3402, 1645.7103, 1979.1604, 2373.0405, & 2836.7806, 3381.0007, 4017.5409, 4764.3911, 5638.7912, & 6660.3412, 7851.2316, 9236.5722, 10866.302, 12783.703, & 15039.303, 17693.0, 20172.2, 21733.4, 22432.1, 22390.9, 21806.1, 21020.7, & 20020.2, 18903.0, 17708.5, 16469.8, 15196.1, 13955.8, 12715.6, 11513.2, 10378.1, & 9286.23, 8250.31, 7287.11, 6379.95, 5522.36, 4715.73, 3957.50, 3236.60, 2549.27, & 1904.52, 1307.47, 761.539, 282.945, 3.90625, 1.61719, 0.00000 / data b72 / & 0.0000000, 0.0000000, 0.0000000, 0.0000000, 0.0000000, & 0.0000000, 0.0000000, 0.0000000, 0.0000000, 0.0000000, & 0.0000000, 0.0000000, 0.0000000, 0.0000000, 0.0000000, & 0.0000000, 0.0000000, 0.0000000, 0.0000000, 0.0000000, & 0.0000000, 0.0000000, 0.0000000, 0.0000000, 0.0000000, & 0.0000000, 0.0000000, 0.0000000, 0.0000000, 0.0000000, & 0.0000000, 0.0000000, 0.0000000, 0.0000000, 0.0000000, & 0.0000000, 0.0000000, 0.0000000, 0.0000000, 0.0000000, & 0.0000000, 8.17541e-09, 0.00767142, 0.0302422, 0.0681799, & 0.112411 , 0.160985 , 0.210491, 0.261081, 0.311757, & 0.362072, 0.411599, 0.460114, 0.506786, 0.551798, & 0.594652, 0.634973, 0.672975, 0.708562, 0.741651, 0.772540, 0.801453, 0.828552, & 0.854026, 0.878108, 0.900895, 0.922270, 0.942063, 0.960042, 0.975797, 0.987755, & 0.994941, 1.000000 / #endif ! 91-Level from ECMWF (7/8/08) ! ---------------------------- data a91 / & 1.000000, & 2.000040, & 3.980832, & 7.387186, & 12.908319, & 21.413612, & 33.952858, & 51.746601, & 76.167656, & 108.715561, & 150.986023, & 204.637451, & 271.356506, & 352.824493, & 450.685791, & 566.519226, & 701.813354, & 857.945801, & 1036.166504, & 1237.585449, & 1463.163940, & 1713.709595, & 1989.874390, & 2292.155518, & 2620.898438, & 2976.302246, & 3358.425781, & 3767.196045, & 4202.416504, & 4663.776367, & 5150.859863, & 5663.156250, & 6199.839355, & 6759.727051, & 7341.469727, & 7942.926270, & 8564.624023, & 9208.305664, & 9873.560547, & 10558.881836, & 11262.484375, & 11982.662109, & 12713.897461, & 13453.225586, & 14192.009766, & 14922.685547, & 15638.053711, & 16329.560547, & 16990.623047, & 17613.281250, & 18191.029297, & 18716.968750, & 19184.544922, & 19587.513672, & 19919.796875, & 20175.394531, & 20348.916016, & 20434.158203, & 20426.218750, & 20319.011719, & 20107.031250, & 19785.357422, & 19348.775391, & 18798.822266, & 18141.296875, & 17385.595703, & 16544.585938, & 15633.566406, & 14665.645508, & 13653.219727, & 12608.383789, & 11543.166992, & 10471.310547, & 9405.222656, & 8356.252930, & 7335.164551, & 6353.920898, & 5422.802734, & 4550.215820, & 3743.464355, & 3010.146973, & 2356.202637, & 1784.854614, & 1297.656128, & 895.193542, & 576.314148, & 336.772369, & 162.043427, & 54.208336, & 6.575628, & 0.003160, & 0.000000 / data b91 / & 0.000000, & 0.000000, & 0.000000, & 0.000000, & 0.000000, & 0.000000, & 0.000000, & 0.000000, & 0.000000, & 0.000000, & 0.000000, & 0.000000, & 0.000000, & 0.000000, & 0.000000, & 0.000000, & 0.000000, & 0.000000, & 0.000000, & 0.000000, & 0.000000, & 0.000000, & 0.000000, & 0.000000, & 0.000000, & 0.000000, & 0.000000, & 0.000000, & 0.000000, & 0.000000, & 0.000000, & 0.000000, & 0.000000, & 0.000000, & 0.000000, & 0.000014, & 0.000055, & 0.000131, & 0.000279, & 0.000548, & 0.001000, & 0.001701, & 0.002765, & 0.004267, & 0.006322, & 0.009035, & 0.012508, & 0.016860, & 0.022189, & 0.028610, & 0.036227, & 0.045146, & 0.055474, & 0.067316, & 0.080777, & 0.095964, & 0.112979, & 0.131935, & 0.152934, & 0.176091, & 0.201520, & 0.229315, & 0.259554, & 0.291993, & 0.326329, & 0.362203, & 0.399205, & 0.436906, & 0.475016, & 0.513280, & 0.551458, & 0.589317, & 0.626559, & 0.662934, & 0.698224, & 0.732224, & 0.764679, & 0.795385, & 0.824185, & 0.850950, & 0.875518, & 0.897767, & 0.917651, & 0.935157, & 0.950274, & 0.963007, & 0.973466, & 0.982238, & 0.989153, & 0.994204, & 0.997630, & 1.000000 / if( km.ne.91 .and. & km.ne.72 .and. & km.ne.55 .and. & km.ne.32 .and. & km.ne.20 ) then print *, km,' levels not supported' print *, 'Supported levels are: 91,72,70,55,32,20' stop endif if( km.eq.91 ) then do k=1,km+1 ak(k) = a91(k) bk(k) = b91(k) enddo endif if( km.eq.72 ) then do k=1,km+1 ak(k) = a72(k) bk(k) = b72(k) enddo endif if( km.eq.55 ) then do k=1,km+1 ak(k) = a55(k) bk(k) = b55(k) enddo endif if( km.eq.32 ) then do k=1,km+1 ak(k) = a32(k) bk(k) = b32(k) enddo endif if( km.eq.20 ) then do k=1,km+1 ak(k) = a20(k) bk(k) = b20(k) enddo endif ptop = ak(1) return end