  A6  V   k820309    M          9.1         H                                                                                                       
       numdiffusion.f90 NUMDIFFUSION          NUMDIFFUSION_INIT XZ_NUMDIFFSCALAR XZ_NUMDIFFKM XZA_NUMDIFFSCALAR PZ_NUMDIFFVELX XR_NUMDIFFVELZ                  @                       
       MESSAGENOTIFY      @                                  
       DIMXMIN DIMXMAX DIMZMIN DIMZMAX SPCNUM DELX DELZ                                          
       DELTIMELONG                  @                       
       XZ_DX_PZ XZ_DZ_XR XR_DX_PR XR_DZ_XZ PZ_DX_XZ PZ_DZ_PR PR_DX_XR PR_DZ_PZ                  @                       
       STOREPOTTEMPDIFF                  @                       
       STOREMIXRTDIFF                                          u #MESSAGENOTIFYC    #MESSAGENOTIFYI    #     @     @                                    #MESSAGENOTIFYC%TRIM    #LEVEL 	   #WHERE 
   #MESSAGE    #I    #R    #D    #L    #N    #C1    #C2    #C3    #CA              @                         TRIM       
   @                     	            1       
   @                     
            1       
   @                                 1       
  @                                           &                             
  @                                 	          &                             
  @                                
          &                             
  @                                           &                             
  @                                           &                             
  @                                 1       
  @                                 1       
  @                                 1 ,      
  @                                           &                       1 #     @     @                                     #MESSAGENOTIFYI%PRESENT    #NUMBER    #WHERE    #MESSAGE    #I    #R    #D    #L    #N    #C1    #C2     #C3 !   #CA "             @                         PRESENT       
   @                                   
   @                                 1       
  @                                 1       
  @                                           &                             
  @                                 	 	         &                             
  @                                
 
         &                             
  @                                           &                             
  @                                           &                             
  @                                 1       
  @                                  1       
  @                     !            1 ,      
  @                     "                      &                       1 (    `                      #                          
#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              (    `                      0                          
#SPCNUM $   #DIMZMAX %   #DIMZMIN &   #DIMXMAX '   #DIMXMIN (   #PR_VAR 1     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                   
   @                     1            
      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                          
