module pepdgrid_init_module use params_module,only: nlat,nlon,nlonp2,nlonp4,nlevp1 use init_module,only: sin_sundec,cos_sundec use addfld_module,only: addfld implicit none ! ------------------------------------------------------------- ! Routine to expand peh(ISZA,nlevp1) to PEPDGRID(nlonp4,nlat,nlevp1) ! Initial array pesza(isza,nlevp1) ! Initial array pdsza(isza,nlevp1) ! Interpolation array pepdgrid(nlonp4,nlat,nlevp1) ! Transfer pepdgrid(nlonp4,nlat,nlevp1) to module data ! ------------------------------------------------------------- ! Array Dimensions for Initial Readin Arrays ! Summary; nlevp1 = 47; isza = 9. integer,parameter :: isza = 9 real :: pesza(isza,nlevp1) real :: pdsza(isza,nlevp1) real :: szah(isza) ! !----------------------------------------- ! For interpolated photoelectron array structure real,dimension(nlonp4,nlat,nlevp1),save :: pepdgrid real,dimension(nlonp4,nlat,nlevp1),save :: pdgrid real,dimension(nlonp4,nlat,nlevp1),save :: pegrid real,dimension(nlevp1,nlonp4,nlat),save :: pepdgrid2 ! ! Final pepdgrid photoelectron frequency array structure ! These are allocated for task subdomains by alloc_q (called from allocdata) ! * pepdgrid(nlonp4,nlat,nlevp1) ! pe on the vtgcm grid (#/cm3/sec) ! * Not used. ! real,allocatable :: pepdgrid(:,:,:) !------------------------------------------ ! contains !----------------------------------------------------------------------- subroutine pepdgrid_init ! ! ---------------------------------------------------------------------- ! DRIVER TO EXPAND PESZA(ISZA,nlevp1) to PEPDGRID(nlonp4,nlat,nlevp1) ! Values are provided by Gerard for smin only ! Conversion from (/s) to (#/cm3/sec) ! ---------------------------------------------------------------------- ! ! Interpolate pesza and pdsza to vtgcm LAT-LON grid ! * Initial array persza(isza,nlevp1) and pdrsza(isza,nlevp1) ! * Interpolation array pe(nlonp4,nlat,nlevp1) ! Add 4-endpoints to generate final array for allocation and subdomain usage ! * Final array pepdgrid(nlonp4,nlat,nlevp1) transfer to module data ! Domain Summary: ! nlevp1 = 47; nlat = 36; nlon = 72; nlonp4 = 76; zbot = -16.0; ztop = +7.0 ! ---------------------------------------------------------------------- ! use params_module,only: nlat,nlon,nlonp2,nlonp4,nlevp1,dz,zibot use cons_module,only: dphi,dlamda,pi use addfld_module,only: addfld implicit none ! ! Local arrays; real :: szp(nlevp1) real :: dchi(isza) ! ! Local: integer :: k,i,il,j real :: glat,zchi,slt,rat,coslat,sinlat #include ! ! data statements ! -------------------------------------------------------------------------- ! One sza column per data statement so that updates can be made later ! -- nlevp1 = 47; isza = 9 ! 12/7/08 btf: nlevp1 = 69 (-16 to 18 by 0.5): pad pesze w/ 0 at the top ! -------------------------------------------------------------------------- ! szah(:) = (/0.,10.,20.,30.,40.,50.,60.,70.,90./) ! #if (NLEV==68) ! ! Photoelectron Data ! pesza(1,:) = (/ | 0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0., | 0.,0.,0.,7.5276e-016,3.29e-015,9.7215e-015, | 2.3132e-12,1.041e-11,1.6999e-11,2.2095e-11, | 2.532e-11,4.47e-11,3.859e-10,1.0334e-09, | 1.1822e-09,9.2761e-10,2.8569e-09,1.0543e-08, | 2.7754e-08,6.5523e-08,1.1626e-07,1.5237e-07, | 1.7631e-07,1.7733e-07,1.6701e-07,1.52e-07, | 1.3811e-07,1.2088e-07,1.0624e-07,9.3903e-08,8.8841e-08, | 8.4078e-08,8.4078e-08,8.4078e-08,8.4078e-08,8.4078e-08, | 8.4078e-08,8.4078e-08,8.4078e-08,8.4078e-08,8.4078e-08, | 8.4078e-08,8.4078e-08,8.4078e-08,8.4078e-08,8.4078e-08, | 8.4078e-08,8.4078e-08,8.4078e-08,8.4078e-08,8.4078e-08, | 8.4078e-08,8.4078e-08,8.4078e-08/) pesza(2,:) = (/ | 0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0., | 0.,0.,0.,0.,0.,2.1325e-15,5.7449e-13, | 2.6895e-12,9.8427e-12,2.0318e-11,3.373e-11, | 6.5811e-11,2.8691e-10,7.4529e-10,5.0709e-10, | 8.2079e-10,2.9329e-09,9.9435e-09,2.7943e-08, | 6.3884e-08,9.8553e-08,1.4894e-07,1.7165e-07, | 1.7313e-07,1.6641e-07,1.5375e-07,1.4384e-07, | 1.2367e-07,1.1435e-07,1.0236e-07,9.6773e-08, | 8.806e-08,8.806e-08,8.806e-08,8.806e-08,8.806e-08, | 8.806e-08,8.806e-08,8.806e-08,8.806e-08,8.806e-08, | 8.806e-08,8.806e-08,8.806e-08,8.806e-08,8.806e-08, | 8.806e-08,8.806e-08,8.806e-08,8.806e-08,8.806e-08, | 8.806e-08,8.806e-08,8.806e-08/) pesza(3,:) = (/ | 0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0., | 0.,0.,1.2976e-24,3.305e-24,1.113e-24, | 2.6165e-15,4.5267e-14,9.204e-13,9.4818e-13, | 1.7603e-12,3.0341e-11,8.6482e-11,3.2037e-10, | 1.3949e-09,9.0132e-10,6.3572e-10,1.9758e-09, | 9.7515e-09,3.209e-08,6.5674e-08,1.1063e-07, | 1.5071e-07,1.7133e-07,1.7426e-07,1.6665e-07, | 1.6161e-07,1.4485e-07,1.2983e-07,1.1371e-07, | 1.0159e-07,9.2586e-08, | 9.0921e-08,9.0921e-08,9.0921e-08,9.0921e-08,9.0921e-08, | 9.0921e-08,9.0921e-08,9.0921e-08,9.0921e-08,9.0921e-08, | 9.0921e-08,9.0921e-08,9.0921e-08,9.0921e-08,9.0921e-08, | 9.0921e-08,9.0921e-08,9.0921e-08,9.0921e-08,9.0921e-08, | 9.0921e-08,9.0921e-08,9.0921e-08/) pesza(4,:) = (/ | 0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0., | 0.,0.,3.2649e-24,1.5778e-23,8.4891e-19, | 1.4985e-18,3.1677e-19,1.7272e-12,5.0202e-12, | 1.3299e-11,6.7626e-11,1.363e-10,1.9119e-10, | 7.122e-10,9.2244e-10,2.9902e-10,2.3486e-09, | 1.1019e-08,3.3333e-08,7.0864e-08,1.161e-07, | 1.529e-07,1.7746e-07,1.745e-07,1.6557e-07, | 1.5372e-07,1.3299e-07,1.1668e-07,1.0549e-07, | 9.5492e-08,9.0883e-08, | 8.6581e-08,8.6581e-08,8.6581e-08,8.6581e-08,8.6581e-08, | 8.6581e-08,8.6581e-08,8.6581e-08,8.6581e-08,8.6581e-08, | 8.6581e-08,8.6581e-08,8.6581e-08,8.6581e-08,8.6581e-08, | 8.6581e-08,8.6581e-08,8.6581e-08,8.6581e-08,8.6581e-08, | 8.6581e-08,8.6581e-08,8.6581e-08/) pesza(5,:) = (/ | 0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0., | 0.,0.,0.,1.5912e-20,1.2598e-15,3.9313e-15, | 9.2726e-14,2.3543e-13,4.3615e-13,1.3321e-11, | 1.9876e-11,1.4403e-11,1.4543e-10,7.7253e-10, | 1.8308e-09,7.027e-10,2.2107e-09,9.0173e-09, | 2.7308e-08,6.5044e-08,1.1288e-07,1.5165e-07, | 1.6353e-07,1.7939e-07,1.7048e-07,1.6051e-07, | 1.4103e-07,1.26e-07,1.1395e-07,1.0962e-07,9.948e-08, | 8.9983e-08,8.9983e-08,8.9983e-08,8.9983e-08,8.9983e-08, | 8.9983e-08,8.9983e-08,8.9983e-08,8.9983e-08,8.9983e-08, | 8.9983e-08,8.9983e-08,8.9983e-08,8.9983e-08,8.9983e-08, | 8.9983e-08,8.9983e-08,8.9983e-08,8.9983e-08,8.9983e-08, | 8.9983e-08,8.9983e-08,8.9983e-08/) pesza(6,:) = (/ | 0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0., | 0.,0.,0.,5.2695e-20,7.944e-18,4.2896e-17, | 1.7472e-14,2.3054e-13,1.2745e-12,9.4234e-12, | 2.0082e-11,2.2314e-11,1.0697e-10,2.7936e-10, | 5.7214e-10,3.6224e-10,1.3e-09,5.0624e-09, | 1.9561e-08,5.1177e-08,9.8231e-08,1.3852e-07, | 1.7109e-07,1.7099e-07,1.6342e-07,1.5386e-07, | 1.3478e-07,1.2169e-07,1.066e-07,9.7504e-08,9.2852e-08, | 8.5535e-08,8.5535e-08,8.5535e-08,8.5535e-08,8.5535e-08, | 8.5535e-08,8.5535e-08,8.5535e-08,8.5535e-08,8.5535e-08, | 8.5535e-08,8.5535e-08,8.5535e-08,8.5535e-08,8.5535e-08, | 8.5535e-08,8.5535e-08,8.5535e-08,8.5535e-08,8.5535e-08, | 8.5535e-08,8.5535e-08,8.5535e-08/) pesza(7,:) = (/ | 0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0., | 0.,0.,0.,0.,2.9036e-019,3.1325e-18,3.9486e-13, | 1.1763e-12,2.1586e-12,4.1632e-12,7.4839e-12, | 1.5982e-11,7.3219e-11,1.8197e-10,2.6244e-10, | 2.6979e-10,6.6081e-10,2.59e-09,1.2465e-08,4e-08, | 7.9733e-08,1.1908e-07,1.4763e-07,1.568e-07, | 1.5486e-07,1.4745e-07,1.3313e-07,1.1907e-07, | 1.0912e-07,9.7606e-08,8.9261e-08, | 8.3998e-08,8.3998e-08,8.3998e-08,8.3998e-08,8.3998e-08, | 8.3998e-08,8.3998e-08,8.3998e-08,8.3998e-08,8.3998e-08, | 8.3998e-08,8.3998e-08,8.3998e-08,8.3998e-08,8.3998e-08, | 8.3998e-08,8.3998e-08,8.3998e-08,8.3998e-08,8.3998e-08, | 8.3998e-08,8.3998e-08,8.3998e-08/) pesza(8,:) = (/ | 0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0., | 0.,0.,0.,0.,0.,1.6486e-16,2.6044e-15,7.3132e-15, | 1.4981e-14,7.2914e-12,1.4341e-11,2.0083e-11, | 2.3526e-11,3.0692e-11,5.1398e-11,5.299e-10, | 1.0515e-09,1.1856e-09,9.5931e-10,1.9591e-09, | 6.3678e-09,2.1258e-08,5.2228e-08,9.9855e-08, | 1.4355e-07,1.7308e-07,1.7712e-07,1.636e-07, | 1.4557e-07,1.2597e-07,1.0775e-07,9.292e-08, | 8.7816e-08,8.7816e-08,8.7816e-08,8.7816e-08,8.7816e-08, | 8.7816e-08,8.7816e-08,8.7816e-08,8.7816e-08,8.7816e-08, | 8.7816e-08,8.7816e-08,8.7816e-08,8.7816e-08,8.7816e-08, | 8.7816e-08,8.7816e-08,8.7816e-08,8.7816e-08,8.7816e-08, | 8.7816e-08,8.7816e-08,8.7816e-08/) pesza(9,:) = (/ | 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., | 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., | 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., | 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., | 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., | 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., | 0.,0.,0., 0., 0., 0., 0., 0., 0./) ! ! Photodissociation data ! pdsza(1,:) = (/ | 0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0., | 0.,0.,0.,0.,1.2005e-32,5.2466e-32,6.0474e-30, | 3.5796e-29,1.1981e-28,3.3718e-28,-7.8015e-20, | -3.0572e-19,-2.2208e-15,1.6298e-12,1.993e-11, | 2.8876e-10,2.1604e-09,8.7102e-09,2.726e-08, | 6.7747e-08,1.5046e-07,2.4749e-07,3.4559e-07, | 4.1978e-07,4.7383e-07,5.0788e-07,5.2965e-07, | 5.4315e-07,5.5042e-07,5.5469e-07,5.5695e-07,5.5821e-07, | 5.5879e-07,5.5879e-07,5.5879e-07,5.5879e-07,5.5879e-07, | 5.5879e-07,5.5879e-07,5.5879e-07,5.5879e-07,5.5879e-07, | 5.5879e-07,5.5879e-07,5.5879e-07,5.5879e-07,5.5879e-07, | 5.5879e-07,5.5879e-07,5.5879e-07,5.5879e-07,5.5879e-07, | 5.5879e-07,5.5879e-07,5.5879e-07/) pdsza(2,:) = (/ | 0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0., | 0.,0.,0.,0.,1.7805e-32,8.1188e-32,7.4145e-30, | 3.9935e-29,1.3247e-28,3.5881e-28,-2.0342e-19, | -8.5646e-19,-2.9447e-15,2.5332e-12,8.3373e-12, | 3.6396e-10,2.4104e-09,9.019e-09,2.711e-08, | 6.9076e-08,1.5105e-07,2.4606e-07,3.407e-07, | 4.1502e-07,4.7064e-07,5.0607e-07,5.2832e-07, | 5.4175e-07,5.4986e-07,5.5423e-07,5.5676e-07,5.5808e-07, | 5.5877e-07,5.5877e-07,5.5877e-07,5.5877e-07,5.5877e-07, | 5.5877e-07,5.5877e-07,5.5877e-07,5.5877e-07,5.5877e-07, | 5.5877e-07,5.5877e-07,5.5877e-07,5.5877e-07,5.5877e-07, | 5.5877e-07,5.5877e-07,5.5877e-07,5.5877e-07,5.5877e-07, | 5.5877e-07,5.5877e-07,5.5877e-07/) pdsza(3,:) = (/ | 0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0., | 0.,0.,0.,6.3741e-36,1.3359e-32,1.8536e-31, | 8.528e-30,2.9519e-29,1.4459e-28,3.5676e-28, | -3.5121e-19,-6.424e-18,1.8138e-16,2.7926e-12, | 2.0881e-11,4.7887e-10,2.4714e-09,6.9385e-09, | 2.6485e-08,7.2694e-08,1.4855e-07,2.4235e-07, | 3.3416e-07,4.1048e-07,4.6577e-07,5.0281e-07, | 5.2696e-07,5.4072e-07,5.4914e-07,5.5391e-07, | 5.5658e-07,5.58e-07, | 5.5872e-07,5.5872e-07,5.5872e-07,5.5872e-07,5.5872e-07, | 5.5872e-07,5.5872e-07,5.5872e-07,5.5872e-07,5.5872e-07, | 5.5872e-07,5.5872e-07,5.5872e-07,5.5872e-07,5.5872e-07, | 5.5872e-07,5.5872e-07,5.5872e-07,5.5872e-07,5.5872e-07, | 5.5872e-07,5.5872e-07,5.5872e-07/) pdsza(4,:) = (/ | 0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0., | 0.,0.,0.,1.9388e-35,9.3694e-35,3.2293e-31, | 5.4738e-30,3.0612e-29,1.3757e-28,3.0892e-28, | 6.6423e-28,-2.5165e-17,9.168e-15,5.9812e-14, | 3.3887e-11,4.006e-10,1.5077e-09,7.2677e-09, | 2.6242e-08,6.8198e-08,1.3802e-07,2.3494e-07, | 3.2724e-07,4.065e-07,4.6181e-07,5.0013e-07, | 5.247e-07,5.3995e-07,5.4859e-07,5.5374e-07, | 5.5642e-07,5.5789e-07, | 5.5869e-07,5.5869e-07,5.5869e-07,5.5869e-07,5.5869e-07, | 5.5869e-07,5.5869e-07,5.5869e-07,5.5869e-07,5.5869e-07, | 5.5869e-07,5.5869e-07,5.5869e-07,5.5869e-07,5.5869e-07, | 5.5869e-07,5.5869e-07,5.5869e-07,5.5869e-07,5.5869e-07, | 5.5869e-07,5.5869e-07,5.5869e-07/) pdsza(5,:) = (/ | 0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0., | 0.,0.,0.,0.,1.1834e-34,2.4109e-31,7.5197e-31, | 2.4901e-29,1.0006e-28,2.4089e-28,5.6792e-28, | -1.816e-17,-6.167e-17,5.4254e-14,2.2683e-11, | 9.5203e-11,1.1772e-09,5.8041e-09,1.9366e-08, | 5.2663e-08,1.1722e-07,2.1439e-07,3.0989e-07, | 3.9515e-07,4.5455e-07,4.9595e-07,5.2226e-07, | 5.3866e-07,5.4797e-07,5.5336e-07,5.5633e-07,5.5757e-07, | 5.5864e-07,5.5864e-07,5.5864e-07,5.5864e-07,5.5864e-07, | 5.5864e-07,5.5864e-07,5.5864e-07,5.5864e-07,5.5864e-07, | 5.5864e-07,5.5864e-07,5.5864e-07,5.5864e-07,5.5864e-07, | 5.5864e-07,5.5864e-07,5.5864e-07,5.5864e-07,5.5864e-07, | 5.5864e-07,5.5864e-07,5.5864e-07/) pdsza(6,:) = (/ | 0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0., | 0.,0.,0.,0.,1.4129e-35,4.7538e-32,2.5955e-31, | 1.1763e-29,4.781e-29,1.6457e-28,4.0461e-28, | -1.3357e-18,-8.3977e-18,2.1155e-17,5.1639e-12, | 2.424e-11,5.7953e-10,3.1045e-09,9.7158e-09, | 3.1874e-08,8.7454e-08,1.7588e-07,2.7495e-07, | 3.6843e-07,4.3868e-07,4.8679e-07,5.1726e-07, | 5.3631e-07,5.468e-07,5.5283e-07,5.5604e-07,5.5768e-07, | 5.5858e-07,5.5858e-07,5.5858e-07,5.5858e-07,5.5858e-07, | 5.5858e-07,5.5858e-07,5.5858e-07,5.5858e-07,5.5858e-07, | 5.5858e-07,5.5858e-07,5.5858e-07,5.5858e-07,5.5858e-07, | 5.5858e-07,5.5858e-07,5.5858e-07,5.5858e-07,5.5858e-07, | 5.5858e-07,5.5858e-07,5.5858e-07/) pdsza(7,:) = (/ | 0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0., | 0.,0.,0.,0.,0.,1.6923e-33,1.8258e-32,2.8384e-30, | 1.1262e-29,8.1597e-29,2.317e-28,-2.7489e-21, | -9.7368e-20,-1.0426e-15,1.9348e-13,3.1363e-12, | 1.651e-10,9.6988e-10,3.6152e-09,1.6361e-08, | 5.1065e-08,1.1757e-07,2.1576e-07,3.1803e-07, | 4.0468e-07,4.6787e-07,5.0675e-07,5.3051e-07, | 5.4427e-07,5.5164e-07,5.5551e-07,5.5743e-07, | 5.5845e-07,5.5845e-07,5.5845e-07,5.5845e-07,5.5845e-07, | 5.5845e-07,5.5845e-07,5.5845e-07,5.5845e-07,5.5845e-07, | 5.5845e-07,5.5845e-07,5.5845e-07,5.5845e-07,5.5845e-07, | 5.5845e-07,5.5845e-07,5.5845e-07,5.5845e-07,5.5845e-07, | 5.5845e-07,5.5845e-07,5.5845e-07/) pdsza(8,:) = (/ | 0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0., | 0.,0.,0.,0.,0.,2.629e-33,4.1533e-32,3.5353e-30, | 1.1534e-29,8.7456e-29,2.2195e-28,4.7087e-28, | -1.7901e-19,-5.7241e-16,-3.0089e-15,2.3432e-12, | 5.2594e-11,2.9495e-10,1.9614e-09,6.5039e-09, | 1.7339e-08,5.2713e-08,1.2135e-07,2.1612e-07, | 3.2068e-07,4.0971e-07,4.7472e-07,5.1366e-07, | 5.3657e-07,5.4849e-07,5.5432e-07,5.5707e-07, | 5.5838e-07,5.5838e-07,5.5838e-07,5.5838e-07,5.5838e-07, | 5.5838e-07,5.5838e-07,5.5838e-07,5.5838e-07,5.5838e-07, | 5.5838e-07,5.5838e-07,5.5838e-07,5.5838e-07,5.5838e-07, | 5.5838e-07,5.5838e-07,5.5838e-07,5.5838e-07,5.5838e-07, | 5.5838e-07,5.5838e-07,5.5838e-07/) pdsza(9,:) = (/ | 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., | 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., | 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., | 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., | 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., | 0., 0., 0., 0., 0., 0., 0., 0., 0., 0., | 0.,0.,0., 0., 0., 0., 0., 0., 0./) ! !************************************************************** ! Double Res. Data !************************************************************** #elif (NLEV==136) ! ! Photoelectron Data ! pesza(1,:) = (/ | 0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0., | 0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0., | 0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,7.53E-16, | 2.02E-15,3.29E-15,5.77E-15,9.72E-15,1.37E-14, | 2.31E-12,6.36E-12,1.04E-11,1.40E-11,1.70E-11, | 2.00E-11,2.21E-11,2.37E-11,2.53E-11,3.35E-11, | 4.47E-11,5.59E-11,3.86E-10,7.25E-10,1.03E-09, | 1.11E-09,1.18E-09,1.11E-09,9.28E-10,7.41E-10, | 2.86E-09,5.25E-09,1.05E-08,1.85E-08,2.78E-08, | 4.66E-08,6.55E-08,9.09E-08,1.16E-07,1.35E-07, | 1.52E-07,1.65E-07,1.76E-07,1.78E-07,1.77E-07, | 1.73E-07,1.67E-07,1.59E-07,1.52E-07,1.45E-07, | 1.38E-07,1.29E-07,1.21E-07,1.13E-07,1.06E-07, | 1.01E-07,9.39E-08,9.07E-08,8.88E-08,8.62E-08, | 8.41E-08,8.39E-08,8.14E-08,7.98E-08,7.84E-08, | 7.65E-08,7.56E-08,7.56E-08,7.56E-08,7.56E-08, | 7.65E-08,7.56E-08,7.56E-08,7.56E-08,7.56E-08, | 7.65E-08,7.56E-08,7.56E-08,7.56E-08,7.56E-08, | 7.65E-08,7.56E-08,7.56E-08,7.56E-08,7.56E-08, | 7.65E-08,7.56E-08,7.56E-08,7.56E-08,7.56E-08, | 7.65E-08,7.56E-08,7.56E-08,7.56E-08,7.56E-08, | 7.65E-08,7.56E-08,7.56E-08,7.56E-08,7.56E-08, | 7.65E-08,7.56E-08,7.56E-08,7.56E-08,7.56E-08/) pesza(2,:) = (/ | 0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0., | 0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0., | 0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0., | 6.21E-16,2.13E-15,3.64E-15,5.74E-13,1.63E-12, | 2.69E-12,5.27E-12,9.84E-12,1.44E-11,2.03E-11, | 2.70E-11,3.37E-11,4.75E-11,6.58E-11,8.42E-11, | 2.87E-10,5.16E-10,7.45E-10,6.29E-10,5.07E-10, | 5.45E-10,8.21E-10,1.10E-09,2.93E-09,5.44E-09, | 9.94E-09,1.89E-08,2.79E-08,4.53E-08,6.39E-08, | 8.14E-08,9.86E-08,1.21E-07,1.49E-07,1.65E-07, | 1.72E-07,1.74E-07,1.73E-07,1.70E-07,1.66E-07, | 1.59E-07,1.54E-07,1.51E-07,1.44E-07,1.34E-07, | 1.24E-07,1.18E-07,1.14E-07,1.08E-07,1.02E-07, | 9.99E-08,9.68E-08,9.22E-08,8.81E-08,8.38E-08, | 8.16E-08,7.99E-08,7.86E-08,7.73E-08,7.66E-08, | 7.64E-08,7.64E-08,7.64E-08,7.64E-08,7.64E-08, | 7.64E-08,7.64E-08,7.64E-08,7.64E-08,7.64E-08, | 7.64E-08,7.64E-08,7.64E-08,7.64E-08,7.64E-08, | 7.64E-08,7.64E-08,7.64E-08,7.64E-08,7.64E-08, | 7.64E-08,7.64E-08,7.64E-08,7.64E-08,7.64E-08, | 7.64E-08,7.64E-08,7.64E-08,7.64E-08,7.64E-08, | 7.64E-08,7.64E-08,7.64E-08,7.64E-08,7.64E-08, | 7.64E-08,7.64E-08,7.64E-08/) pesza(3,:) = (/ | 0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0., | 0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0., | 0.,0.,0.,0.,0.,0.,0.,1.17E-25,1.30E-24, | 2.48E-24,3.30E-24,2.21E-24,1.11E-24,1.71E-26, | 2.62E-15,5.27E-15,4.53E-14,4.83E-13,9.20E-13, | 1.36E-12,9.48E-13,5.02E-13,1.76E-12,1.61E-11, | 3.03E-11,4.63E-11,8.65E-11,1.27E-10,3.20E-10, | 8.58E-10,1.39E-09,1.49E-09,9.01E-10,3.08E-10, | 6.36E-10,1.23E-09,1.98E-09,5.86E-09,9.75E-09, | 1.93E-08,3.21E-08,4.61E-08,6.57E-08,8.54E-08, | 1.11E-07,1.36E-07,1.51E-07,1.64E-07,1.71E-07, | 1.78E-07,1.74E-07,1.70E-07,1.67E-07,1.64E-07, | 1.62E-07,1.53E-07,1.45E-07,1.38E-07,1.30E-07, | 1.21E-07,1.14E-07,1.09E-07,1.02E-07,9.73E-08, | 9.26E-08,8.91E-08,9.09E-08,8.85E-08,8.55E-08, | 8.28E-08,8.10E-08,7.86E-08,7.81E-08,7.71E-08, | 7.71E-08,7.71E-08,7.71E-08,7.71E-08,7.71E-08, | 7.71E-08,7.71E-08,7.71E-08,7.71E-08,7.71E-08, | 7.71E-08,7.71E-08,7.71E-08,7.71E-08,7.71E-08, | 7.71E-08,7.71E-08,7.71E-08,7.71E-08,7.71E-08, | 7.71E-08,7.71E-08,7.71E-08,7.71E-08,7.71E-08, | 7.71E-08,7.71E-08,7.71E-08,7.71E-08,7.71E-08, | 7.71E-08,7.71E-08,7.71E-08,7.71E-08,7.71E-08, | 7.71E-08,7.71E-08/) pesza(4,:) = (/ | 0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0., | 0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0., | 0.,0.,0.,0.,0.,0.,0.,0.,3.26E-24,9.52E-24, | 1.58E-23,3.07E-19,8.49E-19,1.39E-18,1.50E-18, | 9.08E-19,3.17E-19,5.47E-13,1.73E-12,2.91E-12, | 5.02E-12,9.16E-12,1.33E-11,3.18E-11,6.76E-11, | 1.03E-10,1.36E-10,1.64E-10,1.91E-10,3.93E-10, | 7.12E-10,1.03E-09,9.22E-10,5.89E-10,2.99E-10, | 1.32E-09,2.35E-09,4.90E-09,1.10E-08,1.71E-08, | 3.33E-08,5.08E-08,7.09E-08,9.42E-08,1.16E-07, | 1.35E-07,1.53E-07,1.67E-07,1.77E-07,1.77E-07, | 1.74E-07,1.69E-07,1.66E-07,1.64E-07,1.54E-07, | 1.43E-07,1.33E-07,1.25E-07,1.17E-07,1.10E-07, | 1.05E-07,1.01E-07,9.55E-08,9.29E-08,9.09E-08, | 8.85E-08,8.66E-08,8.29E-08,8.13E-08,8.02E-08, | 7.82E-08,7.64E-08,7.39E-08,7.32E-08,7.32E-08, | 7.32E-08,7.32E-08,7.32E-08,7.32E-08,7.32E-08, | 7.32E-08,7.32E-08,7.32E-08,7.32E-08,7.32E-08, | 7.32E-08,7.32E-08,7.32E-08,7.32E-08,7.32E-08, | 7.32E-08,7.32E-08,7.32E-08,7.32E-08,7.32E-08, | 7.32E-08,7.32E-08,7.32E-08,7.32E-08,7.32E-08, | 7.32E-08,7.32E-08,7.32E-08,7.32E-08,7.32E-08, | 7.32E-08,7.32E-08,7.32E-08,7.32E-08,7.32E-08, | 7.32E-08/) pesza(5,:) = (/ | 0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0., | 0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0., | 0.,0.,0.,0.,0.,0.,0.,0.,0.,7.48E-21,1.59E-20, | 2.43E-20,1.26E-15,2.60E-15,3.93E-15,4.25E-14, | 9.27E-14,1.43E-13,2.35E-13,3.36E-13,4.36E-13, | 5.74E-12,1.33E-11,2.09E-11,1.99E-11,1.71E-11, | 1.44E-11,7.08E-11,1.45E-10,2.20E-10,7.73E-10, | 1.33E-09,1.83E-09,1.27E-09,7.03E-10,9.37E-10, | 2.21E-09,3.48E-09,9.02E-09,1.64E-08,2.73E-08, | 4.62E-08,6.50E-08,8.88E-08,1.13E-07,1.33E-07, | 1.52E-07,1.59E-07,1.64E-07,1.72E-07,1.79E-07, | 1.74E-07,1.70E-07,1.69E-07,1.61E-07,1.51E-07, | 1.41E-07,1.33E-07,1.26E-07,1.19E-07,1.14E-07, | 1.09E-07,1.10E-07,1.02E-07,9.95E-08,9.63E-08, | 9.00E-08,8.70E-08,8.44E-08,8.29E-08,8.20E-08, | 8.13E-08,8.00E-08,7.81E-08,7.81E-08,7.81E-08, | 7.81E-08,7.81E-08,7.81E-08,7.81E-08,7.81E-08, | 7.81E-08,7.81E-08,7.81E-08,7.81E-08,7.81E-08, | 7.81E-08,7.81E-08,7.81E-08,7.81E-08,7.81E-08, | 7.81E-08,7.81E-08,7.81E-08,7.81E-08,7.81E-08, | 7.81E-08,7.81E-08,7.81E-08,7.81E-08,7.81E-08, | 7.81E-08,7.81E-08,7.81E-08,7.81E-08,7.81E-08, | 7.81E-08,7.81E-08,7.81E-08,7.81E-08,7.81E-08/) pesza(6,:) = (/ | 0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0., | 0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0., | 0.,0.,0.,0.,0.,0.,0.,0.,0.,1.45E-20,5.27E-20, | 9.09E-20,7.94E-18,2.54E-17,4.29E-17,3.74E-15, | 1.75E-14,3.12E-14,2.31E-13,7.53E-13,1.27E-12, | 3.11E-12,9.42E-12,1.57E-11,2.01E-11,2.12E-11, | 2.23E-11,4.31E-11,1.07E-10,1.71E-10,2.79E-10, | 4.26E-10,5.72E-10,5.29E-10,3.62E-10,1.95E-10, | 1.30E-09,2.45E-09,5.06E-09,1.23E-08,1.96E-08, | 3.47E-08,5.12E-08,7.23E-08,9.82E-08,1.20E-07, | 1.39E-07,1.55E-07,1.71E-07,1.74E-07,1.71E-07, | 1.67E-07,1.63E-07,1.59E-07,1.54E-07,1.46E-07, | 1.35E-07,1.28E-07,1.22E-07,1.14E-07,1.07E-07, | 1.01E-07,9.75E-08,9.52E-08,9.29E-08,8.92E-08, | 8.55E-08,8.28E-08,8.26E-08,8.09E-08,7.94E-08, | 7.79E-08,7.64E-08,7.64E-08,7.64E-08,7.64E-08, | 7.64E-08,7.64E-08,7.64E-08,7.64E-08,7.64E-08, | 7.64E-08,7.64E-08,7.64E-08,7.64E-08,7.64E-08, | 7.64E-08,7.64E-08,7.64E-08,7.64E-08,7.64E-08, | 7.64E-08,7.64E-08,7.64E-08,7.64E-08,7.64E-08, | 7.64E-08,7.64E-08,7.64E-08,7.64E-08,7.64E-08, | 7.64E-08,7.64E-08,7.64E-08,7.64E-08,7.64E-08, | 7.64E-08,7.64E-08,7.64E-08,7.64E-08,7.64E-08/) pesza(7,:) = (/ | 0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0., | 0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0., | 0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0., | 2.90E-19,1.71E-18,3.13E-18,8.87E-15,3.95E-13, | 7.81E-13,1.18E-12,1.67E-12,2.16E-12,2.65E-12, | 4.16E-12,5.75E-12,7.48E-12,1.17E-11,1.60E-11, | 2.02E-11,7.32E-11,1.29E-10,1.82E-10,2.22E-10, | 2.62E-10,2.93E-10,2.70E-10,2.46E-10,6.61E-10, | 1.63E-09,2.59E-09,6.66E-09,1.25E-08,2.09E-08, | 4.00E-08,5.91E-08,7.97E-08,1.01E-07,1.19E-07, | 1.37E-07,1.48E-07,1.57E-07,1.57E-07,1.57E-07, | 1.55E-07,1.52E-07,1.47E-07,1.42E-07,1.33E-07, | 1.26E-07,1.19E-07,1.12E-07,1.09E-07,1.03E-07, | 9.76E-08,9.36E-08,8.93E-08,8.55E-08,8.40E-08, | 8.30E-08,7.99E-08,7.78E-08,7.57E-08,7.45E-08, | 7.37E-08,7.37E-08,7.37E-08,7.37E-08,7.37E-08, | 7.37E-08,7.37E-08,7.37E-08,7.37E-08,7.37E-08, | 7.37E-08,7.37E-08,7.37E-08,7.37E-08,7.37E-08, | 7.37E-08,7.37E-08,7.37E-08,7.37E-08,7.37E-08, | 7.37E-08,7.37E-08,7.37E-08,7.37E-08,7.37E-08, | 7.37E-08,7.37E-08,7.37E-08,7.37E-08,7.37E-08, | 7.37E-08,7.37E-08,7.37E-08,7.37E-08,7.37E-08, | 7.37E-08,7.37E-08,7.37E-08,7.37E-08/) pesza(8,:) = (/ | 0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0., | 0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0., | 0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0., | 1.65E-16,1.38E-15,2.60E-15,3.82E-15,7.31E-15, | 1.11E-14,1.50E-14,3.45E-12,7.29E-12,1.11E-11, | 1.43E-11,1.72E-11,2.01E-11,2.20E-11,2.35E-11, | 2.50E-11,3.07E-11,4.10E-11,5.14E-11,2.22E-10, | 5.30E-10,8.38E-10,1.05E-09,1.12E-09,1.19E-09, | 1.13E-09,9.59E-10,7.93E-10,1.96E-09,4.16E-09, | 6.37E-09,1.38E-08,2.13E-08,3.35E-08,5.22E-08, | 7.33E-08,9.99E-08,1.25E-07,1.44E-07,1.60E-07, | 1.73E-07,1.77E-07,1.77E-07,1.72E-07,1.64E-07, | 1.54E-07,1.46E-07,1.37E-07,1.26E-07,1.17E-07, | 1.08E-07,1.01E-07,9.29E-08,9.01E-08,8.78E-08, | 8.50E-08,8.40E-08,8.20E-08,8.00E-08,7.85E-08, | 7.65E-08,7.56E-08,7.56E-08,7.56E-08,7.56E-08, | 7.56E-08,7.56E-08,7.56E-08,7.56E-08,7.56E-08, | 7.56E-08,7.56E-08,7.56E-08,7.56E-08,7.56E-08, | 7.56E-08,7.56E-08,7.56E-08,7.56E-08,7.56E-08, | 7.56E-08,7.56E-08,7.56E-08,7.56E-08,7.56E-08, | 7.56E-08,7.56E-08,7.56E-08,7.56E-08,7.56E-08, | 7.56E-08,7.56E-08,7.56E-08,7.56E-08,7.56E-08, | 7.56E-08,7.56E-08,7.56E-08,7.56E-08/) pesza(9,:) = (/ | 0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0., | 0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0., | 0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0., | 0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0., | 0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0., | 0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0., | 0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0., | 0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0., | 0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0., | 0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0., | 0.,0.,0.,0.,0.,0.,0./) ! ! Photodissociation data ! pdsza(1,:) = (/ | 0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0., | 0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0., | 0.,0.,0.,0.,0.,0.,0.,0.,0.,0., | 1.20E-32,3.22E-32,5.25E-32,1.93E-30,6.05E-30, | 1.02E-29,3.58E-29,7.78E-29,1.20E-28,2.06E-28, | 3.37E-28,4.69E-28,-7.80E-20,-1.92E-19,-3.06E-19, | -9.02E-16,-2.22E-15,-3.54E-15,1.63E-12,3.31E-12, | 1.99E-11,1.54E-10,2.89E-10,9.89E-10,2.16E-09, | 3.33E-09,8.71E-09,1.46E-08,2.73E-08,4.62E-08, | 6.77E-08,1.09E-07,1.50E-07,1.99E-07,2.47E-07, | 2.96E-07,3.46E-07,3.83E-07,4.20E-07,4.48E-07, | 4.74E-07,4.92E-07,5.08E-07,5.21E-07,5.30E-07, | 5.37E-07,5.43E-07,5.47E-07,5.50E-07,5.53E-07, | 5.55E-07,5.56E-07,5.57E-07,5.58E-07,5.58E-07, | 5.59E-07,5.59E-07,5.59E-07,5.59E-07,5.59E-07, | 5.59E-07,5.59E-07,5.59E-07,5.59E-07,5.59E-07, | 5.59E-07,5.59E-07,5.59E-07,5.59E-07,5.59E-07, | 5.59E-07,5.59E-07,5.59E-07,5.59E-07,5.59E-07, | 5.59E-07,5.59E-07,5.59E-07,5.59E-07,5.59E-07, | 5.59E-07,5.59E-07,5.59E-07,5.59E-07,5.59E-07, | 5.59E-07,5.59E-07,5.59E-07,5.59E-07,5.59E-07, | 5.59E-07,5.59E-07,5.59E-07,5.59E-07,5.59E-07, | 5.59E-07,5.59E-07,5.59E-07,5.59E-07,5.59E-07, | 5.59E-07/) pdsza(2,:) = (/ | 0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0., | 0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0., | 0.,0.,0.,0.,0.,0.,0.,0.,0.,0., | 1.78E-32,4.95E-32,8.12E-32,2.23E-30,7.41E-30, | 1.26E-29,3.99E-29,8.62E-29,1.32E-28,2.19E-28, | 3.59E-28,4.98E-28,-2.03E-19,-5.30E-19,-8.56E-19, | -1.11E-15,-2.94E-15,-4.78E-15,2.53E-12,5.44E-12, | 8.34E-12,1.85E-10,3.64E-10,1.03E-09,2.41E-09, | 3.79E-09,9.02E-09,1.59E-08,2.71E-08,4.81E-08, | 6.91E-08,1.09E-07,1.51E-07,1.98E-07,2.46E-07, | 2.94E-07,3.41E-07,3.80E-07,4.15E-07,4.44E-07, | 4.71E-07,4.89E-07,5.06E-07,5.18E-07,5.28E-07, | 5.36E-07,5.42E-07,5.46E-07,5.50E-07,5.52E-07, | 5.54E-07,5.56E-07,5.57E-07,5.58E-07,5.58E-07, | 5.58E-07,5.59E-07,5.59E-07,5.59E-07,5.59E-07, | 5.59E-07,5.59E-07,5.59E-07,5.59E-07,5.59E-07, | 5.59E-07,5.59E-07,5.59E-07,5.59E-07,5.59E-07, | 5.59E-07,5.59E-07,5.59E-07,5.59E-07,5.59E-07, | 5.59E-07,5.59E-07,5.59E-07,5.59E-07,5.59E-07, | 5.59E-07,5.59E-07,5.59E-07,5.59E-07,5.59E-07, | 5.59E-07,5.59E-07,5.59E-07,5.59E-07,5.59E-07, | 5.59E-07,5.59E-07,5.59E-07,5.59E-07,5.59E-07, | 5.59E-07,5.59E-07,5.59E-07,5.59E-07,5.59E-07, | 5.59E-07/) pdsza(3,:) = (/ | 0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0., | 0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0., | 0.,0.,0.,0.,0.,0.,0.,5.76E-37,6.37E-36, | 1.22E-35,1.34E-32,9.94E-32,1.85E-31,2.71E-31, | 8.53E-30,1.69E-29,2.95E-29,8.71E-29,1.45E-28, | 2.02E-28,3.57E-28,5.16E-28,-3.51E-19,-3.39E-18, | -6.42E-18,2.04E-18,1.81E-16,3.61E-16,2.79E-12, | 1.18E-11,2.09E-11,1.53E-10,4.79E-10,8.05E-10, | 2.47E-09,4.53E-09,6.94E-09,1.67E-08,2.65E-08, | 4.67E-08,7.27E-08,1.02E-07,1.49E-07,1.95E-07, | 2.42E-07,2.90E-07,3.34E-07,3.78E-07,4.10E-07, | 4.42E-07,4.66E-07,4.88E-07,5.03E-07,5.16E-07, | 5.27E-07,5.34E-07,5.41E-07,5.46E-07,5.49E-07, | 5.52E-07,5.54E-07,5.55E-07,5.57E-07,5.57E-07, | 5.58E-07,5.58E-07,5.59E-07,5.59E-07,5.59E-07, | 5.59E-07,5.59E-07,5.59E-07,5.59E-07,5.59E-07, | 5.59E-07,5.59E-07,5.59E-07,5.59E-07,5.59E-07, | 5.59E-07,5.59E-07,5.59E-07,5.59E-07,5.59E-07, | 5.59E-07,5.59E-07,5.59E-07,5.59E-07,5.59E-07, | 5.59E-07,5.59E-07,5.59E-07,5.59E-07,5.59E-07, | 5.59E-07,5.59E-07,5.59E-07,5.59E-07,5.59E-07, | 5.59E-07,5.59E-07,5.59E-07,5.59E-07,5.59E-07, | 5.59E-07,5.59E-07,5.59E-07,5.59E-07,5.59E-07, | 5.59E-07,5.59E-07/) pdsza(4,:) = (/ | 0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0., | 0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0., | 0.,0.,0.,0.,0.,0.,0.,0.,1.94E-35,5.65E-35, | 9.37E-35,1.17E-31,3.23E-31,5.29E-31,5.47E-30, | 1.80E-29,3.06E-29,6.91E-29,1.38E-28,2.06E-28, | 3.09E-28,4.87E-28,6.64E-28,-7.83E-18,-2.52E-17, | -4.25E-17,9.17E-15,3.45E-14,5.98E-14,1.27E-11, | 3.39E-11,5.50E-11,4.01E-10,9.17E-10,1.51E-09, | 4.39E-09,7.27E-09,1.32E-08,2.62E-08,3.93E-08, | 6.82E-08,9.90E-08,1.38E-07,1.87E-07,2.35E-07, | 2.82E-07,3.27E-07,3.69E-07,4.07E-07,4.35E-07, | 4.62E-07,4.83E-07,5.00E-07,5.14E-07,5.25E-07, | 5.33E-07,5.40E-07,5.45E-07,5.49E-07,5.52E-07, | 5.54E-07,5.55E-07,5.56E-07,5.57E-07,5.58E-07, | 5.58E-07,5.59E-07,5.59E-07,5.59E-07,5.59E-07, | 5.59E-07,5.59E-07,5.59E-07,5.59E-07,5.59E-07, | 5.59E-07,5.59E-07,5.59E-07,5.59E-07,5.59E-07, | 5.59E-07,5.59E-07,5.59E-07,5.59E-07,5.59E-07, | 5.59E-07,5.59E-07,5.59E-07,5.59E-07,5.59E-07, | 5.59E-07,5.59E-07,5.59E-07,5.59E-07,5.59E-07, | 5.59E-07,5.59E-07,5.59E-07,5.59E-07,5.59E-07, | 5.59E-07,5.59E-07,5.59E-07,5.59E-07,5.59E-07, | 5.59E-07,5.59E-07,5.59E-07,5.59E-07,5.59E-07, | 5.59E-07/) pdsza(5,:) = (/ | 0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0., | 0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0., | 0.,0.,0.,0.,0.,0.,0.,0.,0.,5.57E-35, | 1.18E-34,1.81E-34,2.41E-31,4.97E-31,7.52E-31, | 1.12E-29,2.49E-29,3.86E-29,1.00E-28,1.70E-28, | 2.41E-28,3.88E-28,5.68E-28,7.48E-28,-1.82E-17, | -3.99E-17,-6.17E-17,2.35E-14,5.43E-14,8.50E-14, | 2.27E-11,4.55E-11,9.52E-11,6.36E-10,1.18E-09, | 2.79E-09,5.80E-09,8.82E-09,1.94E-08,3.31E-08, | 5.27E-08,8.49E-08,1.17E-07,1.65E-07,2.14E-07, | 2.62E-07,3.10E-07,3.53E-07,3.95E-07,4.26E-07, | 4.55E-07,4.77E-07,4.96E-07,5.11E-07,5.22E-07, | 5.32E-07,5.39E-07,5.44E-07,5.48E-07,5.51E-07, | 5.53E-07,5.55E-07,5.56E-07,5.57E-07,5.58E-07, | 5.58E-07,5.59E-07,5.59E-07,5.59E-07,5.59E-07, | 5.59E-07,5.59E-07,5.59E-07,5.59E-07,5.59E-07, | 5.59E-07,5.59E-07,5.59E-07,5.59E-07,5.59E-07, | 5.59E-07,5.59E-07,5.59E-07,5.59E-07,5.59E-07, | 5.59E-07,5.59E-07,5.59E-07,5.59E-07,5.59E-07, | 5.59E-07,5.59E-07,5.59E-07,5.59E-07,5.59E-07, | 5.59E-07,5.59E-07,5.59E-07,5.59E-07,5.59E-07, | 5.59E-07,5.59E-07,5.59E-07,5.59E-07,5.59E-07, | 5.59E-07,5.59E-07,5.59E-07,5.59E-07,5.59E-07, | 5.59E-07/) pdsza(6,:) = (/ | 0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0., | 0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0., | 0.,0.,0.,0.,0.,0.,0.,0.,0.,3.88E-36, | 1.41E-35,2.44E-35,4.75E-32,1.54E-31,2.60E-31, | 2.76E-30,1.18E-29,2.08E-29,4.78E-29,1.06E-28, | 1.65E-28,2.46E-28,4.05E-28,5.63E-28,-1.34E-18, | -4.87E-18,-8.40E-18,-3.18E-18,2.12E-17,4.55E-17, | 5.16E-12,1.47E-11,2.42E-11,2.36E-10,5.80E-10, | 9.23E-10,3.10E-09,5.35E-09,9.72E-09,2.08E-08, | 3.19E-08,5.84E-08,8.75E-08,1.26E-07,1.76E-07, | 2.25E-07,2.75E-07,3.23E-07,3.68E-07,4.05E-07, | 4.39E-07,4.64E-07,4.87E-07,5.04E-07,5.17E-07, | 5.28E-07,5.36E-07,5.42E-07,5.47E-07,5.50E-07, | 5.53E-07,5.55E-07,5.56E-07,5.57E-07,5.58E-07, | 5.58E-07,5.59E-07,5.59E-07,5.59E-07,5.59E-07, | 5.59E-07,5.59E-07,5.59E-07,5.59E-07,5.59E-07, | 5.59E-07,5.59E-07,5.59E-07,5.59E-07,5.59E-07, | 5.59E-07,5.59E-07,5.59E-07,5.59E-07,5.59E-07, | 5.59E-07,5.59E-07,5.59E-07,5.59E-07,5.59E-07, | 5.59E-07,5.59E-07,5.59E-07,5.59E-07,5.59E-07, | 5.59E-07,5.59E-07,5.59E-07,5.59E-07,5.59E-07, | 5.59E-07,5.59E-07,5.59E-07,5.59E-07,5.59E-07, | 5.59E-07,5.59E-07,5.59E-07,5.59E-07,5.59E-07, | 5.59E-07/) pdsza(7,:) = (/ | 0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0., | 0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0., | 0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0., | 1.69E-33,9.97E-33,1.83E-32,8.95E-32,2.84E-30, | 5.59E-30,1.13E-29,4.64E-29,8.16E-29,1.17E-28, | 2.32E-28,3.52E-28,-2.75E-21,-5.01E-20,-9.74E-20, | -1.45E-19,-1.04E-15,-2.14E-15,1.93E-13,1.66E-12, | 3.14E-12,2.49E-11,1.65E-10,3.05E-10,9.70E-10, | 2.29E-09,3.62E-09,8.88E-09,1.64E-08,2.71E-08, | 5.11E-08,7.51E-08,1.18E-07,1.65E-07,2.16E-07, | 2.68E-07,3.18E-07,3.68E-07,4.05E-07,4.39E-07, | 4.68E-07,4.89E-07,5.07E-07,5.21E-07,5.31E-07, | 5.38E-07,5.44E-07,5.49E-07,5.52E-07,5.54E-07, | 5.56E-07,5.57E-07,5.57E-07,5.58E-07,5.58E-07, | 5.59E-07,5.59E-07,5.59E-07,5.59E-07,5.59E-07, | 5.59E-07,5.59E-07,5.59E-07,5.59E-07,5.59E-07, | 5.59E-07,5.59E-07,5.59E-07,5.59E-07,5.59E-07, | 5.59E-07,5.59E-07,5.59E-07,5.59E-07,5.59E-07, | 5.59E-07,5.59E-07,5.59E-07,5.59E-07,5.59E-07, | 5.59E-07,5.59E-07,5.59E-07,5.59E-07,5.59E-07, | 5.59E-07,5.59E-07,5.59E-07,5.59E-07,5.59E-07, | 5.59E-07,5.59E-07,5.59E-07,5.59E-07,5.59E-07, | 5.59E-07,5.59E-07,5.59E-07,5.59E-07/) pdsza(8,:) = (/ | 0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0., | 0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0., | 0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0., | 2.63E-33,2.21E-32,4.15E-32,6.10E-32,3.54E-30, | 7.53E-30,1.15E-29,4.76E-29,8.75E-29,1.27E-28, | 2.22E-28,3.46E-28,4.71E-28,-7.27E-20,-1.79E-19, | -2.85E-19,-5.72E-16,-1.79E-15,-3.01E-15,8.15E-13, | 2.34E-12,3.87E-12,5.26E-11,1.74E-10,2.95E-10, | 9.19E-10,1.96E-09,3.00E-09,6.50E-09,1.19E-08, | 1.73E-08,3.49E-08,5.27E-08,8.03E-08,1.21E-07, | 1.65E-07,2.16E-07,2.67E-07,3.21E-07,3.68E-07, | 4.10E-07,4.45E-07,4.75E-07,4.96E-07,5.14E-07, | 5.27E-07,5.37E-07,5.44E-07,5.48E-07,5.52E-07, | 5.54E-07,5.56E-07,5.57E-07,5.58E-07,5.58E-07, | 5.59E-07,5.59E-07,5.59E-07,5.59E-07,5.59E-07, | 5.59E-07,5.59E-07,5.59E-07,5.59E-07,5.59E-07, | 5.59E-07,5.59E-07,5.59E-07,5.59E-07,5.59E-07, | 5.59E-07,5.59E-07,5.59E-07,5.59E-07,5.59E-07, | 5.59E-07,5.59E-07,5.59E-07,5.59E-07,5.59E-07, | 5.59E-07,5.59E-07,5.59E-07,5.59E-07,5.59E-07, | 5.59E-07,5.59E-07,5.59E-07,5.59E-07,5.59E-07, | 5.59E-07,5.59E-07,5.59E-07,5.59E-07,5.59E-07, | 5.59E-07,5.59E-07,5.59E-07,5.59E-07/) pdsza(9,:) = (/ | 0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0., | 0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0., | 0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0., | 0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0., | 0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0., | 0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0., | 0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0., | 0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0., | 0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0., | 0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0.,0., | 0.,0.,0.,0.,0.,0.,0./) #else write(6,"('ERROR: unsupported NLEV = ',i6)") nlev #endif !************************************************************* ! Check for dres or res; if dres then interpolate ! !#if (NLEV==68) ! do j = 1,9 ! pesza(j,:) = pesza(j,:) ! pdsza(j,:) = pdsza(j,:) ! enddo !#elif (NLEV==136) ! do j = 1,9 ! count = 1 ! pesza(j,1) = pesza(j,1) ! pdsza(j,1) = pdsza(j,1) ! do i = 2,69 ! count = count+1 ! pesza(j,count) = (pesza(j,i-1)+pesza(j,i))/2. ! pesza(j,count+1) = pesza(j,i) ! pdsza(j,count) = (pdsza(j,i-1)+pdsza(j,i))/2. ! pdsza(j,count+1) = pdsza(j,i) ! count = count+1 ! enddo ! enddo !#else ! write(6,"('ERROR: unsupported NLEV = ',i6)") nlev !#endif !************************************************************* ! ! ! ! Initialize photoelectron frequency arrays (final beyond SZA = 90 degrees) ! -- pepdgrid (nlonp4,nlat,nlevp1) pepdgrid (:,:,:) = 0. pegrid (:,:,:) = 0. pdgrid (:,:,:) = 0. ! ! Local vectors : DO IL=1,isza DCHI(IL)=SZAH(IL)*pi/180. ENDDO DO k = 1,nlevp1 szp(k) = zibot + (k-1)*dz ENDDO ! ! --------------------------------------------------------------------- ! Interpolation Loops ! --------------------------------------------------------------------- ! * LATITUDE LOOP (SOUTH TO NORTH) DO J=1,nlat GLAT= -pi*0.5 + (FLOAT(J-1)+0.5)*dphi COSLAT = COS(GLAT) SINLAT = SIN(GLAT) ! * PRESSURE LOOP ! * NEW VTGCM VERTICAL SCALE NEEDED! (ZBOT= -16.0 ZTOP = 7.0) DO K=1,nlevp1 ! * LONGITUDE LOOP (add endpoints later) DO I=1,nlonp4 ! chapman.F REPRESENTATION OF SZA: Any sundec angle slt=amod((float(i-3)*dlamda+pi)*12./pi,24.) zchi=acos(sin_sundec*sinlat+cos_sundec*coslat* | cos(pi*(slt-12.)/12.)) ZCHI=ABS(ZCHI) DO 8 IL=2,isza IF(ZCHI .GE. DCHI(IL)) GO TO 8 RAT=(ZCHI-DCHI(IL-1))/(DCHI(IL)-DCHI(IL-1)) ! --------------------------------------------------------------- ! PHOTOELECTRON FREQUENCY : /sec ! Convert to #/cm3/sec in qjr.F where n2 is available ! PEGRID(I,J,K)=PESZA(IL-1,K)-RAT* ! | (PESZA(IL-1,K)-PESZA(IL,K)) ! Adding the photoelectron and photodissociation data together ! into one variable => pepdgrid PEPDGRID(I,J,K)=(PDSZA(IL-1,K)+PESZA(IL-1,K))-RAT* | ((PDSZA(IL-1,K)+PESZA(IL-1,K))- | (PDSZA(IL,K)+PESZA(IL,K))) PEGRID(I,J,K)=(PESZA(IL-1,K))-RAT* | (PESZA(IL-1,K)-PESZA(IL,K)) PDGRID(I,J,K)=(PDSZA(IL-1,K))-RAT* | (PDSZA(IL-1,K)-PDSZA(IL,K)) ! --------------------------------------------------------------- GO TO 9 8 CONTINUE 9 CONTINUE ENDDO ENDDO ENDDO DO J=1,nlat DO K=1,nlevp1 DO I=1,nlonp4 !re-arranging for plotting pepdgrid2(k,I,J) = pepdgrid(I,J,k) ENDDO ENDDO ENDDO write(6,"('end pepdgrid_init: pepdgrid min,max=',2e12.4)") | minval(pepdgrid),maxval(pepdgrid) write(6,"('end pepdgrid_init: pegrid min,max=',2e12.4)") | minval(pegrid),maxval(pegrid) write(6,"('end pepdgrid_init: pdgrid min,max=',2e12.4)") | minval(pdgrid),maxval(pdgrid) ! stop 'pepdgrid_init' ! ! No periodic points needed since nlon is used in dt.F ! Where is allocation of new pepdgrid 3-D array accomplished so that ! subdomain usage is avalable throughout code? !?? Carry entire domain for pepdgrid to qrj.F end subroutine pepdgrid_init !----------------------------------------------------------------------- ! subroutine alloc_pepdgrid ! ! Allocate pepdgrid for task subdomain: ! Called once per run from init_fields. ! ! use mpi_module,only: lon0,lon1,lat0,lat1 ! integer :: istat ! ! -- pepdgrid (nlonp4,nlat,nlevp1) ! allocate(pepdgrid(lon0:lon1,lat0:lat1,nlevp1),stat=istat) ! if (istat /= 0) write(6,"('>>> alloc_pepdgrid: error allocating', ! | ' pepdgrid: stat=',i3)") istat ! write(6,"('pepdgrid_init sub alloc_pepdgrid returning..')") ! ! end subroutine alloc_pepdgrid ! ------------------------------------------------------------- end module pepdgrid_init_module ! -------------------------------------------------------------