pkg/radtrans/radtrans_sunmod.F

C $Header$
C $Name$

#include "RADTRANS_OPTIONS.h"

CBOP
C     !ROUTINE: RADTRANS_SUNMOD

C     !INTERFACE: ======================================================
      subroutine radtrans_sunmod(rad,iday,imon,iyr,gmt,up,no,ea,
     O                           sunz,rs)

C     !DESCRIPTION:
c MODIFIED only by separating subroutines that were previously included
c
c  Given year, day of year, time in hours (GMT) and latitude and 
c  longitude, returns an accurate solar zenith and azimuth angle.  
c  Based on IAU 1976 Earth ellipsoid.  Method for computing solar 
c  vector and local vertical from Patt and Gregg, 1994, Int. J. 
c  Remote Sensing.  Only outputs solar zenith angle.  This version
c  utilizes a pre-calculation of the local up, north, and east 
c  vectors, since the locations where the solar zenith angle are 
c  calculated in the model are fixed.
c
c  Subroutines required: sun2000
c                        gha2000
c                        jd
c
C     !USES: ===========================================================
      IMPLICIT NONE

C     !INPUT PARAMETERS: ===============================================
      INTEGER iday, imon, iyr
      _RL rad, gmt, up(3), no(3), ea(3)
c     INTEGER myThid

C     !OUTPUT PARAMETERS: ==============================================
      _RL sunz, rs

C     !FUNCTIONS: ======================================================

C     !LOCAL VARIABLES: ================================================
      INTEGER j
      _RL radeg, sec, day, gha, ghar, sunv, sunn, sune
      _RL suni(3),sung(3)
CEOP


      radeg = rad
c  Compute sun vector
c   Compute unit sun vector in geocentric inertial coordinates
      sec = gmt*3600.0D0
      call radtrans_sun2000 (radeg, iyr, imon, iday, sec, suni, rs)

c   Get Greenwich mean sidereal angle
      day = float(iday) + sec/86400.0D0
      call radtrans_gha2000 (radeg, iyr, imon, day, gha)
      ghar = gha/radeg

c   Transform Sun vector into geocentric rotating frame
      sung(1) = suni(1)*cos(ghar) + suni(2)*sin(ghar)
      sung(2) = suni(2)*cos(ghar) - suni(1)*sin(ghar)
      sung(3) = suni(3)
c
c  Compute components of spacecraft and sun vector in the
c  vertical (up), North (no), and East (ea) vectors frame
      sunv = 0.0
      sunn = 0.0
      sune = 0.0
      do j = 1,3
       sunv = sunv + sung(j)*up(j)
       sunn = sunn + sung(j)*no(j)
       sune = sune + sung(j)*ea(j)
      enddo
c
c  Compute the solar zenith 
      sunz = radeg*atan2(sqrt(sunn*sunn+sune*sune),sunv)
c
      return
      end
c
c *****************************************************************
c HERE SUBROUTINE sun2000 removed
c *********************************************************************
c HERE SUBROUTINE gha2000 removed
c HERE SUBROUTINE ephparms removed
c HERE SUBROUTINE nutate removed
c HERE FUNCTION jd removed
c ************************************************************************
1	jahn	1.1	C $Header$
2			C $Name$
3
4			#include "RADTRANS_OPTIONS.h"
5
6			CBOP
7			C !ROUTINE: RADTRANS_SUNMOD
8
9			C !INTERFACE: ======================================================
10			subroutine radtrans_sunmod(rad,iday,imon,iyr,gmt,up,no,ea,
11			O sunz,rs)
12
13			C !DESCRIPTION:
14			c MODIFIED only by separating subroutines that were previously included
15			c
16			c Given year, day of year, time in hours (GMT) and latitude and
17			c longitude, returns an accurate solar zenith and azimuth angle.
18			c Based on IAU 1976 Earth ellipsoid. Method for computing solar
19			c vector and local vertical from Patt and Gregg, 1994, Int. J.
20			c Remote Sensing. Only outputs solar zenith angle. This version
21			c utilizes a pre-calculation of the local up, north, and east
22			c vectors, since the locations where the solar zenith angle are
23			c calculated in the model are fixed.
24			c
25			c Subroutines required: sun2000
26			c gha2000
27			c jd
28			c
29			C !USES: ===========================================================
30			IMPLICIT NONE
31
32			C !INPUT PARAMETERS: ===============================================
33			INTEGER iday, imon, iyr
34			_RL rad, gmt, up(3), no(3), ea(3)
35			c INTEGER myThid
36
37			C !OUTPUT PARAMETERS: ==============================================
38			_RL sunz, rs
39
40			C !FUNCTIONS: ======================================================
41
42			C !LOCAL VARIABLES: ================================================
43			INTEGER j
44			_RL radeg, sec, day, gha, ghar, sunv, sunn, sune
45			_RL suni(3),sung(3)
46			CEOP
47
48
49			radeg = rad
50			c Compute sun vector
51			c Compute unit sun vector in geocentric inertial coordinates
52			sec = gmt*3600.0D0
53			call radtrans_sun2000 (radeg, iyr, imon, iday, sec, suni, rs)
54
55			c Get Greenwich mean sidereal angle
56			day = float(iday) + sec/86400.0D0
57			call radtrans_gha2000 (radeg, iyr, imon, day, gha)
58			ghar = gha/radeg
59
60			c Transform Sun vector into geocentric rotating frame
61			sung(1) = suni(1)cos(ghar) + suni(2)sin(ghar)
62			sung(2) = suni(2)cos(ghar) - suni(1)sin(ghar)
63			sung(3) = suni(3)
64			c
65			c Compute components of spacecraft and sun vector in the
66			c vertical (up), North (no), and East (ea) vectors frame
67			sunv = 0.0
68			sunn = 0.0
69			sune = 0.0
70			do j = 1,3
71			sunv = sunv + sung(j)*up(j)
72			sunn = sunn + sung(j)*no(j)
73			sune = sune + sung(j)*ea(j)
74			enddo
75			c
76			c Compute the solar zenith
77			sunz = radegatan2(sqrt(sunnsunn+sune*sune),sunv)
78			c
79			return
80			end
81			c
82			c *****************************************************************
83			c HERE SUBROUTINE sun2000 removed
84			c *********************************************************************
85			c HERE SUBROUTINE gha2000 removed
86			c HERE SUBROUTINE ephparms removed
87			c HERE SUBROUTINE nutate removed
88			c HERE FUNCTION jd removed
89			c ************************************************************************