Page 1 Source Listing CALDWP2 2014-12-17 20:47 CALDWP2.f 1 SUBROUTINE CALDWP2(P1D,Q1D,TDWP,T1D) 2 C$$$ SUBPROGRAM DOCUMENTATION BLOCK 3 C . . . 4 C SUBPROGRAM: CALDWP2 COMPUTES 5 C PRGRMMR: TREADON ORG: W/NP2 DATE: 92-12-22 6 C 7 C ABSTRACT: COMPUTES ANEMOMETER LEVEL WINDS. AT EACH WIND 8 C POINT IN THE GRID FIT A LOGRITHMIC, THEN LINEAR PROFILE 9 C BETWEEN THE WIND COMPONENTS AT Z0 AND THOSE IN THE FIRST 10 C ATMOSPHERIC ETA LAYER. WE MATCH THE TWO PROFILES AT A 11 C CRITICAL HEIGHT, ZC, MAINTAINING ZEROTH AND FIRST ORDER 12 C CONTINUITY. IF THE DISTANCE BETWEEN THE SURFACE AND 13 C THE FIRST ATMOSPHERIC ETA LAYER EXCEEDS A SWITCH HEIGHT, 14 C ZSWITCH, WE ASSUME AN INVERSE EXPONENTIAL PROFILE BETWEEN 15 C THE WIND COMPONENTS AT Z0 AND THOSE IN THE FIRST ATMOSPHERIC 16 C ETA LAYER. THIS INVERSE EXPONENTIAL PROFILE WAS A TEMPORARY 17 C FIX TO AVOID UNREASONABLY SMALL ANEMOMETER LEVEL WINDS. 18 C . 19 C 20 C PROGRAM HISTORY LOG: 21 C 92-12-22 RUSS TREADON 22 C 93-10-04 RUSS TREADON - ADDED CHECK TO BOUND DEWPOINT 23 C TEMPERATURE TO NOT EXCEED THE 24 C AMBIENT TEMPERATURE. 25 C 98-06-08 T BLACK - CONVERSION FROM 1-D TO 2-D 26 C 00-01-04 JIM TUCCILLO - MPI VERSION 27 C 28 C USAGE: CALL CALDWP2(P1D,Q1D,TDWP,T1D) 29 C INPUT ARGUMENT LIST: 30 C P1D - PRESSURE (PA) 31 C Q1D - SPECIFIC HUMIDITY (KG/KG) 32 C T1D - TEMPERATURE (K) 33 C 34 C OUTPUT ARGUMENT LIST: 35 C TDWP - DEWPOINT TEMPERATURE (K) 36 37 C 38 C OUTPUT FILES: 39 C NONE 40 C 41 C SUBPROGRAMS CALLED: 42 C UTILITIES: 43 C DEWPOINT - COMPUTES DEWPOINT GIVEN VAPOR PRESSURE. 44 C LIBRARY: 45 C NONE 46 C 47 C ATTRIBUTES: 48 C LANGUAGE: FORTRAN 90 49 C MACHINE : CRAY C-90 50 C$$$ 51 C 52 C 53 C SET PARAMETERS. 54 INCLUDE "parmeta" 72 INCLUDE "params" 114 INCLUDE "CTLBLK.comm" 115 C Page 2 Source Listing CALDWP2 2014-12-17 20:47 CALDWP2.f 116 C DECLARE VARIABLES. 117 C 132 REAL EVP(IM,JM),P1D(IM,JM),Q1D(IM,JM),TDWP(IM,JM),T1D(IM,JM) 133 C 134 C 135 C**************************************************************************** 136 C START CALDWP2 HERE. 137 C 138 C COMPUTE VAPOR PRESSURE. CONVERT TO CENITBARS. 139 C 140 DO J=JSTA,JEND 141 DO I=1,IM 142 EVP(I,J) = P1D(I,J)*Q1D(I,J)/(EPS+ONEPS*Q1D(I,J)) 143 EVP(I,J) = EVP(I,J)*D001 144 EVP(I,J) = AMAX1(H1M12,EVP(I,J)) 145 ENDDO 146 ENDDO 147 C 148 C COMPUTE DEWPOINT TEMPERATURE. 149 C 150 CALL DEWPOINT(EVP,TDWP,IM,JM) 151 C 152 C ENSURE DEWPOINT TEMPERATURE DOES NOT EXCEED AMBIENT TEMPERATURE. 153 C 154 DO J=JSTA,JEND 155 DO I=1,IM 156 IF (TDWP(I,J).GT.T1D(I,J)) TDWP(I,J) = T1D(I,J) 157 ENDDO 158 ENDDO 159 C 160 C END OF ROUTINE. 161 C 162 RETURN 163 END Page 3 Source Listing CALDWP2 2014-12-17 20:47 Entry Points CALDWP2.f ENTRY POINTS Name caldwp2_ 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 AMAX1 Func 72 scalar 72 BETA Param 37 R(4) 4 scalar 37 BTG Param 37 R(4) 4 scalar CALDWP2 Subr 1 CAPA Param 28 R(4) 4 scalar 35 CFT0 Param 33 R(4) 4 scalar 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 D001 Param 25 R(4) 4 scalar 71 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 DEWPOINT Subr 78 78 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 ELWV Param 36 R(4) 4 scalar 36 EPS Param 39 R(4) 4 scalar 39,70 EPSQ2 Param 40 R(4) 4 scalar ERAD Param 34 R(4) 4 scalar EVP Local 60 R(4) 4 2 91719 70,71,72,78 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 H100 Param 27 R(4) 4 scalar Page 4 Source Listing CALDWP2 2014-12-17 20:47 Symbol Table CALDWP2.f Name Object Declared Type Bytes Dimen Elements Attributes References 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 72 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 I Local 69 I(4) 4 scalar 69,70,71,72,83,84 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,60,69,78,83 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 68 I(4) 4 scalar 68,70,71,72,82,84 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 68,82 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 4,6,9,12,15,60,78 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 68,82 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 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 Page 5 Source Listing CALDWP2 2014-12-17 20:47 Symbol Table CALDWP2.f Name Object Declared Type Bytes Dimen Elements Attributes References 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 70 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 70 PARALLEL Common 9 8240 PQ0 Param 41 R(4) 4 scalar Q1D Dummy 1 R(4) 4 2 91719 ARG,INOUT 70 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 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 T1D Dummy 1 R(4) 4 2 91719 ARG,INOUT 84 TDWP Dummy 1 R(4) 4 2 91719 ARG,INOUT 78,84 TFRZ Param 39 R(4) 4 scalar Page 6 Source Listing CALDWP2 2014-12-17 20:47 Subprograms/Common Blocks CALDWP2.f SUBPROGRAMS/COMMON BLOCKS Name Object Declared Type Bytes Dimen Elements Attributes References CALDWP2 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 7 Source Listing CALDWP2 2014-12-17 20:47 CALDWP2.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 : CALDWP2.lst -o filename : none COMPILER: Intel(R) Fortran 12.1-2100