/[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.24 by heimbach, Thu May 25 18:32:55 2006 UTC revision 1.25 by mlosch, Thu Jun 29 15:41:58 2006 UTC
# Line 56  c     == end of interface == Line 56  c     == end of interface ==
56        imax = snx+olx        imax = snx+olx
57    
58        do bj = jtlo,jthi        do bj = jtlo,jthi
59          do bi = itlo,ithi         do bi = itlo,ithi
60    
61            do j = jmin,jmax          do j = jmin,jmax
62              do i = imin,imax           do i = imin,imax
63  c  c
64  c             Heat flux.  c     Heat flux.
65                if ( ( hflux(i,j,bi,bj) .GT. 1600. .OR.            if ( ( hflux(i,j,bi,bj) .GT. 1600. .OR.
66       &               hflux(i,j,bi,bj) .LT. -500. ) .AND.       &         hflux(i,j,bi,bj) .LT. -500. ) .AND.
67       &              maskC(i,j,1,bi,bj) .NE. 0. ) then       &         maskC(i,j,1,bi,bj) .NE. 0. ) then
68                   write(standardmessageunit,'(A,3(1X,I6),2X,D22.15)')             write(standardmessageunit,'(A,5(1X,I6),2X,D22.15)')
69       &                'EXF WARNING: hflux out of range for i,j,it= ',       &          'EXF WARNING: hflux out of range for bi,bj,i,j,it= ',
70       &                i, j, myiter, hflux(i,j,bi,bj)       &          bi, bj, i, j, myiter, hflux(i,j,bi,bj)
71                   exferr = 1             exferr = 1
72                endif            endif
73  c  c
74  c             Salt flux.  c     Salt flux.
75                if ( ABS(sflux(i,j,bi,bj)) .GT. 1.E-6 .AND.            if ( ABS(sflux(i,j,bi,bj)) .GT. 1.E-6 .AND.
76       &             maskC(i,j,1,bi,bj) .NE. 0. ) then       &         maskC(i,j,1,bi,bj) .NE. 0. ) then
77                   write(standardmessageunit,'(A,3(1X,I6),2X,D22.15)')             write(standardmessageunit,'(A,5(1X,I6),2X,D22.15)')
78       &                'EXF WARNING: sflux out of range for i,j,it= ',       &          'EXF WARNING: sflux out of range for bi,bj,i,j,it= ',
79       &                i, j, myiter, sflux(i,j,bi,bj)       &          bi, bj, i, j, myiter, sflux(i,j,bi,bj)
80                   exferr = 1             exferr = 1
81                endif            endif
82  c  c
83  c             Zonal wind stress.  c     Zonal wind stress.
84                if ( ABS(ustress(i,j,bi,bj)) .GT. 2.7 .AND.            if ( ABS(ustress(i,j,bi,bj)) .GT. 2.7 .AND.
85       &             maskW(i,j,1,bi,bj) .NE. 0. ) then       &         maskW(i,j,1,bi,bj) .NE. 0. ) then
86                   write(standardmessageunit,'(A,3(1X,I6),2X,D22.15)')             write(standardmessageunit,'(A,5(1X,I6),2X,D22.15)')
87       &                'EXF WARNING: ustress out of range for i,j,it= ',       &          'EXF WARNING: ustress out of range for bi,bj,i,j,it= ',
88       &                i, j, myiter, ustress(i,j,bi,bj)       &          bi, bj, i, j, myiter, ustress(i,j,bi,bj)
89                   exferr = 1             exferr = 1
90                endif            endif
91  c  c
92  c             Meridional wind stress.  c     Meridional wind stress.
93                if ( ABS(vstress(i,j,bi,bj)) .GT. 2.3 .AND.            if ( ABS(vstress(i,j,bi,bj)) .GT. 2.3 .AND.
94       &             maskS(i,j,1,bi,bj) .NE. 0. ) then       &         maskS(i,j,1,bi,bj) .NE. 0. ) then
95                   write(standardmessageunit,'(A,3(1X,I6),2X,D22.15)')             write(standardmessageunit,'(A,5(1X,I6),2X,D22.15)')
96       &                'EXF WARNING: vstress out of range for i,j,it= ',       &          'EXF WARNING: vstress out of range for bi,bj,i,j,it= ',
97       &                i, j, myiter, vstress(i,j,bi,bj)       &          bi, bj, i, j, myiter, vstress(i,j,bi,bj)
98                   exferr = 1             exferr = 1
99                endif            endif
100  c  c
101  #ifdef ALLOW_ATM_WIND  #ifdef ALLOW_ATM_WIND
102  c             zonal wind speed  c     zonal wind speed
103                if ( ABS(uwind(i,j,bi,bj)) .GT. 100. .AND.            if ( ABS(uwind(i,j,bi,bj)) .GT. 100. .AND.
104       &             maskW(i,j,1,bi,bj) .NE. 0. ) then       &         maskW(i,j,1,bi,bj) .NE. 0. ) then
105                   write(standardmessageunit,'(A,3(1X,I6),2X,D22.15)')             write(standardmessageunit,'(A,5(1X,I6),2X,D22.15)')
106       &                'EXF WARNING: uwind out of range for i,j,it= ',       &          'EXF WARNING: uwind out of range for bi,bj,i,j,it= ',
107       &                i, j, myiter, uwind(i,j,bi,bj)       &          bi, bj, i, j, myiter, uwind(i,j,bi,bj)
108                   exferr = 1             exferr = 1
109                endif            endif
110  c  c
111  c             zonal wind speed  c     zonal wind speed
112                if ( ABS(vwind(i,j,bi,bj)) .GT. 100. .AND.            if ( ABS(vwind(i,j,bi,bj)) .GT. 100. .AND.
113       &             maskS(i,j,1,bi,bj) .NE. 0. ) then       &         maskS(i,j,1,bi,bj) .NE. 0. ) then
114                   write(standardmessageunit,'(A,3(1X,I6),2X,D22.15)')             write(standardmessageunit,'(A,5(1X,I6),2X,D22.15)')
115       &                'EXF WARNING: vwind out of range for i,j,it= ',       &          'EXF WARNING: vwind out of range for bi,bj,i,j,it= ',
116       &                i, j, myiter, vwind(i,j,bi,bj)       &          bi, bj, i, j, myiter, vwind(i,j,bi,bj)
117                   exferr = 1             exferr = 1
118                endif            endif
119  #endif  #endif
120  c  c
121  c             wind speed modulus  c     wind speed modulus
122                if ( ( wspeed(i,j,bi,bj) .LT. 0. .OR.            if ( ( wspeed(i,j,bi,bj) .LT. 0. .OR.
123       &              wspeed(i,j,bi,bj) .GT. 100. ) .AND.       &         wspeed(i,j,bi,bj) .GT. 100. ) .AND.
124       &              maskS(i,j,1,bi,bj) .NE. 0. ) then       &         maskS(i,j,1,bi,bj) .NE. 0. ) then
125                   write(standardmessageunit,'(A,3(1X,I6),2X,D22.15)')             write(standardmessageunit,'(A,5(1X,I6),2X,D22.15)')
126       &              'EXF WARNING: wspeed out of range for i,j,it= ',       &          'EXF WARNING: wspeed out of range for bi,bj,i,j,it= ',
127       &              i, j, myiter, wspeed(i,j,bi,bj)       &          bi, bj, i, j, myiter, wspeed(i,j,bi,bj)
128                   exferr = 1             exferr = 1
129                endif            endif
130              
131  #ifdef ALLOW_ATM_TEMP  #ifdef ALLOW_ATM_TEMP
132  c             2-m air temperature  c     2-m air temperature
133                if ( (atemp(i,j,bi,bj) .LT. 183 .OR.            if ( (atemp(i,j,bi,bj) .LT. 183 .OR.
134       &             atemp(i,j,bi,bj) .GT. 343 ) .AND.       &         atemp(i,j,bi,bj) .GT. 343 ) .AND.
135       &             maskC(i,j,1,bi,bj) .NE. 0. ) then       &         maskC(i,j,1,bi,bj) .NE. 0. ) then
136                   write(standardmessageunit,'(2A,3(1X,I6),2X,D22.15)')             write(standardmessageunit,'(2A,5(1X,I6),2X,D22.15)')
137       &                'EXF WARNING: atemp + exf_offset_atemp ',       &          'EXF WARNING: atemp + exf_offset_atemp ',
138       &                'out of range for i,j,it= ',       &          'out of range for bi,bj,i,j,it= ',
139       &                i, j, myiter, atemp(i,j,bi,bj)       &          bi, bj, i, j, myiter, atemp(i,j,bi,bj)
140                   exferr = 1             exferr = 1
141                endif            endif
142  c  c
143  c             2-m specific humidity  c     2-m specific humidity
144                if ( (aqh(i,j,bi,bj) .LT. 0. .OR.            if ( (aqh(i,j,bi,bj) .LT. 0. .OR.
145       &             aqh(i,j,bi,bj) .GT. 0.1 ) .AND.       &         aqh(i,j,bi,bj) .GT. 0.1 ) .AND.
146       &             maskC(i,j,1,bi,bj) .NE. 0. ) then       &         maskC(i,j,1,bi,bj) .NE. 0. ) then
147                   write(standardmessageunit,'(A,3(1X,I6),2X,D22.15)')             write(standardmessageunit,'(A,5(1X,I6),2X,D22.15)')
148       &                'EXF WARNING: aqh out of range for i,j,it= ',       &          'EXF WARNING: aqh out of range for bi,bj,i,j,it= ',
149       &                i, j, myiter, aqh(i,j,bi,bj)       &          bi, bj, i, j, myiter, aqh(i,j,bi,bj)
150                   exferr = 1             exferr = 1
151                endif            endif
152  c  c    
153  c             precipitation rate  c     precipitation rate
154                if ( (precip(i,j,bi,bj) .LT. 0. .OR.            if ( (precip(i,j,bi,bj) .LT. 0. .OR.
155       &             precip(i,j,bi,bj) .GT. 2.E-6 ) .AND.       &         precip(i,j,bi,bj) .GT. 2.E-6 ) .AND.
156       &             maskC(i,j,1,bi,bj) .NE. 0. ) then       &         maskC(i,j,1,bi,bj) .NE. 0. ) then
157                   write(standardmessageunit,'(A,3(1X,I6),2X,D22.15)')             write(standardmessageunit,'(A,5(1X,I6),2X,D22.15)')
158       &                'EXF WARNING: precip out of range for i,j,it= ',       &          'EXF WARNING: precip out of range for bi,bj,i,j,it= ',
159       &                i, j, myiter, precip(i,j,bi,bj)       &          bi, bj, i, j, myiter, precip(i,j,bi,bj)
160                   exferr = 1             exferr = 1
161                endif            endif
162  c  c    
163  c             snow  c     snow
164                if ( (snowprecip(i,j,bi,bj) .LT. 0. .OR.            if ( (snowprecip(i,j,bi,bj) .LT. 0. .OR.
165       &             snowprecip(i,j,bi,bj) .GT. 2.E-6 ) .AND.       &         snowprecip(i,j,bi,bj) .GT. 2.E-6 ) .AND.
166       &             maskC(i,j,1,bi,bj) .NE. 0. ) then       &         maskC(i,j,1,bi,bj) .NE. 0. ) then
167                   write(standardmessageunit,'(A,3(1X,I6),2X,D22.15)')             write(standardmessageunit,'(2A,5(1X,I6),2X,D22.15)')
168       &             'EXF WARNING: snowprecip out of range for i,j,it= ',       &          'EXF WARNING: snowprecip out of range ',
169       &             i, j, myiter, snowprecip(i,j,bi,bj)       &          'for bi,bj,i,j,it= ',
170                   exferr = 1       &          bi, bj, i, j, myiter, snowprecip(i,j,bi,bj)
171                endif             exferr = 1
172              endif
173  #endif  #endif
174    
175  #ifdef SHORTWAVE_HEATING  #ifdef SHORTWAVE_HEATING
176  c             Short wave radiative flux.  c     Short wave radiative flux.
177                if ( (swflux(i,j,bi,bj) .GT. 1. .OR.            if ( (swflux(i,j,bi,bj) .GT. 1. .OR.
178       &             swflux(i,j,bi,bj) .LT. -1000. ) .AND.       &         swflux(i,j,bi,bj) .LT. -1000. ) .AND.
179       &             maskC(i,j,1,bi,bj) .NE. 0. ) then       &         maskC(i,j,1,bi,bj) .NE. 0. ) then
180                   write(standardmessageunit,'(A,3(1X,I6),2X,D22.15)')             write(standardmessageunit,'(A,5(1X,I6),2X,D22.15)')
181       &                'EXF WARNING: swflux out of range for i,j,it= ',       &          'EXF WARNING: swflux out of range for bi,bj,i,j,it= ',
182       &                i, j, myiter, swflux(i,j,bi,bj)       &          bi, bj, i, j, myiter, swflux(i,j,bi,bj)
183                   exferr = 1             exferr = 1
184                endif            endif
185  #endif  #endif
186    
187  #ifdef ALLOW_RUNOFF  #ifdef ALLOW_RUNOFF
188  c             Runoff.  c     Runoff.
189                if ( (runoff(i,j,bi,bj) .LT. 0. .OR.            if ( (runoff(i,j,bi,bj) .LT. 0. .OR.
190       &             runoff(i,j,bi,bj) .GT. 1.E-6 ) .AND.       &         runoff(i,j,bi,bj) .GT. 1.E-6 ) .AND.
191       &             maskC(i,j,1,bi,bj) .NE. 0. ) then       &         maskC(i,j,1,bi,bj) .NE. 0. ) then
192                   write(standardmessageunit,'(A,3(1X,I6),2X,D22.15)')             write(standardmessageunit,'(A,5(1X,I6),2X,D22.15)')
193       &                'EXF WARNING: runoff out of range for i,j,it= ',       &          'EXF WARNING: runoff out of range for bi,bj,i,j,it= ',
194       &                i, j, myiter, runoff(i,j,bi,bj)       &          bi, bj, i, j, myiter, runoff(i,j,bi,bj)
195                   write(standardmessageunit,'(A)')             write(standardmessageunit,'(A)')
196       &                'Please note that input units for runoff are'       &          'Please note that input units for runoff are'
197                   write(standardmessageunit,'(A)')             write(standardmessageunit,'(A)')
198       &                'm/s not m/yr.  If input file is in m/yr, set'       &          'm/s not m/yr.  If input file is in m/yr, set'
199                   write(standardmessageunit,'(A)')             write(standardmessageunit,'(A)')
200       &                'exf_inscal_runoff=3.170979198E-8'       &          'exf_inscal_runoff=3.170979198E-8'
201                   write(standardmessageunit,'(A)')             write(standardmessageunit,'(A)')
202       &                'in the data.exf input file.'       &          'in the data.exf input file.'
203                   exferr = 1             exferr = 1
204                endif            endif
205  #endif  #endif
206             enddo
               if ( exferr .NE. 0 ) then  
                write(standardmessageunit,'(A)')  
      &                'EXF WARNING: If you think these values are OK '  
                write(standardmessageunit,'(A)')  
      &                'EXF WARNING: then set useExfCheckRange=.FALSE.'  
                STOP 'in S/R exf_check_range'  
               endif  
   
             enddo  
           enddo  
 c  
207          enddo          enddo
208    c    
209           enddo
210        enddo        enddo
211    
212          if ( exferr .NE. 0 ) then
213           write(standardmessageunit,'(A)')
214         &      'EXF WARNING: If you think these values are OK '
215           write(standardmessageunit,'(A)')
216         &      'EXF WARNING: then set useExfCheckRange=.FALSE.'
217           STOP 'in S/R exf_check_range'
218          endif
219    
220        end        end
Legend:
Removed from v.1.24  
changed lines
  Added in v.1.25
  ViewVC Help
Powered by ViewVC 1.1.22