Changeset 3875 for palm/trunk/SOURCE/wind_turbine_model_mod.f90
- Timestamp:
- Apr 8, 2019 5:35:12 PM (5 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/wind_turbine_model_mod.f90
r3832 r3875 26 26 ! ----------------- 27 27 ! $Id$ 28 ! Addaped wtm_tendency to fit the module actions interface 29 ! 30 ! 3832 2019-03-28 13:16:58Z raasch 28 31 ! instrumented with openmp directives 29 32 ! … … 463 466 END INTERFACE wtm_parin 464 467 465 INTERFACE wtm_wrd_global466 MODULE PROCEDURE wtm_wrd_global467 END INTERFACE wtm_wrd_global468 469 INTERFACE wtm_rrd_global470 MODULE PROCEDURE wtm_rrd_global471 END INTERFACE wtm_rrd_global472 473 468 INTERFACE wtm_check_parameters 474 469 MODULE PROCEDURE wtm_check_parameters 475 470 END INTERFACE wtm_check_parameters 476 471 477 472 INTERFACE wtm_init_arrays 478 473 MODULE PROCEDURE wtm_init_arrays … … 482 477 MODULE PROCEDURE wtm_init 483 478 END INTERFACE wtm_init 484 485 INTERFACE wtm_read_blade_tables 486 MODULE PROCEDURE wtm_read_blade_tables 487 END INTERFACE wtm_read_blade_tables 488 489 INTERFACE wtm_forces 490 MODULE PROCEDURE wtm_forces 491 END INTERFACE wtm_forces 492 493 INTERFACE wtm_yawcontrol 494 MODULE PROCEDURE wtm_yawcontrol 495 END INTERFACE wtm_yawcontrol 496 497 INTERFACE wtm_rotate_rotor 498 MODULE PROCEDURE wtm_rotate_rotor 499 END INTERFACE wtm_rotate_rotor 500 501 INTERFACE wtm_speed_control 502 MODULE PROCEDURE wtm_init_speed_control 503 MODULE PROCEDURE wtm_speed_control 504 END INTERFACE wtm_speed_control 505 506 INTERFACE wtm_tendencies 507 MODULE PROCEDURE wtm_tendencies 508 MODULE PROCEDURE wtm_tendencies_ij 509 END INTERFACE wtm_tendencies 510 511 512 PUBLIC wtm_check_parameters, wtm_forces, wtm_init, wtm_init_arrays, & 513 wtm_parin, wtm_wrd_global, wtm_rrd_global, wtm_tendencies 479 480 INTERFACE wtm_actions 481 MODULE PROCEDURE wtm_actions 482 MODULE PROCEDURE wtm_actions_ij 483 END INTERFACE wtm_actions 484 485 INTERFACE wtm_rrd_global 486 MODULE PROCEDURE wtm_rrd_global 487 END INTERFACE wtm_rrd_global 488 489 INTERFACE wtm_wrd_global 490 MODULE PROCEDURE wtm_wrd_global 491 END INTERFACE wtm_wrd_global 492 493 494 PUBLIC & 495 wtm_parin, & 496 wtm_check_parameters, & 497 wtm_init_arrays, & 498 wtm_init, & 499 wtm_actions, & 500 wtm_rrd_global, & 501 wtm_wrd_global 514 502 515 503 … … 1085 1073 IF ( speed_control) THEN 1086 1074 1087 CALL wtm_ speed_control1075 CALL wtm_init_speed_control 1088 1076 1089 1077 IF ( TRIM( initializing_actions ) == 'read_restart_data' ) THEN … … 2771 2759 !> Call for all grid points 2772 2760 !------------------------------------------------------------------------------! 2773 SUBROUTINE wtm_tendencies( component ) 2774 2775 2776 IMPLICIT NONE 2777 2778 INTEGER(iwp) :: component !< prognostic variable (u,v,w) 2761 SUBROUTINE wtm_actions( location ) 2762 2763 2764 CHARACTER (LEN=*) :: location !< 2765 2779 2766 INTEGER(iwp) :: i !< running index 2780 2767 INTEGER(iwp) :: j !< running index … … 2782 2769 2783 2770 2784 SELECT CASE ( component ) 2785 2786 CASE ( 1 ) 2771 SELECT CASE ( location ) 2772 2773 CASE ( 'before_timestep' ) 2774 2775 CALL wtm_forces 2776 2777 CASE ( 'u-tendency' ) 2787 2778 ! 2788 2779 !-- Apply the x-component of the force to the u-component of the flow: … … 2807 2798 ENDIF 2808 2799 2809 CASE ( 2)2800 CASE ( 'v-tendency' ) 2810 2801 ! 2811 2802 !-- Apply the y-component of the force to the v-component of the flow: … … 2827 2818 ENDIF 2828 2819 2829 CASE ( 3)2820 CASE ( 'w-tendency' ) 2830 2821 ! 2831 2822 !-- Apply the z-component of the force to the w-component of the flow: … … 2845 2836 CASE DEFAULT 2846 2837 2847 WRITE( message_string, * ) 'unknown prognostic variable: ', component2848 CALL message( 'wtm_ tendencies', 'PA04??', 1, 2, 0, 6, 0 )2838 WRITE( message_string, * ) 'unknown prognostic variable: ', TRIM( location ) 2839 CALL message( 'wtm_actions', 'PA04??', 1, 2, 0, 6, 0 ) 2849 2840 2850 2841 END SELECT 2851 2842 2852 2843 2853 END SUBROUTINE wtm_ tendencies2844 END SUBROUTINE wtm_actions 2854 2845 2855 2846 … … 2861 2852 !> Call for grid point i,j 2862 2853 !------------------------------------------------------------------------------! 2863 SUBROUTINE wtm_tendencies_ij( i, j, component ) 2864 2865 2866 IMPLICIT NONE 2867 2868 INTEGER(iwp) :: component !< prognostic variable (u,v,w) 2854 SUBROUTINE wtm_actions_ij( i, j, location ) 2855 2856 2857 CHARACTER (LEN=*) :: location !< 2869 2858 INTEGER(iwp) :: i !< running index 2870 2859 INTEGER(iwp) :: j !< running index 2871 2860 INTEGER(iwp) :: k !< running index 2872 2861 2873 SELECT CASE ( component ) 2874 2875 CASE ( 1 ) 2862 SELECT CASE ( location ) 2863 2864 CASE ( 'before_timestep' ) 2865 2866 CALL wtm_forces 2867 2868 CASE ( 'u-tendency' ) 2876 2869 ! 2877 2870 !-- Apply the x-component of the force to the u-component of the flow: … … 2891 2884 ENDIF 2892 2885 2893 CASE ( 2)2886 CASE ( 'v-tendency' ) 2894 2887 ! 2895 2888 !-- Apply the y-component of the force to the v-component of the flow: … … 2907 2900 ENDIF 2908 2901 2909 CASE ( 3)2902 CASE ( 'w-tendency' ) 2910 2903 ! 2911 2904 !-- Apply the z-component of the force to the w-component of the flow: … … 2921 2914 CASE DEFAULT 2922 2915 2923 WRITE( message_string, * ) 'unknown prognostic variable: ', component2924 CALL message( 'wtm_ tendencies', 'PA04??', 1, 2, 0, 6, 0 )2916 WRITE( message_string, * ) 'unknown prognostic variable: ', TRIM( location ) 2917 CALL message( 'wtm_actions', 'PA04??', 1, 2, 0, 6, 0 ) 2925 2918 2926 2919 END SELECT 2927 2920 2928 2921 2929 END SUBROUTINE wtm_ tendencies_ij2922 END SUBROUTINE wtm_actions_ij 2930 2923 2931 2924 END MODULE wind_turbine_model_mod
Note: See TracChangeset
for help on using the changeset viewer.