diff --git a/fortran/calc_uh.f90 b/fortran/calc_uh.f90 index 09eafb6..72acfe6 100644 --- a/fortran/calc_uh.f90 +++ b/fortran/calc_uh.f90 @@ -47,7 +47,7 @@ SUBROUTINE DCALCUH(nx, ny, nz, nzp1, zp, mapfct, dx, dy, uhmnhgt, uhmxhgt, us, & ! Misc local variables INTEGER :: i, j, k, kbot, ktop - REAL(KIND=8) :: twodx, twody, wgtlw, sum, wmean, wsum, wavg + REAL(KIND=8) :: twodx, twody, wgtlw, sum, wmean, wsum !,wavg REAL(KIND=8) :: helbot, heltop, wbot, wtop REAL(KIND=8) :: zbot, ztop diff --git a/fortran/omp.f90 b/fortran/omp.f90 index 81805ff..bea1974 100644 --- a/fortran/omp.f90 +++ b/fortran/omp.f90 @@ -17,6 +17,7 @@ SUBROUTINE fomp_set_num_threads(num_threads) !f2py threadsafe INTEGER, INTENT(IN) :: num_threads + IF (.FALSE.) PRINT *, num_threads END SUBROUTINE fomp_set_num_threads @@ -93,6 +94,7 @@ SUBROUTINE fomp_set_dynamic(dynamic_threads) !f2py threadsafe LOGICAL, INTENT(IN) :: dynamic_threads + IF (.FALSE.) PRINT *, dynamic_threads END SUBROUTINE fomp_set_dynamic @@ -117,6 +119,7 @@ SUBROUTINE fomp_set_nested(nested) !f2py threadsafe LOGICAL, INTENT(IN) :: nested + IF (.FALSE.) PRINT *, nested END SUBROUTINE fomp_set_nested @@ -144,6 +147,7 @@ SUBROUTINE fomp_set_schedule(kind, modifier) INTEGER(KIND=fomp_sched_kind), INTENT(IN) :: kind INTEGER, INTENT(IN) :: modifier + IF (.FALSE.) PRINT *, kind, modifier END SUBROUTINE fomp_set_schedule @@ -185,6 +189,7 @@ SUBROUTINE fomp_set_max_active_levels(max_levels) !f2py threadsafe INTEGER, INTENT(IN) :: max_levels + IF (.FALSE.) PRINT *, max_levels END SUBROUTINE fomp_set_max_active_levels @@ -223,6 +228,7 @@ FUNCTION fomp_get_ancestor_thread_num(level) INTEGER, INTENT(IN) :: level INTEGER :: fomp_get_ancestor_thread_num + IF (.FALSE.) PRINT *, level fomp_get_ancestor_thread_num = -1 @@ -237,6 +243,7 @@ FUNCTION fomp_get_team_size(level) INTEGER, INTENT(IN) :: level INTEGER :: fomp_get_team_size + IF (.FALSE.) PRINT *, level fomp_get_team_size = -1 @@ -308,6 +315,7 @@ SUBROUTINE fomp_destroy_lock(svar) !f2py threadsafe INTEGER(KIND=fomp_lock_kind), INTENT(INOUT) :: svar + IF (.FALSE.) PRINT *, svar END SUBROUTINE fomp_destroy_lock @@ -321,6 +329,7 @@ SUBROUTINE fomp_destroy_nest_lock(nvar) !f2py threadsafe INTEGER(KIND=fomp_nest_lock_kind), INTENT(INOUT) :: nvar + IF (.FALSE.) PRINT *, nvar END SUBROUTINE fomp_destroy_nest_lock @@ -334,6 +343,7 @@ SUBROUTINE fomp_set_lock(svar) !f2py threadsafe INTEGER(KIND=fomp_lock_kind), INTENT(INOUT) :: svar + IF (.FALSE.) PRINT *, svar END SUBROUTINE fomp_set_lock @@ -347,6 +357,7 @@ SUBROUTINE fomp_set_nest_lock(nvar) !f2py threadsafe INTEGER(KIND=fomp_nest_lock_kind), INTENT(INOUT) :: nvar + IF (.FALSE.) PRINT *, nvar END SUBROUTINE fomp_set_nest_lock @@ -360,6 +371,7 @@ SUBROUTINE fomp_unset_lock(svar) !f2py threadsafe INTEGER(KIND=fomp_lock_kind), INTENT(INOUT) :: svar + IF (.FALSE.) PRINT *, svar END SUBROUTINE fomp_unset_lock @@ -373,6 +385,7 @@ SUBROUTINE fomp_unset_nest_lock(nvar) !f2py threadsafe INTEGER(KIND=fomp_nest_lock_kind), INTENT(INOUT) :: nvar + IF (.FALSE.) PRINT *, nvar END SUBROUTINE fomp_unset_nest_lock @@ -387,6 +400,7 @@ FUNCTION fomp_test_lock(svar) INTEGER(KIND=fomp_lock_kind), INTENT(INOUT) :: svar LOGICAL :: fomp_test_lock + IF (.FALSE.) PRINT *, svar fomp_test_lock = .FALSE. @@ -403,6 +417,7 @@ FUNCTION fomp_test_nest_lock(nvar) INTEGER(KIND=fomp_nest_lock_kind), INTENT(INOUT) :: nvar INTEGER :: fomp_test_nest_lock + IF (.FALSE.) PRINT *, nvar fomp_test_nest_lock = -1 diff --git a/fortran/ompgen.F90 b/fortran/ompgen.F90 index ed3d9bd..dea4246 100644 --- a/fortran/ompgen.F90 +++ b/fortran/ompgen.F90 @@ -74,8 +74,11 @@ SUBROUTINE fomp_set_num_threads(num_threads) #ifdef _OPENMP CALL omp_set_num_threads(num_threads) +#else + IF (.FALSE.) PRINT *, num_threads #endif + END SUBROUTINE fomp_set_num_threads @@ -192,6 +195,8 @@ SUBROUTINE fomp_set_dynamic(dynamic_threads) #ifdef _OPENMP CALL omp_set_dynamic(dynamic_threads) +#else + IF (.FALSE.) PRINT *, dynamic_threads #endif END SUBROUTINE fomp_set_dynamic @@ -230,6 +235,8 @@ SUBROUTINE fomp_set_nested(nested) #ifdef _OPENMP CALL omp_set_nested(nested) +#else + IF (.FALSE.) PRINT *, nested #endif END SUBROUTINE fomp_set_nested @@ -270,6 +277,8 @@ SUBROUTINE fomp_set_schedule(kind, modifier) #ifdef _OPENMP CALL omp_set_schedule(kind, modifier) +#else + IF (.FALSE.) PRINT *, kind, modifier #endif END SUBROUTINE fomp_set_schedule @@ -327,14 +336,12 @@ SUBROUTINE fomp_set_max_active_levels(max_levels) !f2py threadsafe -#ifdef _OPENMP - INTEGER, INTENT(IN) :: max_levels -#else INTEGER, INTENT(IN) :: max_levels -#endif #ifdef _OPENMP CALL omp_set_max_active_levels(max_levels) +#else + IF (.FALSE.) PRINT *, max_levels #endif END SUBROUTINE fomp_set_max_active_levels @@ -395,6 +402,7 @@ FUNCTION fomp_get_ancestor_thread_num(level) #ifdef _OPENMP fomp_get_ancestor_thread_num = omp_get_ancestor_thread_num(level) #else + IF (.FALSE.) PRINT *, level fomp_get_ancestor_thread_num = -1 #endif @@ -416,6 +424,7 @@ FUNCTION fomp_get_team_size(level) #ifdef _OPENMP fomp_get_team_size = omp_get_team_size(level) #else + IF (.FALSE.) PRINT *, level fomp_get_team_size = -1 #endif @@ -518,8 +527,11 @@ SUBROUTINE fomp_destroy_lock(svar) #ifdef _OPENMP CALL omp_destroy_lock(svar) +#else + IF (.FALSE.) PRINT *, svar #endif + END SUBROUTINE fomp_destroy_lock @@ -537,6 +549,8 @@ SUBROUTINE fomp_destroy_nest_lock(nvar) #ifdef _OPENMP CALL omp_destroy_nest_lock(nvar) +#else + IF (.FALSE.) PRINT *, nvar #endif END SUBROUTINE fomp_destroy_nest_lock @@ -556,6 +570,8 @@ SUBROUTINE fomp_set_lock(svar) #ifdef _OPENMP CALL omp_set_lock(svar) +#else + IF (.FALSE.) PRINT *, svar #endif END SUBROUTINE fomp_set_lock @@ -575,6 +591,8 @@ SUBROUTINE fomp_set_nest_lock(nvar) #ifdef _OPENMP CALL omp_set_nest_lock(nvar) +#else + IF (.FALSE.) PRINT *, nvar #endif END SUBROUTINE fomp_set_nest_lock @@ -594,6 +612,8 @@ SUBROUTINE fomp_unset_lock(svar) #ifdef _OPENMP CALL omp_unset_lock(svar) +#else + IF (.FALSE.) PRINT *, svar #endif END SUBROUTINE fomp_unset_lock @@ -613,6 +633,8 @@ SUBROUTINE fomp_unset_nest_lock(nvar) #ifdef _OPENMP CALL omp_unset_nest_lock(nvar) +#else + IF (.FALSE.) PRINT *, nvar #endif END SUBROUTINE fomp_unset_nest_lock @@ -634,6 +656,7 @@ FUNCTION fomp_test_lock(svar) #ifdef _OPENMP fomp_test_lock = omp_test_lock(svar) #else + IF (.FALSE.) PRINT *, svar fomp_test_lock = .FALSE. #endif @@ -658,6 +681,7 @@ FUNCTION fomp_test_nest_lock(nvar) #ifdef _OPENMP fomp_test_nest_lock = omp_test_nest_lock(nvar) #else + IF (.FALSE.) PRINT *, nvar fomp_test_nest_lock = -1 #endif diff --git a/fortran/rip_cape.f90 b/fortran/rip_cape.f90 index dfa6f20..dc5b79b 100644 --- a/fortran/rip_cape.f90 +++ b/fortran/rip_cape.f90 @@ -274,7 +274,7 @@ END SUBROUTINE DPFCALC SUBROUTINE DCAPECALC3D(prs,tmk,qvp,ght,ter,sfp,cape,cin,& cmsg,mix,mjy,mkzh,ter_follow,& psafile, errstat, errmsg) - USE wrf_constants, ONLY : ALGERR, CELKEL, G, EZERO, ESLCON1, ESLCON2, & + USE wrf_constants, ONLY : CELKEL, G, EZERO, ESLCON1, ESLCON2, & EPS, RD, CP, GAMMA, CPMD, RGASMD, GAMMAMD, TLCLC1, & TLCLC2, TLCLC3, TLCLC4, THTECON1, THTECON2, THTECON3 @@ -312,7 +312,6 @@ SUBROUTINE DCAPECALC3D(prs,tmk,qvp,ght,ter,sfp,cape,cin,& REAL(KIND=8), DIMENSION(150) :: psadithte, psadiprs REAL(KIND=8), DIMENSION(150,150) :: psaditmk LOGICAL :: elfound - REAL :: t1,t2 REAL(KIND=8), DIMENSION(mkzh,mix,mjy) :: prs_new REAL(KIND=8), DIMENSION(mkzh,mix,mjy) :: tmk_new @@ -325,6 +324,7 @@ SUBROUTINE DCAPECALC3D(prs,tmk,qvp,ght,ter,sfp,cape,cin,& klev = 0 klcl = 0 kel = 0 + IF (.FALSE.) PRINT *,ter ! the comments were taken from a mark stoelinga email, 23 apr 2007, @@ -599,7 +599,7 @@ END SUBROUTINE DCAPECALC3D SUBROUTINE DCAPECALC2D(prs,tmk,qvp,ght,ter,sfp,cape,cin,& cmsg,mix,mjy,mkzh,ter_follow,& psafile, errstat, errmsg) - USE wrf_constants, ONLY : ALGERR, CELKEL, G, EZERO, ESLCON1, ESLCON2, & + USE wrf_constants, ONLY : CELKEL, G, EZERO, ESLCON1, ESLCON2, & EPS, RD, CP, GAMMA, CPMD, RGASMD, GAMMAMD, TLCLC1, & TLCLC2, TLCLC3, TLCLC4, THTECON1, THTECON2, THTECON3 @@ -628,10 +628,10 @@ SUBROUTINE DCAPECALC2D(prs,tmk,qvp,ght,ter,sfp,cape,cin,& ! local variables INTEGER :: i, j, k, ilcl, kel, kk, klcl, klev, klfc, kmax, kpar, kpar1, kpar2 - REAL(KIND=8) :: ethmax, q, p, e, tlcl, zlcl + REAL(KIND=8) :: ethmax, p, e, tlcl, zlcl REAL(KIND=8) :: pavg, tvirtual, p1, p2, pp1, pp2, th, totthe, totqvp, totprs REAL(KIND=8) :: cpm, deltap, ethpari, gammam, qvppari, tmkpari - REAL(KIND=8) :: facden, qvplift, tmklift, tvenv, tvlift, ghtlift, fac1, fac2 + REAL(KIND=8) :: facden, qvplift, tmklift, tvenv, tvlift, ghtlift REAL(KIND=8) :: eslift, tmkenv, qvpenv, tonpsadiabat REAL(KIND=8) :: benamin, dz, pup, pdn REAL(KIND=8), DIMENSION(150) :: buoy, zrel, benaccum @@ -639,7 +639,6 @@ SUBROUTINE DCAPECALC2D(prs,tmk,qvp,ght,ter,sfp,cape,cin,& REAL(KIND=8), DIMENSION(150) :: psadithte, psadiprs REAL(KIND=8), DIMENSION(150,150) :: psaditmk LOGICAL :: elfound - INTEGER :: nthreads REAL(KIND=8), DIMENSION(mkzh) :: eth_temp REAL(KIND=8), DIMENSION(mkzh,mix,mjy) :: prs_new REAL(KIND=8), DIMENSION(mkzh,mix,mjy) :: tmk_new diff --git a/fortran/wrf_cloud_fracf.f90 b/fortran/wrf_cloud_fracf.f90 index 9930206..92ae745 100644 --- a/fortran/wrf_cloud_fracf.f90 +++ b/fortran/wrf_cloud_fracf.f90 @@ -78,7 +78,7 @@ SUBROUTINE DCLOUDFRAC2(vert, rh, vert_inc_w_height, low_thresh, mid_thresh, & ! NCLEND - INTEGER i, j, k, kstart, kend + INTEGER i, j, k INTEGER kchi, kcmi, kclo ! Initialize the output diff --git a/fortran/wrf_user.f90 b/fortran/wrf_user.f90 index 1e133a1..1bfc82c 100644 --- a/fortran/wrf_user.f90 +++ b/fortran/wrf_user.f90 @@ -327,11 +327,12 @@ SUBROUTINE DCOMPUTESEAPRS(nx, ny, nz, z, t, p, q, sea_level_pressure, & INTEGER :: i, j, k INTEGER :: klo, khi - INTEGER :: errcnt, bad_i, bad_j, bad_sfp + INTEGER :: errcnt, bad_i, bad_j + REAL(KIND=8) :: bad_sfp REAL(KIND=8) :: plo, phi, tlo, thi, zlo, zhi REAL(KIND=8) :: p_at_pconst, t_at_pconst, z_at_pconst - REAL(KIND=8) :: z_half_lowest + !REAL(KIND=8) :: z_half_lowest LOGICAL :: l1, l2, l3, found diff --git a/fortran/wrf_vinterp.f90 b/fortran/wrf_vinterp.f90 index 2dc7c17..c34efd7 100644 --- a/fortran/wrf_vinterp.f90 +++ b/fortran/wrf_vinterp.f90 @@ -169,7 +169,6 @@ SUBROUTINE wrf_vintrp(datain, dataout, pres, tk, qvp, ght, terrain,& REAL(KIND=8) :: pbot, zbot, tbotextrap, e REAL(KIND=8) :: tlcl, gammam CHARACTER(LEN=1) :: cvcord - INTEGER :: thd ! Removes the warnings for uninitialized variables cvcord = ''