source: CONFIG/publications/ICOLMDZORINCA_CO2_Transport_GMD_2023/INCA/build/ppsrc/INCA_PP/setrxt.f90 @ 6610

Last change on this file since 6610 was 6610, checked in by acosce, 10 months ago

INCA used for ICOLMDZORINCA_CO2_Transport_GMD_2023

File size: 1.9 KB
Line 
1
2
3
4
5
6
7
8
9
10
11
12                                                                       
13      subroutine SETRXT( rate, temp )                                   
14! Stacy Walters, NCAR, 1998.                                           
15                                                                       
16      use INCA_dim                                                     
17      implicit none                                                     
18                                                                       
19!-------------------------------------------------------               
20!       ... Dummy args                                                 
21!-------------------------------------------------------               
22      real, intent(out) ::     rate(PLON,PLEV,12)                   
23      real, intent(in)  ::     temp(PLON,PLEV)                         
24                                                                       
25      real  ::  itemp(PLON,PLEV)                                       
26                                                                       
27      rate(:PLON,:,4) = 2.100e-06                                       
28      rate(:PLON,:,6) = 1.514e-07                                       
29      rate(:PLON,:,7) = 6.7e-11                                         
30      rate(:PLON,:,8) = 4.4e-11                                         
31      rate(:PLON,:,10) = 2.25e-10                                       
32      rate(:PLON,:,11) = 3.75e-11                                       
33      rate(:PLON,:,12) = 4.98e-11                                       
34      itemp(:PLON,:) = 1. / temp(:PLON,:)                               
35      rate(:PLON,:,5) = 1.800e-12 * EXP( -1550. * itemp(:PLON,:) )     
36      rate(:PLON,:,9) = 2.45e-12 * EXP( -1775. * itemp(:PLON,:) )       
37                                                                       
38      end subroutine SETRXT                                             
Note: See TracBrowser for help on using the repository browser.