Page 1 Source Listing READPCP 2025-03-12 18:21 /tmp/ifortiXcTLZ.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 2025-03-12 18:21 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 2025-03-12 18:21 Entry Points READPCP.F ENTRY POINTS Name readpcp_ SYMBOL CROSS REFERENCE Name Object Declared Type Bytes Dimen Elements Attributes References 10 Label 142 139 BAOPENR Subr 129 129 BITMPCP Local 121 L(1) 1 2 91719 137,147,155,186 DSTRB Subr 196 196 GETGB Subr 137 137 GLB_TABLE Common 99 128 SAVE I Local 152 I(4) 4 scalar 152,155,157,158,159,160,161,162,16 4,179,184,185,186 IDIM1 Param 34 I(4) 4 scalar 40,41,43,44 IDIM2 Param 34 I(4) 4 scalar 40,41,43,44 IGSTL Param 29 I(4) 4 scalar 34 IGSTR Param 29 I(4) 4 scalar 34 IHR Local 134 I(4) 4 scalar 134,136,137,139,196 IM Param 18 I(4) 4 scalar 31,34,105,106,107,108,109,113,114, 119,121,137,152,184 INPES Param 26 I(4) 4 scalar 31,34,89,90,91,100,101 IRET Local 129 I(4) 4 scalar 129,130,138,139,145 ITAIL Param 31 I(4) 4 scalar J Local 151 I(4) 4 scalar 151,155,157,158,159,160,161,162,16 4,179,184,185,186 JDIM1 Param 35 I(4) 4 scalar 40,41,43,44 JDIM2 Param 35 I(4) 4 scalar 40,41,43,44 JGDS Local 122 I(4) 4 1 22 137 JGSTL Param 30 I(4) 4 scalar 35 JGSTR Param 30 I(4) 4 scalar 35 JM Param 18 I(4) 4 scalar 32,35,105,106,107,108,109,113,114, 119,121,137,151,184 JNPES Param 26 I(4) 4 scalar 32,35,89,90,91,100,101 JPDS Local 122 I(4) 4 1 25 131,137 JTAIL Param 32 I(4) 4 scalar K Local 137 I(4) 4 scalar 137 KF Local 137 I(4) 4 scalar 137,139 KGDS Local 122 I(4) 4 1 22 137 KPDS Local 122 I(4) 4 1 25 137,139,140,141 LM Param 18 I(4) 4 scalar 40,114 LSM Param 18 I(4) 4 scalar MAPPINGS Common 118 5024 SAVE MOD Func 140 scalar 140,184 MPPCOM Common 63 1464 SAVE PPTASM Common 39 3878292 SAVE READPCP Subr 2 TEMPCOM Common 104 6603768 SAVE TOPO Common 112 17478548 SAVE Page 4 Source Listing READPCP 2025-03-12 18:21 Symbol Table READPCP.F TYPE COMPONENTS/COMMON VARIABLES Name Type Bytes Offset Dimen Elements Attributes References APREC R(4) 4 3519180 2 17955 COM CLDINC R(4) 4 3734652 2 17955 COM CLDINC7 R(4) 4 3806472 2 17955 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 620 1 8 COM IE_GLB_TABLE I(4) 4 32 1 8 COM IE_LOC_TABLE I(4) 4 556 1 8 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 1064 1 100 COM IQUILT_GROUP I(4) 4 1060 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 8 COM IS_LOC_TABLE I(4) 4 492 1 8 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 3591000 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 96 1 8 COM JE_LOC_TABLE I(4) 4 588 1 8 COM JS_GLB_TABLE I(4) 4 64 1 8 COM JS_LOC_TABLE I(4) 4 524 1 8 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 3591004 scalar COM L2GI I(4) 4 956 1 239 COM L2GJ I(4) 4 3468 1 389 COM MPI_COMM_COMP I(4) 4 652 scalar COM MPI_COMM_INTER I(4) 4 656 scalar COM MPI_COMM_INTER_ARRAY I(4) 4 660 1 100 COM MTSTPE I(4) 4 3591008 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:21 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:21 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 128,135 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 17955 COM PPTDAT R(4) 4 3303720 3 53865 COM 196 TEMP1 R(4) 4 0 2 91719 COM 138,146,157,158,159,160,161,162,16 4,179,185,196 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:21 Symbol Table READPCP.F Name Type Bytes Offset Dimen Elements Attributes References TLAT R(4) 4 71820 3 807975 COM TTVG R(4) 4 371884 2 92971 COM VAPINC R(4) 4 3591012 2 17955 COM VAPINC7 R(4) 4 3662832 2 17955 COM Page 8 Source Listing READPCP 2025-03-12 18:21 READPCP.F 202 Page 9 Source Listing READPCP 2025-03-12 18:21 Subprograms/Common Blocks READPCP.F SUBPROGRAMS/COMMON BLOCKS Name Object Declared Type Bytes Dimen Elements Attributes References GLB_TABLE Common 99 128 SAVE MAPPINGS Common 118 5024 SAVE MPPCOM Common 63 1464 SAVE PPTASM Common 39 3878292 SAVE READPCP Subr 2 TEMPCOM Common 104 6603768 SAVE TOPO Common 112 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:21 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:21 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