pro getprf21,prfnam,oz,temp ; ; get prf ; ; author/date ; ; colin seftor / 22-february-96 ; ; purpose: ; ; get a standard profile ; ; category: ; ; radiative transfer routine ; ; calling sequence: ; ; pro getprf,prfnam,oz,prf ; ; input: ; ; prfnam: profile name ; ; output ; ; oz: ozone profile; real array of 11,ground to top of atm ; ; temp: temperature profile; real array of 11,ground to top of atm ; ; ; US standard atmosphere used in version 8 for all profiles ; temp = [269.7,236.4,216.7,216.7,217.9,222.4,226.7,237.0,250.8,265.9,270.6] ; case prfnam of 'l225': begin oz = [ 15.00, 9.00, 4.94, 5.91,25.75,62.60,57.00,29.40,10.90, 3.20, 1.30] end 'l275': begin oz = [ 15.00, 9.00, 5.61,14.66,50.68,78.25,57.00,29.40,10.90, 3.20, 1.30] end 'l325': begin oz = [ 15.00, 9.00,10.98,30.51,71.05,86.66,57.00,29.40,10.90, 3.20, 1.30] end 'm225': begin oz = [ 15.09, 9.03, 8.68,18.55,45.11,53.84,36.26,22.56,10.75, 3.70, 1.43] end 'm275': begin oz = [ 15.09,11.46,15.41,29.67,61.57,63.64,39.72,22.56,10.75, 3.70, 1.43] end 'm325': begin oz = [ 15.09,14.08,25.51,43.64,75.49,70.32,42.43,22.56,10.75, 3.70, 1.43] end 'm375': begin oz = [ 15.09,16.91,39.73,61.73,87.09,73.51,42.50,22.56,10.75, 3.70, 1.43] end 'm425': begin oz = [ 15.09,19.93,56.13,82.00,96.26,74.25,42.90,22.56,10.75, 3.70, 1.43] end 'm475': begin oz = [ 15.09,23.15,73.24,103.71,104.09,74.28,43.00,22.56,10.75, 3.70, 1.43] end 'm525': begin oz = [ 15.09,26.57,91.54,126.67,110.15,73.54,43.00,22.56,10.75, 3.70, 1.43] end 'm575': begin oz = [ 15.09,30.19,111.11,150.11,114.67,72.39,43.00,22.56,10.75, 3.70, 1.43] end 'h125': begin oz = [ 14.11, 8.47,18.56, 5.91, 5.00,23.28,22.94,15.87, 7.37, 2.53, 0.96] end 'h175': begin oz = [ 14.11, 8.73,18.53,16.91,23.61,33.26,30.08,17.91, 8.06, 2.73, 1.07] end 'h225': begin oz = [ 14.11, 9.89,21.62,28.94,41.27,41.10,35.61,19.65, 8.71, 2.93, 1.17] end 'h275': begin oz = [ 14.11,11.85,27.72,42.14,57.09,47.89,39.54,21.02, 9.28, 3.11, 1.25] end 'h325': begin oz = [ 14.11,14.61,36.84,56.58,71.10,53.60,41.88,21.98, 9.72, 3.27, 1.31] end 'h375': begin oz = [ 14.11,17.80,49.19,72.51,83.51,58.34,42.33,22.35,10.09, 3.41, 1.36] end 'h425': begin oz = [ 14.11,22.50,64.50,89.46,93.69,61.65,41.70,22.40,10.08, 3.51, 1.40] end 'h475': begin oz = [ 14.11,29.41,82.32,106.83,101.07,63.21,40.80,22.34, 9.95, 3.54, 1.42] end 'h525': begin oz = [ 14.11,38.08,102.60,124.66,105.87,63.22,39.45,22.17, 9.90, 3.52, 1.42] end 'h575': begin oz = [ 14.11,48.45,125.21,142.90,108.23,61.79,37.65,21.90, 9.88, 3.47, 1.41] end else: begin print,'v8 profile name not recognized' stop end endcase return end