Changeset 1353 for palm/trunk/SOURCE/check_open.f90
- Timestamp:
- Apr 8, 2014 3:21:23 PM (10 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/check_open.f90
r1329 r1353 20 20 ! Current revisions: 21 21 ! ----------------- 22 ! 22 ! REAL constants provided with KIND-attribute, 23 ! declaration for unused variables xkoor, ykoor, zkoor removed 23 24 ! 24 25 ! Former revisions: … … 147 148 LOGICAL :: twoya = .TRUE. !: 148 149 149 REAL(wp) :: ansx = -999.999 !:150 REAL(wp) :: ansy = -999.999 !:151 REAL(wp) :: gwid = 0.1 !:152 REAL(wp) :: rlegfak = 1.5 !:153 REAL(wp) :: sizex = 250.0 !:154 REAL(wp) :: sizey = 40.0 !:155 REAL(wp) :: texfac = 1.5 !:150 REAL(wp) :: ansx = -999.999_wp !: 151 REAL(wp) :: ansy = -999.999_wp !: 152 REAL(wp) :: gwid = 0.1_wp !: 153 REAL(wp) :: rlegfak = 1.5_wp !: 154 REAL(wp) :: sizex = 250.0_wp !: 155 REAL(wp) :: sizey = 40.0_wp !: 156 REAL(wp) :: texfac = 1.5_wp !: 156 157 157 158 REAL(wp), DIMENSION(:), ALLOCATABLE :: eta !: 158 159 REAL(wp), DIMENSION(:), ALLOCATABLE :: ho !: 159 160 REAL(wp), DIMENSION(:), ALLOCATABLE :: hu !: 160 161 REAL(sp), DIMENSION(:), ALLOCATABLE :: xkoor !: 162 REAL(sp), DIMENSION(:), ALLOCATABLE :: ykoor !: 163 REAL(sp), DIMENSION(:), ALLOCATABLE :: zkoor !: 161 164 162 165 163 … … 363 361 !-- Determine and write ISO2D coordiante header 364 362 ALLOCATE( eta(0:ny+1), ho(0:nx+1), hu(0:nx+1) ) 365 hu = 0.0 363 hu = 0.0_wp 366 364 ho = (ny+1) * dy 367 365 DO i = 1, ny 368 eta(i) = REAL( i ) / ( ny + 1.0 )366 eta(i) = REAL( i ) / ( ny + 1.0_wp ) 369 367 ENDDO 370 eta(0) = 0.0 371 eta(ny+1) = 1.0 368 eta(0) = 0.0_wp 369 eta(ny+1) = 1.0_wp 372 370 373 371 WRITE (21) dx,eta,hu,ho … … 396 394 !-- Determine and write ISO2D coordinate header 397 395 ALLOCATE( eta(0:nz+1), ho(0:nx+1), hu(0:nx+1) ) 398 hu = 0.0 396 hu = 0.0_wp 399 397 ho = zu(nz+1) 400 398 DO i = 1, nz 401 399 eta(i) = REAL( zu(i) ) / zu(nz+1) 402 400 ENDDO 403 eta(0) = 0.0 404 eta(nz+1) = 1.0 401 eta(0) = 0.0_wp 402 eta(nz+1) = 1.0_wp 405 403 406 404 WRITE (22) dx,eta,hu,ho … … 429 427 !-- Determine and write ISO2D coordiante header 430 428 ALLOCATE( eta(0:nz+1), ho(0:ny+1), hu(0:ny+1) ) 431 hu = 0.0 429 hu = 0.0_wp 432 430 ho = zu(nz+1) 433 431 DO i = 1, nz 434 432 eta(i) = REAL( zu(i) ) / zu(nz+1) 435 433 ENDDO 436 eta(0) = 0.0 437 eta(nz+1) = 1.0 434 eta(0) = 0.0_wp 435 eta(nz+1) = 1.0_wp 438 436 439 437 WRITE (23) dx,eta,hu,ho
Note: See TracChangeset
for help on using the changeset viewer.