



          *********************************************************
          *********************************************************
          *** PROGRAM TO  SOLVE FOR THE  MESOSCALE AND SYNOPTIC ***
          *** VARIATIONS OF POTENTIAL TEMPERATURE, U AND V WIND ***
          *** COMPONENTS, AND THE MESOSCALE VERTICAL VELOCITY W ***
          *** PRESSURE AND DISTRIBUTION OF POLLUTANTS C  HAVING ***
          *** SOURCES  Q. THE SYNOPTIC  SCALE COMPONENTS ARE IN ***
          *** QUASSI-STEADY STATE BALANCE,  WHILE THE MESOSCALE ***
          *** PRESSURE AND VELOCITY  W ARE FOUND DIAGNOSTICALLY ***
          *********************************************************
          *********************************************************




------------------------------------------------!
  ***  D A T A   F O R   M E S O - R U N  ***   !
------------------------------------------------!
              INTEGER CONSTANTS                 !
------------------------------------------------!
     L E G E N T             !   V A L U E      !
------------------------------------------------!
------------------------------------------------!
             REAL CONSTANTS                     !
------------------------------------------------!
     L E G E N T             !   V A L U E      !
------------------------------------------------!

 FOR LINEAR STABILITY,  UMAX < 15.78  VMAX <   .00

 ** TOTAL SIMULATION PERIOD    86400.0S, = 24.00 HOURS **
 ** HEAT AMPLITUDE  11.00 K,    LAG = 14.00(H),   BIAS A     .1


 ********   LAND SURFACE FORCING  FUNCTION ********* 

 T(X,Y,T)= 11.0{EXP -((T-14.0)/ 4.9)**2 -  .1}

 ********   TIME UNITS ARE IN HOURS        *********

 WORK ARRAY HAS    17826 ELEMENTS.  


 ACTUAL SOURCE VALUE    .1023E+00 AT GRID POINTS I,J,K   15   1   1
 ACTUAL SOURCE VALUE    .1036E+00 AT GRID POINTS I,J,K   16   1   1
 ACTUAL SOURCE VALUE    .1023E+00 AT GRID POINTS I,J,K   17   1   1
 ACTUAL SOURCE VALUE    .1274E+00 AT GRID POINTS I,J,K   15   1   2
 ACTUAL SOURCE VALUE    .1290E+00 AT GRID POINTS I,J,K   16   1   2
 ACTUAL SOURCE VALUE    .1274E+00 AT GRID POINTS I,J,K   17   1   2
 ACTUAL SOURCE VALUE    .1023E+00 AT GRID POINTS I,J,K   15   1   3
 ACTUAL SOURCE VALUE    .1036E+00 AT GRID POINTS I,J,K   16   1   3
 ACTUAL SOURCE VALUE    .1023E+00 AT GRID POINTS I,J,K   17   1   3



 ***********************  MEAN PROFILES ***********************

 CONVERGENCE AFTER     34 ITERATIONS.    NORM=    .0000E+00
 EQUIVALENT TIME FOR THE MEAN PROFILES  IS    2040.00 SEC. 
 Z0=   .0200   ZMIX= 147.91   U*=   .03466   L(4)=      305.58


   K    HEIGHT (M)       UG(Z)     VG(Z)      TM(Z)       KM(Z)       UM(Z)       VM(Z)  

   1        .00          6.88       .00      290.07         .10         .00         .00
   2      29.86          6.88       .00      290.17         .21         .90        1.66
   3      67.15          6.89       .00      290.33         .25        1.76        2.20
   4     113.73          6.89       .00      290.56         .18        2.78        2.54
   5     171.89          6.90       .00      290.86         .12        3.93        2.71
   6     244.53          6.91       .00      291.22         .10        5.09        2.64
   7     335.24          6.92       .00      291.68         .10        6.17        2.29
   8     448.53          6.93       .00      292.24         .10        7.00        1.69
   9     590.00          6.95       .00      292.95         .10        7.46         .96
  10     766.69          6.97       .00      293.83         .10        7.51         .32
  11     987.34          7.00       .00      294.94         .10        7.31        -.05
  12    1262.90          7.03       .00      296.31         .10        7.10        -.11
  13    1607.03          7.07       .00      298.04         .10        7.05        -.04
  14    2036.80          7.12       .00      300.18         .10        7.11         .00
  15    2573.52          7.19       .00      302.87         .10        7.19         .00
  16    3243.79          7.27       .00      306.22         .10        7.27         .00

 **************************************************************


 ****************  X - FILTER  ***************** 

 1122345780000000000000000000000000000000000000000000000875432211
 


            ****   TOPOGRAPHY  MAX =       1.00 MIN =      -1.00  METERS   ****


 LLLLLLLLLLLLLLLLLWWWWWWWWWWWWLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLLL
 






                  SOLVING THE EQUATION 
   DC/DT=U*DC/DX+V*DC/DY+KX*D2C/DX2+KY*DC/DY+D(KZ*DC/DZ)/DZ + Q,
   WHERE C(X,Y,Z,T0)=0 WITH PERIODIC BOUNDARY CONDITIONS IN    X,Y
   AND FLUX VERTICAL BOUNDARY CONDITIONS. 




   X AXIS IS DIVIDED INTO   64  SUBINTERVALS
   Y AXIS IS DIVIDED INTO    1  SUBINTERVALS
   Z AXIS IS DIVIDED INTO   16  SUBINTERVALS
         720  TIME STEPS WERE REQUESTED 


  BATCH MODE = NO       SPECTRAL FILTER       STRECHED Z, Z=F(H)  

   Z=120.0*(EXP(H/  .3)-1),  UPWIND A=  .5


    FOLLOWING NUMBERS ARE IN USER WORLD UNITS


   X  DOMAIN  (       .0, 504000.0),    DX=  8000.00
   Y  DOMAIN  (       .0,       .0),    DY=  5000.00
   Z  DOMAIN  (       .0,    447.9),    DZ=    29.86
   T  DOMAIN  (       .0,  86280.0),    DT=   120.00

   CENTRAL LATITUDE  40.00 DEG,    GRAPH STEP 720


  DIFFUSIVITIES   (KX,KY,KZ) = (   7000.00   1000.00       .10)

 ** TIME SCHEMES USED **  LEAP  FROG 2ND IN X,Y  CRANK-NICKOLSON IN Z

 T(     0)=      .0  FUM=  .1000E+01  E=  .0000E+00  T=  .0000E+00  P=  .0000E+00
        TS= 290.07    U=  .0000E+00  V=  .0000E+00  W=  .0000E+00
 T(   100)= 12000.0  FUM=  .9473E+04  E=  .3018E+01  T= -.7642E+02  P=  .1976E+00
        TS= 289.06    U=  .4021E+01  V=  .3531E+01  W= -.5349E-04
 T(   200)= 24000.0  FUM=  .1692E+05  E=  .1536E+02  T= -.1363E+02  P=  .5779E+00
        TS= 290.11    U=  .1951E+02  V= -.1619E+02  W= -.7112E-04
 T(   300)= 36000.0  FUM=  .2158E+05  E=  .4492E+03  T=  .7785E+03  P= -.6192E+01
        TS= 294.55    U= -.1660E+03  V= -.3455E+02  W=  .1114E-02
 T(   400)= 48000.0  FUM=  .2589E+05  E=  .1399E+04  T=  .2615E+04  P= -.7640E+02
        TS= 299.75    U= -.3548E+03  V=  .2449E+03  W=  .8267E-02
 T(   500)= 60000.0  FUM=  .2765E+05  E=  .2696E+04  T=  .2964E+04  P= -.1223E+03
        TS= 297.21    U= -.1302E+03  V=  .6298E+03  W=  .1075E-01
 T(   600)= 72000.0  FUM=  .3726E+05  E=  .3897E+04  T=  .2503E+04  P= -.1231E+03
        TS= 291.46    U=  .5585E+03  V=  .5339E+03  W=  .5730E-02
 T(   700)= 84000.0  FUM=  .5024E+05  E=  .4078E+04  T=  .2284E+04  P= -.1243E+03
        TS= 289.27    U=  .7945E+03  V= -.2643E+03  W=  .3784E-02
 **      CPU FOR CONCENTRATIONS  C =            .00   **
 **      CPU FOR TEMPERATURE THETA =            .00   **
 **      CPU FOR EXNER PRESSURE PI =            .00   **
 **      CPU FOR X MOMENTUM EQ.  U =            .00   **
 **      CPU FOR Y MOMENTUM EQ.  V =            .00   **
 **      CPU FOR CONTINUITY EQ.  W =            .00   **
 **      CPU FOR PRINTING UNKNOWNS =            .00   **
 **      CPU FOR  BOUNDARY   LAYER =            .00   **
 **      CPU FOR HIGH  PASS FILTER =            .00   **
 **      CPU FOR SYNOPTIC VARIABLS =            .00   **
 **                                                   **
