Page 1 Source Listing SOLARD 2014-12-17 20:47 /tmpnwprd/iforthUcXhu.i 1 # 1 "SOLARD.F" 2 C&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& 3 SUBROUTINE SOLARD(R1) 4 C ****************************************************************** 5 C$$$ SUBPROGRAM DOCUMENTATION BLOCK 6 C . . . . 7 C SUBPROGRAM: SOLARD COMPUTE THE SOLAR-EARTH DISTANCE 8 C PRGRMMR: Q.ZHAO ORG: W/NMC2 DATE: 96-7-23 9 C 10 C ABSTRACT: 11 C SOLARD CALCULATES THE SOLAR-EARTH DISTANCE ON EACH DAY 12 C FOR USE IN SHORT-WAVE RADIATION. 13 C 14 C PROGRAM HISTORY LOG: 15 C 96-07-23 Q.ZHAO - ORIGINATOR 16 C 98-10-09 Q.ZHAO - CHANGED TO USE IW3JDN IN W3LIB TO 17 C CALCULATE JD. 18 C 19 C USAGE: CALL SOLARD FROM SUBROUTINE INIT 20 C 21 C INPUT ARGUMENT LIST: 22 C NONE 23 C 24 C OUTPUT ARGUMENT LIST: 25 C R1 - THE NON-DIMENSIONAL DISTANCE BETWEEN SUN AND THE EARTH 26 C (LESS THAN 1.0 IN SUMMER AND LARGER THAN 1.0 IN WINTER). 27 C 28 C INPUT FILES: 29 C NONE 30 C 31 C OUTPUT FILES: 32 C NONE 33 C 34 C SUBPROGRAMS CALLED: 35 C 36 C UNIQUE: NONE 37 C 38 C LIBRARY: IW3JDN 39 C 40 C COMMON BLOCKS: CTLBLK 41 C 42 C ATTRIBUTES: 43 C LANGUAGE: FORTRAN 90 44 C MACHINE : IBM SP 45 C$$$ 46 # 1 "./sp.h" 1 47 # 4 48 49 # 46 "SOLARD.F" 2 50 C*********************************************************************** 51 P A R A M E T E R 52 & (PI=3.1415926,PI2=2.*PI) 53 C----------------------------------------------------------------------- 54 L O G I C A L 55 & RUN,FIRST,RESTRT,SIGMA,LEAP 56 C----------------------------------------------------------------------- 57 INCLUDE "CTLBLK.comm" Page 2 Source Listing SOLARD 2014-12-17 20:47 SOLARD.F 58 C----------------------------------------------------------------------- 68 D I M E N S I O N 69 1 NDM (12) 70 C----------------------------------------------------------------------- 71 D A T A 72 1 JYR19/1900/, JMN/0/, CCR/1.3E-6/ 73 C 74 D A T A 75 1 NDM/0,31,59,90,120,151,181,212,243,273,304,334/ 76 C 77 C.....TPP = DAYS BETWEEN EPOCH AND PERIHELION PASSAGE OF 1900 78 C.....JDOR1 = JD OF DECEMBER 30, 1899 AT 12 HOURS UT 79 C.....JDOR2 = JD OF EPOCH WHICH IS JANUARY 0, 1990 AT 12 HOURS UT 80 C 81 D A T A 82 1 TPP/1.55/, JDOR1/2415019/, JDOR2/2415020/ 83 C 84 C ******************************************************************* 85 C COMPUTES JULIAN DAY AND FRACTION FROM YEAR, MONTH, DAY AND TIME UT 86 C ACCURATE ONLY BETWEEN MARCH 1, 1900 AND FEBRUARY 28, 2100 87 C BASED ON JULIAN CALENDAR CORRECTED TO CORRESPOND TO GREGORIAN 88 C CALENDAR DURING THIS PERIOD 89 C ******************************************************************* 90 JYR=IDAT(3) 91 JMNTH=IDAT(1) 92 JDAY=IDAT(2) 93 JHR=IHRST 94 C 95 JD=IW3JDN(JYR,JMNTH,JDAY) 96 IF(JHR.GE.12) THEN 97 JD=JD-1 98 FJD=.5E0+.041666667E0*FLOAT(JHR)+.00069444444E0*FLOAT(JMN) 99 ELSE 100 7 FJD=.041666667E0*FLOAT(JHR-12)+.00069444444E0*FLOAT(JMN) 101 END IF 102 DAYINC=JHR/24.0 103 JD=JD+FJD+DAYINC 104 FJD=JD+FJD+DAYINC-JD 105 C*** 106 C*** CALCULATE THE SOLAR-EARTH DISTANCE 107 C*** 108 DAT=FLOAT(JD-JDOR2)-TPP+FJD 109 C*** 110 C COMPUTES TIME IN JULIAN CENTURIES AFTER EPOCH 111 C*** 112 T=FLOAT(JD-JDOR2)/36525.E0 113 C*** 114 C COMPUTES LENGTH OF ANOMALISTIC AND TROPICAL YEARS (MINUS 365 DAYS) 115 C*** 116 YEAR=.25964134E0+.304E-5*T 117 C*** 118 C COMPUTES ORBIT ECCENTRICITY FROM T 119 C*** 120 EC=.01675104E0-(.418E-4+.126E-6*T)*T 121 YEAR=YEAR+365.E0 122 C*** 123 C DATE=DAYS SINCE LAST PERIHELION PASSAGE Page 3 Source Listing SOLARD 2014-12-17 20:47 SOLARD.F 124 C*** 125 DATE = MOD(DAT,YEAR) 126 C*** 127 C SOLVE ORBIT EQUATIONS BY NEWTON'S METHOD 128 C*** 129 EM=PI2*DATE/YEAR 130 E=1.E0 131 ITER = 0 132 31 EP=E-(E-EC*SIN(E)-EM)/(1.E0-EC*COS(E)) 133 CR=ABS(E-EP) 134 E=EP 135 ITER = ITER + 1 136 IF(ITER.GT.10) GOTO 1031 137 IF(CR.GT.CCR) GO TO 31 138 1031 CONTINUE 139 R1=1.E0-EC*COS(E) 140 C 141 WRITE(6,1000)JYR,JMNTH,JDAY,JHR,R1 142 1000 FORMAT('SUN-EARTH DISTANCE CALCULATION FINISHED IN SOLARD'/ 143 & 'YEAR=',I5,' MONTH=',I3,' DAY=',I3,' HOUR=' 144 &, I3,' R1=',F9.4) 145 C*** 146 C RETURN TO RADTN 147 C*** 148 RETURN 149 END ENTRY POINTS Name solard_ Page 4 Source Listing SOLARD 2014-12-17 20:47 Symbol Table SOLARD.F SYMBOL CROSS REFERENCE Name Object Declared Type Bytes Dimen Elements Attributes References 1000 Label 129 128 1031 Label 125 123 31 Label 119 124 7 Label 87 ABS Func 120 scalar 120 CCR Local 59 R(4) 4 scalar 59,124 COS Func 119 scalar 119,126 CR Local 120 R(4) 4 scalar 120,124 CTLBLK Common 4 112 DAT Local 95 R(4) 4 scalar 95,112 DATE Local 112 R(4) 4 scalar 112,116 DAYINC Local 89 R(4) 4 scalar 89,90,91 DT Scalar 6 R(4) 4 scalar COM E Local 117 R(4) 4 scalar 117,119,120,121,126 EC Local 107 R(4) 4 scalar 107,119,126 EM Local 116 R(4) 4 scalar 116,119 EP Local 119 R(4) 4 scalar 119,120,121 FIRST Scalar 51 L(4) 4 scalar COM FJD Local 85 R(4) 4 scalar 85,87,90,91,95 FLOAT Func 85 scalar 85,87,95,99 IDAT Scalar 5 I(4) 4 1 3 COM 77,78,79 IDTAD Scalar 6 I(4) 4 scalar COM IHRST Scalar 5 I(4) 4 scalar COM 80 IOUT Scalar 5 I(4) 4 scalar COM ITER Local 118 I(4) 4 scalar 118,122,123 IW3JDN Func 82 I(4) 4 scalar 82 JD Local 82 I(4) 4 scalar 82,84,90,91,95,99 JDAY Local 79 I(4) 4 scalar 79,82,128 JDOR1 Local 69 I(4) 4 scalar 69 JDOR2 Local 69 I(4) 4 scalar 69,95,99 JHR Local 80 I(4) 4 scalar 80,83,85,87,89,128 JMN Local 59 I(4) 4 scalar 59,85,87 JMNTH Local 78 I(4) 4 scalar 78,82,128 JYR Local 77 I(4) 4 scalar 77,82,128 JYR19 Local 59 I(4) 4 scalar 59 LEAP Local 51 L(4) 4 scalar LIST Scalar 5 I(4) 4 scalar COM MOD Func 112 scalar 112 NBC Scalar 5 I(4) 4 scalar COM NBOCO Scalar 6 I(4) 4 scalar COM NCNVC Scalar 7 I(4) 4 scalar COM NCP Scalar 6 I(4) 4 scalar COM NDM Local 56 I(4) 4 1 12 62 NEST Scalar 7 I(4) 4 scalar COM NFCST Scalar 5 I(4) 4 scalar COM NPHS Scalar 7 I(4) 4 scalar COM NPREC Scalar 6 I(4) 4 scalar COM NRADL Scalar 7 I(4) 4 scalar COM NRADS Scalar 7 I(4) 4 scalar COM NSHDE Scalar 6 I(4) 4 scalar COM NSTART Scalar 6 I(4) 4 scalar COM Page 5 Source Listing SOLARD 2014-12-17 20:47 Symbol Table SOLARD.F Name Object Declared Type Bytes Dimen Elements Attributes References NTDDMP Scalar 6 I(4) 4 scalar COM NTSD Scalar 6 I(4) 4 scalar COM NTSTM Scalar 6 I(4) 4 scalar COM PI Param 48 R(4) 4 scalar 48 PI2 Param 48 R(4) 4 scalar 116 R1 Dummy 2 R(4) 4 scalar ARG,INOUT 126,128 RESTRT Scalar 51 L(4) 4 scalar COM RUN Scalar 51 L(4) 4 scalar COM SIGMA Scalar 51 L(4) 4 scalar COM SIN Func 119 scalar 119 SINGLRST Local 2 L(4) 4 scalar SOLARD Subr 2 SUBPOST Scalar 2 L(4) 4 scalar COM T Local 99 R(4) 4 scalar 99,103,107 TPP Local 69 R(4) 4 scalar 69,95 YEAR Local 103 R(4) 4 scalar 103,108,112,116 Page 6 Source Listing SOLARD 2014-12-17 20:47 Subprograms/Common Blocks SOLARD.F SUBPROGRAMS/COMMON BLOCKS Name Object Declared Type Bytes Dimen Elements Attributes References CTLBLK Common 4 112 SOLARD Subr 2 COMPILER OPTIONS BEING USED -align nocommons -align nodcommons -align noqcommons -align records -align nosequence -align norec1byte -align norec2byte -align norec4byte -align norec8byte -align norec16byte -altparam -assume accuracy_sensitive -assume nobscc -assume nobuffered_io -assume byterecl -assume nocc_omp -assume nocstring -assume nodummy_aliases -assume nofpe_summary -assume noieee_fpe_flags -assume nominus0 -assume noold_boz -assume old_unit_star -assume old_ldout_format -assume noold_logical_ldio -assume old_maxminloc -assume old_xor -assume protect_constants -assume noprotect_parens -assume split_common -assume source_include -assume nostd_intent_in -assume nostd_mod_proc_name -assume norealloc_lhs -assume underscore -assume no2underscores no -auto -auto_scalar no -bintext -ccdefault default -check noargs -check noarg_temp_created -check nobounds -check noformat -check nooutput_conversion -check nooverflow -check nopointers -check power -check noshape -check nounderflow -check nouninitialized -coarray-num-procs 0 no -coarray-config-file -convert big_endian -cross_reference -D __INTEL_COMPILER=1210 -D __unix__ -D __unix -D __linux__ -D __linux -D __gnu_linux__ -D unix -D linux -D __ELF__ -D __x86_64 -D __x86_64__ -D _MT -D __INTEL_COMPILER_BUILD_DATE=20120612 -D __i686 -D __i686__ -D __pentiumpro -D __pentiumpro__ -D __pentium4 -D __pentium4__ -D __tune_pentium4__ -D __SSE2__ -D __SSE__ -D __MMX__ -double_size 64 no -d_lines no -Qdyncom -error_limit 30 no -f66 no -f77rtl no -fast -fpscomp nofilesfromcmd -fpscomp nogeneral -fpscomp noioformat -fpscomp noldio_spacing -fpscomp nologicals -fixed no -fpconstant Page 7 Source Listing SOLARD 2014-12-17 20:47 SOLARD.F -fpe3 -fprm nearest no -ftz -fp_model precise -fp_model nofast -fp_model nostrict -fp_model nosource -fp_model nodouble -fp_model noextended -fp_model novery_fast -fp_model noexcept -fp_model nono_except -fp_modbits nofp_contract -fp_modbits nono_fp_contract -fp_modbits nofenv_access -fp_modbits nono_fenv_access -fp_modbits nocx_limited_range -fp_modbits nono_cx_limited_range -fp_modbits noprec_div -fp_modbits nono_prec_div -fp_modbits noprec_sqrt -fp_modbits nono_prec_sqrt -fp_modbits noftz -fp_modbits no_ftz -fp_modbits nointrin_limited_range -fp_modbits nono_intrin_limited_range -fp_modbits notrunc_compares -fp_modbits nono_trunc_compares -fp_modbits noieee_nan_compares -fp_modbits nono_ieee_nan_compares -fp_modbits nohonor_f32_conversion -fp_modbits nono_honor_f32_conversion -fp_modbits nohonor_f64_conversion -fp_modbits nono_honor_f64_conversion -fp_modbits nono_x87_copy -fp_modbits nono_no_x87_copy -fp_modbits noexception_semantics -fp_modbits nono_exception_semantics -fp_modbits noprecise_libm_functions -fp_modbits nono_precise_libm_functions -heap_arrays 0 no -threadprivate_compat -g2 -iface nomixed_str_len_arg -iface nono_mixed_str_len_arg no -intconstant -integer_size 32 no -mixed_str_len_arg no -module -names lowercase no -noinclude -O2 no -pad_source -real_size 32 no -recursive -reentrancy none no -sharable_localsaves -vec=simd -show noinclude -show map -show options no -syntax_only no -threadcom no -U no -vms -w noall -w nonone -w alignments -w noargument_checking -w nodeclarations -w general -w noignore_bounds -w noignore_loc -w nointerfaces -w notruncated_source -w uncalled -w uninitialized -w nounused -w usage -includepath : /usrx/local/intel/composerxe/tbb/include/,/usr/include/,./,/opt/ibmhpc/pe1308/mpich2/intel/include64/, /opt/ibmhpc/pe1308/base/include64/,/usrx/local/intel/composerxe/mkl/include/,/usrx/local/intel/composerxe/tbb/include/, /gpfs/tp1/usrx/local/intel/composer_xe_2011_sp1.11.339/compiler/include/intel64/,/gpfs/tp1/usrx/local/intel/composer_xe_2011_sp1.11.339/compiler/include/, /usr/local/include/,/usr/lib/gcc/x86_64-redhat-linux/4.4.7/include/,/usr/include/,/usr/include/ -list filename : SOLARD.lst -o filename : none COMPILER: Intel(R) Fortran 12.1-2100