62 |
do i = imin,imax |
do i = imin,imax |
63 |
c |
c |
64 |
c Heat flux. |
c Heat flux. |
65 |
if ( ABS(hflux(i,j,bi,bj)) .GT. 1000. .AND. |
if ( ( hflux(i,j,bi,bj) .GT. 1600. .OR. |
66 |
|
& hflux(i,j,bi,bj) .LT. -500. ) .AND. |
67 |
& hFacC(i,j,1,bi,bj) .NE. 0. ) then |
& hFacC(i,j,1,bi,bj) .NE. 0. ) then |
68 |
print *, 'EXF WARNING: hflux out of range for i,j= ', |
print *, 'EXF WARNING: hflux out of range for i,j= ', |
69 |
& i, j, hflux(i,j,bi,bj) |
& i, j, hflux(i,j,bi,bj) |
134 |
c |
c |
135 |
c precipitation rate |
c precipitation rate |
136 |
if ( (precip(i,j,bi,bj) .LT. 0. .OR. |
if ( (precip(i,j,bi,bj) .LT. 0. .OR. |
137 |
& precip(i,j,bi,bj) .GT. 0.1 ) .AND. |
& precip(i,j,bi,bj) .GT. 2.E-6 ) .AND. |
138 |
& hFacC(i,j,1,bi,bj) .NE. 0. ) then |
& hFacC(i,j,1,bi,bj) .NE. 0. ) then |
139 |
print *, 'EXF WARNING: precip out of range for i,j= ', |
print *, 'EXF WARNING: precip out of range for i,j= ', |
140 |
& i, j, precip(i,j,bi,bj) |
& i, j, precip(i,j,bi,bj) |
145 |
#ifdef SHORTWAVE_HEATING |
#ifdef SHORTWAVE_HEATING |
146 |
c Short wave radiative flux. |
c Short wave radiative flux. |
147 |
if ( (swflux(i,j,bi,bj) .GT. 1. .OR. |
if ( (swflux(i,j,bi,bj) .GT. 1. .OR. |
148 |
& swflux(i,j,bi,bj) .LT. -500. ) .AND. |
& swflux(i,j,bi,bj) .LT. -800. ) .AND. |
149 |
& hFacC(i,j,1,bi,bj) .NE. 0. ) then |
& hFacC(i,j,1,bi,bj) .NE. 0. ) then |
150 |
print *, 'EXF WARNING: swflux out of range for i,j= ', |
print *, 'EXF WARNING: swflux out of range for i,j= ', |
151 |
& i, j, swflux(i,j,bi,bj) |
& i, j, swflux(i,j,bi,bj) |
155 |
|
|
156 |
#ifdef ALLOW_RUNOFF |
#ifdef ALLOW_RUNOFF |
157 |
c Runoff. |
c Runoff. |
158 |
if ( (runoff(i,j,bi,bj) .LT. -5.E-8 .OR. |
if ( (runoff(i,j,bi,bj) .LT. 0. .OR. |
159 |
& runoff(i,j,bi,bj) .GT. 1.E-6 ) .AND. |
& runoff(i,j,bi,bj) .GT. 1.E-6 ) .AND. |
160 |
& hFacC(i,j,1,bi,bj) .NE. 0. ) then |
& hFacC(i,j,1,bi,bj) .NE. 0. ) then |
161 |
print *, 'EXF WARNING: runoff out of range for i,j= ', |
print *, 'EXF WARNING: runoff out of range for i,j= ', |