Page 1 Source Listing CALTHTE 2014-12-17 20:47 CALTHTE.f 1 SUBROUTINE CALTHTE(P1D,T1D,Q1D,THTE) 2 C$$$ SUBPROGRAM DOCUMENTATION BLOCK 3 C . . . 4 C SUBPROGRAM: CALTHTE COMPUTES THETA-E 5 C PRGRMMR: TREADON ORG: W/NP2 DATE: 93-06-18 6 C 7 C ABSTRACT: 8 C THIS ROUTINE COMPUTES THE EQUIVALENT POTENTIAL TEMPERATURE 9 C GIVEN PRESSURE, TEMPERATURE, AND SPECIFIC HUMIDITY. THE 10 C EQUATIONS OF BOLTON (MWR,1980) ARE USED. 11 C . 12 C 13 C PROGRAM HISTORY LOG: 14 C 93-06-18 RUSS TREADON 15 C 98-06-16 T BLACK - CONVERSION FROM 1-D TO 2-D 16 C 00-01-04 JIM TUCCILLO - MPI VERSION 17 C 18 C USAGE: CALL CALTHTE(P1D,T1D,Q1D,THTE) 19 C INPUT ARGUMENT LIST: 20 C P1D - PRESSURE (PA) 21 C T1D - TEMPERATURE (K) 22 C Q1D - SPECIFIC HUMIDITY (KG/KG) 23 C 24 C OUTPUT ARGUMENT LIST: 25 C THTE - THETA-E (K) 26 C 27 C OUTPUT FILES: 28 C NONE 29 C 30 C SUBPROGRAMS CALLED: 31 C UTILITIES: 32 C VAPOR - FUNCTION TO CALCULATE VAPOR PRESSURE. 33 C LIBRARY: 34 C NONE 35 C 36 C ATTRIBUTES: 37 C LANGUAGE: FORTRAN 38 C MACHINE : CRAY C-90 39 C$$$ 40 C 41 C 42 C INCLUDE/SET PARAMETERS. 43 INCLUDE "parmeta" 61 INCLUDE "params" 103 INCLUDE "CTLBLK.comm" 104 C 119 PARAMETER (KG2G=1.E3) 120 PARAMETER (D35=3.5,D4805=4.805,H2840=2840.,H55=55.) 121 PARAMETER (D2845=0.2845,D00028=0.00028,D3376=3.376) 122 PARAMETER (D00254=0.00254,D00081=0.00081,D81=0.81) 123 PARAMETER (D28=0.28,H2675=2675.) 124 C 125 C DECLARE VARIABLES. 126 C 127 REAL P1D(IM,JM),T1D(IM,JM),Q1D(IM,JM),THTE(IM,JM) 128 C 129 C*************************************************************** Page 2 Source Listing CALTHTE 2014-12-17 20:47 CALTHTE.f 130 C START CALTHTE. 131 C 132 C ZERO THETA-E ARRAY 133 DO J=JSTA,JEND 134 DO I=1,IM 135 THTE(I,J)=D00 136 ENDDO 137 ENDDO 138 C 139 C COMPUTE THETA-E. 140 C 141 DO J=JSTA_M,JEND_M 142 DO I=2,IM-1 143 P = P1D(I,J) 144 T = T1D(I,J) 145 Q = Q1D(I,J) 146 EVP = P*Q/(EPS+ONEPS*Q) 147 RMX = EPS*EVP/(P-EVP) 148 CKAPA = D2845*(1.-D28*RMX) 149 RKAPA = 1./CKAPA 150 ARG = EVP*D01 151 ARG = AMAX1(H1M12,ARG) 152 DENOM = D35*ALOG(T) - ALOG(EVP*D01) - D4805 153 TLCL = H2840/DENOM + H55 154 PLCL = P*(TLCL/T)**RKAPA 155 FAC = (P1000/P)**CKAPA 156 ETERM = (D3376/TLCL-D00254)*(RMX*KG2G*(H1+D81*RMX)) 157 THETAE = T*FAC*EXP(ETERM) 158 THTE(I,J)= THETAE 159 ENDDO 160 ENDDO 161 C 162 C END OF ROUTINE. 163 C 164 RETURN 165 END Page 3 Source Listing CALTHTE 2014-12-17 20:47 Entry Points CALTHTE.f ENTRY POINTS Name calthte_ SYMBOL CROSS REFERENCE Name Object Declared Type Bytes Dimen Elements Attributes References A2 Param 41 R(4) 4 scalar A3 Param 41 R(4) 4 scalar A4 Param 41 R(4) 4 scalar AD05 Param 33 R(4) 4 scalar ALOG Func 80 scalar 80 AMAX1 Func 79 scalar 79 ARG Local 78 R(4) 4 scalar 78,79 BETA Param 37 R(4) 4 scalar 37 BTG Param 37 R(4) 4 scalar CALTHTE Subr 1 CAPA Param 28 R(4) 4 scalar 35 CFT0 Param 33 R(4) 4 scalar CKAPA Local 76 R(4) 4 scalar 76,77,83 CM1 Param 23 R(4) 4 scalar CM2 Param 23 R(4) 4 scalar CM3 Param 23 R(4) 4 scalar CP Param 34 R(4) 4 scalar 36,40 CPDR Param 36 R(4) 4 scalar CTLBLK Common 2 96 D00 Param 25 R(4) 4 scalar 63 D00028 Param 49 R(4) 4 scalar D00081 Param 50 R(4) 4 scalar D001 Param 25 R(4) 4 scalar D00254 Param 50 R(4) 4 scalar 84 D0065 Param 29 R(4) 4 scalar D01 Param 25 R(4) 4 scalar 78,80 D05 Param 29 R(4) 4 scalar 33 D1 Param 25 R(4) 4 scalar D115 Param 29 R(4) 4 scalar D125 Param 27 R(4) 4 scalar D25 Param 27 R(4) 4 scalar D28 Param 51 R(4) 4 scalar 76 D2845 Param 49 R(4) 4 scalar 76 D3376 Param 49 R(4) 4 scalar 84 D35 Param 48 R(4) 4 scalar 80 D4805 Param 48 R(4) 4 scalar 80 D50 Param 25 R(4) 4 scalar D608 Param 29 R(4) 4 scalar D81 Param 50 R(4) 4 scalar 84 DENOM Local 80 R(4) 4 scalar 80,81 DP10M Param 35 R(4) 4 scalar DT Scalar 4 R(4) 4 scalar COM DTR Param 34 R(4) 4 scalar 34 ELDR Param 36 R(4) 4 scalar Page 4 Source Listing CALTHTE 2014-12-17 20:47 Symbol Table CALTHTE.f Name Object Declared Type Bytes Dimen Elements Attributes References ELWV Param 36 R(4) 4 scalar 36 EPS Param 39 R(4) 4 scalar 39,74,75 EPSQ2 Param 40 R(4) 4 scalar ERAD Param 34 R(4) 4 scalar ETERM Local 84 R(4) 4 scalar 84,85 EVP Local 74 R(4) 4 scalar 74,75,78,80 EXP Func 85 scalar 85 FAC Local 83 R(4) 4 scalar 83,85 FIRST Scalar 3 R(4) 4 scalar COM FMD Param 39 R(4) 4 scalar 39 FMW Param 39 R(4) 4 scalar 39 G Param 32 R(4) 4 scalar 32 GI Param 32 R(4) 4 scalar 36 H1 Param 26 R(4) 4 scalar 33,84 H100 Param 27 R(4) 4 scalar H1000 Param 27 R(4) 4 scalar H10000 Param 28 R(4) 4 scalar H10E5 Param 28 R(4) 4 scalar H1E2 Param 30 R(4) 4 scalar H1E4 Param 30 R(4) 4 scalar H1M12 Param 26 R(4) 4 scalar 79 H1M2 Param 30 R(4) 4 scalar H1M5 Param 27 R(4) 4 scalar H2 Param 26 R(4) 4 scalar H2675 Param 51 R(4) 4 scalar H2840 Param 48 R(4) 4 scalar 81 H4 Param 26 R(4) 4 scalar H55 Param 48 R(4) 4 scalar 81 H999 Param 25 R(4) 4 scalar H99999 Param 26 R(4) 4 scalar I Local 62 I(4) 4 scalar 62,63,70,71,72,73,86 ICNT Scalar 7 I(4) 4 1 1024 COM IDAT Scalar 3 I(4) 4 1 3 COM IDIM1 Param 14 I(4) 4 scalar IDIM2 Param 14 I(4) 4 scalar IDN Scalar 7 I(4) 4 scalar COM IDSP Scalar 7 I(4) 4 1 1024 COM IDTAD Scalar 4 I(4) 4 scalar COM IGSTL Param 9 I(4) 4 scalar 14 IGSTR Param 9 I(4) 4 scalar 14 IHRST Scalar 3 I(4) 4 scalar COM IM Param 3 I(4) 4 scalar 4,6,7,9,11,14,55,62,70 IM1 Param 6 I(4) 4 scalar IMT Param 4 I(4) 4 scalar INPES Param 6 I(4) 4 scalar 11,14 IOUT Scalar 3 I(4) 4 scalar COM ITAIL Param 11 I(4) 4 scalar IUP Scalar 7 I(4) 4 scalar COM J Local 61 I(4) 4 scalar 61,63,69,71,72,73,86 JAM Param 6 I(4) 4 scalar JDIM1 Param 15 I(4) 4 scalar JDIM2 Param 15 I(4) 4 scalar JEND Scalar 6 I(4) 4 scalar COM 61 JEND_2U Scalar 8 I(4) 4 scalar COM JEND_M Scalar 7 I(4) 4 scalar COM 69 Page 5 Source Listing CALTHTE 2014-12-17 20:47 Symbol Table CALTHTE.f Name Object Declared Type Bytes Dimen Elements Attributes References JEND_M2 Scalar 7 I(4) 4 scalar COM JGSTL Param 10 I(4) 4 scalar 15 JGSTR Param 10 I(4) 4 scalar 15 JM Param 3 I(4) 4 scalar 4,6,9,12,15,55 JM2 Param 6 I(4) 4 scalar JMT Param 4 I(4) 4 scalar JNPES Param 6 I(4) 4 scalar 12,15 JSTA Scalar 6 I(4) 4 scalar COM 61 JSTA_2L Scalar 8 I(4) 4 scalar COM JSTA_M Scalar 7 I(4) 4 scalar COM 69 JSTA_M2 Scalar 7 I(4) 4 scalar COM JTAIL Param 12 I(4) 4 scalar KG2G Param 47 I(4) 4 scalar 84 KSAV Param 14 I(4) 4 scalar KSFI Param 13 I(4) 4 scalar KSLI Param 14 I(4) 4 scalar KSLP Param 14 I(4) 4 scalar KSLPD Param 13 I(4) 4 scalar KSMUD Param 13 I(4) 4 scalar KSRH Param 14 I(4) 4 scalar KSSLP Param 13 I(4) 4 scalar KST Param 13 I(4) 4 scalar KSUV Param 14 I(4) 4 scalar LB Param 9 I(4) 4 scalar LHEAT Param 40 I(4) 4 scalar 40 LHTOCP Param 40 I(4) 4 scalar LIST Scalar 3 I(4) 4 scalar COM LM Param 3 I(4) 4 scalar 8,19 LM1 Param 8 I(4) 4 scalar LP1 Param 8 I(4) 4 scalar LSM Param 3 I(4) 4 scalar ME Scalar 6 I(4) 4 scalar COM NBC Scalar 3 I(4) 4 scalar COM NBOCO Scalar 4 I(4) 4 scalar COM NCP Scalar 4 I(4) 4 scalar COM NDDAMP Scalar 4 I(4) 4 scalar COM NFCST Scalar 3 I(4) 4 scalar COM NINC Param 7 I(4) 4 scalar NPHS Scalar 5 I(4) 4 scalar COM NPREC Scalar 4 I(4) 4 scalar COM NRADL Scalar 5 I(4) 4 scalar COM NRADS Scalar 5 I(4) 4 scalar COM NSHDE Scalar 4 I(4) 4 scalar COM NSTAT Param 19 I(4) 4 scalar NSUMD Param 15 I(4) 4 scalar NTSD Scalar 4 I(4) 4 scalar COM NTSTM Scalar 4 I(4) 4 scalar COM NUM_PROCS Scalar 6 I(4) 4 scalar COM NWORD Param 19 I(4) 4 scalar 19 NWORD3 Param 19 I(4) 4 scalar NWORD6 Param 19 I(4) 4 scalar ONEPS Param 39 R(4) 4 scalar 74 OVERRC Param 33 R(4) 4 scalar 33 P Local 71 R(4) 4 scalar 71,74,75,82,83 P1000 Param 35 R(4) 4 scalar 83 Page 6 Source Listing CALTHTE 2014-12-17 20:47 Symbol Table CALTHTE.f Name Object Declared Type Bytes Dimen Elements Attributes References P1D Dummy 1 R(4) 4 2 91719 ARG,INOUT 71 PARALLEL Common 9 8240 PLCL Local 82 R(4) 4 scalar 82 PQ0 Param 41 R(4) 4 scalar Q Local 73 R(4) 4 scalar 73,74 Q1D Dummy 1 R(4) 4 2 91719 ARG,INOUT 73 RCAPA Param 35 R(4) 4 scalar RD Param 32 R(4) 4 scalar 32,36,40 RDOCP Param 40 R(4) 4 scalar RESTRT Scalar 3 R(4) 4 scalar COM RG Param 36 R(4) 4 scalar RKAPA Local 77 R(4) 4 scalar 77,82 RMX Local 75 R(4) 4 scalar 75,76,84 ROG Param 32 R(4) 4 scalar RTD Param 34 R(4) 4 scalar RUN Scalar 3 R(4) 4 scalar COM SIGMA Scalar 3 R(4) 4 scalar COM SMTHA Param 15 R(4) 4 scalar SMTHB Param 15 R(4) 4 scalar T Local 72 R(4) 4 scalar 72,80,82,85 T1D Dummy 1 R(4) 4 2 91719 ARG,INOUT 72 TFRZ Param 39 R(4) 4 scalar THETAE Local 85 R(4) 4 scalar 85,86 THTE Dummy 1 R(4) 4 2 91719 ARG,INOUT 63,86 TLCL Local 81 R(4) 4 scalar 81,82,84 Page 7 Source Listing CALTHTE 2014-12-17 20:47 Subprograms/Common Blocks CALTHTE.f SUBPROGRAMS/COMMON BLOCKS Name Object Declared Type Bytes Dimen Elements Attributes References CALTHTE Subr 1 CTLBLK Common 2 96 PARALLEL Common 9 8240 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 8 Source Listing CALTHTE 2014-12-17 20:47 CALTHTE.f -fixed no -fpconstant -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 : CALTHTE.lst -o filename : none COMPILER: Intel(R) Fortran 12.1-2100