M o d e l   c o d e   f o r   " M o r e   t h a n   e m i s s i o n s   a n d   c h e m i s t r y :   F i r e   s i z e ,   d i l u t i o n ,   a n d   b a c k g r o u n d   a e r o s o l   a l s o   g r e a t l y   i n f l u e n c e   n e a r - f i e l d   b i o m a s s   b u r n i n g   a e r o s o l   a g i n g " ,   s u b m i t t e d   t o   J o u r n a l   o f   G e o p h y s i c a l   R e s e a r c h :   A t m o s p h e r e s   
 
 
 
 T h e   f o l l o w i n g   c o d e   i s   f o r   t h e   T O M A S - V B S   m o d e l   u s e d   i n   H o d s h i r e   e t   a l .   ( 2 0 1 9 )   ( t i t l e   g i v e n   a b o v e ) .   T h e   m o d e l   i s   w r i t t e n   i n   F O R T R A N   a n d   i s   s e t   u p   t o   b e   r u n   w i t h   a   p y t h o n   w r a p p e r .     A   s a m p l e   p y t h o n   s c r i p t   " r u n A r r a y _ a l l f i r e s i z e s _ s t a b i l i t y c l a s s D _ f o r q s u b 1 6 p r o c s . p y "   t o   r u n   t h e   m o d e l   f o r   t h e   b a s e   c a s e s   p r e s e n t e d   i n   t h e   m a i n   t e x t   ( f o r   a l l   f i r e   s i z e s ,   u s i n g   a   s t a b i l i t y   c l a s s   o f   D   a n d   t h e   " u p p e r - b o u n d   c h e m i s t r y " ) .   
 
 
 
 F o r   a n y   q u e s t i o n s ,   p l e a s e   c o n t a c t   A n n a   H o d s h i r e   a t   h o d s h i r e @ r a m s . c o l o s t a t e . e d u   o r   o r   J e f f r e y   P i e r c e   a t   j e f f r e y . p i e r c e @ c o l o s t a t e . e d u .   
 
 
 
 F i l e   i n v e n t o r y : 
 
 
 
 a d d b a c k g r o u n d . f 
 
 a e r o d e n s . f 
 
 a e r o d i a g . f 
 
 a r r a y . f 
 
 B B 1 9 2 S M 9 . C O M 
 
 b o x . f 
 
 b o x i n i t . f 
 
 B T 2 6 3 b o x . C O M 
 
 c a l c u l a t e _ i n i t _ N _ s i z e p a r a m e t e r s . p y   ( P y t h o n   s c r i p t   f o r   g e t t i n g   i n i t i a l   t o t a l   n u m b e r   f o r   c h o s e n   f i r e   a r e a   a n d   s t a b i l i t y   c l a s s . ) 
 
 d m d t _ i n t . f 
 
 e z c o n d . f 
 
 e z m a s s . f 
 
 e z n h 3 e q m . f 
 
 e z w a t e r e q m . f 
 
 g a s d i f f . f 
 
 g e t C o n d S i n k . f 
 
 g e t l e n . f 
 
 g e t M a s s S i n k . f 
 
 i n i t b o u n d s . f 
 
 l o g i n i t . f 
 
 l o g i n i t _ 2 m o d e . f 
 
 m a k e f i l e 
 
 m n f i x . f 
 
 r e a d m e . t x t   ( R e a d m e   i n f o r m a t i o n   a b o u t   r u n n i n g   m o d e l ) 
 
 r e p o r t . f 
 
 r u n A r r a y _ a l l f i r e s i z e s _ s t a b i l i t y c l a s s D _ f o r q s u b 1 6 p r o c s . p y   ( P y t h o n   s c r i p t   f o r   r u n n i n g   b a s e   c a s e s   f o r   m o d e l . ) 
 
 s c a l e m o m . f 
 
 s i z e c o d e . C O M 
 
 s o a c o n d . f 
 
 s p e c i e s m a p . f 
 
 s t o r e n m . f 
 
 t m c o n d . f 
 
 w a t e r n a c l . f 
 
 w a t e r s o 4 . f 
 
 
 
 T h e   . f   a n d   . C O M   f i l e s   a r e   F o r t r a n   f i l e s .   
 
 
 
 
 
 I n   o r d e r   t o   r u n   T O M A S , f i r s t   t y p e   " m a k e c l e a n " ,   t h e n   " m a k e   b o x . e x e "   i n   t h e   c o m m a n d   l i n e .   O n c e   b o x . e x e   s u c c e s s f u l l y   c o m p i l e s ,   t h e   " r u n A r r a y _ a l l f i r e s i z e s _ s t a b i l i t y c l a s s D _ f o r q s u b 1 6 p r o c s . p y "   f i l e   i s   i n t e n d e d   a s   a   p y t h o n   s c r i p t   t o   r u n   T O M A S .   " r u n A r r a y _ a l l f i r e s i z e s _ s t a b i l i t y c l a s s D _ f o r q s u b 1 6 p r o c s . p y "   h a s   b e e n   s e t   u p   t o   b e   s u b m i t t e d   t h r o u g h   a   s h e l l   ( . s h )   s c r i p t   u s i n g   1 6   p r o c e s s o r s ;   t h i s   c a n   b e   a l t e r e d   t o   o n l y   r u n   o n   o n e   p r o c e s s o r .   
 
 
 
 T h e   f i l e   " c a l c u l a t e _ i n i t _ N _ s i z e p a r a m e t e r s . p y "   a l l o w s   r e s e a r c h e r s   t o   f i n d   t h e   i n i t i a l   n u m b e r   c o n c e n t r a t i o n   a n d   s i g m a _ y   ( S e c t   2 . 2   o f   B i a n   e t   a l . ,   2 0 1 7 )   f o r   a   g i v e n   f i r e   a r e a   a n d   s t a b i l i t y   c l a s s .   F o r   f u r t h e r   d e t a i l s   o n   t h e   i n p u t   a n d   o u t p u t   p a r a m e t e r s   o f   t h i s   f i l e ,   p l e a s e   r e f e r   t o   H o d s h i r e   e t   a l .   ( 2 0 1 9 )   a n d   S e c t   2 . 2   o f   B i a n   e t   a l . ,   2 0 1 7 .   
 
 
 
 D a t e   o f   l a s t   a l t e r a t i o n   b e f o r e   s u b m i s s i o n :   2 / 5 / 1 9 
 
 
 
 C i t a t i o n   o f   t h e   p a p e r   ( a t   t h e   t i m e   o f   r e v i e w   a n d   c r e a t i o n   o f   t h i s   f i l e ) :   A .   L .   H o d s h i r e ,   Q .   B i a n ,   E .   R a m n a r i n e ,   C .   R .   L o n s d a l e ,   M .   J .   A l v a r a d o ,   S .   M .   K r e i d e n w e i s ,   S .   H .   J a t h a r ,   J .   R .   P i e r c e :   M o r e   t h a n   e m i s s i o n s   a n d   c h e m i s t r y :   F i r e   s i z e ,   d i l u t i o n ,   a n d   b a c k g r o u n d   a e r o s o l   a l s o   g r e a t l y   i n f l u e n c e   n e a r - f i e l d   b i o m a s s   b u r n i n g   a e r o s o l   a g i n g ,   s u b m i t t e d   t o   J .   G e o p h y s .   R e s . ,   2 0 1 8 . 
 
 
 
 A d d i t i o n a l   c i t a t i o n s : 
 
 B i a n ,   Q . ,   J a t h a r ,   S .   H . ,   K o d r o s ,   J .   K . ,   B a r s a n t i ,   K .   C . ,   H a t c h ,   L .   E . ,   M a y ,   A .   A . ,   e t   a l .   ( 2 0 1 7 ) .   S e c o n d a r y   o r g a n i c   a e r o s o l   f o r m a t i o n   i n   b i o m a s s - b u r n i n g   p l u m e s :   T h e o r e t i c a l   a n a l y s i s   o f   l a b   s t u d i e s   a n d   a m b i e n t   p l u m e s .   A t m o s p h e r i c   C h e m i s t r y   a n d   P h y s i c s ,   1 7 ( 8 ) ,   5 4 5 9 - 5 4 7 5 .   h t t p s : / / d o i . o r g / 1 0 . 5 1 9 4 / a c p - 1 7 - 5 4 5 9 - 2 0 1 7 
 
 
 
 
 
 