c=========================================================================

      subroutine subsum(nz,z,p,t,wh,wo) 1
      parameter (mz=33)
      dimension z(*),p(*),t(*),wh(*),wo(*)             
      dimension z4(mz),p4(mz),t4(mz),wh4(mz),wo4(mz)
c                                                                       
c     model: subarctic summer mc'clatchey                             
c                                                                       
      data z4/
     &    0.,    1.,    2.,    3.,    4.,    5.,    6.,    7.,    8.,   
     &    9.,   10.,   11.,   12.,   13.,   14.,   15.,   16.,   17.,   
     &   18.,   19.,   20.,   21.,   22.,   23.,   24.,   25.,   30.,   
     &   35.,   40.,   45.,   50.,   70.,  100./                 
      data p4/
     &  1.010e+03,8.960e+02,7.929e+02,7.000e+02,6.160e+02,5.410e+02,      
     &  4.730e+02,4.130e+02,3.590e+02,3.107e+02,2.677e+02,2.300e+02,      
     &  1.977e+02,1.700e+02,1.460e+02,1.250e+02,1.080e+02,9.280e+01,      
     &  7.980e+01,6.860e+01,5.890e+01,5.070e+01,4.360e+01,3.750e+01,      
     &  3.227e+01,2.780e+01,1.340e+01,6.610e+00,3.400e+00,1.810e+00,      
     &  9.870e-01,7.070e-02,3.000e-04/                          
      data t4/
     &  2.870e+02,2.820e+02,2.760e+02,2.710e+02,2.660e+02,2.600e+02,      
     &  2.530e+02,2.460e+02,2.390e+02,2.320e+02,2.250e+02,2.250e+02,      
     &  2.250e+02,2.250e+02,2.250e+02,2.250e+02,2.250e+02,2.250e+02,      
     &  2.250e+02,2.250e+02,2.250e+02,2.250e+02,2.250e+02,2.250e+02,      
     &  2.260e+02,2.280e+02,2.350e+02,2.470e+02,2.620e+02,2.740e+02,      
     &  2.770e+02,2.160e+02,2.100e+02/                          
      data wh4/
     &  9.100e+00,6.000e+00,4.200e+00,2.700e+00,1.700e+00,1.000e+00,      
     &  5.400e-01,2.900e-01,1.300e-01,4.200e-02,1.500e-02,9.400e-03,      
     &  6.000e-03,1.800e-03,1.000e-03,7.600e-04,6.400e-04,5.600e-04,      
     &  5.000e-04,4.900e-04,4.500e-04,5.100e-04,5.100e-04,5.400e-04,      
     &  6.000e-04,6.700e-04,3.600e-04,1.100e-04,4.300e-05,1.900e-05,      
     &  6.300e-06,1.400e-07,1.000e-09/                          
      data wo4/
     &  4.900e-05,5.400e-05,5.600e-05,5.800e-05,6.000e-05,6.400e-05,      
     &  7.100e-05,7.500e-05,7.900e-05,1.100e-04,1.300e-04,1.800e-04,      
     &  2.100e-04,2.600e-04,2.800e-04,3.200e-04,3.400e-04,3.900e-04,      
     &  4.100e-04,4.100e-04,3.900e-04,3.600e-04,3.200e-04,3.000e-04,      
     &  2.800e-04,2.600e-04,1.400e-04,9.200e-05,4.100e-05,1.300e-05,      
     &  4.300e-06,8.600e-08,4.300e-11/                          
      do 1 i=1,mz                                                       
        z(i)=z4(i)                                                        
        p(i)=p4(i)                                                        
        t(i)=t4(i)                                                        
        wh(i)=wh4(i)                                                      
        wo(i)=wo4(i)                                                      
    1 continue             
      nz=mz                                             
      return                                                            
      end