diff --git a/src/HISTORY b/src/HISTORY
index dd1438d..5b5a92e 100644
--- a/src/HISTORY
+++ b/src/HISTORY
@@ -588,4 +588,5 @@ DEBUG about the problem on "MULTI=ON", "init_disp_fault=0" and "multrupt=1" is a
 (contributed by naotaka-chikasada. Thanks!).
 
 --------------------------------------------------------------------------------
-
+2020.04.02
+- Upload again to get a DOI badge
diff --git a/src/mod_fxy.f90 b/src/mod_fxy.f90
index 9e9a7eb..d6d2c31 100644
--- a/src/mod_fxy.f90
+++ b/src/mod_fxy.f90
@@ -721,7 +721,7 @@ subroutine fxynl_rwg(wfld,dfld,ffld,ifz,cfs,cfl,cflag,dt,th0,dth,joff,nlon,nlat,
                         if(hz(i+1,j) + ddx(i,j) > hz(i,j) + ddx(i,j)) discharge = -discharge
                         fx(i,j) = discharge
                         if(brokenx(i,j) == 0) then
-                           btx(i,j) = broken_rate*(btx(i,j) - dhigh) + dhigh
+                           btx(i,j) = broken_rate*(btx(i,j) + dhigh) - dhigh
                            brokenx(i,j) = 1
                         end if
                      end if
@@ -784,7 +784,7 @@ subroutine fxynl_rwg(wfld,dfld,ffld,ifz,cfs,cfl,cflag,dt,th0,dth,joff,nlon,nlat,
                         if(hz(i,j+1) + ddy(i,j) > hz(i,j) + ddy(i,j)) discharge = -discharge
                         fy(i,j) = discharge
                         if(brokeny(i,j) == 0) then
-                           bty(i,j) = broken_rate*(bty(i,j) - dhigh) + dhigh
+                           bty(i,j) = broken_rate*(bty(i,j) + dhigh) - dhigh
                            brokeny(i,j) = 1
                         end if
                      end if
diff --git a/src/mod_fxy_cartesian.f90 b/src/mod_fxy_cartesian.f90
index a9e6b5e..cd97837 100644
--- a/src/mod_fxy_cartesian.f90
+++ b/src/mod_fxy_cartesian.f90
@@ -479,7 +479,7 @@ subroutine fxynl_rwg(wfld,dfld,ffld,ifz,cfs,cfl,dt,dxdy,nlon,nlat,gflag,smallh,b
                         if(hz(i+1,j) + ddx(i,j) > hz(i,j) + ddx(i,j)) discharge = -discharge
                         fx(i,j) = discharge
                         if(brokenx(i,j) == 0) then
-                           btx(i,j) = broken_rate*(btx(i,j) - dhigh) + dhigh
+                           btx(i,j) = broken_rate*(btx(i,j) + dhigh) - dhigh
                            brokenx(i,j) = 1
                         end if
                      end if
@@ -718,7 +718,7 @@ subroutine fxynl_rwg(wfld,dfld,ffld,ifz,cfs,cfl,dt,dxdy,nlon,nlat,gflag,smallh,b
                         if(hz(i,j+1) + ddy(i,j) > hz(i,j) + ddy(i,j)) discharge = -discharge
                         fy(i,j) = discharge
                         if(brokeny(i,j) == 0) then
-                           bty(i,j) = broken_rate*(bty(i,j) - dhigh) + dhigh
+                           bty(i,j) = broken_rate*(bty(i,j) + dhigh) - dhigh
                            brokeny(i,j) = 1
                         end if
                      end if
diff --git a/src/mod_fxy_disp.f90 b/src/mod_fxy_disp.f90
index 4a569f6..2bcd3b2 100644
--- a/src/mod_fxy_disp.f90
+++ b/src/mod_fxy_disp.f90
@@ -983,7 +983,7 @@ subroutine fxynl_rwg_disp(wfld,dfld,ffld,ifz,cfs,cfl,cflag,dt,th0,dth,joff,nlon,
                         if(hz(i+1,j) + ddx(i,j) > hz(i,j) + ddx(i,j)) discharge = -discharge
                         fx(i,j) = discharge
                         if(brokenx(i,j) == 0) then
-                           btx(i,j) = broken_rate*(btx(i,j) - dhigh) + dhigh
+                           btx(i,j) = broken_rate*(btx(i,j) + dhigh) - dhigh
                            brokenx(i,j) = 1
                         end if
                      end if
@@ -1235,7 +1235,7 @@ subroutine fxynl_rwg_disp(wfld,dfld,ffld,ifz,cfs,cfl,cflag,dt,th0,dth,joff,nlon,
                         if(hz(i,j+1) + ddy(i,j) > hz(i,j) + ddy(i,j)) discharge = -discharge
                         fy(i,j) = discharge
                         if(brokeny(i,j) == 0) then
-                           bty(i,j) = broken_rate*(bty(i,j) - dhigh) + dhigh
+                           bty(i,j) = broken_rate*(bty(i,j) + dhigh) - dhigh
                            brokeny(i,j) = 1
                         end if
                      end if
diff --git a/src/mod_fxy_disp_cartesian.f90 b/src/mod_fxy_disp_cartesian.f90
index cdbc052..a3ed9f7 100644
--- a/src/mod_fxy_disp_cartesian.f90
+++ b/src/mod_fxy_disp_cartesian.f90
@@ -900,7 +900,7 @@ subroutine fxynl_rwg_disp(wfld,dfld,ffld,ifz,cfs,cfl,dt,dxdy,nlon,nlat, &
                         if(hz(i+1,j) + ddx(i,j) > hz(i,j) + ddx(i,j)) discharge = -discharge
                         fx(i,j) = discharge
                         if(brokenx(i,j) == 0) then
-                           btx(i,j) = broken_rate*(btx(i,j) - dhigh) + dhigh
+                           btx(i,j) = broken_rate*(btx(i,j) + dhigh) - dhigh
                            brokenx(i,j) = 1
                         end if
                      end if
@@ -1139,7 +1139,7 @@ subroutine fxynl_rwg_disp(wfld,dfld,ffld,ifz,cfs,cfl,dt,dxdy,nlon,nlat, &
                         if(hz(i,j+1) + ddy(i,j) > hz(i,j) + ddy(i,j)) discharge = -discharge
                         fy(i,j) = discharge
                         if(brokeny(i,j) == 0) then
-                           bty(i,j) = broken_rate*(bty(i,j) - dhigh) + dhigh
+                           bty(i,j) = broken_rate*(bty(i,j) + dhigh) - dhigh
                            brokeny(i,j) = 1
                         end if
                      end if