Page 1 Source Listing LL2EG 2014-12-17 20:48 LL2EG.f 1 SUBROUTINE LL2EG (YLAT,WLON,XI,YJ) 2 C$$$ SUBPROGRAM DOCUMENTATION BLOCK 3 C . . . 4 C SUBPROGRAM: LL2EG ETA (LAT,LON) TO (I,J) 5 C PRGRMMR: TREADON ORG: W/NP2 DATE: 92-12-23 6 C 7 C ABSTRACT: 8 C GIVEN (LAT,LON) ON A FILLED E-GRID THIS ROUTINE 9 C RETURNS THE CORRESPONDING FLOATING POINT (I,J) 10 C OF THIS LOCATION. 11 C . 12 C 13 C PROGRAM HISTORY LOG: 14 C ??-??-?? DAVID PLUMMER - SUBROUTINE LL2EG IN ETAPACKC. 15 C 93-02-26 RUSS TREADON - GENERALIZED LL2EG TO HANDLE 16 C VARIABLE SIZE E-GRIDS. 17 C 18 C USAGE: CALL LL2EG(YLAT,WLON,XI,YJ) 19 C INPUT ARGUMENT LIST: 20 C YLAT - LATITUDE OF (I,J) POINT. 21 C WLON - WEST LONGITUDE OF (I,J) POINT. 22 C 23 C OUTPUT ARGUMENT LIST: 24 C XI - FLOATING POINT "I" OF (LAT,LON) ON FILLED E-GRID. 25 C YJ - FLOATING POINT "J" OF (LAT,LON) ON FILLED E-GRID. 26 C 27 C OUTPUT FILES: 28 C NONE 29 C 30 C SUBPROGRAMS CALLED: 31 C UTILITIES: 32 C NONE 33 C LIBRARY: 34 C COMMON - EGRID 35 C 36 C ATTRIBUTES: 37 C LANGUAGE: FORTRAN 38 C MACHINE : CRAY C-90 39 C$$$ 40 C 41 C 42 C DECLARE VARIABLES, COMMON, EQUIVALENCES, AND CONSTANTS. 43 C 44 LOGICAL NORTH 45 REAL ALONVT, DLAM, DPHI, LAM, LAM0, PHI0, POLEI, POLEJ 46 REAL SINPH0, TANPH0, WLONC, XMESHL, YLATC 47 C 48 C INCLUDE EGRID SPECIFICATIONS. 49 INCLUDE "EGRID.comm" 50 C 56 EQUIVALENCE (YLATC,ALONVT), (DPHI, POLEI) 57 EQUIVALENCE (WLONC, POLEJ), (DLAM,XMESHL) 58 C 59 DATA CONV2R/0.017453293/,CONV2D/57.2958/ 60 DATA EARTHR/6371.2/ 61 C 62 C Page 2 Source Listing LL2EG 2014-12-17 20:48 LL2EG.f 63 C******************************************************************** 64 C START LL2EG HERE. 65 C 66 C CONVERT (LAT,LON) TO FLOATING POINT (I,J). 67 PHI = YLAT*CONV2R 68 LAM = WLON*CONV2R 69 COSPHI = COS(PHI) 70 SINPHI = SIN(PHI) 71 72 X1 = COSPHI*COS(LAM-LAM0) 73 X = COSPH0*X1+SINPH0*SINPHI 74 Y = -COSPHI*SIN(LAM-LAM0) 75 Z = -SINPH0*X1+COSPH0*SINPHI 76 BIGPHI = ATAN(Z/(SQRT(X*X+Y*Y)))*CONV2D 77 BIGLAM = ATAN(Y/X)*CONV2D 78 XI = (BIGLAM/DLAM)+(IDIM+1)/2.0 79 YJ = (BIGPHI/DPHI)+(JDIM+1)/2.0 80 C 81 C END OF ROUTINE. 82 C 83 RETURN 84 END ENTRY POINTS Name ll2eg_ Page 3 Source Listing LL2EG 2014-12-17 20:48 Symbol Table LL2EG.f SYMBOL CROSS REFERENCE Name Object Declared Type Bytes Dimen Elements Attributes References ALONVT Scalar 45 R(4) 4 scalar COM ATAN Func 71 scalar 71,72 BIGLAM Local 72 R(4) 4 scalar 72,73 BIGPHI Local 71 R(4) 4 scalar 71,74 CONV2D Local 54 R(4) 4 scalar 54,71,72 CONV2R Local 54 R(4) 4 scalar 54,62,63 COS Func 64 scalar 64,67 COSPH0 Scalar 3 R(4) 4 scalar COM 68,70 COSPHI Local 64 R(4) 4 scalar 64,67,69 DLAM Scalar 45 R(4) 4 scalar COM 73 DPHI Scalar 45 R(4) 4 scalar COM 74 EARTHR Local 55 R(4) 4 scalar 55 EGRID Common 2 48 IDIM Scalar 3 I(4) 4 scalar COM 73 JDIM Scalar 3 I(4) 4 scalar COM 74 LAM Local 45 R(4) 4 scalar 63,67,69 LAM0 Scalar 45 R(4) 4 scalar COM 67,69 LL2EG Subr 1 NORTH Scalar 44 L(4) 4 scalar COM PHI Local 62 R(4) 4 scalar 62,64,65 PHI0 Scalar 45 R(4) 4 scalar COM POLEI Scalar 45 R(4) 4 scalar COM POLEJ Scalar 45 R(4) 4 scalar COM SIN Func 65 scalar 65,69 SINPH0 Scalar 46 R(4) 4 scalar COM 68,70 SINPHI Local 65 R(4) 4 scalar 65,68,70 SQRT Func 71 scalar 71 TANPH0 Scalar 46 R(4) 4 scalar COM WLON Dummy 1 R(4) 4 scalar ARG,INOUT 63 WLONC Scalar 46 R(4) 4 scalar COM X Local 68 R(4) 4 scalar 68,71,72 X1 Local 67 R(4) 4 scalar 67,68,70 XI Dummy 1 R(4) 4 scalar ARG,INOUT 73 XMESHL Scalar 46 R(4) 4 scalar COM Y Local 69 R(4) 4 scalar 69,71,72 YJ Dummy 1 R(4) 4 scalar ARG,INOUT 74 YLAT Dummy 1 R(4) 4 scalar ARG,INOUT 62 YLATC Scalar 46 R(4) 4 scalar COM Z Local 70 R(4) 4 scalar 70,71 Page 4 Source Listing LL2EG 2014-12-17 20:48 Subprograms/Common Blocks LL2EG.f SUBPROGRAMS/COMMON BLOCKS Name Object Declared Type Bytes Dimen Elements Attributes References EGRID Common 2 48 LL2EG Subr 1 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 5 Source Listing LL2EG 2014-12-17 20:48 LL2EG.f -fpe3 -fprm nearest no -ftz -fp_model noprecise -fp_model fast -fp_model nostrict -fp_model nosource -fp_model nodouble -fp_model noextended -fp_model novery_fast -fp_model noexcept -fp_model nono_except -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 : /gpfs/tp1/usrx/local/intel/composer_xe_2011_sp1.11.339/compiler/include/,.f,./.f,/usrx/local/intel/composerxe/mkl/include/.f, /usrx/local/intel/composerxe/tbb/include/.f,/gpfs/tp1/usrx/local/intel/composer_xe_2011_sp1.11.339/compiler/include/intel64/.f, /gpfs/tp1/usrx/local/intel/composer_xe_2011_sp1.11.339/compiler/include/.f,/usr/local/include/.f,/usr/lib/gcc/x86_64-redhat-linux/4.4.7/include/.f, /usr/include/.f,/usr/include/.f -list filename : LL2EG.lst -o filename : none COMPILER: Intel(R) Fortran 12.1-2100