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