Changeset 3045 for palm/trunk/SOURCE/large_scale_forcing_nudging_mod.f90
- Timestamp:
- May 28, 2018 7:55:41 AM (6 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/large_scale_forcing_nudging_mod.f90
r3026 r3045 25 25 ! ----------------- 26 26 ! $Id$ 27 ! Error messages revised 28 ! 29 ! 3026 2018-05-22 10:30:53Z schwenkel 27 30 ! Changed the name specific humidity to mixing ratio, since we are computing 28 31 ! mixing ratios. … … 732 735 !-- Check nudging and large scale forcing from external file 733 736 IF ( nudging .AND. ( .NOT. large_scale_forcing ) ) THEN 734 message_string = 'Nudging requires large_scale_forcing = .T.. &'//&735 'Surface fluxes and geostrophic wind should be &'//&737 message_string = 'Nudging requires large_scale_forcing = .T.. '// & 738 'Surface fluxes and geostrophic wind should be '// & 736 739 'prescribed in file LSF_DATA' 737 740 CALL message( 'check_parameters', 'PA0374', 1, 2, 0, 6, 0 ) … … 740 743 IF ( large_scale_forcing .AND. ( bc_lr /= 'cyclic' .OR. & 741 744 bc_ns /= 'cyclic' ) ) THEN 742 message_string = 'Non-cyclic lateral boundaries do not allow for &'//&745 message_string = 'Non-cyclic lateral boundaries do not allow for '//& 743 746 'the usage of large scale forcing from external file.' 744 747 CALL message( 'check_parameters', 'PA0375', 1, 2, 0, 6, 0 ) … … 746 749 747 750 IF ( large_scale_forcing .AND. ( .NOT. humidity ) ) THEN 748 message_string = 'The usage of large scale forcing from external &'//&751 message_string = 'The usage of large scale forcing from external '//& 749 752 'file LSF_DATA requires humidity = .T..' 750 753 CALL message( 'check_parameters', 'PA0376', 1, 2, 0, 6, 0 ) … … 752 755 753 756 IF ( large_scale_forcing .AND. passive_scalar ) THEN 754 message_string = 'The usage of large scale forcing from external &'//&757 message_string = 'The usage of large scale forcing from external '// & 755 758 'file LSF_DATA is not implemented for passive scalars' 756 759 CALL message( 'check_parameters', 'PA0440', 1, 2, 0, 6, 0 ) … … 759 762 IF ( large_scale_forcing .AND. topography /= 'flat' & 760 763 .AND. .NOT. lsf_exception ) THEN 761 message_string = 'The usage of large scale forcing from external &'//&764 message_string = 'The usage of large scale forcing from external '//& 762 765 'file LSF_DATA is not implemented for non-flat topography' 763 766 CALL message( 'check_parameters', 'PA0377', 1, 2, 0, 6, 0 ) … … 765 768 766 769 IF ( large_scale_forcing .AND. ocean ) THEN 767 message_string = 'The usage of large scale forcing from external &'//&770 message_string = 'The usage of large scale forcing from external '// & 768 771 'file LSF_DATA is not implemented for ocean runs' 769 772 CALL message( 'check_parameters', 'PA0378', 1, 2, 0, 6, 0 ) … … 858 861 TRIM( data_output_pr(var_count) ) // & 859 862 ' is not implemented for ' // & 860 ' large_scale_forcing = .FALSE.'861 CALL message( 'lsf_nudging_check_data_output_pr', 'PA039 3', &863 'nudging = .FALSE.' 864 CALL message( 'lsf_nudging_check_data_output_pr', 'PA0394', & 862 865 1, 2, 0, 6, 0 ) 863 866 ELSE … … 873 876 TRIM( data_output_pr(var_count) ) // & 874 877 ' is not implemented for ' // & 875 ' large_scale_forcing = .FALSE.'876 CALL message( 'lsf_nudging_check_data_output_pr', 'PA039 3', &878 'nudging = .FALSE.' 879 CALL message( 'lsf_nudging_check_data_output_pr', 'PA0394', & 877 880 1, 2, 0, 6, 0 ) 878 881 ELSE … … 888 891 TRIM( data_output_pr(var_count) ) // & 889 892 ' is not implemented for ' // & 890 ' large_scale_forcing = .FALSE.'891 CALL message( 'lsf_nudging_check_data_output_pr', 'PA039 3', &893 'nudging = .FALSE.' 894 CALL message( 'lsf_nudging_check_data_output_pr', 'PA0394', & 892 895 1, 2, 0, 6, 0 ) 893 896 ELSE … … 903 906 TRIM( data_output_pr(var_count) ) // & 904 907 ' is not implemented for ' // & 905 ' large_scale_forcing = .FALSE.'906 CALL message( 'lsf_nudging_check_data_output_pr', 'PA039 3', &908 'nudging = .FALSE.' 909 CALL message( 'lsf_nudging_check_data_output_pr', 'PA0394', & 907 910 1, 2, 0, 6, 0 ) 908 911 ELSE … … 1280 1283 1281 1284 IF ( ierrn /= 0 ) THEN 1282 WRITE ( message_string, * ) 'No time dependent surface variables ' //&1283 ' in&LSF_DATA for end of run found'1285 WRITE ( message_string, * ) 'No time dependent surface ' // & 1286 'variables in LSF_DATA for end of run found' 1284 1287 1285 1288 CALL message( 'ls_forcing', 'PA0363', 1, 2, 0, 6, 0 ) … … 1289 1292 IF ( time_surf(1) > end_time ) THEN 1290 1293 WRITE ( message_string, * ) 'Time dependent surface variables in ' // & 1291 ' &LSF_DATA set in after end of ' ,&1294 'LSF_DATA set in after end of ' , & 1292 1295 'simulation - lsf_surf is set to FALSE' 1293 1296 CALL message( 'ls_forcing', 'PA0371', 0, 0, 0, 6, 0 ) … … 1316 1319 IF ( ierrn < 0 ) THEN 1317 1320 WRITE( message_string, * ) 'No time dependent vertical profiles',& 1318 ' in &LSF_DATA for end of run found'1321 ' in LSF_DATA for end of run found' 1319 1322 CALL message( 'ls_forcing', 'PA0372', 1, 2, 0, 6, 0 ) 1320 1323 ENDIF … … 1362 1365 1363 1366 IF ( ierrn /= 0 ) THEN 1364 WRITE( message_string, * ) 'zu(',k,') = ', zu(k), 'm ', &1365 'is higher than the maximum height in LSF_DATA which ',&1366 ' is ', lowheight, 'm. Interpolation on PALM ',&1367 WRITE( message_string, * ) 'zu(',k,') = ', zu(k), 'm ', & 1368 'is higher than the maximum height in LSF_DATA ', & 1369 'which is ', lowheight, 'm. Interpolation on PALM ',& 1367 1370 'grid is not possible.' 1368 1371 CALL message( 'ls_forcing', 'PA0395', 1, 2, 0, 6, 0 ) … … 1832 1835 1833 1836 ierrn = 0 1834 READ ( finput, *, IOSTAT=ierrn ) highheight , hightnudge , &1835 highunudge , highvnudge , &1836 highwnudge , highptnudge, &1837 READ ( finput, *, IOSTAT=ierrn ) highheight , hightnudge , & 1838 highunudge , highvnudge , & 1839 highwnudge , highptnudge, & 1837 1840 highqnudge 1838 1841 IF (ierrn /= 0 ) THEN 1839 WRITE( message_string, * ) 'zu(',k,') = ', zu(k), 'm is ', &1840 'higher than the maximum height in NUDING_DATA which ', 1841 'is ', lowheight, 'm. Interpolation on PALM ', 1842 WRITE( message_string, * ) 'zu(',k,') = ', zu(k), 'm is ', & 1843 'higher than the maximum height in NUDING_DATA which ',& 1844 'is ', lowheight, 'm. Interpolation on PALM ', & 1842 1845 'grid is not possible.' 1843 1846 CALL message( 'nudging', 'PA0364', 1, 2, 0, 6, 0 ) … … 1872 1875 ENDIF 1873 1876 1874 WRITE( message_string, * ) 'Initial profiles of u, v and ',&1875 ' scalarsfrom NUDGING_DATA are used.'1877 WRITE( message_string, * ) 'Initial profiles of u, v, pt and q ', & 1878 'from NUDGING_DATA are used.' 1876 1879 CALL message( 'large_scale_forcing_nudging', 'PA0370', 0, 0, 0, 6, 0 ) 1877 1880 ENDIF
Note: See TracChangeset
for help on using the changeset viewer.
