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) |
87 |
endif |
endif |
88 |
c |
c |
89 |
c Meridional wind stress. |
c Meridional wind stress. |
90 |
if ( ABS(vstress(i,j,bi,bj)) .GT. 2. .AND. |
if ( ABS(vstress(i,j,bi,bj)) .GT. 2.2 .AND. |
91 |
& hFacS(i,j,1,bi,bj) .NE. 0. ) then |
& hFacS(i,j,1,bi,bj) .NE. 0. ) then |
92 |
print *, 'EXF WARNING: vstress out of range for i,j= ', |
print *, 'EXF WARNING: vstress out of range for i,j= ', |
93 |
& i, j, vstress(i,j,bi,bj) |
& i, j, vstress(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. -1000. ) .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) |