GF deep convection divide by zero

[dudhia]

Will just post email thread here to make sure we have a record in Issues. On Thu, Oct 1, 2020 at 9:33 AM Pedro Jimenez Munoz [email protected] wrote: the line is highlighted in bold below. We identified this issue before and we have a thread opened with Georg but we never completed it. I will forward the thread.

        do ki=jmin(i)  ,1,-1
         dzo=zo_cup(i,ki+1)-zo_cup(i,ki)
         h_entr=.5*(heo(i,ki)+.5*(hco(i,ki)+hco(i,ki+1)))
         **ucd(i,ki)=(ucd(i,ki+1)*zdo(i,ki+1)                                   &
                     -.5*dd_massdetru(i,ki)*ucd(i,ki+1)+                      &
                    dd_massentru(i,ki)*us(i,ki)                               &
                    -pgcon*zdo(i,ki+1)*(us(i,ki+1)-us(i,ki)))   /             &
                    (zdo(i,ki+1)-.5*dd_massdetru(i,ki)+dd_massentru(i,ki))**
         vcd(i,ki)=(vcd(i,ki+1)*zdo(i,ki+1)                                   &
                     -.5*dd_massdetru(i,ki)*vcd(i,ki+1)+                      &
                    dd_massentru(i,ki)*vs(i,ki)                               &
                    -pgcon*zdo(i,ki+1)*(vs(i,ki+1)-vs(i,ki)))   /             &
                    (zdo(i,ki+1)-.5*dd_massdetru(i,ki)+dd_massentru(i,ki))
         hcdo(i,ki)=(hcdo(i,ki+1)*zdo(i,ki+1)                                 &
                     -.5*dd_massdetro(i,ki)*hcdo(i,ki+1)+                     &
                    dd_massentro(i,ki)*h_entr)   /                            &
                    (zdo(i,ki+1)-.5*dd_massdetro(i,ki)+dd_massentro(i,ki))
         dbydo(i,ki)=hcdo(i,ki)-heso_cup(i,ki)
         bud(i)=bud(i)+dbydo(i,ki)*dzo
        enddo

On Thu, Oct 1, 2020 at 9:24 AM Jimy Dudhia [email protected] wrote: Pedro, What line is this? Jimy

module_cu_gf_deep.f90:1007

---

[pedro-jm]

@dudhia This is the line: ucd(i,ki)=(ucd(i,ki+1)zdo(i,ki+1) & **-.5dd_massdetru(i,ki)ucd(i,ki+1)+ &* dd_massentru(i,ki)us(i,ki) & -pgconzdo(i,ki+1)(us(i,ki+1)-us(i,ki))) / & (zdo(i,ki+1)-.5dd_massdetru(i,ki)+dd_massentru(i,ki))

[pedro-jm]

@dudhia It seems github instead of writing in bold, marks the lines with **

[dudhia]

Yes, I tried editing the Issue text and bold does that.

On Thu, Nov 5, 2020 at 10:44 AM pedro-jm [email protected] wrote:

@dudhia https://github.com/dudhia It seems github instead of writing in bold, marks the lines with **

— You are receiving this because you were mentioned. Reply to this email directly, view it on GitHub https://github.com/wrf-model/WRF/issues/1318#issuecomment-722532781, or unsubscribe https://github.com/notifications/unsubscribe-auth/AEIZ77DRFE7MUCXD4GWJBQLSOLQA5ANCNFSM4TLU5WAQ .