Page 1 Source Listing CALPW 2014-12-17 20:47 CALPW.f 1 SUBROUTINE CALPW(PW) 2 C$$$ SUBPROGRAM DOCUMENTATION BLOCK 3 C . . . 4 C SUBPROGRAM: CALPW COMPUTES 5 C PRGRMMR: TREADON ORG: W/NP2 DATE: 92-12-24 6 C 7 C ABSTRACT: 8 C THIS ROUTINE COMPUTES PRECIPITABLE WATER IN A COLUMN 9 C EXTENDING FROM THE FIRST ATMOSPHERIC ETA LAYER TO THE 10 C MODEL TOP. THE DEFINITION USED IS 11 C TOP 12 C PRECIPITABLE WATER = SUM (Q+CLDW) DP*HTM/G 13 C BOT 14 C WHERE, 15 C BOT IS THE FIRST ETA LAYER, 16 C TOP IS THE MODEL TOP, 17 C Q IS THE SPECIFIC HUMIDITY (KG/KG) IN THE LAYER 18 C CLDW IS THE CLOUD WATER (KG/KG) IN THE LAYER 19 C DP (Pa) IS THE LAYER THICKNESS. 20 C HTM IS THE HEIGHT MASK AT THAT LAYER (=0 IF BELOW GROUND) 21 C G IS THE GRAVITATIONAL CONSTANT 22 C 23 C PROGRAM HISTORY LOG: 24 C 92-12-24 RUSS TREADON 25 C 96-03-04 MIKE BALDWIN - ADD CLOUD WATER AND SPEED UP CODE 26 C 98-06-15 T BLACK - CONVERSION FROM 1-D TO 2-D 27 C 00-01-04 JIM TUCCILLO - MPI VERSION 28 C 29 C USAGE: CALL CALPW(PW) 30 C INPUT ARGUMENT LIST: 31 C PW - ARRAY OF PRECIPITABLE WATER. 32 C 33 C OUTPUT ARGUMENT LIST: 34 C NONE 35 C 36 C OUTPUT FILES: 37 C NONE 38 C 39 C SUBPROGRAMS CALLED: 40 C UTILITIES: 41 C NONE 42 C LIBRARY: 43 C COMMON - LOOPS 44 C EXTRA 45 C VRBLS 46 C MASKS 47 C CLDWTR 48 C 49 C ATTRIBUTES: 50 C LANGUAGE: FORTRAN 51 C MACHINE : CRAY C-90 52 C$$$ 53 C 54 C 55 C INCLUDE/SET PARAMETERS. 56 INCLUDE "parmeta" 74 INCLUDE "params" Page 2 Source Listing CALPW 2014-12-17 20:47 CALPW.f 75 C 76 C SET DENSITY OF WATER AT 1 ATMOSPHERE PRESSURE, 0C. 77 C UNITS ARE KG/M**3. 119 PARAMETER (RHOWAT=1.E3) 120 C 121 C DECLARE VARIABLES. 122 C 123 INTEGER LLMH 124 REAL ALPM,DZ,PM,PWSUM,RHOAIR 125 REAL PW(IM,JM) 126 C 127 C INCLUDE COMMON BLOCKS. 128 INCLUDE "VRBLS.comm" 135 INCLUDE "EXTRA.comm" 147 INCLUDE "LOOPS.comm" 155 INCLUDE "MASKS.comm" 163 INCLUDE "CLDWTR.comm" 169 INCLUDE "CTLBLK.comm" 170 C 171 C*************************************************************** 172 C START CALPW HERE. 173 C 174 C INITIALIZE PW TO 0. 175 C 190 PW = 0. 191 C 192 C OUTER LOOP OVER VERTICAL DIMENSION. 193 C INNER LOOP OVER HORIZONTAL GRID. 194 C 195 !$omp parallel do 196 !$omp& private(dp) 197 DO L = 1,LM 198 DO J=JSTA,JEND 199 DO I=1,IM 200 DP =PINT(I,J,L+1)-PINT(I,J,L) 201 PW(I,J)=PW(I,J)+(Q(I,J,L)+CWM(I,J,L))*DP*GI*HTM(I,J,L) 202 ENDDO 203 ENDDO 204 ENDDO 205 C 206 C AT ONE TIME THE SUM WAS DIVIDED BY THE DENSITY OF 207 C WATER. (TO GET PW IN M) THIS IS NO LONGER DONE. 208 C PW(I,J) = PWSUM/RHOWAT 209 C 210 C END OF ROUTINE. 211 C 212 RETURN 213 END Page 3 Source Listing CALPW 2014-12-17 20:47 Entry Points CALPW.f ENTRY POINTS Name calpw_ 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 ALPINT Scalar 4 R(4) 4 3 4219074 COM ALPM Local 66 R(4) 4 scalar BETA Param 37 R(4) 4 scalar 37 BTG Param 37 R(4) 4 scalar CALPW Subr 1 CAPA Param 28 R(4) 4 scalar 35 CFT0 Param 33 R(4) 4 scalar CLDWTR Common 2 17610408 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 CWM Scalar 3 R(4) 4 3 4127355 COM 93 D00 Param 25 R(4) 4 scalar D001 Param 25 R(4) 4 scalar D0065 Param 29 R(4) 4 scalar D01 Param 25 R(4) 4 scalar 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 D50 Param 25 R(4) 4 scalar D608 Param 29 R(4) 4 scalar DP Local 92 R(4) 4 scalar 92,93 DP10M Param 35 R(4) 4 scalar DT Scalar 4 R(4) 4 scalar COM DTR Param 34 R(4) 4 scalar 34 DZ Local 66 R(4) 4 scalar ELDR Param 36 R(4) 4 scalar ELWV Param 36 R(4) 4 scalar 36 EPS Param 39 R(4) 4 scalar 39 EPSQ2 Param 40 R(4) 4 scalar ERAD Param 34 R(4) 4 scalar EXTRA Common 3 71506268 FIRST Scalar 3 R(4) 4 scalar COM FIS Scalar 3 R(4) 4 2 91719 COM FMD Param 39 R(4) 4 scalar 39 Page 4 Source Listing CALPW 2014-12-17 20:47 Symbol Table CALPW.f Name Object Declared Type Bytes Dimen Elements Attributes References 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,93 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 H1M2 Param 30 R(4) 4 scalar H1M5 Param 27 R(4) 4 scalar H2 Param 26 R(4) 4 scalar H4 Param 26 R(4) 4 scalar H999 Param 25 R(4) 4 scalar H99999 Param 26 R(4) 4 scalar HBM2 Scalar 3 R(4) 4 2 91719 COM HTM Scalar 5 R(4) 4 3 4127355 COM 93 I Local 91 I(4) 4 scalar 91,92,93 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 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,67,91 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 IOUT Scalar 3 I(4) 4 scalar COM 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 90 I(4) 4 scalar 90,92,93 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 90 JEND_2U Scalar 8 I(4) 4 scalar COM JEND_M Scalar 7 I(4) 4 scalar COM 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,67 JM2 Param 6 I(4) 4 scalar Page 5 Source Listing CALPW 2014-12-17 20:47 Symbol Table CALPW.f Name Object Declared Type Bytes Dimen Elements Attributes References 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 90 JSTA_2L Scalar 8 I(4) 4 scalar COM JSTA_M Scalar 7 I(4) 4 scalar COM 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 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 89 I(4) 4 scalar 89,92,93 LB Param 9 I(4) 4 scalar LC Scalar 3 I(4) 4 2 91719 COM 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 65 I(4) 4 scalar LM Param 3 I(4) 4 scalar 3,5,8,19,89 LM1 Param 8 I(4) 4 scalar LMH Scalar 5 I(4) 4 2 91719 COM 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 MASKS Common 2 34853220 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 Page 6 Source Listing CALPW 2014-12-17 20:47 Symbol Table CALPW.f Name Object Declared Type Bytes Dimen Elements Attributes References OLDRD Scalar 10 R(4) 4 scalar COM ONEPS Param 39 R(4) 4 scalar OVERRC Param 33 R(4) 4 scalar 33 P1000 Param 35 R(4) 4 scalar 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 92 PM Local 66 R(4) 4 scalar 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 PW Dummy 1 R(4) 4 2 91719 ARG,INOUT 82,93 PWSUM Local 66 R(4) 4 scalar Q Scalar 5 R(4) 4 3 4127355 COM 93 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 RHOAIR Local 66 R(4) 4 scalar RHOWAT Param 61 R(4) 4 scalar ROG Param 32 R(4) 4 scalar RTD Param 34 R(4) 4 scalar RUN Scalar 3 R(4) 4 scalar COM SICE Scalar 3 R(4) 4 2 91719 COM SIGMA Scalar 3 R(4) 4 scalar COM SLP Scalar 5 R(4) 4 2 91719 COM SM Scalar 3 R(4) 4 2 91719 COM SMTHA Param 15 R(4) 4 scalar SMTHB Param 15 R(4) 4 scalar SR Scalar 3 R(4) 4 2 91719 COM STDRD Scalar 10 R(4) 4 scalar COM T Scalar 5 R(4) 4 3 4127355 COM T500 Scalar 6 R(4) 4 2 91719 COM TFRZ Param 39 R(4) 4 scalar 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 U00 Scalar 3 R(4) 4 2 91719 COM UL Scalar 3 R(4) 4 1 90 COM V Scalar 5 R(4) 4 3 4127355 COM VBM2 Scalar 3 R(4) 4 2 91719 COM VBM3 Scalar 3 R(4) 4 2 91719 COM VRBLS Common 2 67138308 VTM Scalar 5 R(4) 4 3 4127355 COM Z1000 Scalar 9 R(4) 4 2 91719 COM ZINT Scalar 4 R(4) 4 3 4219074 COM Page 7 Source Listing CALPW 2014-12-17 20:47 Subprograms/Common Blocks CALPW.f SUBPROGRAMS/COMMON BLOCKS Name Object Declared Type Bytes Dimen Elements Attributes References CALPW Subr 1 CLDWTR Common 2 17610408 CTLBLK Common 2 96 EXTRA Common 3 71506268 LOOPS Common 2 748952 MASKS Common 2 34853220 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 Page 8 Source Listing CALPW 2014-12-17 20:47 CALPW.f 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 -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 : CALPW.lst -o filename : none COMPILER: Intel(R) Fortran 12.1-2100