| 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 |