From 16cbaf50ce4be4994dbf7cee496ab440e50050d7 Mon Sep 17 00:00:00 2001 From: Bill Ladwig Date: Wed, 10 Jan 2018 12:16:54 -0700 Subject: [PATCH] Fixed a segfault issue when writing the LCL and LFC values --- fortran/rip_cape.f90 | 6 +++--- 1 file changed, 3 insertions(+), 3 deletions(-) diff --git a/fortran/rip_cape.f90 b/fortran/rip_cape.f90 index be05b79..c92be3f 100644 --- a/fortran/rip_cape.f90 +++ b/fortran/rip_cape.f90 @@ -710,7 +710,7 @@ SUBROUTINE DCAPECALC2D(prs,tmk,qvp,ght,ter,sfp,cape,cin,& !$OMP facden, tmkenv, qvpenv, eslift, qvplift, buoy, benamin, & !$OMP benaccum, zrel, kmax, dz, elfound, & !$OMP kel, klfc, pavg, p2, p1, totthe, totqvp, totprs, & - !$OMP i, j, k, kpar, kpar1, kpar2, qvppari, tmkpari,p, pup, pdn, th, & + !$OMP i, j, k, kpar, kpar1, kpar2, qvppari, tmkpari, p, pup, pdn, th, & !$OMP pp1, pp2, ethmax, eth_temp, klev) SCHEDULE(runtime) DO j = 1,mjy DO i = 1,mix @@ -964,9 +964,9 @@ SUBROUTINE DCAPECALC2D(prs,tmk,qvp,ght,ter,sfp,cape,cin,& cape(i,j,mkzh) = cape(i,j,kpar1) cin(i,j,mkzh) = cin(i,j,kpar1) ! meters agl - cin(i,j,mkzh-1) = zrel(klcl) + ght_new(kpar,i,j) - ter(i,j) + cin(i,j,mkzh-1) = zrel(klcl) + ght_new(kpar1,i,j) - ter(i,j) ! meters agl - cin(i,j,mkzh-2) = zrel(klfc) + ght_new(kpar,i,j) - ter(i,j) + cin(i,j,mkzh-2) = zrel(klfc) + ght_new(kpar1,i,j) - ter(i,j) END DO END DO