Page 1 Source Listing READPCP 2014-12-17 20:47 /tmpnwprd/ifort9jlJOn.i 1 # 1 "READPCP.F" 2 SUBROUTINE READPCP 3 C **************************************************************** 4 C * * 5 C * PRECIPITATION ASSIMILATION INITIALIZATION. * 6 C * READ IN PRECIP ANALYSIS AND DATA MASK AND SET UP ALL * 7 C * APPROPRIATE VARIABLES. * 8 C * MIKE BALDWIN, MARCH 1994 * 9 C * Adapted to 2-D code, Ying Lin, Mar 1996 * 10 C * If bad precip read, set precip to undefined 1/2012 * 11 C * * 12 C **************************************************************** 13 C----------------------------------------------------------------------- 14 INCLUDE "parmeta" 15 C----------------------------------------------------------------------- 38 INCLUDE "PPTASM.comm" 62 INCLUDE "mpp.h" 121 LOGICAL*1 BITMPCP(IM,JM) 122 INTEGER JPDS(25), JGDS(22), KGDS(22), KPDS(25) 123 124 C----------------------------------------------------------------------- 125 C 126 C READ ANALYSES AND DATA MASK 127 C 128 IF(MYPE.EQ.0) THEN 129 call baopenr(41,'fort.41',iret) 130 write(0,*) 'baopenr on unit41', ' iret=', iret 131 JPDS=-1 132 ENDIF 133 C 134 DO IHR=1,3 135 IF(MYPE.EQ.0) THEN 136 print*,'ihr=',ihr 137 CALL GETGB(41,0,IM*JM,IHR-1,JPDS,JGDS,KF,K,KPDS,KGDS,BITMPCP, 138 & temp1,IRET) 139 write(0,10) ihr, iret, kf, kpds(5), 140 & (kpds(21)*100+kpds(8))/100-1, mod(kpds(8),100),kpds(9), 141 & kpds(10), kpds(11) 142 10 format('ihr=',i1,' iret=',i3,' kf=', i6,' fld=', i3, 143 & 2x,'1h accum from ', 5i2.2) 144 C 1/2012: WNE if bad read 145 if (iret.ne.0) then 146 temp1 = 999.0 147 BITMPCP = .true. 148 endif 149 150 C 151 DO J=1,JM 152 DO I=1,IM 153 c if(temp1(i,j).gt.1.0) print*,'i,j,pcp=',i,j,temp1(i,j) 154 c print*,'i,j,pcp=',i,j,temp1(i,j) 155 IF (BITMPCP(I,J)) THEN 156 C Convert data from mm to m: 157 if (temp1(i,j).lt.0.0) then 158 write(0,*) " precip < 0. will be set to 999.", I,J, temp1(i,j) 159 temp1(i,j) = 999.0 160 else if (temp1(i,j).gt.24.0) then Page 2 Source Listing READPCP 2014-12-17 20:47 READPCP.F 161 write(0,*) " precip >24. will be set to 999.", I,J, temp1(i,j) 162 temp1(i,j) = 999.0 163 else 164 temp1(I,J)=temp1(I,J)*0.001 165 endif 166 c if(temp1(i,j).gt.0.04) then 167 c if(temp1(i,j).gt.0.001) then 168 c if(i.ge. 99.and.i.le.105.and.(j.ge.172.and.j.le.174).or 169 c * .(j.ge.180.and.j.le.184)) then 170 c temp1(i,j)=999. 171 c bitmpcp(i,j)=.f. 172 c endif 173 c if(i.ge. 99.and.i.le.105.and.j.ge.113.and.j.le.119) then 174 c print*,'i,j,pcp=',i,j,temp1(i,j) 175 c temp1(i,j)=999. 176 c bitmpcp(i,j)=.f. 177 c endif 178 ELSE 179 TEMP1(I,J)=999. 180 ENDIF 181 c 182 c Mask out outer 2 rows and outer 2 colums - 19 March 2004 PS and DJ 183 c 184 if(i.le.1.or.i.ge.im-1+mod(j,2) .or. j.le.2.or.j.ge.jm-1) then 185 temp1(i,j)=999. 186 bitmpcp(i,j)=.FALSE. 187 endif 188 ENDDO 189 ENDDO 190 WRITE(0,*) ' FINISHED READING PRECIP ANALYSES AND DATA MASK' 191 ENDIF 192 c 193 c special test 194 c 195 C 196 CALL DSTRB(TEMP1,PPTDAT,1,3,IHR) 197 C 198 ENDDO 199 C 200 RETURN 201 END Page 3 Source Listing READPCP 2014-12-17 20:47 Entry Points READPCP.F ENTRY POINTS Name readpcp_ SYMBOL CROSS REFERENCE Name Object Declared Type Bytes Dimen Elements Attributes References 10 Label 38 35 APREC Scalar 3 R(4) 4 2 17955 COM BAOPENR Subr 25 25 BITMPCP Local 17 L(1) 1 2 91719 33,43,51,82 CLDINC Scalar 6 R(4) 4 2 17955 COM CLDINC7 Scalar 6 R(4) 4 2 17955 COM DSTRB Subr 92 92 G2LI Scalar 55 I(4) 4 1 239 COM G2LJ Scalar 55 I(4) 4 1 389 COM GETGB Subr 33 33 GLB_TABLE Common 37 128 HTMG Scalar 52 R(4) 4 3 4183695 COM I Local 48 I(4) 4 scalar 48,51,53,54,55,56,57,58,60,75,80,8 1,82 IBROW Scalar 22 I(4) 4 scalar COM ICHUNKTAB Scalar 29 I(4) 4 1 8 COM IDIM1 Param 19 I(4) 4 scalar 2,3,5,6 IDIM2 Param 19 I(4) 4 scalar 2,3,5,6 IE_GLB_TABLE Scalar 38 I(4) 4 1 8 COM IE_LOC_TABLE Scalar 28 I(4) 4 1 8 COM IGSTL Param 14 I(4) 4 scalar 19 IGSTR Param 14 I(4) 4 scalar 19 IHR Local 30 I(4) 4 scalar 30,32,33,35,92 ILCOL Scalar 22 I(4) 4 scalar COM ILPAD1 Scalar 23 I(4) 4 scalar COM ILPAD2 Scalar 23 I(4) 4 scalar COM ILPAD3 Scalar 23 I(4) 4 scalar COM ILPAD4 Scalar 23 I(4) 4 scalar COM ILPAD5 Scalar 23 I(4) 4 scalar COM IM Param 3 I(4) 4 scalar 16,17,19,33,43,44,45,46,47,48,51,5 2,57,80 INPES Param 11 I(4) 4 scalar 16,19,27,28,29,38,39 INUMQ Scalar 34 I(4) 4 1 100 COM IQUILT_GROUP Scalar 33 I(4) 4 scalar COM IRCOL Scalar 22 I(4) 4 scalar COM IRET Local 25 I(4) 4 scalar 25,26,34,35,41 IRPAD1 Scalar 24 I(4) 4 scalar COM IRPAD2 Scalar 24 I(4) 4 scalar COM IRPAD3 Scalar 24 I(4) 4 scalar COM IRPAD4 Scalar 24 I(4) 4 scalar COM IRPAD5 Scalar 24 I(4) 4 scalar COM IS_GLB_TABLE Scalar 38 I(4) 4 1 8 COM IS_LOC_TABLE Scalar 27 I(4) 4 1 8 COM ITAIL Param 16 I(4) 4 scalar Page 4 Source Listing READPCP 2014-12-17 20:47 Symbol Table READPCP.F Name Object Declared Type Bytes Dimen Elements Attributes References ITEMP Scalar 47 I(4) 4 2 91719 COM ITEMP2 Scalar 47 I(4) 4 2 91719 COM ITROW Scalar 22 I(4) 4 scalar COM ITSTLOC Scalar 4 I(4) 4 scalar COM J Local 47 I(4) 4 scalar 47,51,53,54,55,56,57,58,60,75,80,8 1,82 JBPAD1 Scalar 25 I(4) 4 scalar COM JBPAD2 Scalar 25 I(4) 4 scalar COM JBPAD3 Scalar 25 I(4) 4 scalar COM JBPAD4 Scalar 25 I(4) 4 scalar COM JBPAD5 Scalar 25 I(4) 4 scalar COM JDIM1 Param 20 I(4) 4 scalar 2,3,5,6 JDIM2 Param 20 I(4) 4 scalar 2,3,5,6 JE_GLB_TABLE Scalar 39 I(4) 4 1 8 COM JE_LOC_TABLE Scalar 28 I(4) 4 1 8 COM JGDS Local 18 I(4) 4 1 22 33 JGSTL Param 15 I(4) 4 scalar 20 JGSTR Param 15 I(4) 4 scalar 20 JM Param 3 I(4) 4 scalar 17,20,33,43,44,45,46,47,51,52,57,8 0 JNPES Param 11 I(4) 4 scalar 17,20,27,28,29,38,39 JPDS Local 18 I(4) 4 1 25 27,33 JS_GLB_TABLE Scalar 39 I(4) 4 1 8 COM JS_LOC_TABLE Scalar 27 I(4) 4 1 8 COM JTAIL Param 17 I(4) 4 scalar JTPAD1 Scalar 26 I(4) 4 scalar COM JTPAD2 Scalar 26 I(4) 4 scalar COM JTPAD3 Scalar 26 I(4) 4 scalar COM JTPAD4 Scalar 26 I(4) 4 scalar COM JTPAD5 Scalar 26 I(4) 4 scalar COM JTSTLOC Scalar 4 I(4) 4 scalar COM K Local 33 I(4) 4 scalar 33 KF Local 33 I(4) 4 scalar 33,35 KGDS Local 18 I(4) 4 1 22 33 KPDS Local 18 I(4) 4 1 25 33,35,36,37 L2GI Scalar 57 I(4) 4 1 239 COM L2GJ Scalar 57 I(4) 4 1 389 COM LM Param 3 I(4) 4 scalar 2,52 LSM Param 3 I(4) 4 scalar MAPPINGS Common 56 5024 MOD Func 36 scalar 36,80 MPI_COMM_COMP Scalar 30 I(4) 4 scalar COM MPI_COMM_INTER Scalar 31 I(4) 4 scalar COM MPI_COMM_INTER_ARRAY Scalar 32 I(4) 4 1 100 COM MPPCOM Common 1 1464 MTSTPE Scalar 4 I(4) 4 scalar COM MYIE Scalar 4 I(4) 4 scalar COM MYIE1 Scalar 4 I(4) 4 scalar COM MYIE1_P1 Scalar 8 I(4) 4 scalar COM MYIE1_P2 Scalar 8 I(4) 4 scalar COM MYIE1_P3 Scalar 8 I(4) 4 scalar COM MYIE1_P4 Scalar 8 I(4) 4 scalar COM MYIE2 Scalar 4 I(4) 4 scalar COM MYIE2_P1 Scalar 9 I(4) 4 scalar COM MYIE_P1 Scalar 7 I(4) 4 scalar COM Page 5 Source Listing READPCP 2014-12-17 20:47 Symbol Table READPCP.F Name Object Declared Type Bytes Dimen Elements Attributes References MYIE_P2 Scalar 7 I(4) 4 scalar COM MYIE_P3 Scalar 7 I(4) 4 scalar COM MYIE_P4 Scalar 7 I(4) 4 scalar COM MYIE_P5 Scalar 7 I(4) 4 scalar COM MYIS Scalar 4 I(4) 4 scalar COM MYIS1 Scalar 4 I(4) 4 scalar COM MYIS1_P1 Scalar 6 I(4) 4 scalar COM MYIS1_P2 Scalar 6 I(4) 4 scalar COM MYIS1_P3 Scalar 6 I(4) 4 scalar COM MYIS1_P4 Scalar 6 I(4) 4 scalar COM MYIS2 Scalar 4 I(4) 4 scalar COM MYIS_P1 Scalar 5 I(4) 4 scalar COM MYIS_P2 Scalar 5 I(4) 4 scalar COM MYIS_P3 Scalar 5 I(4) 4 scalar COM MYIS_P4 Scalar 5 I(4) 4 scalar COM MYIS_P5 Scalar 5 I(4) 4 scalar COM MYJE Scalar 15 I(4) 4 scalar COM MYJE1 Scalar 15 I(4) 4 scalar COM MYJE1_P1 Scalar 17 I(4) 4 scalar COM MYJE1_P2 Scalar 17 I(4) 4 scalar COM MYJE1_P3 Scalar 17 I(4) 4 scalar COM MYJE1_P4 Scalar 17 I(4) 4 scalar COM MYJE2 Scalar 15 I(4) 4 scalar COM MYJE2_P1 Scalar 18 I(4) 4 scalar COM MYJE2_P2 Scalar 18 I(4) 4 scalar COM MYJE2_P3 Scalar 18 I(4) 4 scalar COM MYJE2_P4 Scalar 18 I(4) 4 scalar COM MYJE3 Scalar 15 I(4) 4 scalar COM MYJE3_P4 Scalar 18 I(4) 4 scalar COM MYJE4 Scalar 15 I(4) 4 scalar COM MYJE4_P1 Scalar 19 I(4) 4 scalar COM MYJE4_P4 Scalar 19 I(4) 4 scalar COM MYJE5 Scalar 15 I(4) 4 scalar COM MYJE5_P1 Scalar 19 I(4) 4 scalar COM MYJE5_P2 Scalar 19 I(4) 4 scalar COM MYJE_P1 Scalar 16 I(4) 4 scalar COM MYJE_P2 Scalar 16 I(4) 4 scalar COM MYJE_P3 Scalar 16 I(4) 4 scalar COM MYJE_P4 Scalar 16 I(4) 4 scalar COM MYJE_P5 Scalar 16 I(4) 4 scalar COM MYJS Scalar 10 I(4) 4 scalar COM MYJS1 Scalar 10 I(4) 4 scalar COM MYJS1_P1 Scalar 12 I(4) 4 scalar COM MYJS1_P2 Scalar 12 I(4) 4 scalar COM MYJS1_P3 Scalar 12 I(4) 4 scalar COM MYJS1_P4 Scalar 12 I(4) 4 scalar COM MYJS2 Scalar 10 I(4) 4 scalar COM MYJS2_P1 Scalar 13 I(4) 4 scalar COM MYJS2_P2 Scalar 13 I(4) 4 scalar COM MYJS2_P3 Scalar 13 I(4) 4 scalar COM MYJS2_P4 Scalar 13 I(4) 4 scalar COM MYJS3 Scalar 10 I(4) 4 scalar COM MYJS3_P4 Scalar 13 I(4) 4 scalar COM MYJS4 Scalar 10 I(4) 4 scalar COM MYJS4_P1 Scalar 14 I(4) 4 scalar COM Page 6 Source Listing READPCP 2014-12-17 20:47 Symbol Table READPCP.F Name Object Declared Type Bytes Dimen Elements Attributes References MYJS4_P4 Scalar 14 I(4) 4 scalar COM MYJS5 Scalar 10 I(4) 4 scalar COM MYJS5_P1 Scalar 14 I(4) 4 scalar COM MYJS5_P2 Scalar 14 I(4) 4 scalar COM MYJS_P1 Scalar 11 I(4) 4 scalar COM MYJS_P2 Scalar 11 I(4) 4 scalar COM MYJS_P3 Scalar 11 I(4) 4 scalar COM MYJS_P4 Scalar 11 I(4) 4 scalar COM MYJS_P5 Scalar 11 I(4) 4 scalar COM MYPE Scalar 2 I(4) 4 scalar COM 24,31 MY_E Scalar 20 I(4) 4 scalar COM MY_IE_GLB Scalar 2 I(4) 4 scalar COM MY_IE_LOC Scalar 3 I(4) 4 scalar COM MY_IS_GLB Scalar 2 I(4) 4 scalar COM MY_IS_LOC Scalar 3 I(4) 4 scalar COM MY_JE_GLB Scalar 2 I(4) 4 scalar COM MY_JE_LOC Scalar 3 I(4) 4 scalar COM MY_JS_GLB Scalar 2 I(4) 4 scalar COM MY_JS_LOC Scalar 3 I(4) 4 scalar COM MY_N Scalar 20 I(4) 4 scalar COM MY_NE Scalar 21 I(4) 4 scalar COM MY_NEB Scalar 21 I(4) 4 1 8 COM MY_NW Scalar 21 I(4) 4 scalar COM MY_S Scalar 20 I(4) 4 scalar COM MY_SE Scalar 21 I(4) 4 scalar COM MY_SW Scalar 21 I(4) 4 scalar COM MY_W Scalar 20 I(4) 4 scalar COM NPES Scalar 2 I(4) 4 scalar COM PHOUR Scalar 2 R(4) 4 2 17955 COM PPTASM Common 1 3878292 PPTDAT Scalar 3 R(4) 4 3 53865 COM 92 READPCP Subr 1 TEMP1 Scalar 43 R(4) 4 2 91719 COM 34,42,53,54,55,56,57,58,60,75,81,9 2 TEMP10 Scalar 45 R(4) 4 2 91719 COM TEMP11 Scalar 45 R(4) 4 2 91719 COM TEMP12 Scalar 45 R(4) 4 2 91719 COM TEMP13 Scalar 46 R(4) 4 2 91719 COM TEMP14 Scalar 46 R(4) 4 2 91719 COM TEMP15 Scalar 46 R(4) 4 2 91719 COM TEMP16 Scalar 46 R(4) 4 2 91719 COM TEMP2 Scalar 43 R(4) 4 2 91719 COM TEMP2X Scalar 51 R(4) 4 2 92971 COM TEMP3 Scalar 43 R(4) 4 2 91719 COM TEMP4 Scalar 43 R(4) 4 2 91719 COM TEMP5 Scalar 44 R(4) 4 2 91719 COM TEMP6 Scalar 44 R(4) 4 2 91719 COM TEMP7 Scalar 44 R(4) 4 2 91719 COM TEMP8 Scalar 44 R(4) 4 2 91719 COM TEMP9 Scalar 45 R(4) 4 2 91719 COM TEMPCOM Common 42 6603768 TLAT Scalar 2 R(4) 4 3 807975 COM TOPO Common 50 17478548 TTVG Scalar 51 R(4) 4 2 92971 COM VAPINC Scalar 5 R(4) 4 2 17955 COM Page 7 Source Listing READPCP 2014-12-17 20:47 Symbol Table READPCP.F Name Object Declared Type Bytes Dimen Elements Attributes References VAPINC7 Scalar 5 R(4) 4 2 17955 COM Page 8 Source Listing READPCP 2014-12-17 20:47 READPCP.F 202 Page 9 Source Listing READPCP 2014-12-17 20:47 Subprograms/Common Blocks READPCP.F SUBPROGRAMS/COMMON BLOCKS Name Object Declared Type Bytes Dimen Elements Attributes References GLB_TABLE Common 37 128 MAPPINGS Common 56 5024 MPPCOM Common 1 1464 PPTASM Common 1 3878292 READPCP Subr 1 TEMPCOM Common 42 6603768 TOPO Common 50 17478548 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 Page 10 Source Listing READPCP 2014-12-17 20:47 READPCP.F 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 precise -fp_model nofast -fp_model nostrict -fp_model nosource -fp_model nodouble -fp_model noextended -fp_model novery_fast -fp_model noexcept -fp_model nono_except -fp_modbits nofp_contract -fp_modbits nono_fp_contract -fp_modbits nofenv_access -fp_modbits nono_fenv_access -fp_modbits nocx_limited_range -fp_modbits nono_cx_limited_range -fp_modbits noprec_div -fp_modbits nono_prec_div -fp_modbits noprec_sqrt -fp_modbits nono_prec_sqrt -fp_modbits noftz -fp_modbits no_ftz -fp_modbits nointrin_limited_range -fp_modbits nono_intrin_limited_range -fp_modbits notrunc_compares -fp_modbits nono_trunc_compares -fp_modbits noieee_nan_compares -fp_modbits nono_ieee_nan_compares -fp_modbits nohonor_f32_conversion -fp_modbits nono_honor_f32_conversion -fp_modbits nohonor_f64_conversion -fp_modbits nono_honor_f64_conversion -fp_modbits nono_x87_copy -fp_modbits nono_no_x87_copy -fp_modbits noexception_semantics -fp_modbits nono_exception_semantics -fp_modbits noprecise_libm_functions -fp_modbits nono_precise_libm_functions -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 : /usrx/local/intel/composerxe/tbb/include/,/usr/include/,./,/opt/ibmhpc/pe1308/mpich2/intel/include64/, /opt/ibmhpc/pe1308/base/include64/,/usrx/local/intel/composerxe/mkl/include/,/usrx/local/intel/composerxe/tbb/include/, /gpfs/tp1/usrx/local/intel/composer_xe_2011_sp1.11.339/compiler/include/intel64/,/gpfs/tp1/usrx/local/intel/composer_xe_2011_sp1.11.339/compiler/include/, /usr/local/include/,/usr/lib/gcc/x86_64-redhat-linux/4.4.7/include/,/usr/include/,/usr/include/ -list filename : READPCP.lst -o filename : none COMPILER: Intel(R) Fortran 12.1-2100