Page 1 Source Listing CALVIS 2014-12-17 20:47 CALVIS.f 1 c**********************************************************************c 2 SUBROUTINE CALVIS(QV,QC,QR,QI,QS,TT,PP,VIS) 3 c 4 c This routine computes horizontal visibility at the 5 c surface or lowest model layer, from qc, qr, qi, and qs. 6 c qv--water vapor mixing ratio (kg/kg) 7 c qc--cloud water mixing ratio (kg/kg) 8 c qr--rain water mixing ratio (kg/kg) 9 c qi--cloud ice mixing ratio (kg/kg) 10 c qs--snow mixing ratio (kg/kg) 11 c tt--temperature (k) 12 c pp--pressure (Pa) 13 c 14 c If iice=0: 15 c qprc=qr qrain=qr and qclw=qc if T>0C 16 c qcld=qc =0 =0 if T<0C 17 c qsnow=qs and qclice=qc if T<0C 18 c =0 =0 if T>0C 19 c If iice=1: 20 c qprc=qr+qs qrain=qr and qclw=qc 21 c qcld=qc+qi qsnow=qs and qclice=qc 22 c 23 c Independent of the above definitions, the scheme can use different 24 c assumptions of the state of hydrometeors: 25 c meth='d': qprc is all frozen if T<0, liquid if T>0 26 c meth='b': Bocchieri scheme used to determine whether qprc 27 c is rain or snow. A temperature assumption is used to 28 c determine whether qcld is liquid or frozen. 29 c meth='r': Uses the four mixing ratios qrain, qsnow, qclw, 30 c and qclice 31 c 32 c The routine uses the following 33 c expressions for extinction coefficient, beta (in km**-1), 34 c with C being the mass concentration (in g/m**3): 35 c 36 c cloud water: beta = 144.7 * C ** (0.8800) 37 c rain water: beta = 2.24 * C ** (0.7500) 38 c cloud ice: beta = 327.8 * C ** (1.0000) 39 c snow: beta = 10.36 * C ** (0.7776) 40 c 41 c These expressions were obtained from the following sources: 42 c 43 c for cloud water: from Kunkel (1984) 44 c for rainwater: from M-P dist'n, with No=8e6 m**-4 and 45 c rho_w=1000 kg/m**3 46 c for cloud ice: assume randomly oriented plates which follow 47 c mass-diameter relationship from Rutledge and Hobbs (1983) 48 c for snow: from Stallabrass (1985), assuming beta = -ln(.02)/vis 49 c 50 c The extinction coefficient for each water species present is 51 c calculated, and then all applicable betas are summed to yield 52 c a single beta. Then the following relationship is used to 53 c determine visibility (in km), where epsilon is the threshhold 54 c of contrast, usually taken to be .02: 55 c 56 c vis = -ln(epsilon)/beta [found in Kunkel (1984)] 57 c Page 2 Source Listing CALVIS 2014-12-17 20:47 CALVIS.f 58 C------------------------------------------------------------------ 59 INCLUDE "parmeta" 77 INCLUDE "params" 119 INCLUDE "CTLBLK.comm" 120 C 135 REAL QV(IM,JM),QC(IM,JM),QI(IM,JM),QR(IM,JM),QS(IM,JM) 136 1, TT(IM,JM),PP(IM,JM),VIS(IM,JM) 137 CHARACTER METH*1 138 C------------------------------------------------------------------ 139 C------------------------------------------------------------------ 140 CELKEL=273.15 141 TICE=CELKEL-10. 142 COEFLC=144.7 143 COEFLP=2.24 144 COEFFC=327.8 145 COEFFP=10.36 146 EXPONLC=0.8800 147 EXPONLP=0.7500 148 EXPONFC=1.0000 149 EXPONFP=0.7776 150 CONST1=-LOG(.02) 151 RHOICE=917. 152 RHOWAT=1000. 153 C 154 DO J=JSTA,JEND 155 DO I=1,IM 156 c IF(IICE.EQ.0)THEN 157 c QPRC=QR 158 c QCLD=QC 159 c IF(TT.LT.CELKEL)THEN 160 c QRAIN=0. 161 c QSNOW=QPRC 162 c QCLW=0. 163 c QCLICE=QCLD 164 c ELSE 165 c QRAIN=QPRC 166 c QSNOW=0. 167 c QCLW=QCLD 168 c QCLICE=0. 169 c ENDIF 170 c ELSE 171 QPRC=QR(I,J)+QS(I,J) 172 QCLD=QC(I,J)+QI(I,J) 173 QRAIN=QR(I,J) 174 QSNOW=QS(I,J) 175 QCLW=QC(I,J) 176 QCLICE=QI(I,J) 177 c ENDIF 178 c TV=VIRTUAL(TT,QV) 179 TV=TT(I,J)*(H1+D608*QV(I,J)) 180 RHOAIR=PP(I,J)/(RD*TV) 181 c IF(METH.EQ.'D')THEN 182 c IF(TT.LT.CELKEL)THEN 183 c VOVERMD=(1.+QV)/RHOAIR+(QPRC+QCLD)/RHOICE 184 c CONCLC = 0. 185 c CONCLP = 0. 186 c CONCFC = QCLD/VOVERMD*1000. Page 3 Source Listing CALVIS 2014-12-17 20:47 CALVIS.f 187 c CONCFP = QPRC/VOVERMD*1000. 188 c ELSE 189 c VOVERMD=(1.+QV)/RHOAIR+(QPRC+QCLD)/RHOWAT 190 c CONCLC = QCLD/VOVERMD*1000. 191 c CONCLP = QPRC/VOVERMD*1000. 192 c CONCFC = 0. 193 c CONCFP = 0. 194 c ENDIF 195 c ELSEIF(METH.EQ.'B')THEN 196 c IF(TT.LT.TICE)THEN 197 c VOVERMD=(1.+QV)/RHOAIR+(QPRC+QCLD)/RHOICE 198 c CONCLC = 0. 199 c CONCLP = 0. 200 c CONCFC = QCLD/VOVERMD*1000. 201 c CONCFP = QPRC/VOVERMD*1000. 202 c ELSEIF(PRSNOW.GE.50.)THEN 203 c VOVERMD=(1.+QV)/RHOAIR+QPRC/RHOICE+QCLD/RHOWAT 204 c CONCLC = QCLD/VOVERMD*1000. 205 c CONCLP = 0. 206 c CONCFC = 0. 207 c CONCFP = QPRC/VOVERMD*1000. 208 c ELSE 209 c VOVERMD=(1.+QV)/RHOAIR+(QPRC+QCLD)/RHOWAT 210 c CONCLC = QCLD/VOVERMD*1000. 211 c CONCLP = QPRC/VOVERMD*1000. 212 c CONCFC = 0. 213 c CONCFP = 0. 214 c ENDIF 215 c ELSEIF(METH.EQ.'R')THEN 216 VOVERMD=(1.+QV(I,J))/RHOAIR+(QCLW+QRAIN)/RHOWAT+ 217 1 (QCLICE+QSNOW)/RHOICE 218 CONCLC=QCLW/VOVERMD*1000. 219 CONCLP=QRAIN/VOVERMD*1000. 220 CONCFC=QCLICE/VOVERMD*1000. 221 CONCFP=QSNOW/VOVERMD*1000. 222 c ENDIF 223 BETAV=COEFFC*CONCFC**EXPONFC+COEFFP*CONCFP**EXPONFP 224 1 +COEFLC*CONCLC**EXPONLC+COEFLP*CONCLP**EXPONLP 225 2 +1.E-10 226 c CHANGED GSM 3-10-00 --> no point in distinguishing values 227 c above 20 km, so make that value the max (prev max was 80) 228 VIS(I,J)=1.E3*MIN(20.,CONST1/BETAV) ! max of 20km 229 ENDDO 230 ENDDO 231 C 232 RETURN 233 END Page 4 Source Listing CALVIS 2014-12-17 20:47 Entry Points CALVIS.f ENTRY POINTS Name calvis_ 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 BETA Param 37 R(4) 4 scalar 37 BETAV Local 151 R(4) 4 scalar 151,156 BTG Param 37 R(4) 4 scalar CALVIS Subr 2 CAPA Param 28 R(4) 4 scalar 35 CELKEL Local 68 R(4) 4 scalar 68,69 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 COEFFC Local 72 R(4) 4 scalar 72,151 COEFFP Local 73 R(4) 4 scalar 73,151 COEFLC Local 70 R(4) 4 scalar 70,152 COEFLP Local 71 R(4) 4 scalar 71,152 CONCFC Local 148 R(4) 4 scalar 148,151 CONCFP Local 149 R(4) 4 scalar 149,151 CONCLC Local 146 R(4) 4 scalar 146,152 CONCLP Local 147 R(4) 4 scalar 147,152 CONST1 Local 78 R(4) 4 scalar 78,156 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 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 107 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 EPSQ2 Param 40 R(4) 4 scalar Page 5 Source Listing CALVIS 2014-12-17 20:47 Symbol Table CALVIS.f Name Object Declared Type Bytes Dimen Elements Attributes References ERAD Param 34 R(4) 4 scalar EXPONFC Local 76 R(4) 4 scalar 76,151 EXPONFP Local 77 R(4) 4 scalar 77,151 EXPONLC Local 74 R(4) 4 scalar 74,152 EXPONLP Local 75 R(4) 4 scalar 75,152 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,107 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 I Local 83 I(4) 4 scalar 83,99,100,101,102,103,104,107,108, 144,156 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,63,64,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 82 I(4) 4 scalar 82,99,100,101,102,103,104,107,108, 144,156 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 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,63,64 Page 6 Source Listing CALVIS 2014-12-17 20:47 Symbol Table CALVIS.f Name Object Declared Type Bytes Dimen Elements Attributes References 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 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 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 LOG Func 78 scalar 78 LP1 Param 8 I(4) 4 scalar LSM Param 3 I(4) 4 scalar ME Scalar 6 I(4) 4 scalar COM METH Local 65 CHAR 1 scalar MIN Func 156 scalar 156 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 OVERRC Param 33 R(4) 4 scalar 33 P1000 Param 35 R(4) 4 scalar PARALLEL Common 9 8240 PP Dummy 2 R(4) 4 2 91719 ARG,INOUT 108 PQ0 Param 41 R(4) 4 scalar Page 7 Source Listing CALVIS 2014-12-17 20:47 Symbol Table CALVIS.f Name Object Declared Type Bytes Dimen Elements Attributes References QC Dummy 2 R(4) 4 2 91719 ARG,INOUT 100,103 QCLD Local 100 R(4) 4 scalar 100 QCLICE Local 104 R(4) 4 scalar 104,145,148 QCLW Local 103 R(4) 4 scalar 103,144,146 QI Dummy 2 R(4) 4 2 91719 ARG,INOUT 100,104 QPRC Local 99 R(4) 4 scalar 99 QR Dummy 2 R(4) 4 2 91719 ARG,INOUT 99,101 QRAIN Local 101 R(4) 4 scalar 101,144,147 QS Dummy 2 R(4) 4 2 91719 ARG,INOUT 99,102 QSNOW Local 102 R(4) 4 scalar 102,145,149 QV Dummy 2 R(4) 4 2 91719 ARG,INOUT 107,144 RCAPA Param 35 R(4) 4 scalar RD Param 32 R(4) 4 scalar 32,36,40,108 RDOCP Param 40 R(4) 4 scalar RESTRT Scalar 3 R(4) 4 scalar COM RG Param 36 R(4) 4 scalar RHOAIR Local 108 R(4) 4 scalar 108,144 RHOICE Local 79 R(4) 4 scalar 79,145 RHOWAT Local 80 R(4) 4 scalar 80,144 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 TFRZ Param 39 R(4) 4 scalar TICE Local 69 R(4) 4 scalar 69 TT Dummy 2 R(4) 4 2 91719 ARG,INOUT 107 TV Local 107 R(4) 4 scalar 107,108 VIS Dummy 2 R(4) 4 2 91719 ARG,INOUT 156 VOVERMD Local 144 R(4) 4 scalar 144,146,147,148,149 Page 8 Source Listing CALVIS 2014-12-17 20:47 Subprograms/Common Blocks CALVIS.f SUBPROGRAMS/COMMON BLOCKS Name Object Declared Type Bytes Dimen Elements Attributes References CALVIS Subr 2 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 9 Source Listing CALVIS 2014-12-17 20:47 CALVIS.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 : CALVIS.lst -o filename : none COMPILER: Intel(R) Fortran 12.1-2100