- Timestamp:
- Jul 29, 2014, 12:29:13 PM (10 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/WRFV3/dyn_em/solve_em.F
r6 r143 859 859 im2 = ims + (ime - ims) / 2 860 860 jm2 = jms + (jme - jms) / 2 861 862 IF (grid%t_2(im2,iz,jm2) /= grid%t_2(im2,iz,jm2) .OR. ABS(grid%t_2(im2,iz,jm2)) > 10000. ) THEN 863 PRINT errmsg 864 WRITE(wrF_err_message,*)'solve_em: wrong T value=', & 865 grid%t_2(im2,iz,jm2),' at: ', im2,', ', iz,', ', jm2,' !!!' 861 DO iz = kms, kme 862 863 IF (grid%t_2(im2,iz,jm2) /= grid%t_2(im2,iz,jm2) .OR. ABS(grid%t_2(im2,iz,jm2)) > 10000. ) THEN 864 PRINT *,TRIM(errmsg) 865 WRITE(wrF_err_message,*)'solve_em: wrong T value=', & 866 grid%t_2(im2,iz,jm2),' at: ', im2,', ', iz,', ', jm2,' !!!' 866 867 #ifdef DM_PARALLEL 867 CALL wrf_error_fatal(TRIM(wrf_err_message))868 CALL wrf_error_fatal(TRIM(wrf_err_message)) 868 869 #else 869 PRINT *,TRIM(wrf_err_message) 870 STOP 871 #endif 872 END IF 870 PRINT *,TRIM(wrf_err_message) 871 STOP 872 #endif 873 END IF 874 END DO 873 875 ! Checking for NaNs (should not be necessary but....) 874 876 IF (grid%psfc(im2,jm2) /= grid%psfc(im2,jm2) .OR. ABS(grid%psfc(im2,jm2)) > 1000000. ) THEN 875 PRINT errmsg877 PRINT *,errmsg 876 878 WRITE(wrF_err_message,*)'solve_em: wrong PSFC value=', & 877 879 grid%psfc(im2,jm2),' at: ', im2 ,', ', jm2, ' !!!'
Note: See TracChangeset
for help on using the changeset viewer.