/[MITgcm]/MITgcm/pkg/exf/exf_check_range.F
ViewVC logotype

Diff of /MITgcm/pkg/exf/exf_check_range.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.1 by heimbach, Thu Feb 26 19:34:34 2004 UTC revision 1.2 by heimbach, Thu Feb 26 22:30:27 2004 UTC
# Line 58  c     == end of interface == Line 58  c     == end of interface ==
58              do i = imin,imax              do i = imin,imax
59  c  c
60  c             Heat flux.  c             Heat flux.
61                if ( ABS(hflux(i,j,bi,bj)) .GT. 1000. )                if ( ABS(hflux(i,j,bi,bj)) .GT. 1000. .AND.
62         &           hFacC(i,j,1,bi,bj) .NE. 0. )
63       &           print *, 'EXF WARNING: hflux out of range for i,j= ',       &           print *, 'EXF WARNING: hflux out of range for i,j= ',
64       &           hflux(i,j,bi,bj)       &           i, j, hflux(i,j,bi,bj)
65  c  c
66  c             Salt flux.  c             Salt flux.
67                if ( ABS(sflux(i,j,bi,bj)) .GT. 1.E-6 )                if ( ABS(sflux(i,j,bi,bj)) .GT. 1.E-6 .AND.
68         &           hFacC(i,j,1,bi,bj) .NE. 0. )
69       &           print *, 'EXF WARNING: sflux out of range for i,j= ',       &           print *, 'EXF WARNING: sflux out of range for i,j= ',
70       &           sflux(i,j,bi,bj)       &           i, j, sflux(i,j,bi,bj)
71  c  c
72  c             Zonal wind stress.  c             Zonal wind stress.
73                if ( ABS(ustress(i,j,bi,bj)) .GT. 2. )                if ( ABS(ustress(i,j,bi,bj)) .GT. 2. .AND.
74         &           hFacW(i,j,1,bi,bj) .NE. 0. )
75       &           print *, 'EXF WARNING: ustress out of range for i,j= ',       &           print *, 'EXF WARNING: ustress out of range for i,j= ',
76       &           ustress(i,j,bi,bj)       &           i, j, ustress(i,j,bi,bj)
77  c  c
78  c             Meridional wind stress.  c             Meridional wind stress.
79                if ( ABS(vstress(i,j,bi,bj)) .GT. 2. )                if ( ABS(vstress(i,j,bi,bj)) .GT. 2. .AND.
80         &           hFacS(i,j,1,bi,bj) .NE. 0. )
81       &           print *, 'EXF WARNING: vstress out of range for i,j= ',       &           print *, 'EXF WARNING: vstress out of range for i,j= ',
82       &           vstress(i,j,bi,bj)       &           i, j, vstress(i,j,bi,bj)
83  c  c
84  #ifdef ALLOW_ATM_WIND  #ifdef ALLOW_ATM_WIND
85  c             zonal wind speed  c             zonal wind speed
86                if ( ABS(uwind(i,j,bi,bj)) .GT. 40. )                if ( ABS(uwind(i,j,bi,bj)) .GT. 40. .AND.
87         &           hFacW(i,j,1,bi,bj) .NE. 0. )
88       &           print *, 'EXF WARNING: uwind out of range for i,j= ',       &           print *, 'EXF WARNING: uwind out of range for i,j= ',
89       &           uwind(i,j,bi,bj)       &           i, j, uwind(i,j,bi,bj)
90  c  c
91  c             zonal wind speed  c             zonal wind speed
92                if ( ABS(vwind(i,j,bi,bj)) .GT. 40. )                if ( ABS(vwind(i,j,bi,bj)) .GT. 40. .AND.
93         &           hFacS(i,j,1,bi,bj) .NE. 0. )
94       &           print *, 'EXF WARNING: vwind out of range for i,j= ',       &           print *, 'EXF WARNING: vwind out of range for i,j= ',
95       &           vwind(i,j,bi,bj)       &           i, j, vwind(i,j,bi,bj)
96  #endif  #endif
97    
98  #ifdef ALLOW_ATM_TEMP  #ifdef ALLOW_ATM_TEMP
99  c             2-m air temperature  c             2-m air temperature
100                if ( atemp(i,j,bi,bj) .LT. 223 .OR.                if ( (atemp(i,j,bi,bj) .LT. 223 .OR.
101       &             atemp(i,j,bi,bj) .GT. 323 )       &              atemp(i,j,bi,bj) .GT. 323 ) .AND.
102         &           hFacC(i,j,1,bi,bj) .NE. 0. )
103       &           print *, 'EXF WARNING: atemp out of range for i,j= ',       &           print *, 'EXF WARNING: atemp out of range for i,j= ',
104       &           atemp(i,j,bi,bj)       &           i, j, atemp(i,j,bi,bj)
105  c  c
106  c             2-m specific humidity  c             2-m specific humidity
107                if ( aqh(i,j,bi,bj) .LT. 0. .OR.                if ( (aqh(i,j,bi,bj) .LT. 0. .OR.
108       &             aqh(i,j,bi,bj) .GT. 1.E-6 )       &              aqh(i,j,bi,bj) .GT. 1.E-6 ) .AND.
109         &           hFacC(i,j,1,bi,bj) .NE. 0. )
110       &           print *, 'EXF WARNING: aqh out of range for i,j= ',       &           print *, 'EXF WARNING: aqh out of range for i,j= ',
111       &           aqh(i,j,bi,bj)       &           i, j, aqh(i,j,bi,bj)
112  c  c
113  c             precipitation rate  c             precipitation rate
114                if ( precip(i,j,bi,bj) .LT. 0. .OR.                if ( (precip(i,j,bi,bj) .LT. 0. .OR.
115       &             precip(i,j,bi,bj) .GT. 0.1 )       &              precip(i,j,bi,bj) .GT. 0.1 ) .AND.
116         &           hFacC(i,j,1,bi,bj) .NE. 0. )
117       &           print *, 'EXF WARNING: precip out of range for i,j= ',       &           print *, 'EXF WARNING: precip out of range for i,j= ',
118       &           precip(i,j,bi,bj)       &           i, j, precip(i,j,bi,bj)
119  #endif  #endif
120    
121  #ifdef SHORTWAVE_HEATING  #ifdef SHORTWAVE_HEATING
122  c             Short wave radiative flux.  c             Short wave radiative flux.
123                if ( swflux(i,j,bi,bj) .GT. 0. .OR.                if ( (swflux(i,j,bi,bj) .GT. 0. .OR.
124       &             swflux(i,j,bi,bj) .LT. -500. )       &              swflux(i,j,bi,bj) .LT. -500. ) .AND.
125         &           hFacC(i,j,1,bi,bj) .NE. 0. )
126       &           print *, 'EXF WARNING: swflux out of range for i,j= ',       &           print *, 'EXF WARNING: swflux out of range for i,j= ',
127       &           swflux(i,j,bi,bj)       &           i, j, swflux(i,j,bi,bj)
128  #endif  #endif
129    
130              enddo              enddo
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22