Page 1 Source Listing CALLCL 2014-12-17 20:47 CALLCL.f 1 SUBROUTINE CALLCL(P1D,T1D,Q1D,PLCL,ZLCL) 2 C$$$ SUBPROGRAM DOCUMENTATION BLOCK 3 C . . . 4 C SUBPROGRAM: CALLCL COMPUTES LCL HEIGHTS AND PRESSURE 5 C PRGRMMR: TREADON ORG: W/NP2 DATE: 93-03-15 6 C 7 C ABSTRACT: 8 C THIS ROUTINE COMPUTES THE LIFTING CONDENSATION LEVEL 9 C PRESSURE AND HEIGHT IN EACH COLUMN AT MASS POINTS. 10 C THE HEIGHT IS ABOVE GROUND LEVEL. THE EQUATION USED 11 C TO FIND THE LCL PRESSURE IS FROM BOLTAN (1980,MWR) 12 C AND IS THE SAME AS THAT USED IN SUBROUTINE CALCAPE. 13 C 14 C THIS ROUTINE IS A TEST VERSION. STILL TO BE RESOLVED 15 C IS THE "BEST" PARCEL TO LIFT. 16 C . 17 C 18 C PROGRAM HISTORY LOG: 19 C 93-03-15 RUSS TREADON 20 C 98-06-16 T BLACK - CONVERSION FROM 1-D TO 2-D 21 C 00-01-04 JIM TUCCILLO - MPI VERSION 22 C 23 C USAGE: CALL CALLCL(P1D,T1D,Q1D,PLCL,ZLCL) 24 C INPUT ARGUMENT LIST: 25 C P1D - ARRAY OF PARCEL PRESSURES (PA) 26 C T1D - ARRAY OF PARCEL TEMPERATURES (K) 27 C Q1D - ARRAY OF PARCEL SPECIFIC HUMIDITIES (KG/KG) 28 C 29 C OUTPUT ARGUMENT LIST: 30 C PLCL - PARCEL PRESSURE AT LCL (PA) 31 C ZLCL - PARCEL AGL HEIGHT AT LCL (M) 32 C 33 C OUTPUT FILES: 34 C NONE 35 C 36 C SUBPROGRAMS CALLED: 37 C UTILITIES: 38 C NONE 39 C LIBRARY: 40 C COMMON - LOOPS 41 C VRBLS 42 C EXTRA 43 C OPTIONS 44 C 45 C ATTRIBUTES: 46 C LANGUAGE: FORTRAN 90 47 C MACHINE : CRAY C-90 48 C$$$ 49 C 50 C 51 C 52 C INCLUDE/SET PARAMETERS. 53 INCLUDE "parmeta" 71 INCLUDE "params" 72 C 114 PARAMETER (D35=3.5,D4805=4.805,H2840=2840.,H55=55.) 115 PARAMETER (D2845=0.2845,D28=0.28) Page 2 Source Listing CALLCL 2014-12-17 20:47 CALLCL.f 116 C 117 C DECLARE VARIABLES. 118 C 119 LOGICAL OLDRD,STDRD 120 REAL P1D(IM,JM),T1D(IM,JM),Q1D(IM,JM) 121 REAL PLCL(IM,JM),TLCL(IM,JM),ZLCL(IM,JM) 122 CX REAL PFAL(IM,JM),TFAL(IM,JM),QFAL(IM,JM) 123 CX REAL PSAL(IM,JM),TSAL(IM,JM),QSAL(IM,JM) 124 CX REAL PBAR(IM,JM),TBAR(IM,JM),QBAR(IM,JM) 125 C 126 C INCLUDE COMMON BLOCKS. 127 INCLUDE "EXTRA.comm" 139 INCLUDE "LOOPS.comm" 147 INCLUDE "VRBLS.comm" 154 INCLUDE "OPTIONS.comm" 158 INCLUDE "CTLBLK.comm" 159 C 160 C 161 C********************************************************************** 162 C START CALLCL HERE. 163 C 164 C LOAD OUTPUT ARRAYS WITH SPECIAL VALUE. 165 C 180 DO J=JSTA,JEND 181 DO I=1,IM 182 PLCL(I,J)=SPVAL 183 TLCL(I,J)=SPVAL 184 ZLCL(I,J)=SPVAL 185 ENDDO 186 ENDDO 187 188 CX DO J=JSTA,JEND 189 CX DO I=1,IM 190 CX LLMH = LMH(I,J) 191 CX ALPFAL = D50*(ALPINT(I,J,LLMH)+ALPINT(I,J,LLMH+1)) 192 CX ALPSAL = D50*(ALPINT(I,J,LLMH-1)+ALPINT(I,J,LLMH)) 193 CX PFAL(I,J) = EXP(ALPFAL) 194 CX PSAL(I,J) = EXP(ALPSAL) 195 CX TFAL(I,J) = T(I,J,LLMH) 196 CX TSAL(I,J) = T(I,J,LLMH-1) 197 CX QFAL(I,J) = Q(I,J,LLMH) 198 CX QSAL(I,J) = Q(I,J,LLMH-1) 199 CX PBAR(I,J) = D50*(PFAL(I,J)+PSAL(I,J)) 200 CX TBAR(I,J) = D50*(TFAL(I,J)+TSAL(I,J)) 201 CX QBAR(I,J) = D50*(QFAL(I,J)+QSAL(I,J)) 202 CX ENDDO 203 CX ENDDO 204 205 C 206 C COMPUTE PRESSURE, TEMPERATURE AND AGL HEIGHT AT LCL. 207 C 208 DO 30 J=JSTA_M,JEND_M 209 DO 30 I=2,IM-1 210 EVP =P1D(I,J)*Q1D(I,J)/(EPS+ONEPS*Q1D(I,J)) 211 RMX =EPS*EVP/(P1D(I,J)-EVP) 212 CKAPA =D2845*(1.-D28*RMX) 213 RKAPA =1./CKAPA Page 3 Source Listing CALLCL 2014-12-17 20:47 CALLCL.f 214 ARG =EVP*D01 215 ARG =AMAX1(H1M12,ARG) 216 DENOM =D35*ALOG(T1D(I,J))-ALOG(ARG)-D4805 217 TLCL(I,J)=H2840/DENOM+H55 218 PLCL(I,J)=P1D(I,J)*(TLCL(I,J)/T1D(I,J))**RKAPA 219 ALPLCL =ALOG(PLCL(I,J)) 220 LLMH =LMH(I,J) 221 C 222 DO 20 L=LLMH,1,-1 223 IF(ALPINT(I,J,L).LT.ALPLCL)THEN 224 DLPLCL =ALPLCL-ALPINT(I,J,L+1) 225 DALP =ALPINT(I,J,L)-ALPINT(I,J,L+1) 226 DZ =ZINT(I,J,L)-ZINT(I,J,L+1) 227 ZLCL(I,J)=ZINT(I,J,L+1)+DZ*DLPLCL/DALP 228 ZSFC =FIS(I,J)*GI 229 ZLCL(I,J)=ZLCL(I,J)-ZSFC 230 ZLCL(I,J)=AMAX1(D00,ZLCL(I,J)) 231 CX PSFC=PD(I,J)+50.E2 232 CX WRITE(81,1234)I,J,LLMH,L,PSFC*D01,T1D(I,J),ZSFC 233 CX WRITE(81,1234)I,J,LLMH,L,PLCL(I,J)*D01,TLCL(I,J),ZLCL(I,J) 234 CX WRITE(81,*)' ' 235 CX 1234 FORMAT(I3,1X,I3,1X,2(I2,1X),4(G12.6,1X)) 236 237 GOTO 30 238 ENDIF 239 20 CONTINUE 240 30 CONTINUE 241 C 242 C END OF ROUTINE. 243 C 244 RETURN 245 END Page 4 Source Listing CALLCL 2014-12-17 20:47 Entry Points CALLCL.f ENTRY POINTS Name callcl_ SYMBOL CROSS REFERENCE Name Object Declared Type Bytes Dimen Elements Attributes References 20 Label 140 123 30 Label 141 109,110,138 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 117 scalar 117,120 ALPINT Scalar 4 R(4) 4 3 4219074 COM 124,125,126 ALPLCL Local 120 R(4) 4 scalar 120,124,125 AMAX1 Func 116 scalar 116,131 ARG Local 115 R(4) 4 scalar 115,116,117 BETA Param 37 R(4) 4 scalar 37 BTG Param 37 R(4) 4 scalar CALLCL Subr 1 CAPA Param 28 R(4) 4 scalar 35 CFT0 Param 33 R(4) 4 scalar CKAPA Local 113 R(4) 4 scalar 113,114 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 131 D001 Param 25 R(4) 4 scalar D0065 Param 29 R(4) 4 scalar D01 Param 25 R(4) 4 scalar 115 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 57 R(4) 4 scalar 113 D2845 Param 57 R(4) 4 scalar 113 D35 Param 56 R(4) 4 scalar 117 D4805 Param 56 R(4) 4 scalar 117 D50 Param 25 R(4) 4 scalar D608 Param 29 R(4) 4 scalar DALP Local 126 R(4) 4 scalar 126,128 DENOM Local 117 R(4) 4 scalar 117,118 DLPLCL Local 125 R(4) 4 scalar 125,128 DP10M Param 35 R(4) 4 scalar DT Scalar 4 R(4) 4 scalar COM DTR Param 34 R(4) 4 scalar 34 Page 5 Source Listing CALLCL 2014-12-17 20:47 Symbol Table CALLCL.f Name Object Declared Type Bytes Dimen Elements Attributes References DZ Local 127 R(4) 4 scalar 127,128 ELDR Param 36 R(4) 4 scalar ELWV Param 36 R(4) 4 scalar 36 EPS Param 39 R(4) 4 scalar 39,111,112 EPSQ2 Param 40 R(4) 4 scalar ERAD Param 34 R(4) 4 scalar EVP Local 111 R(4) 4 scalar 111,112,115 EXTRA Common 3 71506268 FIRST Scalar 3 R(4) 4 scalar COM FIS Scalar 3 R(4) 4 2 91719 COM 129 FMD Param 39 R(4) 4 scalar 39 FMW Param 39 R(4) 4 scalar 39 FSL Scalar 8 R(4) 4 2 91719 COM G Param 32 R(4) 4 scalar 32 GI Param 32 R(4) 4 scalar 36,129 H1 Param 26 R(4) 4 scalar 33 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 116 H1M2 Param 30 R(4) 4 scalar H1M5 Param 27 R(4) 4 scalar H2 Param 26 R(4) 4 scalar H2840 Param 56 R(4) 4 scalar 118 H4 Param 26 R(4) 4 scalar H55 Param 56 R(4) 4 scalar 118 H999 Param 25 R(4) 4 scalar H99999 Param 26 R(4) 4 scalar I Local 82 I(4) 4 scalar 82,83,84,85,110,111,112,117,118,11 9,120,121,124,125,126,127,128,129, 130,131 IBESSL Scalar 2 I(4) 4 scalar COM 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 IFLAG Scalar 2 I(4) 4 scalar COM IGSTL Param 9 I(4) 4 scalar 14 IGSTR Param 9 I(4) 4 scalar 14 IHHA Scalar 3 I(4) 4 1 760 COM IHLA Scalar 3 I(4) 4 1 760 COM IHRST Scalar 3 I(4) 4 scalar COM IM Param 3 I(4) 4 scalar 2,3,4,5,6,7,8,9,11,14,62,63,82,110 IM1 Param 6 I(4) 4 scalar IMJM Param 2 I(4) 4 scalar 5 IMT Param 4 I(4) 4 scalar INPES Param 6 I(4) 4 scalar 11,14 IOFFS Scalar 2 I(4) 4 scalar COM IOUT Scalar 3 I(4) 4 scalar COM Page 6 Source Listing CALLCL 2014-12-17 20:47 Symbol Table CALLCL.f Name Object Declared Type Bytes Dimen Elements Attributes References ITAIL Param 11 I(4) 4 scalar IUP Scalar 7 I(4) 4 scalar COM IVHA Scalar 3 I(4) 4 1 760 COM IVLA Scalar 3 I(4) 4 1 760 COM J Local 81 I(4) 4 scalar 81,83,84,85,109,111,112,117,118,11 9,120,121,124,125,126,127,128,129, 130,131 JAM Param 6 I(4) 4 scalar 3 JDIM1 Param 15 I(4) 4 scalar JDIM2 Param 15 I(4) 4 scalar JEND Scalar 6 I(4) 4 scalar COM 81 JEND_2U Scalar 8 I(4) 4 scalar COM JEND_M Scalar 7 I(4) 4 scalar COM 109 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 2,3,4,5,6,7,8,9,12,15,62,63 JM2 Param 6 I(4) 4 scalar JMT Param 4 I(4) 4 scalar JNPES Param 6 I(4) 4 scalar 12,15 JRA Scalar 3 I(4) 4 1 760 COM JSTA Scalar 6 I(4) 4 scalar COM 81 JSTA_2L Scalar 8 I(4) 4 scalar COM JSTA_M Scalar 7 I(4) 4 scalar COM 109 JSTA_M2 Scalar 7 I(4) 4 scalar COM JTAIL Param 12 I(4) 4 scalar KMNT Scalar 5 I(4) 4 2 4118670 COM KMNTM Scalar 5 I(4) 4 1 45 COM KSAV Param 14 I(4) 4 scalar KSB Scalar 2 I(4) 4 scalar COM 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 L Local 123 I(4) 4 scalar 123,124,125,126,127,128 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 LLMH Local 121 I(4) 4 scalar 121,123 LM Param 3 I(4) 4 scalar 5,8,19 LM1 Param 8 I(4) 4 scalar LMH Scalar 5 I(4) 4 2 91719 COM 121 LMV Scalar 5 I(4) 4 2 91719 COM LOOPS Common 2 748952 LP1 Param 8 I(4) 4 scalar 4 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 Page 7 Source Listing CALLCL 2014-12-17 20:47 Symbol Table CALLCL.f Name Object Declared Type Bytes Dimen Elements Attributes References 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 OLDRD Scalar 61 L(4) 4 scalar COM ONEPS Param 39 R(4) 4 scalar 111 OPTIONS Common 2 24 OVERRC Param 33 R(4) 4 scalar 33 P1000 Param 35 R(4) 4 scalar P1D Dummy 1 R(4) 4 2 91719 ARG,INOUT 111,112,119 PARALLEL Common 9 8240 PD Scalar 3 R(4) 4 2 91719 COM PDSL Scalar 6 R(4) 4 2 91719 COM PDVP1 Scalar 6 R(4) 4 2 91719 COM PFSL Scalar 7 R(4) 4 2 91719 COM PINT Scalar 4 R(4) 4 3 4219074 COM PLCL Dummy 1 R(4) 4 2 91719 ARG,INOUT 83,119,120 PQ0 Param 41 R(4) 4 scalar PSLP Scalar 6 R(4) 4 2 91719 COM PTSL Scalar 7 R(4) 4 2 91719 COM Q Scalar 5 R(4) 4 3 4127355 COM Q1D Dummy 1 R(4) 4 2 91719 ARG,INOUT 111 QSL Scalar 9 R(4) 4 2 91719 COM RCAPA Param 35 R(4) 4 scalar RD Param 32 R(4) 4 scalar 32,36,40 RDOCP Param 40 R(4) 4 scalar RES Scalar 3 R(4) 4 2 91719 COM RESTRT Scalar 3 R(4) 4 scalar COM RG Param 36 R(4) 4 scalar RKAPA Local 114 R(4) 4 scalar 114,119 RMX Local 112 R(4) 4 scalar 112,113 ROG Param 32 R(4) 4 scalar RTD Param 34 R(4) 4 scalar RUN Scalar 3 R(4) 4 scalar COM SATDEL Scalar 2 R(4) 4 scalar COM SIGMA Scalar 3 R(4) 4 scalar COM SLP Scalar 5 R(4) 4 2 91719 COM SMTHA Param 15 R(4) 4 scalar SMTHB Param 15 R(4) 4 scalar SPVAL Scalar 2 R(4) 4 scalar COM 83,84,85 STDRD Scalar 61 L(4) 4 scalar COM T Scalar 5 R(4) 4 3 4127355 COM Page 8 Source Listing CALLCL 2014-12-17 20:47 Symbol Table CALLCL.f Name Object Declared Type Bytes Dimen Elements Attributes References T1D Dummy 1 R(4) 4 2 91719 ARG,INOUT 117,119 T500 Scalar 6 R(4) 4 2 91719 COM TFRZ Param 39 R(4) 4 scalar TLCL Local 63 R(4) 4 2 91719 84,118,119 TSL Scalar 8 R(4) 4 2 91719 COM TTV Scalar 8 R(4) 4 2 91719 COM U Scalar 5 R(4) 4 3 4127355 COM V Scalar 5 R(4) 4 3 4127355 COM VRBLS Common 2 67138308 Z1000 Scalar 9 R(4) 4 2 91719 COM ZINT Scalar 4 R(4) 4 3 4219074 COM 127,128 ZLCL Dummy 1 R(4) 4 2 91719 ARG,INOUT 85,128,130,131 ZSFC Local 129 R(4) 4 scalar 129,130 Page 9 Source Listing CALLCL 2014-12-17 20:47 Subprograms/Common Blocks CALLCL.f SUBPROGRAMS/COMMON BLOCKS Name Object Declared Type Bytes Dimen Elements Attributes References CALLCL Subr 1 CTLBLK Common 2 96 EXTRA Common 3 71506268 LOOPS Common 2 748952 OPTIONS Common 2 24 PARALLEL Common 9 8240 VRBLS Common 2 67138308 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 Page 10 Source Listing CALLCL 2014-12-17 20:47 CALLCL.f no -f66 no -f77rtl no -fast -fpscomp nofilesfromcmd -fpscomp nogeneral -fpscomp noioformat -fpscomp noldio_spacing -fpscomp nologicals -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 : CALLCL.lst -o filename : none COMPILER: Intel(R) Fortran 12.1-2100