#SPCNUM $   #DIMZMAX %   #DIMZMIN &   #DIMXMAX '   #DIMXMIN (   #XZ_VAR 3     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                   
   @                     3            
      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              (    `                      4           	               
#SPCNUM $   #DIMZMAX %   #DIMZMIN &   #DIMXMAX '   #DIMXMIN (   #XZ_VAR 5     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            
      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              (    `                      6                          
#SPCNUM $   #DIMZMAX %   #DIMZMIN &   #DIMXMAX '   #DIMXMIN (   #PR_VAR 7     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                   
   @                     7            
 
     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              (    `                      8                          
#SPCNUM $   #DIMZMAX %   #DIMZMIN &   #DIMXMAX '   #DIMXMIN (   #XR_VAR 9     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                   
   @                     9            
      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 (   #WORK =        
   @                     =            
      5 r >     5 r ?   5 r @   p    5 r @     & 5 r @   5 r ?     & 5 r >   5 r A         5 r ?   5 r @   p        5 r A   5 r >   p              #     @                         B                  #SPCNUM $   #DIMZMAX %   #DIMZMIN &   #DIMXMAX '   #DIMXMIN (   #WORK C        
   @                     C            
        p      5 r D   5 r E   p    5 r E     5 r F   5 r G   p    5 r G     & 5 r G   5 r F     & 5 r E   5 r D     5 r H         5 r F   5 r G   p        5 r D   5 r E   p      5 r H             #     @                          I                    (    `  @                    J                           
#XZ_SCALAR K     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            
      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                           
#XZ_SCALAR 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                   
  @@                     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           	                
#XZA_SCALAR O       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 $                  
   @                     O            
        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 $             (    `                       P                           
#PZ_VARX 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                   
  @@                     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                           
#XR_VARZ S     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            
      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                 @                     E                 @                     D                 @                     G                 @                     F                 @                     >                 @                     A                 @                     @                 @                     ?                 @                     *                 @                     -                 @                     ,                 @                     +              &      fn#fn "      l   b   uapp(NUMDIFFUSION    .  B   J  DC_MESSAGE    p  e   j  GRIDSET      @   J  TIMESET &     |   J  DIFFERENTIATE_CENTER4      E   J  STOREPOTTEMP      C   J  STOREMIXRT -     \       gen@MESSAGENOTIFY+DC_MESSAGE 9   u        MESSAGENOTIFYC+DC_MESSAGE=MESSAGENOTIFYC 4   0  1      MESSAGENOTIFYC%TRIM+DC_MESSAGE=TRIM 0   a  8   e   MESSAGENOTIFYC%LEVEL+DC_MESSAGE 0     8   e   MESSAGENOTIFYC%WHERE+DC_MESSAGE 2     8   e   MESSAGENOTIFYC%MESSAGE+DC_MESSAGE ,   	  `   e   MESSAGENOTIFYC%I+DC_MESSAGE ,   i  `   e   MESSAGENOTIFYC%R+DC_MESSAGE ,     `   e   MESSAGENOTIFYC%D+DC_MESSAGE ,   )  `   e   MESSAGENOTIFYC%L+DC_MESSAGE ,     `   e   MESSAGENOTIFYC%N+DC_MESSAGE -     8   e   MESSAGENOTIFYC%C1+DC_MESSAGE -   !  8   e   MESSAGENOTIFYC%C2+DC_MESSAGE -   Y  8   e   MESSAGENOTIFYC%C3+DC_MESSAGE -     d   e   MESSAGENOTIFYC%CA+DC_MESSAGE 9           MESSAGENOTIFYI+DC_MESSAGE=MESSAGENOTIFYI :     4      MESSAGENOTIFYI%PRESENT+DC_MESSAGE=PRESENT 1     0   e   MESSAGENOTIFYI%NUMBER+DC_MESSAGE 0   	  8   e   MESSAGENOTIFYI%WHERE+DC_MESSAGE 2   P	  8   e   MESSAGENOTIFYI%MESSAGE+DC_MESSAGE ,   	  `   e   MESSAGENOTIFYI%I+DC_MESSAGE ,   	  `   e   MESSAGENOTIFYI%R+DC_MESSAGE ,   H
  `   e   MESSAGENOTIFYI%D+DC_MESSAGE ,   
  `   e   MESSAGENOTIFYI%L+DC_MESSAGE ,     `   e   MESSAGENOTIFYI%N+DC_MESSAGE -   h  8   e   MESSAGENOTIFYI%C1+DC_MESSAGE -     8   e   MESSAGENOTIFYI%C2+DC_MESSAGE -     8   e   MESSAGENOTIFYI%C3+DC_MESSAGE -     d   e   MESSAGENOTIFYI%CA+DC_MESSAGE /   t        XZ_DX_PZ+DIFFERENTIATE_CENTER4      0     SPCNUM+GRIDSET     (  0     DIMZMAX+GRIDSET     X  0     DIMZMIN+GRIDSET       0     DIMXMAX+GRIDSET       0     DIMXMIN+GRIDSET 6     ,  e   XZ_DX_PZ%PZ_VAR+DIFFERENTIATE_CENTER4 /           XZ_DZ_XR+DIFFERENTIATE_CENTER4 6     ,  e   XZ_DZ_XR%XR_VAR+DIFFERENTIATE_CENTER4 /           XR_DX_PR+DIFFERENTIATE_CENTER4 6   H  ,  e   XR_DX_PR%PR_VAR+DIFFERENTIATE_CENTER4 /   t        XR_DZ_XZ+DIFFERENTIATE_CENTER4 6     ,  e   XR_DZ_XZ%XZ_VAR+DIFFERENTIATE_CENTER4 /   $        PZ_DX_XZ+DIFFERENTIATE_CENTER4 6     ,  e   PZ_DX_XZ%XZ_VAR+DIFFERENTIATE_CENTER4 /           PZ_DZ_PR+DIFFERENTIATE_CENTER4 6   X  ,  e   PZ_DZ_PR%PR_VAR+DIFFERENTIATE_CENTER4 /           PR_DX_XR+DIFFERENTIATE_CENTER4 6     ,  e   PR_DX_XR%XR_VAR+DIFFERENTIATE_CENTER4 /   4         PR_DZ_PZ+DIFFERENTIATE_CENTER4 6   !  ,  e   PR_DZ_PZ%PZ_VAR+DIFFERENTIATE_CENTER4 .   "         STOREPOTTEMPDIFF+STOREPOTTEMP 3   j#  ,  e   STOREPOTTEMPDIFF%WORK+STOREPOTTEMP *   $         STOREMIXRTDIFF+STOREMIXRT /   %    e   STOREMIXRTDIFF%WORK+STOREMIXRT "   &  <       NUMDIFFUSION_INIT !   &  K      XZ_NUMDIFFSCALAR +   '(  ,  a   XZ_NUMDIFFSCALAR%XZ_SCALAR    S)  G      XZ_NUMDIFFKM '   *  ,  a   XZ_NUMDIFFKM%XZ_SCALAR "   +        XZA_NUMDIFFSCALAR -   k-    a   XZA_NUMDIFFSCALAR%XZA_SCALAR    .  E      PZ_NUMDIFFVELX '   40  ,  a   PZ_NUMDIFFVELX%PZ_VARX    `1  E      XR_NUMDIFFVELZ '   2  ,  a   XR_NUMDIFFVELZ%XR_VARZ &   3  0      SPCNUM+GRIDSET=SPCNUM (   4  0      DIMZMIN+GRIDSET=DIMZMIN (   14  0      DIMZMAX+GRIDSET=DIMZMAX (   a4  0      DIMXMIN+GRIDSET=DIMXMIN (   4  0      DIMXMAX+GRIDSET=DIMXMAX (   4  0      DIMZMIN+GRIDSET=DIMZMIN (   4  0      DIMZMAX+GRIDSET=DIMZMAX (   !5  0      DIMXMIN+GRIDSET=DIMXMIN (   Q5  0      DIMXMAX+GRIDSET=DIMXMAX (   5  0      DIMZMIN+GRIDSET=DIMZMIN (   5  0      DIMZMAX+GRIDSET=DIMZMAX (   5  0      DIMXMIN+GRIDSET=DIMXMIN (   6  0      DIMXMAX+GRIDSET=DIMXMAX 