Page 1 Source Listing READPCP 2025-03-12 18:22 /tmp/ifort92HPoY.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----------------------------------------------------------------------- 39 INCLUDE "PPTASM.comm" 63 INCLUDE "mpp.h" 122 LOGICAL*1 BITMPCP(IM,JM) 123 INTEGER JPDS(25), JGDS(22), KGDS(22), KPDS(25) 124 125 C----------------------------------------------------------------------- 126 C 127 C READ ANALYSES AND DATA MASK 128 C 129 IF(MYPE.EQ.0) THEN 130 call baopenr(41,'fort.41',iret) 131 write(0,*) 'baopenr on unit41', ' iret=', iret 132 JPDS=-1 133 ENDIF 134 C 135 DO IHR=1,3 136 IF(MYPE.EQ.0) THEN 137 print*,'ihr=',ihr 138 CALL GETGB(41,0,IM*JM,IHR-1,JPDS,JGDS,KF,K,KPDS,KGDS,BITMPCP, 139 & temp1,IRET) 140 write(0,10) ihr, iret, kf, kpds(5), 141 & (kpds(21)*100+kpds(8))/100-1, mod(kpds(8),100),kpds(9), 142 & kpds(10), kpds(11) 143 10 format('ihr=',i1,' iret=',i3,' kf=', i6,' fld=', i3, 144 & 2x,'1h accum from ', 5i2.2) 145 C 1/2012: WNE if bad read 146 if (iret.ne.0) then 147 temp1 = 999.0 148 BITMPCP = .true. 149 endif 150 151 C 152 DO J=1,JM 153 DO I=1,IM 154 c if(temp1(i,j).gt.1.0) print*,'i,j,pcp=',i,j,temp1(i,j) 155 c print*,'i,j,pcp=',i,j,temp1(i,j) 156 IF (BITMPCP(I,J)) THEN 157 C Convert data from mm to m: 158 if (temp1(i,j).lt.0.0) then 159 write(0,*) " precip < 0. will be set to 999.", I,J, temp1(i,j) 160 temp1(i,j) = 999.0 161 else if (temp1(i,j).gt.24.0) then Page 2 Source Listing READPCP 2025-03-12 18:22 READPCP.F 162 write(0,*) " precip >24. will be set to 999.", I,J, temp1(i,j) 163 temp1(i,j) = 999.0 164 else 165 temp1(I,J)=temp1(I,J)*0.001 166 endif 167 c if(temp1(i,j).gt.0.04) then 168 c if(temp1(i,j).gt.0.001) then 169 c if(i.ge. 99.and.i.le.105.and.(j.ge.172.and.j.le.174).or 170 c * .(j.ge.180.and.j.le.184)) then 171 c temp1(i,j)=999. 172 c bitmpcp(i,j)=.f. 173 c endif 174 c if(i.ge. 99.and.i.le.105.and.j.ge.113.and.j.le.119) then 175 c print*,'i,j,pcp=',i,j,temp1(i,j) 176 c temp1(i,j)=999. 177 c bitmpcp(i,j)=.f. 178 c endif 179 ELSE 180 TEMP1(I,J)=999. 181 ENDIF 182 c 183 c Mask out outer 2 rows and outer 2 colums - 19 March 2004 PS and DJ 184 c 185 if(i.le.1.or.i.ge.im-1+mod(j,2) .or. j.le.2.or.j.ge.jm-1) then 186 temp1(i,j)=999. 187 bitmpcp(i,j)=.FALSE. 188 endif 189 ENDDO 190 ENDDO 191 WRITE(0,*) ' FINISHED READING PRECIP ANALYSES AND DATA MASK' 192 ENDIF 193 c 194 c special test 195 c 196 C 197 CALL DSTRB(TEMP1,PPTDAT,1,3,IHR) 198 C 199 ENDDO 200 C 201 RETURN 202 END Page 3 Source Listing READPCP 2025-03-12 18:22 Entry Points READPCP.F ENTRY POINTS Name readpcp_ SYMBOL CROSS REFERENCE Name Object Declared Type Bytes Dimen Elements Attributes References 10 Label 143 140 BAOPENR Subr 130 130 BITMPCP Local 122 L(1) 1 2 91719 138,148,156,187 DSTRB Subr 197 197 GETGB Subr 138 138 GLB_TABLE Common 100 336 SAVE I Local 153 I(4) 4 scalar 153,156,158,159,160,161,162,163,16 5,180,185,186,187 IDIM1 Param 35 I(4) 4 scalar 41,42,44,45 IDIM2 Param 35 I(4) 4 scalar 41,42,44,45 IGSTL Param 30 I(4) 4 scalar 35 IGSTR Param 30 I(4) 4 scalar 35 IHR Local 135 I(4) 4 scalar 135,137,138,140,197 IM Param 18 I(4) 4 scalar 32,35,106,107,108,109,110,114,115, 120,122,138,153,185 INPES Param 27 I(4) 4 scalar 32,35,90,91,92,101,102 IRET Local 130 I(4) 4 scalar 130,131,139,140,146 ITAIL Param 32 I(4) 4 scalar J Local 152 I(4) 4 scalar 152,156,158,159,160,161,162,163,16 5,180,185,186,187 JDIM1 Param 36 I(4) 4 scalar 41,42,44,45 JDIM2 Param 36 I(4) 4 scalar 41,42,44,45 JGDS Local 123 I(4) 4 1 22 138 JGSTL Param 31 I(4) 4 scalar 36 JGSTR Param 31 I(4) 4 scalar 36 JM Param 18 I(4) 4 scalar 33,36,106,107,108,109,110,114,115, 120,122,138,152,185 JNPES Param 27 I(4) 4 scalar 33,36,90,91,92,101,102 JPDS Local 123 I(4) 4 1 25 132,138 JTAIL Param 33 I(4) 4 scalar K Local 138 I(4) 4 scalar 138 KF Local 138 I(4) 4 scalar 138,140 KGDS Local 123 I(4) 4 1 22 138 KPDS Local 123 I(4) 4 1 25 138,140,141,142 LM Param 18 I(4) 4 scalar 41,115 LSM Param 18 I(4) 4 scalar MAPPINGS Common 119 5024 SAVE MOD Func 141 scalar 141,185 MPPCOM Common 64 1724 SAVE PPTASM Common 40 1364916 SAVE READPCP Subr 2 TEMPCOM Common 105 6603768 SAVE TOPO Common 113 17478548 SAVE Page 4 Source Listing READPCP 2025-03-12 18:22 Symbol Table READPCP.F TYPE COMPONENTS/COMMON VARIABLES Name Type Bytes Offset Dimen Elements Attributes References APREC R(4) 4 1238524 2 6319 COM CLDINC R(4) 4 1314364 2 6319 COM CLDINC7 R(4) 4 1339640 2 6319 COM G2LI I(4) 4 0 1 239 COM G2LJ I(4) 4 1912 1 389 COM HTMG R(4) 4 743768 3 4183695 COM IBROW I(4) 4 404 scalar COM ICHUNKTAB I(4) 4 828 1 21 COM IE_GLB_TABLE I(4) 4 84 1 21 COM IE_LOC_TABLE I(4) 4 660 1 21 COM ILCOL I(4) 4 396 scalar COM ILPAD1 I(4) 4 412 scalar COM ILPAD2 I(4) 4 416 scalar COM ILPAD3 I(4) 4 420 scalar COM ILPAD4 I(4) 4 424 scalar COM ILPAD5 I(4) 4 428 scalar COM INUMQ I(4) 4 1324 1 100 COM IQUILT_GROUP I(4) 4 1320 scalar COM IRCOL I(4) 4 400 scalar COM IRPAD1 I(4) 4 432 scalar COM IRPAD2 I(4) 4 436 scalar COM IRPAD3 I(4) 4 440 scalar COM IRPAD4 I(4) 4 444 scalar COM IRPAD5 I(4) 4 448 scalar COM IS_GLB_TABLE I(4) 4 0 1 21 COM IS_LOC_TABLE I(4) 4 492 1 21 COM ITEMP I(4) 4 5870016 2 91719 COM ITEMP2 I(4) 4 6236892 2 91719 COM ITROW I(4) 4 408 scalar COM ITSTLOC I(4) 4 1263800 scalar COM JBPAD1 I(4) 4 452 scalar COM JBPAD2 I(4) 4 456 scalar COM JBPAD3 I(4) 4 460 scalar COM JBPAD4 I(4) 4 464 scalar COM JBPAD5 I(4) 4 468 scalar COM JE_GLB_TABLE I(4) 4 252 1 21 COM JE_LOC_TABLE I(4) 4 744 1 21 COM JS_GLB_TABLE I(4) 4 168 1 21 COM JS_LOC_TABLE I(4) 4 576 1 21 COM JTPAD1 I(4) 4 472 scalar COM JTPAD2 I(4) 4 476 scalar COM JTPAD3 I(4) 4 480 scalar COM JTPAD4 I(4) 4 484 scalar COM JTPAD5 I(4) 4 488 scalar COM JTSTLOC I(4) 4 1263804 scalar COM L2GI I(4) 4 956 1 239 COM L2GJ I(4) 4 3468 1 389 COM MPI_COMM_COMP I(4) 4 912 scalar COM MPI_COMM_INTER I(4) 4 916 scalar COM MPI_COMM_INTER_ARRAY I(4) 4 920 1 100 COM MTSTPE I(4) 4 1263808 scalar COM MYIE I(4) 4 52 scalar COM MYIE1 I(4) 4 56 scalar COM Page 5 Source Listing READPCP 2025-03-12 18:22 Symbol Table READPCP.F Name Type Bytes Offset Dimen Elements Attributes References MYIE1_P1 I(4) 4 120 scalar COM MYIE1_P2 I(4) 4 124 scalar COM MYIE1_P3 I(4) 4 128 scalar COM MYIE1_P4 I(4) 4 132 scalar COM MYIE2 I(4) 4 60 scalar COM MYIE2_P1 I(4) 4 136 scalar COM MYIE_P1 I(4) 4 100 scalar COM MYIE_P2 I(4) 4 104 scalar COM MYIE_P3 I(4) 4 108 scalar COM MYIE_P4 I(4) 4 112 scalar COM MYIE_P5 I(4) 4 116 scalar COM MYIS I(4) 4 40 scalar COM MYIS1 I(4) 4 44 scalar COM MYIS1_P1 I(4) 4 84 scalar COM MYIS1_P2 I(4) 4 88 scalar COM MYIS1_P3 I(4) 4 92 scalar COM MYIS1_P4 I(4) 4 96 scalar COM MYIS2 I(4) 4 48 scalar COM MYIS_P1 I(4) 4 64 scalar COM MYIS_P2 I(4) 4 68 scalar COM MYIS_P3 I(4) 4 72 scalar COM MYIS_P4 I(4) 4 76 scalar COM MYIS_P5 I(4) 4 80 scalar COM MYJE I(4) 4 236 scalar COM MYJE1 I(4) 4 240 scalar COM MYJE1_P1 I(4) 4 280 scalar COM MYJE1_P2 I(4) 4 284 scalar COM MYJE1_P3 I(4) 4 288 scalar COM MYJE1_P4 I(4) 4 292 scalar COM MYJE2 I(4) 4 244 scalar COM MYJE2_P1 I(4) 4 296 scalar COM MYJE2_P2 I(4) 4 300 scalar COM MYJE2_P3 I(4) 4 304 scalar COM MYJE2_P4 I(4) 4 308 scalar COM MYJE3 I(4) 4 248 scalar COM MYJE3_P4 I(4) 4 312 scalar COM MYJE4 I(4) 4 252 scalar COM MYJE4_P1 I(4) 4 316 scalar COM MYJE4_P4 I(4) 4 320 scalar COM MYJE5 I(4) 4 256 scalar COM MYJE5_P1 I(4) 4 324 scalar COM MYJE5_P2 I(4) 4 328 scalar COM MYJE_P1 I(4) 4 260 scalar COM MYJE_P2 I(4) 4 264 scalar COM MYJE_P3 I(4) 4 268 scalar COM MYJE_P4 I(4) 4 272 scalar COM MYJE_P5 I(4) 4 276 scalar COM MYJS I(4) 4 140 scalar COM MYJS1 I(4) 4 144 scalar COM MYJS1_P1 I(4) 4 184 scalar COM MYJS1_P2 I(4) 4 188 scalar COM MYJS1_P3 I(4) 4 192 scalar COM MYJS1_P4 I(4) 4 196 scalar COM MYJS2 I(4) 4 148 scalar COM MYJS2_P1 I(4) 4 200 scalar COM Page 6 Source Listing READPCP 2025-03-12 18:22 Symbol Table READPCP.F Name Type Bytes Offset Dimen Elements Attributes References MYJS2_P2 I(4) 4 204 scalar COM MYJS2_P3 I(4) 4 208 scalar COM MYJS2_P4 I(4) 4 212 scalar COM MYJS3 I(4) 4 152 scalar COM MYJS3_P4 I(4) 4 216 scalar COM MYJS4 I(4) 4 156 scalar COM MYJS4_P1 I(4) 4 220 scalar COM MYJS4_P4 I(4) 4 224 scalar COM MYJS5 I(4) 4 160 scalar COM MYJS5_P1 I(4) 4 228 scalar COM MYJS5_P2 I(4) 4 232 scalar COM MYJS_P1 I(4) 4 164 scalar COM MYJS_P2 I(4) 4 168 scalar COM MYJS_P3 I(4) 4 172 scalar COM MYJS_P4 I(4) 4 176 scalar COM MYJS_P5 I(4) 4 180 scalar COM MYPE I(4) 4 0 scalar COM 129,136 MY_E I(4) 4 336 scalar COM MY_IE_GLB I(4) 4 12 scalar COM MY_IE_LOC I(4) 4 28 scalar COM MY_IS_GLB I(4) 4 8 scalar COM MY_IS_LOC I(4) 4 24 scalar COM MY_JE_GLB I(4) 4 20 scalar COM MY_JE_LOC I(4) 4 36 scalar COM MY_JS_GLB I(4) 4 16 scalar COM MY_JS_LOC I(4) 4 32 scalar COM MY_N I(4) 4 332 scalar COM MY_NE I(4) 4 348 scalar COM MY_NEB I(4) 4 364 1 8 COM MY_NW I(4) 4 360 scalar COM MY_S I(4) 4 340 scalar COM MY_SE I(4) 4 352 scalar COM MY_SW I(4) 4 356 scalar COM MY_W I(4) 4 344 scalar COM NPES I(4) 4 4 scalar COM PHOUR R(4) 4 0 2 6319 COM PPTDAT R(4) 4 1162696 3 18957 COM 197 TEMP1 R(4) 4 0 2 91719 COM 139,147,158,159,160,161,162,163,16 5,180,186,197 TEMP10 R(4) 4 3301884 2 91719 COM TEMP11 R(4) 4 3668760 2 91719 COM TEMP12 R(4) 4 4035636 2 91719 COM TEMP13 R(4) 4 4402512 2 91719 COM TEMP14 R(4) 4 4769388 2 91719 COM TEMP15 R(4) 4 5136264 2 91719 COM TEMP16 R(4) 4 5503140 2 91719 COM TEMP2 R(4) 4 366876 2 91719 COM TEMP2X R(4) 4 0 2 92971 COM TEMP3 R(4) 4 733752 2 91719 COM TEMP4 R(4) 4 1100628 2 91719 COM TEMP5 R(4) 4 1467504 2 91719 COM TEMP6 R(4) 4 1834380 2 91719 COM TEMP7 R(4) 4 2201256 2 91719 COM TEMP8 R(4) 4 2568132 2 91719 COM TEMP9 R(4) 4 2935008 2 91719 COM Page 7 Source Listing READPCP 2025-03-12 18:22 Symbol Table READPCP.F Name Type Bytes Offset Dimen Elements Attributes References TLAT R(4) 4 25276 3 284355 COM TTVG R(4) 4 371884 2 92971 COM VAPINC R(4) 4 1263812 2 6319 COM VAPINC7 R(4) 4 1289088 2 6319 COM Page 8 Source Listing READPCP 2025-03-12 18:22 READPCP.F 203 Page 9 Source Listing READPCP 2025-03-12 18:22 Subprograms/Common Blocks READPCP.F SUBPROGRAMS/COMMON BLOCKS Name Object Declared Type Bytes Dimen Elements Attributes References GLB_TABLE Common 100 336 SAVE MAPPINGS Common 119 5024 SAVE MPPCOM Common 64 1724 SAVE PPTASM Common 40 1364916 SAVE READPCP Subr 2 TEMPCOM Common 105 6603768 SAVE TOPO Common 113 17478548 SAVE COMPILER OPTIONS BEING USED -align noall -align nonone -align nocommons -align nodcommons -align noqcommons -align nozcommons -align records -align nosequence -align norec1byte -align norec2byte -align norec4byte -align norec8byte -align norec16byte -align norec32byte -align norec64byte -align noarray8byte -align noarray16byte -align noarray32byte -align noarray64byte -align noarray128byte -align noarray256byte -altparam -assume accuracy_sensitive -assume nobscc -assume nobuffered_io -assume nobuffered_stdout -assume byterecl -assume nocontiguous_assumed_shape -assume nocontiguous_pointer -assume nocc_omp -assume nocstring -assume nodummy_aliases -assume nofpe_summary -assume noieee_fpe_flags -assume nominus0 -assume noold_boz -assume old_complex_align -assume old_unit_star -assume old_inquire_recl -assume old_ldout_format -assume old_ldout_zero -assume noold_logical_assign -assume noold_logical_ldio -assume old_maxminloc -assume old_xor -assume noprotect_allocates -assume protect_constants -assume noprotect_parens -assume split_common -assume source_include -assume nostd_intent_in -assume std_minus0_rounding -assume nostd_mod_proc_name -assume std_value -assume realloc_lhs -assume underscore -assume no2underscores -assume norecursion no -auto -auto_scalar no -bintext -ccdefault default -check noarg_temp_created -check noassume -check nobounds -check nocontiguous -check noformat -check nooutput_conversion -check nooverflow -check nopointers -check noshape -check nostack -check nouninitialized -check noudio_iostat -coarray-num-procs 0 no -coarray-config-file -convert big_endian -cross_reference -D __INTEL_COMPILER=1910 -D __INTEL_COMPILER_UPDATE=3 -D __unix__ -D __unix Page 10 Source Listing READPCP 2025-03-12 18:22 READPCP.F -D __linux__ -D __linux -D __gnu_linux__ -D unix -D linux -D __ELF__ -D __x86_64 -D __x86_64__ -D __amd64 -D __amd64__ -D __INTEL_COMPILER_BUILD_DATE=20200925 -D __INTEL_OFFLOAD -D __MMX__ -D __SSE__ -D __SSE_MATH__ -D __SSE2__ -D __SSE2_MATH__ -D __SSE3__ -D __SSSE3__ -D __SSE4_1__ -D __SSE4_2__ -D __POPCNT__ -D __PCLMUL__ -D __AES__ -D __AVX__ -D __F16C__ -D __AVX_I__ -D __RDRND__ -D __FMA__ -D __FP_FAST_FMA -D __FP_FAST_FMAF -D __BMI__ -D __LZCNT__ -D __AVX2__ -D __haswell -D __haswell__ -D __tune_haswell__ -D __core_avx2 -D __core_avx2__ -D __tune_core_avx2__ -D __CRAY_X86_ROME -D __CRAYXT_COMPUTE_LINUX_TARGET -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 -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 -init noarrays -init nohuge -init noinfinity -init nominus_huge -init nominus_infinity -init nominus_tiny -init nonan -init nosnan -init notiny -init nozero no -intconstant -integer_size 32 no -mixed_str_len_arg Page 11 Source Listing READPCP 2025-03-12 18:22 READPCP.F no -module -names lowercase no -noinclude no -o -offload-build=host -openmp-simd -O2 no -pad_source -real_size 32 no -recursive -reentrancy threaded -vec=simd -show nofullpath -show noinclude -show map -show options no -syntax_only no -threadcom no -U no -vms -w noall -w nonone -w alignments -w nodeclarations -w noexternals -w general -w noignore_bounds -w noignore_loc -w nointerfaces -w noshape -w notruncated_source -w uncalled -w uninitialized -w nounused -w usage no -wrap-margins -includepath : /pe/intel/compilers_and_libraries_2020.4.304/linux/pstl/include/,/pe/intel/compilers_and_libraries_2020.4.304/linux/compiler/include/, .f90,./.f90,/opt/cray/pe/mpich/8.1.12/ofi/intel/19.0/include/.f90,/pe/intel/compilers_and_libraries_2020.4.304/linux/ipp/include/.f90, /pe/intel/compilers_and_libraries_2020.4.304/linux/mkl/include/.f90,/pe/intel/compilers_and_libraries_2020.4.304/linux/pstl/include/.f90, /pe/intel/compilers_and_libraries_2020.4.304/linux/pstl/stdlib/.f90,/pe/intel/compilers_and_libraries_2020.4.304/linux/tbb/include/.f90, /pe/intel/compilers_and_libraries_2020.4.304/linux/compiler/include/intel64/.f90,/pe/intel/compilers_and_libraries_2020.4.304/linux/compiler/include/icc/.f90, /pe/intel/compilers_and_libraries_2020.4.304/linux/compiler/include/.f90,/usr/lib64/gcc/x86_64-suse-linux/7/include/.f90, /usr/lib64/gcc/x86_64-suse-linux/7/include-fixed/.f90,/usr/include/.f90,/usr/include/.f90,/usr/include/.f90 -list filename : READPCP.lst no -o COMPILER: Intel(R) Fortran 19.1-1655