source: palm/trunk/LIB/rrtmg/rrlw_ncpar.f90 @ 4874

Last change on this file since 4874 was 4858, checked in by raasch, 4 years ago

rrtmg input file names changed, default value for unscheduled calls of the radiation model set to FALSE

File size: 1.8 KB
RevLine 
[1585]1module rrlw_ncpar
2        use parkind ,only : im => kind_im, rb => kind_rb
3
4        implicit none
5        save
6       
[4858]7        character(len=100) ::  rrtmg_lw_input_file_name = 'RRTMG_LW'
8
[1585]9        real(kind=rb), parameter :: cpdair = 1003.5  ! Specific heat capacity of dry air
10                                                         ! at constant pressure at 273 K
11                                                         ! (J kg-1 K-1)
12
13       
14        integer(kind=im), parameter :: maxAbsorberNameLength = 5, &
15                          Absorber              = 12
16    character(len = maxAbsorberNameLength), dimension(Absorber), parameter :: &
17    AbsorberNames = (/        &
18                                'N2   ',  &
19                                'CCL4 ',  &
20                                'CFC11',  &
21                                'CFC12',  &
22                                'CFC22',  &
23                                'H2O  ',  &
24                                'CO2  ',  &
25                                'O3   ',  &
26                                'N2O  ',  & 
27                                'CO   ',  &
28                                'CH4  ',  &
29                                'O2   '  /)
30       
31        integer(kind=im), dimension(40) :: status
32        integer(kind=im) :: i
33        integer(kind=im), parameter :: keylower  = 9,   &
34                                                  keyupper  = 5,   &
35                                                  Tdiff     = 5,   &
36                                                  ps        = 59,  &
37                                                  plower    = 13,  &
38                                                  pupper    = 47,  &
39                                                  Tself     = 10,  &
40                                                  Tforeign  = 4,   &
41                                                  pforeign  = 4,   &
42                                                  T         = 19,  &
43                                                  Tplanck   = 181, &
44                                                  band      = 16,  &
45                                                  GPoint    = 16,  &
46                                                  GPointSet = 2
47                                                 
48        contains
49       
50        subroutine getAbsorberIndex(AbsorberName,AbsorberIndex)
51                character(len = *), intent(in) :: AbsorberName
52                integer(kind=im), intent(out)           :: AbsorberIndex
53               
54                integer(kind=im) :: m
55       
56                AbsorberIndex = -1
57                do m = 1, Absorber
58                        if (trim(AbsorberNames(m)) == trim(AbsorberName)) then
59                                AbsorberIndex = m
60                        end if
61                end do
62               
63                if (AbsorberIndex == -1) then
64                        print*, "Absorber name index lookup failed."
65                end if
66        end subroutine getAbsorberIndex
67
68end module rrlw_ncpar
Note: See TracBrowser for help on using the repository browser.