Page 1 Source Listing DEWPOINT 2014-12-17 20:47 DEWPOINT.f 1 SUBROUTINE DEWPOINT( VP, TD, IM,JM ) 2 C$$$ SUBPROGRAM DOCUMENTATION BLOCK 3 C 4 C SUBPROGRAM: DEWPOINT COMPUTES DEWPOINTS FROM VAPOR PRESSURE 5 C PRGMMR: J TUCCILLO ORG: W/NP2 DATE: 90-05-19 6 C 7 C ABSTRACT: COMPUTES THE DEWPOINTS FOR THE N VALUES 8 C OF VAPOR PRESSURE IN ARRAY VP. 9 C . 10 C . THE FORMULA 11 C . 12 C . VP = 0.611 * (X**A) * EXP( (A+B)*(1-X) ) 13 C . 14 C IS USED TO GET DEWPOINT TEMPERATURE T, WHERE 15 C . 16 C X = T3/T, T3=TRIPLE PT TEMPERATURE, 17 C VP=VAPOR PRESSURE IN CBS, 0.611=VP AT T3, 18 C A=(SPEC. HT. OF WATER-CSUBP OF VAPOR)/GAS CONST OF VAPOR 19 C . AND 20 C B=LATENT HEAT AT T3/(GAS CONST OF VAPOR TIMES T3). 21 C . 22 C ON THE FIRST CALL, A TABLE TDP IS CONSTRUCTED GIVING 23 C DEWPOINT AS A FUNCTION OF VAPOR PRESSURE. 24 C . 25 C VALUES OF VP LESS THAN THE FIRST TABLE ENTRY 26 C (RVP1 IN THE CODE) WILL BE GIVEN DEWPOINTS FOR 27 C THAT BEGINNING VALUE. SIMILARLY , VP VALUES THAT 28 C EXCEED THE MAXIMUM TABLE VALUE (RVP2 IN THE CODE) 29 C WILL BE ASSIGNED DEWPOINTS FOR THAT MAXIMUM VALUE. 30 C . 31 C THE VALUES 0.02 AND 8.0 FOR RVP1 AND RVP2 YIELD 32 C DEWPOINTS OF 233.6K AND 314.7K,RESPECTIVELY. 33 C . 34 C 35 C PROGRAM HISTORY LOG: 36 C 90-05-19 J TUCCILLO 37 C 93-05-12 R TREADON - EXPANDED TABLE SIZE AND RESET 38 C RANGE OF PRESSURES COVERED BY 39 C TABLE. 40 C 98-06-12 T BLACK - CONVERSION FROM 1-D TO 2-D 41 C 00-01-04 JIM TUCCILLO - MPI VERSION 42 C 43 C USAGE: CALL DEWPOINT( VP, TD, N ) 44 C INPUT ARGUMENT LIST: 45 C VP - ARRAY OF N VAPOR PRESSURES(CENTIBARS) 46 C IM,JM - DIMENSIONS OF THE INPUT ARRAY VP. 47 C 48 C OUTPUT ARGUMENT LIST: 49 C TD - DEWPOINT IN DEGREES ABSOLUTE 50 C 51 C SUBPROGRAMS CALLED: 52 C LIBRARY: 53 C 54 C ATTRIBUTES: 55 C LANGUAGE: FORTRAN 56 C MACHINE: CRAY C-90 57 C Page 2 Source Listing DEWPOINT 2014-12-17 20:47 DEWPOINT.f 58 C$$$ 59 INCLUDE "CTLBLK.comm" 60 C 61 C NT IS THE TABLE SIZE 76 PARAMETER (NT=2000) 77 C...TRANSLATED BY FPP 3.00Z36 11/09/90 14:48:53 78 C...SWITCHES: OPTON=I47,OPTOFF=VAE0 79 DIMENSION TDP(NT),VP(IM,JM),TD(IM,JM) 80 CX LOGICAL LFRST 81 CX DATA LFRST / .TRUE. / 82 C 83 C 84 CX IF ( LFRST ) THEN 85 C 86 CX LFRST = .FALSE. 87 C PREPARE THE TABLE (TDP=DEWPT AS FCN OF VAPOR PRESS). 88 C RANGE IN CENTIBARS IS FROM RVP1 THRU RVP2 89 coff RVP1=0.02E0 90 coff RVP2=8.E0 91 rvp1=0.0001E0 92 rvp2=10.E0 93 C THE TRIPLE POINT 94 RT3=273.16E0 95 C VAPOR PRESS AT THE TRIPLE POINT 96 RVP3=0.611E0 97 RLOG3=LOG(RVP3) 98 C (SPEC HT OF WATER -CSUBP OF VAPOR)/GAS CONST OF VAPOR. 99 RA=5.0065E0 100 C LATENT HEAT AT T3/(GAS CONST OF VAPOR * TRIPLE PT TEMP). 101 RB=19.83923E0 102 RAPB=RA+RB 103 C CRITERION FOR CONVERGENCE OF NEWTON ITERATION 104 RTEST=1.E-6 105 CMEB RTEST=1.E-8 ! PROBABLY WON'T CONVERGE WITH 32-BIT AT THIS CRITERION 106 C 107 RNT=FLOAT(NT) 108 C TABLE INCREMENT IN VAPOR PRESS 109 RDVP=(RVP2-RVP1)/(RNT-1.E0) 110 C RGS WILL BE THE GUESSED VALUE OF (T3 / DEWPOINT) 111 RGS=1.E0 112 RVP=RVP1-RDVP 113 C 114 DO 20 NN=1,NT 115 RVP=RVP+RDVP 116 RLVP=LOG(RVP)-RLOG3-RAPB 117 C ***** ENTER NEWTON ITERATION LOOP 118 10 RN=RA*LOG(RGS)-RAPB*RGS-RLVP 119 C THAT WAS VALUE OF FUNCTION 120 C NOW GET ITS DERIVATIVE 121 RD=(RA/RGS)-RAPB 122 C THE DESIRED CHANGE IN THE GUESS 123 RCH=RN/RD 124 IF( ABS(RCH) .LT. RTEST ) GO TO 15 125 C NEED MORE ITERATIONS 126 RGS=RGS-RCH 127 GO TO 10 128 C ***** Page 3 Source Listing DEWPOINT 2014-12-17 20:47 DEWPOINT.f 129 C HAVE ACCURATE ENUF VALUE OF RGS=T3/DEWPOINT. 130 15 RT=RT3/RGS 131 TDP(NN)=RT 132 C 133 20 CONTINUE 134 C PRINT 25,RVP1,RVP2,TDP(1),TDP(NT) 135 C 25 FORMAT(/'0', 'IN SUBROUTINE DEWPOINT, THE DEWPT TABLE ', 136 C 1 'HAS RVP1=', 1PE13.6, ', RVP2=', 1PE13.6, 137 C 2 ', TDP(1)=', 1PE13.6, ', AND TDP(NT)=', 138 C 3 1PE13.6, '.'/) 139 C CONSTANTS FOR USING THE TABLE 140 A=1./RDVP 141 B=1.-A*RVP1 142 DNTM1=FLOAT(NT) -.01 143 C 144 CX END IF 145 C 146 C *********** ENTER TO USE THE TABLE. ************ 147 C 148 DO J=JSTA,JEND 149 DO I=1,IM 150 W1=AMIN1(AMAX1((A*VP(I,J)+B),1.0),DNTM1) 151 W2=AINT(W1) 152 JNT=INT(W2) 153 TD(I,J) = TDP(JNT)+(W1-W2)*(TDP(JNT+1)-TDP(JNT)) 154 ENDDO 155 ENDDO 156 C 157 C 158 RETURN 159 END Page 4 Source Listing DEWPOINT 2014-12-17 20:47 Entry Points DEWPOINT.f ENTRY POINTS Name dewpoint_ SYMBOL CROSS REFERENCE Name Object Declared Type Bytes Dimen Elements Attributes References 10 Label 104 113 15 Label 116 110 20 Label 119 100 A Local 126 R(4) 4 scalar 126,127,136 ABS Func 110 scalar 110 AINT Func 137 scalar 137 AMAX1 Func 136 scalar 136 AMIN1 Func 136 scalar 136 B Local 127 R(4) 4 scalar 127,136 CTLBLK Common 2 96 DEWPOINT Subr 1 DNTM1 Local 128 R(4) 4 scalar 128,136 DT Scalar 4 R(4) 4 scalar COM FIRST Scalar 3 R(4) 4 scalar COM FLOAT Func 93 scalar 93,128 I Local 135 I(4) 4 scalar 135,136,139 ICNT Scalar 7 I(4) 4 1 1024 COM IDAT Scalar 3 I(4) 4 1 3 COM 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 IHRST Scalar 3 I(4) 4 scalar COM IM Dummy 1 I(4) 4 scalar ARG,INOUT 65,135 INT Func 138 scalar 138 IOUT Scalar 3 I(4) 4 scalar COM IUP Scalar 7 I(4) 4 scalar COM J Local 134 I(4) 4 scalar 134,136,139 JEND Scalar 6 I(4) 4 scalar COM 134 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 JM Dummy 1 I(4) 4 scalar ARG,INOUT 65 JNT Local 138 I(4) 4 scalar 138,139 JSTA Scalar 6 I(4) 4 scalar COM 134 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 LIST Scalar 3 I(4) 4 scalar COM LOG Func 83 scalar 83,102,104 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 Page 5 Source Listing DEWPOINT 2014-12-17 20:47 Symbol Table DEWPOINT.f Name Object Declared Type Bytes Dimen Elements Attributes References NFCST Scalar 3 I(4) 4 scalar COM NN Local 100 I(4) 4 scalar 100,117 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 NT Param 62 I(4) 4 scalar 65,93,100,128 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 PARALLEL Common 9 8240 RA Local 85 R(4) 4 scalar 85,88,104,107 RAPB Local 88 R(4) 4 scalar 88,102,104,107 RB Local 87 R(4) 4 scalar 87,88 RCH Local 109 R(4) 4 scalar 109,110,112 RD Local 107 R(4) 4 scalar 107,109 RDVP Local 95 R(4) 4 scalar 95,98,101,126 RESTRT Scalar 3 R(4) 4 scalar COM RGS Local 97 R(4) 4 scalar 97,104,107,112,116 RLOG3 Local 83 R(4) 4 scalar 83,102 RLVP Local 102 R(4) 4 scalar 102,104 RN Local 104 R(4) 4 scalar 104,109 RNT Local 93 R(4) 4 scalar 93,95 RT Local 116 R(4) 4 scalar 116,117 RT3 Local 80 R(4) 4 scalar 80,116 RTEST Local 90 R(4) 4 scalar 90,110 RUN Scalar 3 R(4) 4 scalar COM RVP Local 98 R(4) 4 scalar 98,101,102 RVP1 Local 77 R(4) 4 scalar 77,95,98,127 RVP2 Local 78 R(4) 4 scalar 78,95 RVP3 Local 82 R(4) 4 scalar 82,83 SIGMA Scalar 3 R(4) 4 scalar COM TD Dummy 1 R(4) 4 2 0 ARG,INOUT 139 TDP Local 65 R(4) 4 1 2000 117,139 VP Dummy 1 R(4) 4 2 0 ARG,INOUT 136 W1 Local 136 R(4) 4 scalar 136,137,139 W2 Local 137 R(4) 4 scalar 137,138,139 Page 6 Source Listing DEWPOINT 2014-12-17 20:47 Subprograms/Common Blocks DEWPOINT.f SUBPROGRAMS/COMMON BLOCKS Name Object Declared Type Bytes Dimen Elements Attributes References CTLBLK Common 2 96 DEWPOINT Subr 1 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 DEWPOINT 2014-12-17 20:47 DEWPOINT.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 : DEWPOINT.lst -o filename : none COMPILER: Intel(R) Fortran 12.1-2100