From ad00aa904f8eada15aca3b493a7edeafa9a03c25 Mon Sep 17 00:00:00 2001 From: Nick Golding Date: Fri, 18 Aug 2023 12:50:47 +0800 Subject: [PATCH] patch rainydays divide by 0 error in micro_global() --- R/micro_global.R | 12 ++++++++++-- 1 file changed, 10 insertions(+), 2 deletions(-) diff --git a/R/micro_global.R b/R/micro_global.R index d5e119f5..d3ef0319 100644 --- a/R/micro_global.R +++ b/R/micro_global.R @@ -912,11 +912,19 @@ micro_global <- function( if(k==1){ RAINFALL1[m]<-RAINFALL[i]*rainfrac*rainmult # if first day of month, make user-specified fraction of monthly rainfall fall on first day }else{ - RAINFALL1[m]<-(RAINFALL[i]*(1-rainfrac)*rainmult)/RAINYDAYS[i] # make remaining rain fall evenly over the remaining number of rainy days for the month, starting at the beginning of the month + if(RAINYDAYS[i]>0){ + RAINFALL1[m]<-(RAINFALL[i]*(1-rainfrac)*rainmult)/RAINYDAYS[i] # make remaining rain fall evenly over the remaining number of rainy days for the month, starting at the beginning of the month + }else{ + RAINFALL1[m]<-0 + } } }else{ if(rainfrac==0){ - RAINFALL1[m]<-(RAINFALL[i]*rainmult)/RAINYDAYS[i] + if(RAINYDAYS[i]>0){ + RAINFALL1[m]<-(RAINFALL[i]*rainmult)/RAINYDAYS[i] + }else{ + RAINFALL1[m]<-0 + } }else{ RAINFALL1[m]<-0 }