- Timestamp:
- Jul 23, 2014, 1:07:57 PM (10 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/LMDZ_WRFmeas/WRFV3/lmdz/diagphy_mod.F90
r35 r36 79 79 !C 80 80 integer pas 81 82 save pas 83 data pas/0/ 84 !$OMP THREADPRIVATE(pas) 85 !C 86 81 87 ! L. Fita, LMD July 2014 82 88 CHARACTER(LEN=50) :: errmsg, fname, varname … … 88 94 largest = 10000. 89 95 90 save pas91 data pas/0/92 !$OMP THREADPRIVATE(pas)93 !C94 96 pas=pas+1 95 97 stops=0. … … 219 221 varname = 'evap' 220 222 CALL check_var(fname, varn, evap, klon, largest, .FALSE.) 223 END IF 221 224 PRINT *,' rain_fall_tot= ',rain_fall_tot 222 225 IF (rain_fall_tot .NE. rain_fall_tot .OR. ABS(rain_fall_tot) > largest) THEN 223 226 varname = 'rain_fall' 224 227 CALL check_var(fname, varn, rain_fall, klon, largest, .FALSE.) 228 END IF 225 229 PRINT *,' snow_fall_tot= ',snow_fall_tot 226 230 IF (snow_fall_tot .NE. snow_fall_tot .OR. ABS(snow_fall_tot) > largest) THEN 227 231 varname = 'snow_fall' 228 232 CALL check_var(fname, varn, snow_fall, klon, largest, .FALSE.) 233 END IF 229 234 STOP 230 235 END IF … … 241 246 SUBROUTINE check_var(funcn, varn, var, sizev, bigvalue, stoprun) 242 247 ! Subroutine to check the consistency of a variable 248 ! * NaN value: by definition is variable /= variable 249 ! * bigvalue: allowd threshold for variable 243 250 244 251 IMPLICIT NONE … … 270 277 wrongi = 0 271 278 DO i=1,sizev 272 IF (var(i) /= var(i) .OR. var(i) > bigvalule ) THEN279 IF (var(i) /= var(i) .OR. ABS(var(i)) > bigvalule ) THEN 273 280 IF (wrongi == 0) found = .TRUE. 274 281 wrongi = wrongi + 1
Note: See TracChangeset
for help on using the changeset viewer.