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 |