Skip to content

Commit

Permalink
bug fix for the lon-lat boxes of topographical scaling to use 1.2, 3,…
Browse files Browse the repository at this point in the history 
… 6, 9 degree boxes, respectively
  • Loading branch information
@huitang-earth
huitang-earth committed Oct 28, 2022
1 parent 1f48970 commit 516857a
Showing 1 changed file with 8 additions and 8 deletions.
16 changes: 8 additions & 8 deletions src/FLEXDUST.f90
View file
Original file line number Diff line number Diff line change
Expand Up @@ -354,36 +354,36 @@ program FLEXDUST
erodibility_total(ix,iy)=0

!get lower left corner (ix_ll, iy_ll) of erosion area
call getGridPointWind(lat_out-2., lon_out-2,dummy_int, ix_ll, iy_ll, dummy_int, dummy_int)
call getGridPointWind(lat_out-0.6, lon_out-0.6,dummy_int, ix_ll, iy_ll, dummy_int, dummy_int)
!get upper right corner (ix_ur, iy_ur) of erosion area
call getGridPointWind(lat_out+2., lon_out+2.,dummy_int, ix_ur, iy_ur,dummy_int,dummy_int)
call getGridPointWind(lat_out+0.6, lon_out+0.6,dummy_int, ix_ur, iy_ur,dummy_int,dummy_int)
!scale erodibility in this area
call getErodibility(erodibility(ix,iy), ix_wind(ix), iy_wind(iy), ix_ll, ix_ur, iy_ll, iy_ur)

erodibility_total(ix,iy)=erodibility_total(ix,iy)+erodibility(ix,iy)

!get lower left corner (ix_ll, iy_ll) of erosion area
call getGridPointWind(lat_out-5., lon_out-5,dummy_int, ix_ll, iy_ll, dummy_int, dummy_int)
call getGridPointWind(lat_out-1.5, lon_out-1.5,dummy_int, ix_ll, iy_ll, dummy_int, dummy_int)
!get upper right corner (ix_ur, iy_ur) of erosion area
call getGridPointWind(lat_out+5., lon_out+5.,dummy_int, ix_ur, iy_ur,dummy_int,dummy_int)
call getGridPointWind(lat_out+1.5, lon_out+1.5,dummy_int, ix_ur, iy_ur,dummy_int,dummy_int)
!scale erodibility in this area
call getErodibility(erodibility(ix,iy), ix_wind(ix), iy_wind(iy), ix_ll, ix_ur, iy_ll, iy_ur)

erodibility_total(ix,iy)=erodibility_total(ix,iy)+erodibility(ix,iy)

!get lower left corner (ix_ll, iy_ll) of erosion area
call getGridPointWind(lat_out-10., lon_out-10,dummy_int, ix_ll, iy_ll, dummy_int, dummy_int)
call getGridPointWind(lat_out-3, lon_out-3,dummy_int, ix_ll, iy_ll, dummy_int, dummy_int)
!get upper right corner (ix_ur, iy_ur) of erosion area
call getGridPointWind(lat_out+10., lon_out+10.,dummy_int, ix_ur, iy_ur,dummy_int,dummy_int)
call getGridPointWind(lat_out+3, lon_out+3,dummy_int, ix_ur, iy_ur,dummy_int,dummy_int)
!scale erodibility in this area
call getErodibility(erodibility(ix,iy), ix_wind(ix), iy_wind(iy), ix_ll, ix_ur, iy_ll, iy_ur)

erodibility_total(ix,iy)=erodibility_total(ix,iy)+erodibility(ix,iy)

!get lower left corner (ix_ll, iy_ll) of erosion area
call getGridPointWind(lat_out-15., lon_out-15,dummy_int, ix_ll, iy_ll, dummy_int, dummy_int)
call getGridPointWind(lat_out-4.5, lon_out-4.5,dummy_int, ix_ll, iy_ll, dummy_int, dummy_int)
!get upper right corner (ix_ur, iy_ur) of erosion area
call getGridPointWind(lat_out+15., lon_out+15.,dummy_int, ix_ur, iy_ur,dummy_int,dummy_int)
call getGridPointWind(lat_out+4.5, lon_out+4.5,dummy_int, ix_ur, iy_ur,dummy_int,dummy_int)
!scale erodibility in this area
call getErodibility(erodibility(ix,iy), ix_wind(ix), iy_wind(iy), ix_ll, ix_ur, iy_ll, iy_ur)

Expand Down

0 comments on commit 516857a

Please sign in to comment.