Page 1 Source Listing CALDRG 2014-12-17 20:47 CALDRG.f 1 SUBROUTINE CALDRG(DRAGCO) 2 C$$$ SUBPROGRAM DOCUMENTATION BLOCK 3 C . . . 4 C SUBPROGRAM: CALDRG COMPUTE DRAG COEFFICIENT 5 C PRGRMMR: TREADON ORG: W/NP2 DATE: 93-09-01 6 C 7 C ABSTRACT: THIS ROUTINE COMPUTES A SURFACE LAYER DRAG 8 C COEFFICIENT USING EQUATION (7.4.1A) IN "AN INTRODUCTION 9 C TO BOUNDARY LAYER METEOROLOGY" BY STULL (1988, KLUWER 10 C ACADEMIC PUBLISHERS). 11 C . 12 C 13 C PROGRAM HISTORY LOG: 14 C 93-09-01 RUSS TREADON 15 C 98-06-15 T BLACK - CONVERSION FROM 1-D TO 2-D 16 C 00-01-04 JIM TUCCILLO - MPI VERSION 17 C 18 C USAGE: CALL CALDRG(DRAGCO) 19 C INPUT ARGUMENT LIST: 20 C NONE 21 C 22 C OUTPUT ARGUMENT LIST: 23 C DRAGCO - SURFACE LAYER DRAG COEFFICIENT 24 C 25 C OUTPUT FILES: 26 C NONE 27 C 28 C SUBPROGRAMS CALLED: 29 C UTILITIES: 30 C NONE 31 C LIBRARY: 32 C COMMON - LOOPS 33 C VRBLS 34 C SRFDSP 35 C PVRBLS 36 C 37 C ATTRIBUTES: 38 C LANGUAGE: FORTRAN 90 39 C MACHINE : CRAY C-90 40 C$$$ 41 C 42 C 43 C INCLUDE/SET PARAMETERS. 44 C 45 INCLUDE "parmeta" 63 INCLUDE "params" 64 C 65 C DECLARE VARIABLES. 107 REAL DRAGCO(IM,JM) 108 C 109 C INCLUDE COMMON BLOCKS. 110 C 111 INCLUDE "LOOPS.comm" 119 INCLUDE "VRBLS.comm" 126 INCLUDE "SRFDSP.comm" 132 INCLUDE "PVRBLS.comm" 144 INCLUDE "INDX.comm" Page 2 Source Listing CALDRG 2014-12-17 20:47 CALDRG.f 150 INCLUDE "CTLBLK.comm" 151 C 152 C 153 C******************************************************************** 154 C START CALDRG HERE. 155 C 156 C INITIALIZE DRAG COEFFICIENT ARRAY TO ZERO. 157 C 172 DO J=JSTA,JEND 173 DO I=1,IM 174 DRAGCO(I,J) = D00 175 ENDDO 176 ENDDO 177 C 178 C COMPUTE A SURFACE LAYER DRAG COEFFICIENT. 179 C 180 CALL EXCH(UZ0) 181 CALL EXCH(VZ0) 182 C 183 C ASSUMET THAT U AND V HAVE HALOS 184 C 185 DO 30 J=JSTA_M,JEND_M 186 DO 30 I=2,IM-1 187 C 188 C COMPUTE A MEAN MASS POINT WIND IN THE 189 C FIRST ATMOSPHERIC ETA LAYER. 190 C 191 LMHK=LMH(I,J) 192 IE=I+IHE(J) 193 IW=I+IHW(J) 194 SUMU=U(IE,J,LMV(IE,J))+U(IW,J,LMV(IW,J))+U(I,J-1,LMV(I,J-1)) 195 1 +U(I,J+1,LMV(I,J+1)) 196 SUMV=V(IE,J,LMV(IE,J))+V(IW,J,LMV(IW,J))+V(I,J-1,LMV(I,J-1)) 197 1 +V(I,J+1,LMV(I,J+1)) 198 ULMH=D25*SUMU 199 VLMH=D25*SUMV 200 C 201 C COMPUTE A MEAN MASS POINT WIND AT HEIGHT Z0. 202 C 203 UZ0H=D25*(UZ0(IE,J)+UZ0(IW,J)+UZ0(I,J-1)+UZ0(I,J+1)) 204 VZ0H=D25*(VZ0(IE,J)+VZ0(IW,J)+VZ0(I,J-1)+VZ0(I,J+1)) 205 C 206 C COMPUTE A MEAN MASS POINT WIND SPEED BETWEEN THE 207 C FIRST ATMOSPHERIC ETA LAYER AND Z0. 208 C 209 UBAR=D50*(ULMH+UZ0H) 210 VBAR=D50*(VLMH+VZ0H) 211 WSPDSQ=UBAR*UBAR+VBAR*VBAR 212 C 213 C COMPUTE A DRAG COEFFICIENT. 214 C 215 USTRSQ=USTAR(I,J)*USTAR(I,J) 216 DRAGCO(I,J)=USTRSQ/WSPDSQ 217 C 218 30 CONTINUE 219 C 220 C END OF ROUTINE. Page 3 Source Listing CALDRG 2014-12-17 20:47 CALDRG.f 221 C 222 RETURN 223 END ENTRY POINTS Name caldrg_ SYMBOL CROSS REFERENCE Name Object Declared Type Bytes Dimen Elements Attributes References 30 Label 112 79,80 A2 Param 41 R(4) 4 scalar A3 Param 41 R(4) 4 scalar A4 Param 41 R(4) 4 scalar ACCLIQ Scalar 7 R(4) 4 2 91719 COM ACPREC Scalar 7 R(4) 4 2 91719 COM AD05 Param 33 R(4) 4 scalar AKHS Scalar 3 R(4) 4 2 91719 COM AKMS Scalar 3 R(4) 4 2 91719 COM BETA Param 37 R(4) 4 scalar 37 BTG Param 37 R(4) 4 scalar CALDRG Subr 1 CAPA Param 28 R(4) 4 scalar 35 CFT0 Param 33 R(4) 4 scalar CLDEFI Scalar 6 R(4) 4 2 91719 COM 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 CUPREC Scalar 7 R(4) 4 2 91719 COM D00 Param 25 R(4) 4 scalar 68 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 92,93,97,98 D50 Param 25 R(4) 4 scalar 103,104 D608 Param 29 R(4) 4 scalar DP10M Param 35 R(4) 4 scalar DRAGCO Dummy 1 R(4) 4 2 91719 ARG,INOUT 68,110 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 4 Source Listing CALDRG 2014-12-17 20:47 Symbol Table CALDRG.f Name Object Declared Type Bytes Dimen Elements Attributes References ERAD Param 34 R(4) 4 scalar EXCH Subr 74 74,75 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 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 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 67 I(4) 4 scalar 67,68,80,85,86,87,88,89,90,91,97,9 8,109,110 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 IE Local 86 I(4) 4 scalar 86,88,90,97,98 IGSTL Param 9 I(4) 4 scalar 14 IGSTR Param 9 I(4) 4 scalar 14 IHE Scalar 3 I(4) 4 1 387 COM 86 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 IHW Scalar 3 I(4) 4 1 387 COM 87 IM Param 3 I(4) 4 scalar 3,4,5,6,7,9,11,14,49,67,80 IM1 Param 6 I(4) 4 scalar IMT Param 4 I(4) 4 scalar INDX Common 2 6192 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 IVE Scalar 3 I(4) 4 1 387 COM IVHA Scalar 3 I(4) 4 1 760 COM IVLA Scalar 3 I(4) 4 1 760 COM IVW Scalar 3 I(4) 4 1 387 COM IW Local 87 I(4) 4 scalar 87,88,90,97,98 J Local 66 I(4) 4 scalar 66,68,79,85,86,87,88,89,90,91,97,9 8,109,110 JAM Param 6 I(4) 4 scalar 3 Page 5 Source Listing CALDRG 2014-12-17 20:47 Symbol Table CALDRG.f Name Object Declared Type Bytes Dimen Elements Attributes References JDIM1 Param 15 I(4) 4 scalar JDIM2 Param 15 I(4) 4 scalar JEND Scalar 6 I(4) 4 scalar COM 66 JEND_2U Scalar 8 I(4) 4 scalar COM JEND_M Scalar 7 I(4) 4 scalar COM 79 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 3,4,5,6,7,9,12,15,49 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 66 JSTA_2L Scalar 8 I(4) 4 scalar COM JSTA_M Scalar 7 I(4) 4 scalar COM 79 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 5,8,9,19 LM1 Param 8 I(4) 4 scalar LMH Scalar 5 I(4) 4 2 91719 COM 85 LMHK Local 85 I(4) 4 scalar 85 LMV Scalar 5 I(4) 4 2 91719 COM 88,89,90,91 LOOPS Common 2 748952 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 Page 6 Source Listing CALDRG 2014-12-17 20:47 Symbol Table CALDRG.f Name Object Declared Type Bytes Dimen Elements Attributes References 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 PD Scalar 3 R(4) 4 2 91719 COM PQ0 Param 41 R(4) 4 scalar PREC Scalar 7 R(4) 4 2 91719 COM PSHLTR Scalar 4 R(4) 4 2 91719 COM PVRBLS Common 2 23480064 Q Scalar 5 R(4) 4 3 4127355 COM Q10 Scalar 3 R(4) 4 2 91719 COM Q2 Scalar 9 R(4) 4 3 4127355 COM Q30 Scalar 5 R(4) 4 2 91719 COM QS Scalar 3 R(4) 4 2 91719 COM QSHLTR Scalar 4 R(4) 4 2 91719 COM QWBS Scalar 5 R(4) 4 2 91719 COM QZ0 Scalar 4 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 RF Scalar 5 R(4) 4 2 91719 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 SI Scalar 6 R(4) 4 2 91719 COM SIGMA Scalar 3 R(4) 4 scalar COM SMTHA Param 15 R(4) 4 scalar SMTHB Param 15 R(4) 4 scalar SNO Scalar 6 R(4) 4 2 91719 COM SRFDSP Common 2 4402512 SUMU Local 88 R(4) 4 scalar 88,92 SUMV Local 90 R(4) 4 scalar 90,93 T Scalar 5 R(4) 4 3 4127355 COM TFRZ Param 39 R(4) 4 scalar TH10 Scalar 3 R(4) 4 2 91719 COM TH30 Scalar 5 R(4) 4 2 91719 COM THS Scalar 3 R(4) 4 2 91719 COM THZ0 Scalar 4 R(4) 4 2 91719 COM TSHLTR Scalar 4 R(4) 4 2 91719 COM TWBS Scalar 5 R(4) 4 2 91719 COM U Scalar 5 R(4) 4 3 4127355 COM 88,89 U10 Scalar 3 R(4) 4 2 91719 COM U30 Scalar 5 R(4) 4 2 91719 COM UBAR Local 103 R(4) 4 scalar 103,105 ULMH Local 92 R(4) 4 scalar 92,103 USTAR Scalar 3 R(4) 4 2 91719 COM 109 USTRSQ Local 109 R(4) 4 scalar 109,110 Page 7 Source Listing CALDRG 2014-12-17 20:47 Symbol Table CALDRG.f Name Object Declared Type Bytes Dimen Elements Attributes References UZ0 Scalar 4 R(4) 4 2 91719 COM 74,97 UZ0H Local 97 R(4) 4 scalar 97,103 V Scalar 5 R(4) 4 3 4127355 COM 90,91 V10 Scalar 3 R(4) 4 2 91719 COM V30 Scalar 5 R(4) 4 2 91719 COM VBAR Local 104 R(4) 4 scalar 104,105 VLMH Local 93 R(4) 4 scalar 93,104 VRBLS Common 2 67138308 VZ0 Scalar 4 R(4) 4 2 91719 COM 75,98 VZ0H Local 98 R(4) 4 scalar 98,104 WSPDSQ Local 105 R(4) 4 scalar 105,110 Z0 Scalar 3 R(4) 4 2 91719 COM Page 8 Source Listing CALDRG 2014-12-17 20:47 Subprograms/Common Blocks CALDRG.f SUBPROGRAMS/COMMON BLOCKS Name Object Declared Type Bytes Dimen Elements Attributes References CALDRG Subr 1 CTLBLK Common 2 96 INDX Common 2 6192 LOOPS Common 2 748952 PARALLEL Common 9 8240 PVRBLS Common 2 23480064 SRFDSP Common 2 4402512 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 9 Source Listing CALDRG 2014-12-17 20:47 CALDRG.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 : CALDRG.lst -o filename : none COMPILER: Intel(R) Fortran 12.1-2100