Changeset 214 for palm/trunk
- Timestamp:
- Nov 14, 2008 10:06:34 AM (16 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/check_parameters.f90
r213 r214 403 403 IF ( momentum_advec /= 'pw-scheme' .AND. momentum_advec /= 'ups-scheme' ) & 404 404 THEN 405 IF ( myid == 0 ) PRINT*, '+++ check_parameters: unknown advection ',&406 'scheme: momentum_advec=', momentum_advec407 CALL local_stop405 message_string = 'unknown advection scheme: momentum_advec = "' // & 406 TRIM( momentum_advec ) // '"' 407 CALL handle_palm_message( 'check_parameters', 'PA0022', 1, 2, 0, 6, 0 ) 408 408 ENDIF 409 409 IF ( ( momentum_advec == 'ups-scheme' .OR. scalar_advec == 'ups-scheme' )& 410 410 .AND. timestep_scheme /= 'euler' ) THEN 411 IF ( myid == 0 ) PRINT*, '+++ check_parameters: momentum_advec=',&412 momentum_advec, ' is not allowed with ',&413 'timestep_scheme=', timestep_scheme414 CALL local_stop411 message_string = 'momentum_advec = "' // TRIM( momentum_advec ) // & 412 '" is not allowed with timestep_scheme = "' // & 413 TRIM( timestep_scheme ) // '"' 414 CALL handle_palm_message( 'check_parameters', 'PA0023', 1, 2, 0, 6, 0 ) 415 415 ENDIF 416 416 417 417 IF ( scalar_advec /= 'pw-scheme' .AND. scalar_advec /= 'bc-scheme' .AND.& 418 418 scalar_advec /= 'ups-scheme' ) THEN 419 IF ( myid == 0 ) PRINT*, '+++ check_parameters: unknown advection ',&420 'scheme: scalar_advec=', scalar_advec421 CALL local_stop419 message_string = 'unknown advection scheme: scalar_advec = "' // & 420 TRIM( scalar_advec ) // '"' 421 CALL handle_palm_message( 'check_parameters', 'PA0024', 1, 2, 0, 6, 0 ) 422 422 ENDIF 423 423 424 424 IF ( use_sgs_for_particles .AND. .NOT. use_upstream_for_tke ) THEN 425 425 use_upstream_for_tke = .TRUE. 426 IF ( myid == 0 ) THEN 427 PRINT*, '+++ WARNING: check_parameters: use_upstream_for_tke set ', & 428 '.TRUE. because use_sgs_for_particles = .TRUE.' 429 ENDIF 426 message_string = 'use_upstream_for_tke set .TRUE. because ' // & 427 'use_sgs_for_particles = .TRUE.' 428 CALL handle_palm_message( 'check_parameters', 'PA0025', 0, 1, 0, 6, 0 ) 430 429 ENDIF 431 430 432 431 IF ( use_upstream_for_tke .AND. timestep_scheme(1:8) == 'leapfrog' ) THEN 433 IF ( myid == 0 ) THEN 434 PRINT*, '+++ check_parameters: use_upstream_for_tke = .TRUE. ', & 435 'not allowed with timestep_scheme=', timestep_scheme 436 ENDIF 437 CALL local_stop 432 message_string = 'use_upstream_for_tke = .TRUE. not allowed with ' // & 433 'timestep_scheme = "' // TRIM( timestep_scheme ) // '"' 434 CALL handle_palm_message( 'check_parameters', 'PA0026', 1, 2, 0, 6, 0 ) 438 435 ENDIF 439 436 … … 458 455 459 456 CASE DEFAULT 460 IF ( myid == 0 ) PRINT*, '+++ check_parameters: unknown timestep ',& 461 'scheme: timestep_scheme=', timestep_scheme 462 CALL local_stop 457 message_string = 'unknown timestep scheme: timestep_scheme = "' // & 458 TRIM( timestep_scheme ) // '"' 459 CALL handle_palm_message( 'check_parameters', 'PA0027', 1, 2, 0, 6, & 460 0 ) 463 461 464 462 END SELECT … … 466 464 IF ( scalar_advec == 'ups-scheme' .AND. timestep_scheme(1:5) == 'runge' )& 467 465 THEN 468 IF ( myid == 0 ) THEN 469 PRINT*, '+++ check_parameters: scalar advection scheme "', & 470 TRIM( scalar_advec ), '"' 471 PRINT*, ' does not work with timestep_scheme "', & 472 TRIM( timestep_scheme ), '"' 473 ENDIF 474 CALL local_stop 466 message_string = 'scalar advection scheme "' // TRIM( scalar_advec ) // & 467 '" & does not work with timestep_scheme "' // & 468 TRIM( timestep_scheme ) // '"' 469 CALL handle_palm_message( 'check_parameters', 'PA0028', 1, 2, 0, 6, 0 ) 475 470 ENDIF 476 471 477 472 IF ( momentum_advec /= 'pw-scheme' .AND. timestep_scheme(1:5) == 'runge' ) & 478 473 THEN 479 IF ( myid == 0 ) THEN 480 PRINT*, '+++ check_parameters: momentum advection scheme "', & 481 TRIM( momentum_advec ), '"' 482 PRINT*, ' does not work with timestep_scheme "', & 483 TRIM( timestep_scheme ), '"' 484 ENDIF 485 CALL local_stop 486 ENDIF 487 488 489 IF ( initializing_actions == ' ' ) THEN 490 IF ( myid == 0 ) THEN 491 PRINT*, '+++ check parameters:' 492 PRINT*, ' no value found for initializing_actions' 493 ENDIF 494 CALL local_stop 474 message_string = 'momentum advection scheme "' // & 475 TRIM( momentum_advec ) // '" & does not work with ' // & 476 'timestep_scheme "' // TRIM( timestep_scheme ) // '"' 477 CALL handle_palm_message( 'check_parameters', 'PA0029', 1, 2, 0, 6, 0 ) 495 478 ENDIF 496 479 … … 498 481 TRIM( initializing_actions ) /= 'read_data_for_recycling' ) THEN 499 482 ! 500 !-- No model continuation run;several initialising actions are possible483 !-- No restart run: several initialising actions are possible 501 484 action = initializing_actions 502 485 DO WHILE ( TRIM( action ) /= '' ) … … 509 492 510 493 CASE DEFAULT 511 IF ( myid == 0 ) PRINT*, '+++ check_parameters: initializi',&512 'ng_action unkown or not allowed: action = "', &513 TRIM(action), '"'514 CALL local_stop494 message_string = 'initializing_action = "' // & 495 TRIM( action ) // '" unkown or not allowed' 496 CALL handle_palm_message( 'check_parameters', 'PA0030', 1, 2, & 497 0, 6, 0 ) 515 498 516 499 END SELECT 517 500 ENDDO 518 501 ENDIF 502 519 503 IF ( INDEX( initializing_actions, 'set_constant_profiles' ) /= 0 .AND. & 520 504 INDEX( initializing_actions, 'set_1d-model_profiles' ) /= 0 ) THEN 521 IF ( myid == 0 ) PRINT*, '+++ check_parameters: initializing_actions', & 522 '"set_constant_profiles" and "set_1d-model_profiles" are not', & 523 ' allowed simultaneously' 524 CALL local_stop 525 ENDIF 505 message_string = 'initializing_actions = "set_constant_profiles"' // & 506 ' and "set_1d-model_profiles" are not allowed ' // & 507 'simultaneously' 508 CALL handle_palm_message( 'check_parameters', 'PA0031', 1, 2, 0, 6, 0 ) 509 ENDIF 510 526 511 IF ( INDEX( initializing_actions, 'set_constant_profiles' ) /= 0 .AND. & 527 512 INDEX( initializing_actions, 'by_user' ) /= 0 ) THEN 528 IF ( myid == 0 ) PRINT*, '+++ check_parameters: initializing_actions',&529 '"set_constant_profiles" and "by_user" are not', &530 ' allowed simultaneously'531 CALL local_stop532 ENDIF 513 message_string = 'initializing_actions = "set_constant_profiles"' // & 514 ' and "by_user" are not allowed simultaneously' 515 CALL handle_palm_message( 'check_parameters', 'PA0032', 1, 2, 0, 6, 0 ) 516 ENDIF 517 533 518 IF ( INDEX( initializing_actions, 'by_user' ) /= 0 .AND. & 534 519 INDEX( initializing_actions, 'set_1d-model_profiles' ) /= 0 ) THEN 535 IF ( myid == 0 ) PRINT*, '+++ check_parameters: initializing_actions', & 536 '"by_user" and "set_1d-model_profiles" are not', & 537 ' allowed simultaneously' 538 CALL local_stop 520 message_string = 'initializing_actions = "by_user" and ' // & 521 '"set_1d-model_profiles" are not allowed simultaneously' 522 CALL handle_palm_message( 'check_parameters', 'PA0033', 1, 2, 0, 6, 0 ) 539 523 ENDIF 540 524 541 525 IF ( cloud_physics .AND. .NOT. humidity ) THEN 542 IF ( myid == 0 ) PRINT*, '+++ check_parameters: cloud_physics =', & 543 cloud_physics, ' is not allowed with ', & 544 'humidity =', humidity 545 CALL local_stop 526 WRITE( message_string, * ) 'cloud_physics = ', cloud_physics, ' is ', & 527 'not allowed with humidity = ', humidity 528 CALL handle_palm_message( 'check_parameters', 'PA0034', 1, 2, 0, 6, 0 ) 546 529 ENDIF 547 530 548 531 IF ( precipitation .AND. .NOT. cloud_physics ) THEN 549 IF ( myid == 0 ) PRINT*, '+++ check_parameters: precipitation =', & 550 precipitation, ' is not allowed with ', & 551 'cloud_physics =', cloud_physics 552 CALL local_stop 532 WRITE( message_string, * ) 'precipitation = ', precipitation, ' is ', & 533 'not allowed with cloud_physics = ', cloud_physics 534 CALL handle_palm_message( 'check_parameters', 'PA0035', 1, 2, 0, 6, 0 ) 553 535 ENDIF 554 536 555 537 IF ( humidity .AND. sloping_surface ) THEN 556 IF ( myid == 0 ) PRINT*, '+++ check_parameters: humidity = TRUE', & 557 'and hang = TRUE are not', & 558 ' allowed simultaneously' 559 CALL local_stop 538 message_string = 'humidity = .TRUE. and sloping_surface = .TRUE. ' // & 539 'are not allowed simultaneously' 540 CALL handle_palm_message( 'check_parameters', 'PA0036', 1, 2, 0, 6, 0 ) 560 541 ENDIF 561 542 562 543 IF ( humidity .AND. scalar_advec == 'ups-scheme' ) THEN 563 IF ( myid == 0 ) PRINT*, '+++ check_parameters: UPS-scheme', & 564 'is not implemented for humidity' 565 CALL local_stop 544 message_string = 'UPS-scheme is not implemented for humidity = .TRUE.' 545 CALL handle_palm_message( 'check_parameters', 'PA0037', 1, 2, 0, 6, 0 ) 566 546 ENDIF 567 547 568 548 IF ( passive_scalar .AND. humidity ) THEN 569 IF ( myid == 0 ) PRINT*, '+++ check_parameters: humidity = TRUE and', & 570 'passive_scalar = TRUE is not allowed ', & 571 'simultaneously' 572 CALL local_stop 549 message_string = 'humidity = .TRUE. and passive_scalar = .TRUE. ' // & 550 'is not allowed simultaneously' 551 CALL handle_palm_message( 'check_parameters', 'PA0038', 1, 2, 0, 6, 0 ) 573 552 ENDIF 574 553 575 554 IF ( passive_scalar .AND. scalar_advec == 'ups-scheme' ) THEN 576 IF ( myid == 0 ) PRINT*, '+++ check_parameters: UPS-scheme',&577 'is not implemented for passive_scalar'578 CALL local_stop555 message_string = 'UPS-scheme is not implemented for passive_scalar' // & 556 ' = .TRUE.' 557 CALL handle_palm_message( 'check_parameters', 'PA0039', 1, 2, 0, 6, 0 ) 579 558 ENDIF 580 559 581 560 IF ( grid_matching /= 'strict' .AND. grid_matching /= 'match' ) THEN 582 IF ( myid == 0 ) PRINT*, '+++ check_parameters: illegal value "', & 583 TRIM( grid_matching ), & 584 '" found for parameter grid_matching' 585 CALL local_stop 561 message_string = 'illegal value "' // TRIM( grid_matching ) // & 562 '" found for parameter grid_matching' 563 CALL handle_palm_message( 'check_parameters', 'PA0040', 1, 2, 0, 6, 0 ) 586 564 ENDIF 587 565 588 566 IF ( plant_canopy .AND. ( drag_coefficient == 0.0 ) ) THEN 589 IF ( myid == 0 ) PRINT*, '+++ check_parameters: plant_canopy = TRUE',&590 'requires a non-zero drag_coefficient'591 CALL local_stop567 message_string = 'plant_canopy = .TRUE. requires a non-zero drag ' // & 568 'coefficient & given value is drag_coefficient = 0.0' 569 CALL handle_palm_message( 'check_parameters', 'PA0041', 1, 2, 0, 6, 0 ) 592 570 ENDIF 593 571
Note: See TracChangeset
for help on using the changeset viewer.