51 |
|
|
52 |
if ( runofffile .NE. ' ' ) then |
if ( runofffile .NE. ' ' ) then |
53 |
#ifdef USE_EXF_INTERPOLATION |
#ifdef USE_EXF_INTERPOLATION |
54 |
interp_method = 2 |
interp_method = 1 |
55 |
call exf_interp( |
call exf_interp( |
56 |
& runofffile, exf_iprec |
& runofffile, exf_iprec |
57 |
& , runoff, 1, xC, yC |
& , runoff, 1, xC, yC |
74 |
do bi = mybxlo(mythid), mybxhi(mythid) |
do bi = mybxlo(mythid), mybxhi(mythid) |
75 |
do j = 1, sny |
do j = 1, sny |
76 |
do i = 1, snx |
do i = 1, snx |
77 |
if (ABS(runoff(i,j,bi,bj)).gt.20.) runoff(i,j,bi,bj)=0. |
if(ABS(runoff(i,j,bi,bj)).gt.998.) runoff(i,j,bi,bj)=0. |
78 |
runoff(i,j,bi,bj)=exf_inscal_runoff*runoff(i,j,bi,bj) |
runoff(i,j,bi,bj)=exf_inscal_runoff*runoff(i,j,bi,bj) |
79 |
enddo |
enddo |
80 |
enddo |
enddo |
84 |
|
|
85 |
_EXCH_XY_R4( runoff, mythid ) |
_EXCH_XY_R4( runoff, mythid ) |
86 |
|
|
|
print *, 'WARNING: runoff raw file cut off at values > 20!' |
|
87 |
CALL PLOT_FIELD_XYRL( runoff, |
CALL PLOT_FIELD_XYRL( runoff, |
88 |
& 'Runoff after rescaling and exf_inscal_runoff' , |
& 'Runoff after rescaling and exf_inscal_runoff' , |
89 |
& 0, myThid ) |
& 0, myThid ) |