1 |
jahn |
1.1 |
C $Header$ |
2 |
|
|
C $Name$ |
3 |
|
|
|
4 |
|
|
#include "RADTRANS_OPTIONS.h" |
5 |
|
|
|
6 |
|
|
CBOP |
7 |
|
|
C !ROUTINE: RADTRANS_SFCSOLZ |
8 |
|
|
|
9 |
|
|
C !INTERFACE: ====================================================== |
10 |
|
|
subroutine radtrans_sfcsolz(rad,iyr,imon,iday,isec,xlon,ylat, |
11 |
|
|
O solz) |
12 |
|
|
|
13 |
|
|
C !DESCRIPTION: |
14 |
|
|
c SLIGHTLY MODIFIED FORM WG sfsolz |
15 |
|
|
c original version contained subroutine code for |
16 |
|
|
c localvec |
17 |
|
|
c which is here separated into separate file |
18 |
|
|
c |
19 |
|
|
c Computes solar zenith angle at surface given location |
20 |
|
|
c and time. |
21 |
|
|
c |
22 |
|
|
|
23 |
|
|
C !USES: =========================================================== |
24 |
|
|
IMPLICIT NONE |
25 |
|
|
|
26 |
|
|
C !INPUT PARAMETERS: =============================================== |
27 |
|
|
INTEGER iyr, imon, iday, isec |
28 |
|
|
_RL rad, xlon, ylat |
29 |
|
|
c INTEGER myThid |
30 |
|
|
|
31 |
|
|
C !OUTPUT PARAMETERS: ============================================== |
32 |
|
|
_RL solz |
33 |
|
|
|
34 |
|
|
C !FUNCTIONS: ====================================================== |
35 |
|
|
|
36 |
|
|
C !LOCAL VARIABLES: ================================================ |
37 |
|
|
|
38 |
|
|
c |
39 |
|
|
c ifst is not used, so don't need this ??? |
40 |
|
|
c save |
41 |
|
|
_RL up(3),no(3),ea(3),gmt,rs |
42 |
|
|
c data ifst /0/ |
43 |
|
|
CEOP |
44 |
|
|
|
45 |
|
|
|
46 |
|
|
c Local direction vector |
47 |
|
|
call radtrans_localvec(rad,xlon,ylat,up,no,ea) |
48 |
|
|
c |
49 |
|
|
c Solar zenith angle |
50 |
|
|
solz = -1.0 |
51 |
|
|
gmt = float(isec)/3600. |
52 |
|
|
call radtrans_sunmod(rad,iday,imon,iyr,gmt,up,no,ea,solz,rs) |
53 |
|
|
c |
54 |
|
|
return |
55 |
|
|
end |
56 |
|
|
c |
57 |
|
|
c ********************************************************************** |
58 |
|
|
c HERE SUBROUTINE localvec REMOVED AND PUT IN SEP FILE |