  @X  b   k820309    M          9.1         ÂèH                                                                                                       
       turbulence.f90 TURBULENCE       
   TURBULENCE_INIT XZ_TURBSCALAR XZA_TURBSCALAR PZ_TURBVELX XR_TURBVELZ XZ_SHEARKM XZ_DISPKM XZ_DISPHEAT XZ_BUOYKM EDDYVISCOSITY      @       €                           
       DIMXMIN DIMXMAX DIMZMIN DIMZMAX SPCNUM DELX DELZ                                          
       CPDRY GRAV XZ_POTTEMPBASICZ XZ_EFFMOLWTBASICZ XZ_EXNERBASICZ                  @                       
       XZ_AVR_PZ XZ_AVR_XR XZ_AVR_PR PZ_AVR_XZ PR_AVR_XZ XR_AVR_XZ                  @                       
       XZ_DX_PZ XZ_DZ_XR XR_DX_PR XR_DZ_XZ PZ_DX_XZ PZ_DZ_PR PR_DX_XR PR_DZ_PZ                  @                       
       STOREPOTTEMPDISP STOREPOTTEMPTURB                  @                       
       STOREMIXRTTURB (    `                                                
#SPCNUM    #DIMZMAX 	   #DIMZMIN 
   #DIMXMAX    #DIMXMIN    #PZ_VAR      5 r      5 r    5 r    p    5 r      & 5 r    5 r      & 5 r    5 r          5 r    5 r    p        5 r    5 r    p                    @  @                                      @                     	                 @                     
                 @                                      @                                  
   @                                 
      5 r      5 r    5 r    p    5 r      & 5 r    5 r      & 5 r    5 r          5 r    5 r    p        5 r    5 r    p              (    `                                                
#SPCNUM    #DIMZMAX 	   #DIMZMIN 
   #DIMXMAX    #DIMXMIN    #XR_VAR      5 r      5 r    5 r    p    5 r      & 5 r    5 r      & 5 r    5 r          5 r    5 r    p        5 r    5 r    p                   
   @                                 
      5 r      5 r    5 r    p    5 r      & 5 r    5 r      & 5 r    5 r          5 r    5 r    p        5 r    5 r    p              (    `                                 	               
#SPCNUM    #DIMZMAX 	   #DIMZMIN 
   #DIMXMAX    #DIMXMIN    #PR_VAR      5 r      5 r    5 r    p    5 r      & 5 r    5 r      & 5 r    5 r          5 r    5 r    p        5 r    5 r    p                   
   @                                 
      5 r      5 r    5 r    p    5 r      & 5 r    5 r      & 5 r    5 r          5 r    5 r    p        5 r    5 r    p              (    `                                                
#SPCNUM    #DIMZMAX 	   #DIMZMIN 
   #DIMXMAX    #DIMXMIN    #XZ_VAR      5 r      5 r    5 r    p    5 r      & 5 r    5 r      & 5 r    5 r          5 r    5 r    p        5 r    5 r    p                   
   @                                 
 
     5 r      5 r    5 r    p    5 r      & 5 r    5 r      & 5 r    5 r          5 r    5 r    p        5 r    5 r    p              (    `                                                
#SPCNUM    #DIMZMAX 	   #DIMZMIN 
   #DIMXMAX    #DIMXMIN    #XZ_VAR      5 r      5 r    5 r    p    5 r      & 5 r    5 r      & 5 r    5 r          5 r    5 r    p        5 r    5 r    p                   
   @                                 
      5 r      5 r    5 r    p    5 r      & 5 r    5 r      & 5 r    5 r          5 r    5 r    p        5 r    5 r    p              (    `                                                
#SPCNUM    #DIMZMAX 	   #DIMZMIN 
   #DIMXMAX    #DIMXMIN    #XZ_VAR      5 r      5 r    5 r    p    5 r      & 5 r    5 r      & 5 r    5 r          5 r    5 r    p        5 r    5 r    p                   
   @                                 
      5 r      5 r    5 r    p    5 r      & 5 r    5 r      & 5 r    5 r          5 r    5 r    p        5 r    5 r    p              (    `                                                
#SPCNUM    #DIMZMAX 	   #DIMZMIN 
   #DIMXMAX    #DIMXMIN    #PZ_VAR      5 r      5 r    5 r     p    5 r       & 5 r     5 r      & 5 r    5 r !         5 r    5 r     p        5 r !   5 r    p                   
   @                                 
      5 r      5 r    5 r     p    5 r       & 5 r     5 r      & 5 r    5 r !         5 r    5 r     p        5 r !   5 r    p              (    `                      "                          
#SPCNUM    #DIMZMAX 	   #DIMZMIN 
   #DIMXMAX    #DIMXMIN    #XR_VAR #     5 r      5 r    5 r     p    5 r       & 5 r     5 r      & 5 r    5 r !         5 r    5 r     p        5 r !   5 r    p                   
   @                     #            
      5 r      5 r    5 r     p    5 r       & 5 r     5 r      & 5 r    5 r !         5 r    5 r     p        5 r !   5 r    p              (    `                      $                          
#SPCNUM    #DIMZMAX 	   #DIMZMIN 
   #DIMXMAX    #DIMXMIN    #PR_VAR %     5 r      5 r    5 r     p    5 r       & 5 r     5 r      & 5 r    5 r !         5 r    5 r     p        5 r !   5 r    p                   
   @                     %            
      5 r      5 r    5 r     p    5 r       & 5 r     5 r      & 5 r    5 r !         5 r    5 r     p        5 r !   5 r    p              (    `                      &                          
#SPCNUM    #DIMZMAX 	   #DIMZMIN 
   #DIMXMAX    #DIMXMIN    #XZ_VAR '     5 r      5 r    5 r     p    5 r       & 5 r     5 r      & 5 r    5 r !         5 r    5 r     p        5 r !   5 r    p                   
   @                     '            
      5 r      5 r    5 r     p    5 r       & 5 r     5 r      & 5 r    5 r !         5 r    5 r     p        5 r !   5 r    p              (    `                      (           	               
#SPCNUM    #DIMZMAX 	   #DIMZMIN 
   #DIMXMAX    #DIMXMIN    #XZ_VAR )     5 r      5 r    5 r     p    5 r       & 5 r     5 r      & 5 r    5 r !         5 r    5 r     p        5 r !   5 r    p                   
   @                     )            
      5 r      5 r    5 r     p    5 r       & 5 r     5 r      & 5 r    5 r !         5 r    5 r     p        5 r !   5 r    p              (    `                      *                          
#SPCNUM    #DIMZMAX 	   #DIMZMIN 
   #DIMXMAX    #DIMXMIN    #PR_VAR +     5 r      5 r    5 r     p    5 r       & 5 r     5 r      & 5 r    5 r !         5 r    5 r     p        5 r !   5 r    p                   
   @                     +            
 
     5 r      5 r    5 r     p    5 r       & 5 r     5 r      & 5 r    5 r !         5 r    5 r     p        5 r !   5 r    p              (    `                      ,                          
#SPCNUM    #DIMZMAX 	   #DIMZMIN 
   #DIMXMAX    #DIMXMIN    #XR_VAR -     5 r      5 r    5 r     p    5 r       & 5 r     5 r      & 5 r    5 r !         5 r    5 r     p        5 r !   5 r    p                   
   @                     -            
      5 r      5 r    5 r     p    5 r       & 5 r     5 r      & 5 r    5 r !         5 r    5 r     p        5 r !   5 r    p              (    `                      .                          
#SPCNUM    #DIMZMAX 	   #DIMZMIN 
   #DIMXMAX    #DIMXMIN    #PZ_VAR /     5 r      5 r    5 r     p    5 r       & 5 r     5 r      & 5 r    5 r !         5 r    5 r     p        5 r !   5 r    p                   
   @                     /            
      5 r      5 r    5 r     p    5 r       & 5 r     5 r      & 5 r    5 r !         5 r    5 r     p        5 r !   5 r    p              #     @                         0                  #SPCNUM    #DIMZMAX 	   #DIMZMIN 
   #DIMXMAX    #DIMXMIN    #WORK 1        
   @                     1            
      5 r 2     5 r 3   5 r 4   p    5 r 4     & 5 r 4   5 r 3     & 5 r 2   5 r 5         5 r 3   5 r 4   p        5 r 5   5 r 2   p              #     @                         6                  #SPCNUM    #DIMZMAX 	   #DIMZMIN 
   #DIMXMAX    #DIMXMIN    #WORK 7        
   @                     7            
      5 r 2     5 r 3   5 r 4   p    5 r 4     & 5 r 4   5 r 3     & 5 r 2   5 r 5         5 r 3   5 r 4   p        5 r 5   5 r 2   p              #     @                         8                  #SPCNUM    #DIMZMAX 	   #DIMZMIN 
   #DIMXMAX    #DIMXMIN    #WORK 9        
   @                     9            
        p      5 r :   5 r ;   p    5 r ;     5 r <   5 r =   p    5 r =     & 5 r =   5 r <     & 5 r ;   5 r :     5 r >         5 r <   5 r =   p        5 r :   5 r ;   p      5 r >             #     @                          ?                   #TURBULENCE_INIT%MIN @             @                    @     MIN (    `  @                    A                           
#XZ_VAR B   #XZ_KH C     5 r 
     5 r    5 r    p    5 r      & 5 r    5 r      & 5 r 
   5 r 	         5 r    5 r    p        5 r 	   5 r 
   p                   
  @@                     B            
      5 r 
     5 r    5 r    p    5 r      & 5 r    5 r      & 5 r 
   5 r 	         5 r    5 r    p        5 r 	   5 r 
   p                   
  @@                     C            
      5 r 
     5 r    5 r    p    5 r      & 5 r    5 r      & 5 r 
   5 r 	         5 r    5 r    p        5 r 	   5 r 
   p              (    `  @                    D                           
#XZA_VAR E   #XZ_KH F       p      5 r 	   5 r 
   p    5 r 
     5 r    5 r    p    5 r      & 5 r    5 r      & 5 r 
   5 r 	     5 r          5 r    5 r    p        5 r 	   5 r 
   p      5 r                   
   @                     E            
        p      5 r 	   5 r 
   p    5 r 
     5 r    5 r    p    5 r      & 5 r    5 r      & 5 r 
   5 r 	     5 r          5 r    5 r    p        5 r 	   5 r 
   p      5 r                   
  @@                     F            
      5 r 
     5 r    5 r    p    5 r      & 5 r    5 r      & 5 r 
   5 r 	         5 r    5 r    p        5 r 	   5 r 
   p              (    `                       G                           
#XZ_KM H   #PZ_VELX I   #XR_VELZ J     5 r 
     5 r    5 r    p    5 r      & 5 r    5 r      & 5 r 
   5 r 	         5 r    5 r    p        5 r 	   5 r 
   p                   
  @@                     H            
      5 r 
     5 r    5 r    p    5 r      & 5 r    5 r      & 5 r 
   5 r 	         5 r    5 r    p        5 r 	   5 r 
   p                   
  @@                     I            
 	     5 r 
     5 r    5 r    p    5 r      & 5 r    5 r      & 5 r 
   5 r 	         5 r    5 r    p        5 r 	   5 r 
   p                   
  @@                     J            
 
     5 r 
     5 r    5 r    p    5 r      & 5 r    5 r      & 5 r 
   5 r 	         5 r    5 r    p        5 r 	   5 r 
   p              (    `                       K                           
#XZ_KM L   #PZ_VELX M   #XR_VELZ N     5 r 
     5 r    5 r    p    5 r      & 5 r    5 r      & 5 r 
   5 r 	         5 r    5 r    p        5 r 	   5 r 
   p                   
  @@                     L            
      5 r 
     5 r    5 r    p    5 r      & 5 r    5 r      & 5 r 
   5 r 	         5 r    5 r    p        5 r 	   5 r 
   p                   
  @@                     M            
      5 r 
     5 r    5 r    p    5 r      & 5 r    5 r      & 5 r 
   5 r 	         5 r    5 r    p        5 r 	   5 r 
   p                   
  @@                     N            
      5 r 
     5 r    5 r    p    5 r      & 5 r    5 r      & 5 r 
   5 r 	         5 r    5 r    p        5 r 	   5 r 
   p              (    `                      O                           
#XZ_KM P   #PZ_VELX Q   #XR_VELZ R     5 r 
     5 r    5 r    p    5 r      & 5 r    5 r      & 5 r 
   5 r 	         5 r    5 r    p        5 r 	   5 r 
   p                   
   @                     P            
      5 r 
     5 r    5 r    p    5 r      & 5 r    5 r      & 5 r 
   5 r 	         5 r    5 r    p        5 r 	   5 r 
   p                   
  @@                     Q            
      5 r 
     5 r    5 r    p    5 r      & 5 r    5 r      & 5 r 
   5 r 	         5 r    5 r    p        5 r 	   5 r 
   p                   
  @@                     R            
      5 r 
     5 r    5 r    p    5 r      & 5 r    5 r      & 5 r 
   5 r 	         5 r    5 r    p        5 r 	   5 r 
   p              (    `                       S                           
#XZ_KM T     5 r 
     5 r    5 r    p    5 r      & 5 r    5 r      & 5 r 
   5 r 	         5 r    5 r    p        5 r 	   5 r 
   p                   
   @                     T            
      5 r 
     5 r    5 r    p    5 r      & 5 r    5 r      & 5 r 
   5 r 	         5 r    5 r    p        5 r 	   5 r 
   p              (    `  @                    U                            
#XZ_KM V     5 r 
     5 r    5 r    p    5 r      & 5 r    5 r      & 5 r 
   5 r 	         5 r    5 r    p        5 r 	   5 r 
   p                   
   @                     V            
      5 r 
     5 r    5 r    p    5 r      & 5 r    5 r      & 5 r 
   5 r 	         5 r    5 r    p        5 r 	   5 r 
   p              (    `                      W           #                
#XZ_POTTEMP X     5 r 
     5 r    5 r    p    5 r      & 5 r    5 r      & 5 r 
   5 r 	         5 r    5 r    p        5 r 	   5 r 
   p                   
   @                     X            
 !     5 r 
     5 r    5 r    p    5 r      & 5 r    5 r      & 5 r 
   5 r 	         5 r    5 r    p        5 r 	   5 r 
   p              #     @                          Y                  #EDDYVISCOSITY%MAX Z   #EDDYVISCOSITY%SQRT [   #PZ_VELX \   #XR_VELZ ]   #XZ_POTTEMP ^   #XZ_KM _   #XZ_KMA `             @                    Z     MAX           @                    [     SQRT      
  @@                     \            
 $     5 r 
     5 r    5 r    p    5 r      & 5 r    5 r      & 5 r 
   5 r 	         5 r    5 r    p        5 r 	   5 r 
   p                   
  @@                     ]            
 %     5 r 
     5 r    5 r    p    5 r      & 5 r    5 r      & 5 r 
   5 r 	         5 r    5 r    p        5 r 	   5 r 
   p                   
  @@                     ^            
 &     5 r 
     5 r    5 r    p    5 r      & 5 r    5 r      & 5 r 
   5 r 	         5 r    5 r    p        5 r 	   5 r 
   p                   
  @@                     _            
 '     5 r 
     5 r    5 r    p    5 r      & 5 r    5 r      & 5 r 
   5 r 	         5 r    5 r    p        5 r 	   5 r 
   p                   D  @                     `            
 (      5 r 
     5 r    5 r    p    5 r      & 5 r    5 r      & 5 r 
   5 r 	         5 r    5 r    p        5 r 	   5 r 
   p                    @  @                     >                 @                     ;                 @                     :                 @                     =                 @                     <                 @                     2                 @                     5                 @                     4                 @                     3                 @                                      @                     !                 @                                       @                                      @                                      @                                      @                                      @                                œ   "      fn#fn     ¾   Š   b   uapp(TURBULENCE    H  e   j  GRIDSET    ­  q   J  BASICSET      p   J  AVERAGE &   Ž  |   J  DIFFERENTIATE_CENTER4    
  V   J  STOREPOTTEMP    `  C   J  STOREMIXRT "   £  „      XZ_AVR_PZ+AVERAGE    '  0     SPCNUM+GRIDSET     W  0     DIMZMAX+GRIDSET     ‡  0     DIMZMIN+GRIDSET     ·  0     DIMXMAX+GRIDSET     ç  0     DIMXMIN+GRIDSET )     ,  e   XZ_AVR_PZ%PZ_VAR+AVERAGE "   C  „      XZ_AVR_XR+AVERAGE )   Ç  ,  e   XZ_AVR_XR%XR_VAR+AVERAGE "   ó	  „      XZ_AVR_PR+AVERAGE )   w  ,  e   XZ_AVR_PR%PR_VAR+AVERAGE "   £  „      PZ_AVR_XZ+AVERAGE )   '  ,  e   PZ_AVR_XZ%XZ_VAR+AVERAGE "   S  „      PR_AVR_XZ+AVERAGE )   ×  ,  e   PR_AVR_XZ%XZ_VAR+AVERAGE "     „      XR_AVR_XZ+AVERAGE )   ‡  ,  e   XR_AVR_XZ%XZ_VAR+AVERAGE /   ³  „      XZ_DX_PZ+DIFFERENTIATE_CENTER4 6   7  ,  e   XZ_DX_PZ%PZ_VAR+DIFFERENTIATE_CENTER4 /   c  „      XZ_DZ_XR+DIFFERENTIATE_CENTER4 6   ç  ,  e   XZ_DZ_XR%XR_VAR+DIFFERENTIATE_CENTER4 /     „      XR_DX_PR+DIFFERENTIATE_CENTER4 6   —  ,  e   XR_DX_PR%PR_VAR+DIFFERENTIATE_CENTER4 /   Ã  „      XR_DZ_XZ+DIFFERENTIATE_CENTER4 6   G  ,  e   XR_DZ_XZ%XZ_VAR+DIFFERENTIATE_CENTER4 /   s  „      PZ_DX_XZ+DIFFERENTIATE_CENTER4 6   ÷   ,  e   PZ_DX_XZ%XZ_VAR+DIFFERENTIATE_CENTER4 /   #"  „      PZ_DZ_PR+DIFFERENTIATE_CENTER4 6   §#  ,  e   PZ_DZ_PR%PR_VAR+DIFFERENTIATE_CENTER4 /   Ó$  „      PR_DX_XR+DIFFERENTIATE_CENTER4 6   W&  ,  e   PR_DX_XR%XR_VAR+DIFFERENTIATE_CENTER4 /   ƒ'  „      PR_DZ_PZ+DIFFERENTIATE_CENTER4 6   )  ,  e   PR_DZ_PZ%PZ_VAR+DIFFERENTIATE_CENTER4 .   3*  †       STOREPOTTEMPDISP+STOREPOTTEMP 3   ¹*  ,  e   STOREPOTTEMPDISP%WORK+STOREPOTTEMP .   å+  †       STOREPOTTEMPTURB+STOREPOTTEMP 3   k,  ,  e   STOREPOTTEMPTURB%WORK+STOREPOTTEMP *   —-  †       STOREMIXRTTURB+STOREMIXRT /   .  „  e   STOREMIXRTTURB%WORK+STOREMIXRT     ¡/  U       TURBULENCE_INIT $   ö/  0      TURBULENCE_INIT%MIN    &0  S      XZ_TURBSCALAR %   y1  ,  a   XZ_TURBSCALAR%XZ_VAR $   ¥2  ,  a   XZ_TURBSCALAR%XZ_KH    Ñ3  ­      XZA_TURBSCALAR '   ~5  „  a   XZA_TURBSCALAR%XZA_VAR %   7  ,  a   XZA_TURBSCALAR%XZ_KH    .8  ]      PZ_TURBVELX "   ‹9  ,  a   PZ_TURBVELX%XZ_KM $   ·:  ,  a   PZ_TURBVELX%PZ_VELX $   ã;  ,  a   PZ_TURBVELX%XR_VELZ    =  ]      XR_TURBVELZ "   l>  ,  a   XR_TURBVELZ%XZ_KM $   ˜?  ,  a   XR_TURBVELZ%PZ_VELX $   Ä@  ,  a   XR_TURBVELZ%XR_VELZ    ðA  ]      XZ_SHEARKM !   MC  ,  a   XZ_SHEARKM%XZ_KM #   yD  ,  a   XZ_SHEARKM%PZ_VELX #   ¥E  ,  a   XZ_SHEARKM%XR_VELZ    ÑF  C      XZ_DISPKM     H  ,  a   XZ_DISPKM%XZ_KM    @I  G      XZ_DISPHEAT "   ‡J  ,  a   XZ_DISPHEAT%XZ_KM    ³K  H      XZ_BUOYKM %   ûL  ,  a   XZ_BUOYKM%XZ_POTTEMP    'N  ¬       EDDYVISCOSITY "   ÓN  0      EDDYVISCOSITY%MAX #   O  1      EDDYVISCOSITY%SQRT &   4O  ,  a   EDDYVISCOSITY%PZ_VELX &   `P  ,  a   EDDYVISCOSITY%XR_VELZ )   ŒQ  ,  a   EDDYVISCOSITY%XZ_POTTEMP $   ¸R  ,  a   EDDYVISCOSITY%XZ_KM %   äS  ,  a   EDDYVISCOSITY%XZ_KMA &   U  0      SPCNUM+GRIDSET=SPCNUM (   @U  0      DIMZMIN+GRIDSET=DIMZMIN (   pU  0      DIMZMAX+GRIDSET=DIMZMAX (    U  0      DIMXMIN+GRIDSET=DIMXMIN (   ÐU  0      DIMXMAX+GRIDSET=DIMXMAX (    V  0      DIMZMIN+GRIDSET=DIMZMIN (   0V  0      DIMZMAX+GRIDSET=DIMZMAX (   `V  0      DIMXMIN+GRIDSET=DIMXMIN (   V  0      DIMXMAX+GRIDSET=DIMXMAX (   ÀV  0      DIMZMIN+GRIDSET=DIMZMIN (   ðV  0      DIMZMAX+GRIDSET=DIMZMAX (    W  0      DIMXMIN+GRIDSET=DIMXMIN (   PW  0      DIMXMAX+GRIDSET=DIMXMAX (   €W  0      DIMZMIN+GRIDSET=DIMZMIN (   °W  0      DIMZMAX+GRIDSET=DIMZMAX (   àW  0      DIMXMIN+GRIDSET=DIMXMIN (   X  0      DIMXMAX+GRIDSET=DIMXMAX 