/[MITgcm]/MITgcm_contrib/darwin/pkg/radtrans/radtrans_ephparms.F
ViewVC logotype

Annotation of /MITgcm_contrib/darwin/pkg/radtrans/radtrans_ephparms.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph

Revision 1.1 - (hide annotations) (download)
Wed Jun 9 15:59:37 2010 UTC (15 years, 1 month ago) by jahn
Branch: MAIN
CVS Tags: ctrb_darwin_ckpt63j_20120217, ctrb_darwin_ckpt62w_20110426, ctrb_darwin_ckpt62k_20101005, ctrb_darwin_ckpt63m_20120506, ctrb_darwin_ckpt66l_20171025, ctrb_darwin_ckpt66a_20161020, ctrb_darwin_ckpt62l_20101006, ctrb_darwin_ckpt65n_20150729, ctrb_darwin_ckpt64v_20140411, ctrb_darwin_ckpt62q_20110124, ctrb_darwin_ckpt65b_20140812, ctrb_darwin_ckpt62m_20101024, ctrb_darwin_ckpt63s_20120908, ctrb_darwin_ckpt62z_20110622, ctrb_darwin_ckpt65y_20160801, ctrb_darwin_ckpt66o_20180209, ctrb_darwin_ckpt63q_20120731, ctrb_darwin_ckpt62s_20110217, ctrb_darwin_ckpt63k_20120317, ctrb_darwin_ckpt64y_20140622, ctrb_darwin_ckpt64b_20121224, ctrb_darwin_ckpt64g_20130503, ctrb_darwin_ckpt65c_20140830, ctrb_darwin_ckpt62o_20101125, ctrb_darwin_ckpt64r_20131210, ctrb_darwin_ckpt65x_20160612, ctrb_darwin_ckpt66j_20170815, ctrb_darwin_ckpt64j_20130704, ctrb_darwin_ckpt65u_20160315, ctrb_darwin_ckpt62i_20100729, ctrb_darwin_ckpt62t_20110305, ctrb_darwin_ckpt62h_20100707, ctrb_darwin_ckpt64t_20140202, ctrb_darwin_ckpt64_20121012, ctrb_darwin_ckpt64s_20140105, ctrb_darwin_ckpt66k_20171025, ctrb_darwin_ckpt63i_20120124, ctrb_darwin_ckpt64h_20130528, ctrb_darwin_ckpt64u_20140308, ctrb_darwin_ckpt64k_20130723, ctrb_darwin_ckpt63g_20111220, ctrb_darwin_ckpt63b_20110830, ctrb_darwin_ckpt65s_20160114, ctrb_darwin_ckpt64a_20121116, ctrb_darwin_ckpt62p_20101217, ctrb_darwin_ckpt63_20110728, ctrb_darwin_ckpt62n_20101109, ctrb_darwin_ckpt62r_20110124, ctrb_darwin_ckpt62v_20110412, ctrb_darwin_ckpt64i_20130622, ctrb_darwin_ckpt65h_20141217, ctrb_darwin_ckpt65a_20140728, ctrb_darwin_ckpt63o_20120629, ctrb_darwin_ckpt63a_20110804, ctrb_darwin_ckpt66i_20170718, ctrb_darwin_ckpt62u_20110405, ctrb_darwin_ckpt66n_20180118, ctrb_darwin_ckpt65l_20150504, ctrb_darwin_ckpt63c_20111011, ctrb_darwin_ckpt65o_20150914, ctrb_darwin_ckpt62v_20110406, ctrb_darwin_ckpt65f_20141014, ctrb_darwin_ckpt63p_20120707, ctrb_darwin_ckpt64q_20131118, ctrb_darwin_ckpt65i_20150123, ctrb_darwin_ckpt63d_20111107, ctrb_darwin_ckpt66b_20161219, ctrb_darwin_ckpt65t_20160221, ctrb_darwin_ckpt65k_20150402, ctrb_darwin_ckpt64p_20131024, ctrb_darwin_ckpt65e_20140929, ctrb_darwin_ckpt64c_20130120, ctrb_darwin_ckpt64q_20131024, ctrb_darwin_ckpt65m_20150615, ctrb_darwin_ckpt64m_20130820, ctrb_darwin_ckpt66f_20170407, ctrb_darwin_ckpt63r_20120817, ctrb_darwin_ckpt65w_20160512, ctrb_darwin_ckpt62x_20110513, ctrb_darwin_ckpt65d_20140915, ctrb_darwin_ckpt63f_20111201, ctrb_darwin_ckpt64o_20131024, ctrb_darwin_ckpt64f_20130405, ctrb_darwin_ckpt65r_20151221, ctrb_darwin_ckpt65q_20151118, ctrb_darwin_ckpt66g_20170424, ctrb_darwin_ckpt65j_20150225, ctrb_darwin_ckpt64w_20140502, ctrb_darwin_ckpt62j_20100901, ctrb_darwin_ckpt63e_20111107, ctrb_darwin_ckpt64d_20130219, ctrb_darwin_ckpt63l_20120405, ctrb_darwin_ckpt66c_20170121, ctrb_darwin_ckpt64x_20140524, ctrb_darwin_ckpt66h_20170602, ctrb_darwin_ckpt66d_20170214, ctrb_darwin_ckpt63h_20111230, ctrb_darwin_ckpt64z_20140711, ctrb_darwin_ckpt65v_20160409, ctrb_darwin_ckpt65p_20151023, ctrb_darwin_ckpt65g_20141120, ctrb_darwin_ckpt63n_20120604, ctrb_darwin_ckpt64e_20130305, ctrb_darwin_ckpt66m_20171213, ctrb_darwin_ckpt66e_20170314, ctrb_darwin_ckpt64l_20130806, ctrb_darwin_ckpt65z_20160929, ctrb_darwin_ckpt62y_20110526, ctrb_darwin_ckpt65_20140718, ctrb_darwin_ckpt64n_20130826, HEAD 
add radtrans package
1 jahn 1.1 C $Header$
2     C $Name$
3    
4     #include "RADTRANS_OPTIONS.h"
5    
6     CBOP
7     C !ROUTINE: RADTRANS_EPHPARMS
8    
9     C !INTERFACE: ==========================================================
10     subroutine radtrans_ephparms (t,
11     O xls, gs, xlm, omega)
12    
13     C !DESCRIPTION:
14     c This subroutine computes ephemeris parameters used by other Mission
15     c Operations routines: the solar mean longitude and mean anomaly, and
16     c the lunar mean longitude and mean ascending node. It uses the model
17     c referenced in The Astronomical Almanac for 1984, Section S
18     c (Supplement) and documented and documented in Exact closed-form
19     c geolocation algorithm for Earth survey sensors, by F.S. Patt and
20     c W.W. Gregg, Int. Journal of Remote Sensing, 1993. These parameters
21     c are used to compute the solar longitude and the nutation in
22     c longitude and obliquity.
23     c
24     c Program written by: Frederick S. Patt
25     c General Sciences Corporation
26     c November 2, 1992
27    
28     c Calling Arguments
29     C !USES: ===============================================================
30     IMPLICIT NONE
31    
32     C !INPUT PARAMETERS: ===================================================
33     c t :: Time in days since January 1, 2000 at 12 hours UT
34     _RL t
35     c INTEGER myThid
36    
37     C !OUTPUT PARAMETERS: ==================================================
38     c xls :: Mean solar longitude (degrees)
39     c gs :: Mean solar anomaly (degrees)
40     c xlm :: Mean lunar longitude (degrees)
41     c omega :: Ascending node of mean lunar orbit (degrees)
42     _RL xls, gs, xlm, omega
43     CEOP
44    
45     C !LOCAL VARIABLES: ====================================================
46     c
47    
48     c Sun Mean Longitude
49     xls = 280.46592D0 + 0.9856473516D0*t
50     xls = mod(xls,360.0)
51    
52     c Sun Mean Anomaly
53     gs = 357.52772D0 + 0.9856002831D0*t
54     gs = mod(gs,360.0)
55    
56     c Moon Mean Longitude
57     xlm = 218.31643D0 + 13.17639648D0*t
58     xlm = mod(xlm,360.0)
59    
60     c Ascending Node of Moons Mean Orbit
61     omega = 125.04452D0 - 0.0529537648D0*t
62     omega = mod(omega,360.0)
63    
64     return
65     end
66    
  ViewVC Help
Powered by ViewVC 1.1.22