From 5b50e8d5cf38898fc6aac9c382960f5a1a3dd1eb Mon Sep 17 00:00:00 2001 From: jagurs-admin Date: Thu, 2 Apr 2020 18:13:10 +0900 Subject: [PATCH] release of JAGURS-D0515.a Reupload to get a DOI budge --- src/HISTORY | 3 ++- src/mod_fxy.f90 | 4 ++-- src/mod_fxy_cartesian.f90 | 4 ++-- src/mod_fxy_disp.f90 | 4 ++-- src/mod_fxy_disp_cartesian.f90 | 4 ++-- 5 files changed, 10 insertions(+), 9 deletions(-) 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