1 |
jahn |
1.1 |
C $Header$ |
2 |
|
|
C $Name$ |
3 |
|
|
|
4 |
|
|
#include "RADTRANS_OPTIONS.h" |
5 |
|
|
|
6 |
|
|
CBOP |
7 |
|
|
C !ROUTINE: RADTRANS_SFCRMUD |
8 |
|
|
|
9 |
|
|
C !INTERFACE: ====================================================== |
10 |
|
|
subroutine radtrans_sfcrmud(rad,sunz, |
11 |
|
|
O rmud) |
12 |
|
|
|
13 |
|
|
C !DESCRIPTION: |
14 |
|
|
c Computes average cosine for direct irradiance. |
15 |
|
|
c |
16 |
|
|
C !USES: =========================================================== |
17 |
|
|
IMPLICIT NONE |
18 |
|
|
|
19 |
|
|
C !INPUT PARAMETERS: =============================================== |
20 |
|
|
_RL rad |
21 |
|
|
_RL sunz |
22 |
|
|
c INTEGER myThid |
23 |
|
|
|
24 |
|
|
C !OUTPUT PARAMETERS: ============================================== |
25 |
|
|
_RL rmud |
26 |
|
|
|
27 |
|
|
C !FUNCTIONS: ====================================================== |
28 |
|
|
|
29 |
|
|
C !LOCAL VARIABLES: ================================================ |
30 |
|
|
_RL rsza, sinszaw, szaw, rmudl, rn |
31 |
|
|
data rn /1.341 _d 0/ !refractive index of seawater |
32 |
|
|
CEOP |
33 |
|
|
|
34 |
|
|
c |
35 |
|
|
c Compute average cosine for direct irradiance in the water |
36 |
|
|
c column given solar zenith angle (in degrees) at surface. |
37 |
|
|
c |
38 |
|
|
rsza = sunz/rad |
39 |
|
|
sinszaw = sin(rsza)/rn |
40 |
|
|
szaw = asin(sinszaw) |
41 |
|
|
rmudl = 1.0/cos(szaw) !avg cosine direct (1 over) |
42 |
|
|
rmud = min(rmudl,1.5) |
43 |
|
|
rmud = max(rmud,0.0) |
44 |
|
|
c |
45 |
|
|
return |
46 |
|
|
end |