Changeset 37 in lmdz_wrf
- Timestamp:
- Jul 23, 2014, 1:11:39 PM (10 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/LMDZ_WRFmeas/WRFV3/lmdz/diagphy_mod.F90
r36 r37 158 158 159 159 ! L. Fita, LMD July 2014. Checking for consistency 160 IF (fs_bound .NE. fs_bound .OR. ABS(fs_bound) > lar agest) THEN160 IF (fs_bound .NE. fs_bound .OR. ABS(fs_bound) > largest) THEN 161 161 PRINT *,TRIM(errmsg) 162 162 PRINT *,' ' + TRIM(fname) + ': Wrong fs_bound= ',fs_bound,' !!!' … … 166 166 IF (airetot .NE. airetot .OR. ABS(airetot) > largest) THEN 167 167 varname = 'airephy' 168 CALL check_var(fname, varn , airephy, klon, largest, .FALSE.)168 CALL check_var(fname, varname, airephy, klon, largest, .FALSE.) 169 169 END IF 170 170 PRINT *,' stops= ',stops 171 171 IF (stops .NE. stops .OR. ABS(stops) > largest) THEN 172 172 varname = 'tops' 173 CALL check_var(fname, varn , tops, klon, largest, .FALSE.)173 CALL check_var(fname, varname, tops, klon, largest, .FALSE.) 174 174 END IF 175 175 PRINT *,' stopl= ',stopl 176 176 IF (stopl .NE. stopl .OR. ABS(stopl) > largest) THEN 177 177 varname = 'topl' 178 CALL check_var(fname, varn , topl, klon, largest, .FALSE.)178 CALL check_var(fname, varname, topl, klon, largest, .FALSE.) 179 179 END IF 180 180 PRINT *,' ssols= ',ssols 181 181 IF (ssols .NE. ssols .OR. ABS(ssols) > largest) THEN 182 182 varname = 'sols' 183 CALL check_var(fname, varn , sols, klon, largest, .FALSE.)183 CALL check_var(fname, varname, sols, klon, largest, .FALSE.) 184 184 END IF 185 185 PRINT *,' ssoll= ',ssoll 186 186 IF (ssoll .NE. ssoll .OR. ABS(ssoll) > largest) THEN 187 187 varname = 'soll' 188 CALL check_var(fname, varn , soll, klon, largest, .FALSE.)188 CALL check_var(fname, varname, soll, klon, largest, .FALSE.) 189 189 END IF 190 190 PRINT *,' ssens= ',ssens 191 191 IF (ssens .NE. ssens .OR. ABS(ssens) > largest) THEN 192 192 varname = 'sens' 193 CALL check_var(fname, varn , sens, klon, largest, .FALSE.)193 CALL check_var(fname, varname, sens, klon, largest, .FALSE.) 194 194 END IF 195 195 PRINT *,' sfront= ',sfront 196 196 IF (sfront .NE. sfront .OR. ABS(sfront) > largest) THEN 197 197 varname = 'evap' 198 CALL check_var(fname, varn , evap, klon, largest, .FALSE.)198 CALL check_var(fname, varname, evap, klon, largest, .FALSE.) 199 199 varname = 'rain_fall' 200 CALL check_var(fname, varn , rain_fall, klon, largest, .FALSE.)200 CALL check_var(fname, varname, rain_fall, klon, largest, .FALSE.) 201 201 varname = 'snow_fall' 202 CALL check_var(fname, varn , snow_fall, klon, largest, .FALSE.)202 CALL check_var(fname, varname, snow_fall, klon, largest, .FALSE.) 203 203 varname = 'ts' 204 CALL check_var(fname, varn , ts, klon, largest, .FALSE.)204 CALL check_var(fname, varname, ts, klon, largest, .FALSE.) 205 205 END IF 206 206 STOP … … 216 216 IF (airetot .NE. airetot .OR. ABS(airetot) > largest) THEN 217 217 varname = 'airephy' 218 CALL check_var(fname, varn , airephy, klon, largest, .FALSE.)218 CALL check_var(fname, varname, airephy, klon, largest, .FALSE.) 219 219 END IF 220 220 IF (evap_tot .NE. evap_tot .OR. ABS(evap_tot) > largest) THEN 221 221 varname = 'evap' 222 CALL check_var(fname, varn , evap, klon, largest, .FALSE.)222 CALL check_var(fname, varname, evap, klon, largest, .FALSE.) 223 223 END IF 224 224 PRINT *,' rain_fall_tot= ',rain_fall_tot 225 225 IF (rain_fall_tot .NE. rain_fall_tot .OR. ABS(rain_fall_tot) > largest) THEN 226 226 varname = 'rain_fall' 227 CALL check_var(fname, varn , rain_fall, klon, largest, .FALSE.)227 CALL check_var(fname, varname, rain_fall, klon, largest, .FALSE.) 228 228 END IF 229 229 PRINT *,' snow_fall_tot= ',snow_fall_tot 230 230 IF (snow_fall_tot .NE. snow_fall_tot .OR. ABS(snow_fall_tot) > largest) THEN 231 231 varname = 'snow_fall' 232 CALL check_var(fname, varn , snow_fall, klon, largest, .FALSE.)232 CALL check_var(fname, varname, snow_fall, klon, largest, .FALSE.) 233 233 END IF 234 234 STOP … … 277 277 wrongi = 0 278 278 DO i=1,sizev 279 IF (var(i) /= var(i) .OR. ABS(var(i)) > bigvalu le ) THEN279 IF (var(i) /= var(i) .OR. ABS(var(i)) > bigvalue ) THEN 280 280 IF (wrongi == 0) found = .TRUE. 281 281 wrongi = wrongi + 1 … … 288 288 PRINT *,TRIM(errmsg) 289 289 PRINT *," at '" // TRIM(funcn) // "' variable '" //TRIM(varn)// & 290 "' is wrong at i value___ '290 "' is wrong at i value___" 291 291 DO i=1,wrongi 292 292 PRINT *,wronggridpt(i), wrongvalues(i)
Note: See TracChangeset
for help on using the changeset viewer.