2 integer :: ncol_
,nlayers_
,nbndlw_
,ngptlw_
3 ! changed to arguments
for thread safety
7 integer :: ngsd(nbndlw
)
10 real :: taucmcd(ncol_
, ngptlw_
, nlayers_
+1)
12 real
, dimension(ncol_
, 0:nlayers_
+1) :: pzd
! level (interface
) pressures (hPa
, mb
)
13 ! Dimensions
: (ncol
,0:nlayers
)
14 real
, dimension(ncol_
) :: pwvcmd
! precipitable water
vapor (cm
)
16 real
, dimension(ncol_
,nbndlw_
) :: semissd
! lw surface emissivity
17 ! Dimensions
: (ncol
,nbndlw
)
18 real
, dimension(ncol_
,nlayers_
+1,nbndlw_
) :: planklayd
!
19 ! Dimensions
: (ncol
,nlayers
+1,nbndlw
)
20 real
, dimension(ncol_
,0:nlayers_
+1,nbndlw_
) :: planklevd
!
21 ! Dimensions
: (ncol
,0:nlayers
+1,nbndlw
)
22 real
, dimension(ncol_
,nbndlw_
) :: plankbndd
!
23 ! Dimensions
: (ncol
,nbndlw
)
25 real :: gurad(ncol_
,ngptlw_
,0:nlayers_
+1) ! upward longwave
flux (w
/m2
)
26 real :: gdrad(ncol_
,ngptlw_
,0:nlayers_
+1) ! downward longwave
flux (w
/m2
)
27 real :: gclrurad(ncol_
,ngptlw_
,0:nlayers_
+1) ! clear sky upward longwave
flux (w
/m2
)
28 real :: gclrdrad(ncol_
,ngptlw_
,0:nlayers_
+1) ! clear sky downward longwave
flux (w
/m2
)
30 real :: gdtotuflux_dtd(ncol_
, ngptlw_
, 0:nlayers_
+1) ! change in upward longwave
flux (w
/m1
/k
)
31 ! with respect to surface temperature
33 real :: gdtotuclfl_dtd(ncol_
, ngptlw_
, 0:nlayers_
+1) ! change in clear sky upward longwave
flux (w
/m2
/k
)
34 ! with respect to surface temperature
37 integer :: idrvd
! flag
for calculation of dF
/dt from
38 ! Planck derivative
[0=off
, 1=on
]
42 real :: a0d(nbndlw_
), a1d(nbndlw_
), a2d(nbndlw_
)
43 real :: delwaved(nbndlw_
)
44 real :: totufluxd(ncol_
, 0:nlayers_
+1) ! upward longwave
flux (w
/m2
)
45 real :: totdfluxd(ncol_
, 0:nlayers_
+1) ! downward longwave
flux (w
/m2
)
46 real :: fnetd(ncol_
, 0:nlayers_
+1) ! net longwave
flux (w
/m2
)
47 real :: htrd(ncol_
, 0:nlayers_
+1) ! longwave heating
rate (k
/day
)
48 real :: totuclfld(ncol_
, 0:nlayers_
+1) ! clear sky upward longwave
flux (w
/m2
)
49 real :: totdclfld(ncol_
, 0:nlayers_
+1) ! clear sky downward longwave
flux (w
/m2
)
50 real :: fnetcd(ncol_
, 0:nlayers_
+1) ! clear sky net longwave
flux (w
/m2
)
51 real :: htrcd(ncol_
, 0:nlayers_
+1) ! clear sky longwave heating
rate (k
/day
)
52 real :: dtotuflux_dtd(ncol_
, 0:nlayers_
+1) ! change in upward longwave
flux (w
/m2
/k
)
53 ! with respect to surface temperature
54 real :: dtotuclfl_dtd(ncol_
, 0:nlayers_
+1) ! change in clear sky upward longwave
flux (w
/m2
/k
)
55 ! with respect to surface temperature
56 real :: dplankbnd_dtd(ncol_
,nbndlw_
)