Changeset 2563
- Timestamp:
- Oct 19, 2017 3:36:10 PM (7 years ago)
- Location:
- palm/trunk
- Files:
-
- 11 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SCRIPTS/NCL/spectra.ncl
r2030 r2563 685 685 res@trXMaxF = max(max_x) 686 686 res@trYMinF = min(min_y) 687 res@trYMaxF = max(max_y )687 res@trYMaxF = max(max_y*10) 688 688 689 689 plot_h(q) = gsn_csm_xy(wks,x_axis(q,:),\ -
palm/trunk/SOURCE/Makefile
r2544 r2563 25 25 # ----------------- 26 26 # $Id$ 27 # wind_turbine_model_mod.o and synthetic_turbulence_generator_mod.o were added to 28 # write_var_list.o and virtual_flight_mod.o was deleted from read_var_list.o 29 # 30 # 2544 2017-10-13 18:09:32Z maronga 27 31 # Added date_and_time_mod 28 32 # … … 626 630 spectra_mod.o surface_mod.o urban_surface_mod.o 627 631 read_var_list.o: modules.o mod_kinds.o model_1d_mod.o netcdf_interface_mod.o plant_canopy_model_mod.o \ 628 spectra_mod.o microphysics_mod.o urban_surface_mod.o vertical_nesting_mod.o virtual_flight_mod.o632 spectra_mod.o microphysics_mod.o urban_surface_mod.o vertical_nesting_mod.o 629 633 run_control.o: modules.o cpulog_mod.o mod_kinds.o 630 634 set_slicer_attributes_dvrp.o: modules.o mod_kinds.o … … 697 701 spectra_mod.o surface_mod.o 698 702 write_var_list.o: modules.o mod_kinds.o model_1d_mod.o netcdf_interface_mod.o plant_canopy_model_mod.o\ 699 spectra_mod.o microphysics_mod.o urban_surface_mod.o vertical_nesting_mod.o virtual_flight_mod.o 703 spectra_mod.o microphysics_mod.o urban_surface_mod.o vertical_nesting_mod.o virtual_flight_mod.o\ 704 wind_turbine_model_mod.o synthetic_turbulence_generator_mod.o -
palm/trunk/SOURCE/modules.f90
r2550 r2563 25 25 ! ----------------- 26 26 ! $Id$ 27 ! Variable wind_turbine was added to control_parameters 28 ! 29 ! 2550 2017-10-16 17:12:01Z boeske 27 30 ! complex_terrain namelist parameter added 28 31 ! … … 1242 1245 LOGICAL :: virtual_flight = .FALSE. !< use virtual flight model? 1243 1246 LOGICAL :: wall_adjustment = .TRUE. !< namelist parameter 1247 LOGICAL :: wind_turbine = .FALSE. !< flag for use of wind turbine model 1244 1248 LOGICAL :: write_binary = .FALSE. !< ENVPAR namelist parameter to steer restart I/O (ENVPAR is created by mrun) 1245 1249 LOGICAL :: ws_scheme_sca = .FALSE. !< use Wicker-Skamarock scheme (scalar advection)? -
palm/trunk/SOURCE/parin.f90
r2550 r2563 25 25 ! ----------------- 26 26 ! $Id$ 27 ! Changed position where restart files are closed. 28 ! 29 ! 2550 2017-10-16 17:12:01Z boeske 27 30 ! Added complex_terrain 28 31 ! … … 559 562 12 IF ( TRIM( initializing_actions ) == 'read_restart_data' ) THEN 560 563 CALL read_var_list 561 !562 !-- The restart file will be reopened when reading the subdomain data563 CALL close_file( 13 )564 564 565 565 ! … … 698 698 699 699 ! 700 !-- The restart file will be reopened when reading the subdomain data 701 IF ( TRIM( initializing_actions ) == 'read_restart_data' ) THEN 702 CALL close_file( 13 ) 703 ENDIF 704 705 ! 700 706 !-- Check in case of initial run, if the grid point numbers are well 701 707 !-- defined and allocate some arrays which are already needed in -
palm/trunk/SOURCE/prognostic_equations.f90
r2320 r2563 25 25 ! ----------------- 26 26 ! $Id$ 27 ! Variable wind_turbine moved to module control_parameters 28 ! 29 ! 2320 2017-07-21 12:47:43Z suehring 27 30 ! Modularize large-scale forcing and nudging 28 31 ! … … 260 263 scalar_advec, scalar_advec, simulated_time, sloping_surface, & 261 264 timestep_scheme, tsc, use_subsidence_tendencies, & 262 use_upstream_for_tke, ws_scheme_mom, ws_scheme_sca 265 use_upstream_for_tke, wind_turbine, ws_scheme_mom, & 266 ws_scheme_sca 263 267 264 268 USE cpulog, & … … 358 362 359 363 USE wind_turbine_model_mod, & 360 ONLY: w ind_turbine, wtm_tendencies364 ONLY: wtm_tendencies 361 365 362 366 -
palm/trunk/SOURCE/read_var_list.f90
r2372 r2563 25 25 ! ----------------- 26 26 ! $Id$ 27 ! CALL stg_read_restart_data moved to synthetic_turbulence_generator_mod and 28 ! CALL flight_read_restart_data moved to virtual_flights_mod. Furthermore 29 ! *** end default *** marks the end of the standard parameter list of restart 30 ! files and *** end *** marks the end of all parameter including module 31 ! parameter. Therefore the call of flight_skip_var_list becomes unnecessary. 32 ! 33 ! 2372 2017-08-25 12:37:32Z sward 27 34 ! y_shift added to vars, version no. increased 28 35 ! … … 220 227 USE control_parameters 221 228 222 USE flight_mod, &223 ONLY: flight_read_restart_data224 225 229 USE grid_variables, & 226 230 ONLY: dx, dy … … 255 259 v_max, v_max_ijk, w_max, w_max_ijk 256 260 257 USE synthetic_turbulence_generator_mod, &258 ONLY: stg_read_restart_data259 260 261 USE vertical_nesting_mod, & 261 262 ONLY: vnest_init … … 272 273 !-- Make version number check first 273 274 READ ( 13 ) version_on_file 274 binary_version = '4. 2'275 binary_version = '4.3' 275 276 IF ( TRIM( version_on_file ) /= TRIM( binary_version ) ) THEN 276 277 WRITE( message_string, * ) 'version mismatch concerning control ', & … … 349 350 !-- increased. The same changes must also be done in write_var_list. 350 351 READ ( 13 ) variable_chr 351 DO WHILE ( TRIM( variable_chr ) /= '*** end ***' )352 DO WHILE ( TRIM( variable_chr ) /= '*** end default ***' ) 352 353 353 354 SELECT CASE ( TRIM( variable_chr ) ) … … 842 843 ENDDO 843 844 844 IF ( virtual_flight ) CALL flight_read_restart_data845 846 IF ( synthetic_turbulence_generator ) CALL stg_read_restart_data847 845 848 846 END SUBROUTINE read_var_list … … 984 982 READ ( 13 ) variable_chr 985 983 986 DO WHILE ( TRIM( variable_chr ) /= '*** end ***' )984 DO WHILE ( TRIM( variable_chr ) /= '*** end default ***' ) 987 985 988 986 SELECT CASE ( TRIM( variable_chr ) ) … … 1072 1070 SUBROUTINE skip_var_list 1073 1071 1074 USE control_parameters, &1075 ONLY: virtual_flight1076 1077 USE flight_mod, &1078 ONLY: flight_skip_var_list1079 1072 1080 1073 IMPLICIT NONE … … 1096 1089 1097 1090 ENDDO 1098 !1099 !-- In case of virtual flights, skip also variables related to1100 !-- this module.1101 IF ( virtual_flight ) CALL flight_skip_var_list1102 1091 1103 1092 -
palm/trunk/SOURCE/synthetic_turbulence_generator_mod.f90
r2259 r2563 25 25 ! ----------------- 26 26 ! $Id$ 27 ! stg_read_restart_data is called in stg_parin in the case of a restart run 28 ! 29 ! 2259 2017-06-08 09:09:11Z gronemeier 27 30 ! Initial revision 28 31 ! … … 194 197 !-- Public interfaces 195 198 PUBLIC stg_check_parameters, stg_header, stg_init, stg_main, stg_parin, & 196 stg_ read_restart_data, stg_write_restart_data199 stg_write_restart_data 197 200 198 201 ! … … 659 662 !-- on 660 663 synthetic_turbulence_generator = .TRUE. 664 665 IF ( TRIM( initializing_actions ) == 'read_restart_data' ) THEN 666 CALL stg_read_restart_data 667 ENDIF 661 668 662 669 10 CONTINUE -
palm/trunk/SOURCE/time_integration.f90
r2365 r2563 25 25 ! ----------------- 26 26 ! $Id$ 27 ! Variable wind_turbine moved to module control_parameters 28 ! 29 ! 2365 2017-08-21 14:59:59Z kanani 27 30 ! Vertical grid nesting implemented (SadiqHuq) 28 31 ! … … 323 326 use_initial_profile_as_reference, & 324 327 use_single_reference_value, u_gtrans, v_gtrans, virtual_flight, & 325 w s_scheme_mom, ws_scheme_sca328 wind_turbine, ws_scheme_mom, ws_scheme_sca 326 329 327 330 USE cpulog, & … … 400 403 401 404 USE wind_turbine_model_mod, & 402 ONLY: w ind_turbine, wtm_forces405 ONLY: wtm_forces 403 406 404 407 USE vertical_nesting_mod, & -
palm/trunk/SOURCE/virtual_flight_mod.f90
r2271 r2563 25 25 ! ----------------- 26 26 ! $Id$ 27 ! flight_read_restart_data is called in flight_parin in the case of a restart 28 ! run. flight_skip_var_list is not necessary anymore due to marker changes in 29 ! restart files. 30 ! 31 ! 2271 2017-06-09 12:34:55Z sward 27 32 ! Todo added 28 33 ! … … 115 120 END INTERFACE flight_measurement 116 121 117 INTERFACE flight_skip_var_list118 MODULE PROCEDURE flight_skip_var_list119 END INTERFACE flight_skip_var_list120 121 122 INTERFACE flight_read_restart_data 122 123 MODULE PROCEDURE flight_read_restart_data … … 133 134 !-- Public interfaces 134 135 PUBLIC flight_init, flight_header, flight_parin, flight_measurement, & 135 flight_skip_var_list, flight_read_restart_data, & 136 flight_write_restart_data 136 flight_write_restart_data 137 137 ! 138 138 !-- Public variables … … 192 192 SUBROUTINE flight_parin 193 193 194 194 USE control_parameters, & 195 ONLY: initializing_actions 196 195 197 IMPLICIT NONE 196 198 … … 216 218 !-- Set switch that virtual flights shall be carried out 217 219 virtual_flight = .TRUE. 220 221 IF ( TRIM( initializing_actions ) == 'read_restart_data' ) THEN 222 CALL flight_read_restart_data 223 ENDIF 218 224 219 225 10 CONTINUE … … 854 860 855 861 END SUBROUTINE flight_check_parameters 856 857 !------------------------------------------------------------------------------!858 ! Description:859 ! ------------860 !> Skipping the flight-module variables from restart-file (binary format).861 !------------------------------------------------------------------------------!862 SUBROUTINE flight_skip_var_list863 864 IMPLICIT NONE865 866 CHARACTER (LEN=1) :: cdum867 CHARACTER (LEN=30) :: variable_chr868 869 READ ( 13 ) variable_chr870 DO WHILE ( TRIM( variable_chr ) /= '*** end flight ***' )871 READ ( 13 ) cdum872 READ ( 13 ) variable_chr873 ENDDO874 875 END SUBROUTINE flight_skip_var_list876 862 877 863 !------------------------------------------------------------------------------! -
palm/trunk/SOURCE/wind_turbine_model_mod.f90
r2553 r2563 26 26 ! ----------------- 27 27 ! $Id$ 28 ! Restart runs with wind turbine model are possible now. For this purpose, two 29 ! new subroutines wtm_write_restart_data and wtm_read_restart_data had to be 30 ! defined 31 ! 32 ! 2553 2017-10-18 08:03:45Z Giersch 28 33 ! Bugfix of vertical loop in wtm_tendencies to account for different turbine 29 34 ! heights, bugfix of the interpolation of the u-component concerning the … … 104 109 105 110 USE control_parameters, & 106 ONLY: dt_3d, dz, message_string, simulated_time 111 ONLY: dt_3d, dz, message_string, simulated_time, wind_turbine, & 112 initializing_actions 107 113 108 114 USE cpulog, & … … 124 130 125 131 PRIVATE 126 127 LOGICAL :: wind_turbine=.FALSE. !< switch for use of wind turbine model128 132 129 133 ! … … 265 269 REAL(wp), DIMENSION(:), ALLOCATABLE :: alpha_attack !< 266 270 REAL(wp), DIMENSION(:), ALLOCATABLE :: chord !< 267 REAL(wp), DIMENSION(:), ALLOCATABLE :: omega_gen !< curr. generator speed268 271 REAL(wp), DIMENSION(:), ALLOCATABLE :: phi_rel !< 269 REAL(wp), DIMENSION(:), ALLOCATABLE :: pitch_add_old!<270 272 REAL(wp), DIMENSION(:), ALLOCATABLE :: torque_total !< 271 273 REAL(wp), DIMENSION(:), ALLOCATABLE :: thrust_rotor !< … … 337 339 REAL(wp) :: vs_sysp !< 338 340 REAL(wp) :: lp_coeff !< coeff for the controller low pass filter 339 340 REAL(wp), DIMENSION(:), ALLOCATABLE :: omega_gen_old !< last gen. speed341 REAL(wp), DIMENSION(:), ALLOCATABLE :: omega_gen_f !< filtered gen. sp342 REAL(wp), DIMENSION(:), ALLOCATABLE :: omega_gen_f_old !< last filt. gen. sp343 REAL(wp), DIMENSION(:), ALLOCATABLE :: torque_gen !< generator torque344 REAL(wp), DIMENSION(:), ALLOCATABLE :: torque_gen_old !< last gen. torque345 341 346 342 REAL(wp), DIMENSION(100) :: omega_rot_l = 0.0_wp !< local rot speed [rad/s] 343 347 344 ! 348 345 !-- Fixed variables for the yaw controller … … 362 359 INTEGER(iwp) :: WDSHO !< 363 360 361 ! 362 !-- Variables that have to be saved in the binary file for restarts 363 REAL(wp), DIMENSION(1:100) :: pitch_add_old = 0.0_wp !< old constant pitch angle 364 REAL(wp), DIMENSION(1:100) :: omega_gen = 0.0_wp !< curr. generator speed 365 REAL(wp), DIMENSION(1:100) :: omega_gen_f = 0.0_wp !< filtered generator speed 366 REAL(wp), DIMENSION(1:100) :: omega_gen_old = 0.0_wp !< last generator speed 367 REAL(wp), DIMENSION(1:100) :: omega_gen_f_old = 0.0_wp !< last filtered generator speed 368 REAL(wp), DIMENSION(1:100) :: torque_gen = 0.0_wp !< generator torque 369 REAL(wp), DIMENSION(1:100) :: torque_gen_old = 0.0_wp !< last generator torque 370 364 371 365 372 SAVE … … 369 376 MODULE PROCEDURE wtm_parin 370 377 END INTERFACE wtm_parin 378 379 INTERFACE wtm_write_restart_data 380 MODULE PROCEDURE wtm_write_restart_data 381 END INTERFACE wtm_write_restart_data 382 383 INTERFACE wtm_read_restart_data 384 MODULE PROCEDURE wtm_read_restart_data 385 END INTERFACE wtm_read_restart_data 371 386 372 387 INTERFACE wtm_check_parameters … … 410 425 411 426 PUBLIC wtm_check_parameters, wtm_forces, wtm_init, wtm_init_arrays, & 412 wtm_parin, wtm_tendencies, wind_turbine 427 wtm_parin, wtm_write_restart_data, wtm_tendencies 428 413 429 414 430 CONTAINS … … 467 483 wind_turbine = .TRUE. 468 484 485 IF ( TRIM( initializing_actions ) == 'read_restart_data' ) THEN 486 CALL wtm_read_restart_data 487 ENDIF 488 469 489 10 CONTINUE ! TBD Change from continue, mit ierrn machen 470 490 … … 472 492 END SUBROUTINE wtm_parin 473 493 494 495 !------------------------------------------------------------------------------! 496 ! Description: 497 ! ------------ 498 !> This routine writes the respective restart data. 499 !------------------------------------------------------------------------------! 500 SUBROUTINE wtm_write_restart_data 501 502 IMPLICIT NONE 503 504 WRITE ( 14 ) 'omega_gen ' 505 WRITE ( 14 ) omega_gen 506 WRITE ( 14 ) 'omega_gen_f ' 507 WRITE ( 14 ) omega_gen_f 508 WRITE ( 14 ) 'omega_gen_f_old ' 509 WRITE ( 14 ) omega_gen_f_old 510 WRITE ( 14 ) 'omega_gen_old ' 511 WRITE ( 14 ) omega_gen_old 512 WRITE ( 14 ) 'omega_rot ' 513 WRITE ( 14 ) omega_rot 514 WRITE ( 14 ) 'phi_yaw ' 515 WRITE ( 14 ) phi_yaw(:) 516 WRITE ( 14 ) 'pitch_add ' 517 WRITE ( 14 ) pitch_add 518 WRITE ( 14 ) 'pitch_add_old ' 519 WRITE ( 14 ) pitch_add_old 520 WRITE ( 14 ) 'torque_gen ' 521 WRITE ( 14 ) torque_gen 522 WRITE ( 14 ) 'torque_gen_old ' 523 WRITE ( 14 ) torque_gen_old 524 525 WRITE ( 14 ) '*** end wtm *** ' 526 527 END SUBROUTINE wtm_write_restart_data 528 529 530 !------------------------------------------------------------------------------! 531 ! Description: 532 ! ------------ 533 !> This routine reads the respective restart data. 534 !------------------------------------------------------------------------------! 535 SUBROUTINE wtm_read_restart_data 536 537 538 IMPLICIT NONE 539 540 CHARACTER (LEN=30) :: variable_chr !< dummy variable to read string 541 542 543 READ ( 13 ) variable_chr 544 DO WHILE ( TRIM( variable_chr ) /= '*** end wtm ***' ) 545 546 SELECT CASE ( TRIM( variable_chr ) ) 547 548 CASE ( 'omega_gen' ) 549 READ ( 13 ) omega_gen 550 CASE ( 'omega_gen_f' ) 551 READ ( 13 ) omega_gen_f 552 CASE ( 'omega_gen_old' ) 553 READ ( 13 ) omega_gen_old 554 CASE ( 'omega_gen_f_old' ) 555 READ ( 13 ) omega_gen_f_old 556 CASE ( 'omega_rot' ) 557 READ ( 13 ) omega_rot 558 CASE ( 'phi_yaw' ) 559 READ ( 13 ) phi_yaw 560 CASE ( 'pitch_add' ) 561 READ ( 13 ) pitch_add 562 CASE ( 'pitch_add_old' ) 563 READ ( 13 ) pitch_add_old 564 CASE ( 'torque_gen' ) 565 READ ( 13 ) torque_gen 566 CASE ( 'torque_gen_old' ) 567 READ ( 13 ) torque_gen_old 568 569 END SELECT 570 571 READ ( 13 ) variable_chr 572 573 ENDDO 574 575 END SUBROUTINE wtm_read_restart_data 576 577 578 !------------------------------------------------------------------------------! 579 ! Description: 580 ! ------------ 581 !> Check namelist parameter 582 !------------------------------------------------------------------------------! 474 583 SUBROUTINE wtm_check_parameters 475 584 … … 594 703 595 704 ! 596 !-- Allocation of the 1D arrays for speed pitch_control597 ALLOCATE( omega_gen(1:nturbines) )598 ALLOCATE( omega_gen_old(1:nturbines) )599 ALLOCATE( omega_gen_f(1:nturbines) )600 ALLOCATE( omega_gen_f_old(1:nturbines) )601 ALLOCATE( pitch_add_old(1:nturbines) )602 ALLOCATE( torque_gen(1:nturbines) )603 ALLOCATE( torque_gen_old(1:nturbines) )604 605 !606 705 !-- Allocation of the 1D arrays for yaw control 607 706 ALLOCATE( yawdir(1:nturbines) ) … … 672 771 thrust_rotor(:) = 0.0_wp 673 772 674 omega_gen(:) = 0.0_wp 675 omega_gen_old(:) = 0.0_wp 676 omega_gen_f(:) = 0.0_wp 677 omega_gen_f_old(:) = 0.0_wp 678 pitch_add_old(:) = 0.0_wp 679 torque_gen(:) = 0.0_wp 680 torque_gen_old(:) = 0.0_wp 681 773 IF ( TRIM( initializing_actions ) /= 'read_restart_data' ) THEN 774 omega_gen(:) = 0.0_wp 775 omega_gen_old(:) = 0.0_wp 776 omega_gen_f(:) = 0.0_wp 777 omega_gen_f_old(:) = 0.0_wp 778 pitch_add_old(:) = 0.0_wp 779 torque_gen(:) = 0.0_wp 780 torque_gen_old(:) = 0.0_wp 781 ENDIF 782 682 783 yawdir(:) = 0.0_wp 683 784 wdir(:) = 0.0_wp … … 776 877 CALL wtm_speed_control 777 878 879 IF ( TRIM( initializing_actions ) == 'read_restart_data' ) THEN 880 881 DO inot = 1, nturbines 882 883 IF ( nxl > i_hub(inot) ) THEN 884 torque_gen(inot) = 0.0_wp 885 omega_gen_f(inot) = 0.0_wp 886 omega_rot_l(inot) = 0.0_wp 887 ENDIF 888 889 IF ( nxr < i_hub(inot) ) THEN 890 torque_gen(inot) = 0.0_wp 891 omega_gen_f(inot) = 0.0_wp 892 omega_rot_l(inot) = 0.0_wp 893 ENDIF 894 895 IF ( nys > j_hub(inot) ) THEN 896 torque_gen(inot) = 0.0_wp 897 omega_gen_f(inot) = 0.0_wp 898 omega_rot_l(inot) = 0.0_wp 899 ENDIF 900 901 IF ( nyn < j_hub(inot) ) THEN 902 torque_gen(inot) = 0.0_wp 903 omega_gen_f(inot) = 0.0_wp 904 omega_rot_l(inot) = 0.0_wp 905 ENDIF 906 907 IF ( ( nxl <= i_hub(inot) ) .AND. ( nxr >= i_hub(inot) ) ) THEN 908 IF ( ( nys <= j_hub(inot) ) .AND. ( nyn >= j_hub(inot) ) ) THEN 909 910 omega_rot_l(inot) = omega_gen(inot) / gear_ratio 911 912 ENDIF 913 ENDIF 914 915 END DO 916 917 ENDIF 918 778 919 ENDIF 779 920 … … 814 955 ! 815 956 !-- Change yaw angle to rad: 816 phi_yaw(:) = phi_yaw(:) * pi / 180.0_wp 957 IF ( TRIM( initializing_actions ) /= 'read_restart_data' ) THEN 958 phi_yaw(:) = phi_yaw(:) * pi / 180.0_wp 959 ENDIF 817 960 818 961 … … 2366 2509 !-- At the first timestep the torque is set to its maximum to prevent 2367 2510 !-- an overspeeding of the rotor 2368 torque_gen_old(:) = max_torque_gen 2511 IF ( TRIM( initializing_actions ) /= 'read_restart_data' ) THEN 2512 torque_gen_old(:) = max_torque_gen 2513 ENDIF 2369 2514 2370 2515 END SUBROUTINE wtm_init_speed_control -
palm/trunk/SOURCE/write_var_list.f90
r2372 r2563 25 25 ! ----------------- 26 26 ! $Id$ 27 ! Function call wtm_write_restart_data was added and the end of the standard 28 ! parameter list is now marked with *** end default ***. The end of the whole 29 ! parameter list including module parameter is marked with *** end *** 30 ! 31 ! 2372 2017-08-25 12:37:32Z sward 27 32 ! y_shift added to vars, version no. increased 28 33 ! … … 233 238 ONLY: vnest_init 234 239 240 USE wind_turbine_model_mod, & 241 ONLY: wtm_write_restart_data 242 235 243 IMPLICIT NONE 236 244 … … 238 246 239 247 240 binary_version = '4. 2'248 binary_version = '4.3' 241 249 242 250 WRITE ( 14 ) binary_version … … 739 747 740 748 ! 741 !-- Set the end-of-file mark 742 WRITE ( 14 ) '*** end *** ' 749 !-- Set the end-of-file mark for default parameter 750 WRITE ( 14 ) '*** end default *** ' 751 752 ! 753 !-- If required, write restart data for wind turbine model. 754 IF ( wind_turbine ) CALL wtm_write_restart_data 755 743 756 ! 744 757 !-- If required, write restart data for virtual measurements. … … 749 762 IF ( synthetic_turbulence_generator ) CALL stg_write_restart_data 750 763 764 ! 765 !-- Set the end-of-file mark for default parameter 766 WRITE ( 14 ) '*** end *** ' 767 751 768 752 769 END SUBROUTINE write_var_list
Note: See TracChangeset
for help on using the changeset viewer.