Page 1 Source Listing GRADFS 2014-12-17 20:46 /tmpnwprd/ifortYj0Mnf.i 1 # 1 "GRADFS.F" 2 C&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&&& 3 SUBROUTINE GRADFS 4 C ****************************************************************** 5 C * * 6 C * THIS SUBROUTINE WAS MODIFIED BY Q. ZHAO TO BE USED IN THE * 7 C * ETA MODEL. 1993. 11. 18. * 8 C * * 9 C ****************************************************************** 10 C*** 11 C*** REQUIRED INPUT 12 C*** 13 1 (SIGL,KCCO2,NFILE) 14 C************** 15 C* SIGL(LP1): MIDLAYER PRESSURES IN PA (LP1=LM+1) 16 C* KCCO2: =0 (NOT USED) 17 C* NFILE: THE FILE NUMBER FOR O3 DATA 18 C************** 19 INCLUDE "parmeta" 42 INCLUDE "rdparm" 43 # 1 "./sp.h" 1 44 # 4 45 46 # 21 "GRADFS.F" 2 47 C*********************************************************************** 79 PARAMETER (LNGTH=37*L) 80 COMMON /RDFSAV/ EMISP,EMIST,XLATT,XLATP,Q19001,HP98,H3M6, 81 * HP75,H6M2,HP537,H74E1,H15E1,Q14330,HP2,TWENTY,HNINE, 82 * DEGRAD,HSIGMA,DAYSEC,RCO2, 83 * CAO3SW(5),CAH2SW(5),CBSW(5) 84 C====> ************************************************************** 85 C-- SEASONAL CLIMATOLOGIES OF O3 (OBTAINED FROM A PREVIOUSLY RUN 86 C CODE WHICH INTERPOLATES O3 TO USER VERTICAL COORDINATE). 87 C DEFINED AS 5 DEG LAT MEANS N.P.->S.P. 88 COMMON /SAVMEM/ 89 C- ...WINTER.... ...SPRING.... ...SUMMER.... ....FALL..... 90 1 DDUO3N(37,L), DDO3N2(37,L), DDO3N3(37,L), DDO3N4(37,L) 91 DIMENSION RAD1(LNGTH), RAD2(LNGTH), RAD3(LNGTH), RAD4(LNGTH) 92 EQUIVALENCE (RAD1(1),DDUO3N(1,1)),(RAD2(1),DDO3N2(1,1)) 93 EQUIVALENCE (RAD3(1),DDO3N3(1,1)),(RAD4(1),DDO3N4(1,1)) 94 C====> ************************************************************** 95 DIMENSION SIGL(LP1) 96 CZHAO DIMENSION SIGL(L) 97 DIMENSION XAO3SW(5),XAH2SW(5),XBSW(5) 98 DATA 99 1 XAO3SW / 0., .690, .480, .210, 0./ , 100 2 XAH2SW / 0., .690, .480, .210, 0./ , 101 3 XBSW / 0., .035, .020, .005, 0./ 102 CCCCCCMIC$ GUARD 0004 103 C..************************** 104 C..************************** 105 C ******************************************************** 106 C * ONE TIME COMPUTATION OF NECESSARY QUANTITIES * 107 C ******************************************************** 108 C.... INITIALIZE ARRAYS,GET CONSTANTS,ETC... 109 PI=3.1415927 110 Q19001=19.001 Page 2 Source Listing GRADFS 2014-12-17 20:46 GRADFS.F 111 HP98=0.98 112 H3M6=3.0E-6 113 HP537=0.537 114 H74E1=74.0 115 H15E1=15.0 116 Q14330=1.43306E-6 117 HP2=0.2 118 TWENTY=20.0 119 HNINE=9.0 120 DEGRAD=180.0/PI 121 HSIGMA=5.673E-5 122 DAYSEC=1.1574E-5 123 C ATMOSPERIC CARBON DIOXIDE CONCENTRATION IS NOW READ BY CONRAD, 124 C BUT IT DEFAULTS TO 330 PPM FOR BACKWARD COMPATIBILITY. 125 RCO2=3.3E-4 126 CALL HCONST 127 C.... INTERPOLATE CLIMO O3 TO THE CURRENT VERTICAL COORDINATE... 128 C.. NEED LAYER SIGMA, GET FROM PSFC AND LAYER P FOR I=1..... 129 DO 3 I = 1 , 5 130 CAO3SW(I) = XAO3SW(I) 131 CAH2SW(I) = XAH2SW(I) 132 CBSW(I) = XBSW(I) 133 3 CONTINUE 134 C*** 135 C*** CONVERT SIGL FROM PA TO MB TO BE USED IN O3INT 136 C*** 137 DO 100 LV=1,LP1 138 SIGL(LV)=0.01*SIGL(LV) 139 100 CONTINUE 140 CALL O3INT(SIGL) 141 CALL CONRAD(NFILE) 142 C.... AVERAGE CLIMATOLOGICAL VALUS OF O3 FROM 5 DEG LAT MEANS, SO THAT 143 C TIME AND SPACE INTERPOLATION WILL WORK (DONE ELSEWHERE IN RADFS) 144 DO 5 I=1,LNGTH 145 AVG=.25E0*(RAD1(I)+RAD2(I)+RAD3(I)+RAD4(I)) 146 A1=.5E0*(RAD2(I)-RAD4(I)) 147 B1=.5E0*(RAD1(I)-RAD3(I)) 148 B2=.25E0*((RAD1(I)+RAD3(I))-(RAD2(I)+RAD4(I))) 149 RAD1(I)=AVG 150 RAD2(I)=A1 151 RAD3(I)=B1 152 RAD4(I)=B2 153 5 CONTINUE 154 EMIST = .6E0 155 EMISP = .3E0 156 XLATP = 60.E0 157 XLATT = 30.E0 158 C 159 RETURN 160 END Page 3 Source Listing GRADFS 2014-12-17 20:46 Entry Points GRADFS.F ENTRY POINTS Name gradfs_ SYMBOL CROSS REFERENCE Name Object Declared Type Bytes Dimen Elements Attributes References 100 Label 82 80 3 Label 76 72 5 Label 96 87 A1 Local 89 R(4) 4 scalar 89,93 AVG Local 88 R(4) 4 scalar 88,92 B1 Local 90 R(4) 4 scalar 90,94 B2 Local 91 R(4) 4 scalar 91,95 CAH2SW Scalar 26 R(4) 4 1 5 COM 74 CAO3SW Scalar 26 R(4) 4 1 5 COM 73 CBSW Scalar 26 R(4) 4 1 5 COM 75 CONRAD Subr 84 84 DAYSEC Scalar 25 R(4) 4 scalar COM 65 DDO3N2 Scalar 33 R(4) 4 2 1665 COM DDO3N3 Scalar 33 R(4) 4 2 1665 COM DDO3N4 Scalar 33 R(4) 4 2 1665 COM DDUO3N Scalar 33 R(4) 4 2 1665 COM DEGRAD Scalar 25 R(4) 4 scalar COM 63 EMISP Scalar 23 R(4) 4 scalar COM 98 EMIST Scalar 23 R(4) 4 scalar COM 97 GRADFS Subr 2 H15E1 Scalar 24 R(4) 4 scalar COM 58 H3M6 Scalar 23 R(4) 4 scalar COM 55 H6M2 Scalar 24 R(4) 4 scalar COM H74E1 Scalar 24 R(4) 4 scalar COM 57 HCONST Subr 69 69 HNINE Scalar 24 R(4) 4 scalar COM 62 HP2 Scalar 24 R(4) 4 scalar COM 60 HP537 Scalar 24 R(4) 4 scalar COM 56 HP75 Scalar 24 R(4) 4 scalar COM HP98 Scalar 23 R(4) 4 scalar COM 54 HSIGMA Scalar 25 R(4) 4 scalar COM 64 I Local 72 I(4) 4 scalar 72,73,74,75,87,88,89,90,91,92,93,9 4,95 IDIM1 Param 19 I(4) 4 scalar IDIM2 Param 19 I(4) 4 scalar IGSTL Param 14 I(4) 4 scalar 19 IGSTR Param 14 I(4) 4 scalar 19 IM Param 3 I(4) 4 scalar 15,16,19 IMAX Param 15 I(4) 4 scalar 28 INLTE Param 27 I(4) 4 scalar 27 INLTEP Param 27 I(4) 4 scalar INPES Param 11 I(4) 4 scalar 16,19 ITAIL Param 16 I(4) 4 scalar JDIM1 Param 20 I(4) 4 scalar Page 4 Source Listing GRADFS 2014-12-17 20:46 Symbol Table GRADFS.F Name Object Declared Type Bytes Dimen Elements Attributes References JDIM2 Param 20 I(4) 4 scalar JGSTL Param 15 I(4) 4 scalar 20 JGSTR Param 15 I(4) 4 scalar 20 JM Param 3 I(4) 4 scalar 17,20 JNPES Param 11 I(4) 4 scalar 17,20 JTAIL Param 17 I(4) 4 scalar KCCO2 Dummy 12 I(4) 4 scalar ARG,INOUT KO2 Param 30 I(4) 4 scalar 31 KO21 Param 31 I(4) 4 scalar KO2M Param 31 I(4) 4 scalar L Param 14 I(4) 4 scalar 18,19,20,22,23,25,33 LL Param 20 I(4) 4 scalar 20,21 LL3P Param 25 I(4) 4 scalar 28 LL3PI Param 28 I(4) 4 scalar LLM1 Param 21 I(4) 4 scalar LLM2 Param 21 I(4) 4 scalar LLM3 Param 21 I(4) 4 scalar LLP1 Param 20 I(4) 4 scalar 28 LLP1I Param 28 I(4) 4 scalar LLP2 Param 20 I(4) 4 scalar LLP3 Param 20 I(4) 4 scalar LM Param 3 I(4) 4 scalar 14 LM1 Param 19 I(4) 4 scalar LM2 Param 19 I(4) 4 scalar LM3 Param 19 I(4) 4 scalar LNGTH Param 22 I(4) 4 scalar 34,87 LP1 Param 18 I(4) 4 scalar 22,23,24,28,38,80 LP121 Param 24 I(4) 4 scalar LP1I Param 28 I(4) 4 scalar LP1M Param 22 I(4) 4 scalar 22 LP1M1 Param 22 I(4) 4 scalar LP1V Param 23 I(4) 4 scalar LP2 Param 18 I(4) 4 scalar LP3 Param 18 I(4) 4 scalar LSM Param 3 I(4) 4 scalar LV Local 80 I(4) 4 scalar 80,81 NB Param 26 I(4) 4 scalar 29 NB1 Param 29 I(4) 4 scalar NBLM Param 17 I(4) 4 scalar NBLW Param 16 I(4) 4 scalar NBLX Param 16 I(4) 4 scalar NBLY Param 16 I(4) 4 scalar 17,24 NFILE Dummy 12 I(4) 4 scalar ARG,INOUT 84 NNLTE Param 27 I(4) 4 scalar O3INT Subr 83 83 PI Local 52 R(4) 4 scalar 52,63 Q14330 Scalar 24 R(4) 4 scalar COM 59 Q19001 Scalar 23 R(4) 4 scalar COM 53 RAD1 Scalar 34 R(4) 4 1 1665 COM 88,90,91,92 RAD2 Scalar 34 R(4) 4 1 1665 COM 88,89,91,93 RAD3 Scalar 34 R(4) 4 1 1665 COM 88,90,91,94 RAD4 Scalar 34 R(4) 4 1 1665 COM 88,89,91,95 RCO2 Scalar 25 R(4) 4 scalar COM 68 RDFSAV Common 23 140 SAVMEM Common 31 26640 Page 5 Source Listing GRADFS 2014-12-17 20:46 Symbol Table GRADFS.F Name Object Declared Type Bytes Dimen Elements Attributes References SIGL Dummy 12 R(4) 4 1 46 ARG,INOUT 81,83 TWENTY Scalar 24 R(4) 4 scalar COM 61 XAH2SW Local 40 R(4) 4 1 5 43,74 XAO3SW Local 40 R(4) 4 1 5 42,73 XBSW Local 40 R(4) 4 1 5 44,75 XLATP Scalar 23 R(4) 4 scalar COM 99 XLATT Scalar 23 R(4) 4 scalar COM 100 Page 6 Source Listing GRADFS 2014-12-17 20:46 Subprograms/Common Blocks GRADFS.F SUBPROGRAMS/COMMON BLOCKS Name Object Declared Type Bytes Dimen Elements Attributes References GRADFS Subr 2 RDFSAV Common 23 140 SAVMEM Common 31 26640 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 Page 7 Source Listing GRADFS 2014-12-17 20:46 GRADFS.F -fixed no -fpconstant -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 : GRADFS.lst -o filename : none COMPILER: Intel(R) Fortran 12.1-2100