Changeset 1240 for palm/trunk/SOURCE/init_3d_model.f90
- Timestamp:
- Oct 29, 2013 1:21:31 PM (11 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/init_3d_model.f90
r1239 r1240 23 23 ! Current revisions: 24 24 ! ------------------ 25 ! Overwrite initial profiles in case of nudging 26 ! Inititialize shf and qsws in case of large_scale_forcing 25 ! 27 26 ! 28 27 ! Former revisions: … … 261 260 USE indices 262 261 USE interfaces 263 USE ls_forcing_mod264 262 USE model_1d 265 263 USE netcdf_control … … 782 780 ELSEIF ( INDEX(initializing_actions, 'set_constant_profiles') /= 0 ) & 783 781 THEN 784 785 !786 !-- Overwrite initial profiles in case of nudging787 IF ( nudging ) THEN788 pt_init = ptnudge(:,1)789 u_init = unudge(:,1)790 v_init = vnudge(:,1)791 IF ( humidity .OR. passive_scalar ) THEN792 q_init = qnudge(:,1)793 ENDIF794 795 WRITE( message_string, * ) 'Initial profiles of u, v and ', &796 'scalars from NUDGING_DATA are used.'797 CALL message( 'init_3d_model', 'PA0370', 0, 0, 0, 6, 0 )798 ENDIF799 800 782 ! 801 783 !-- Use constructed initial profiles (velocity constant with height, … … 891 873 IF ( sloping_surface ) CALL init_slope 892 874 893 IF(myid==0) THEN894 PRINT*,"initial profiles from init_3d_model"895 DO k=nzb, nzt+1896 PRINT*, k, u(k,0,0), v(k,0,0), pt(k,0,0), q(k,0,0)897 ENDDO898 ENDIF899 900 875 ELSEIF ( INDEX(initializing_actions, 'by_user') /= 0 ) & 901 876 THEN … … 973 948 ELSE 974 949 shf = surface_heatflux 975 !976 !-- Initialize shf with data from external file LSF_DATA977 IF ( large_scale_forcing .AND. lsf_surf ) THEN978 CALL ls_forcing_surf ( simulated_time )979 ENDIF980 981 950 ! 982 951 !-- Over topography surface_heatflux is replaced by wall_heatflux(0)
Note: See TracChangeset
for help on using the changeset viewer.