Changeset 1353 for palm/trunk/SOURCE/surface_coupler.f90
- Timestamp:
- Apr 8, 2014 3:21:23 PM (10 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/surface_coupler.f90
r1325 r1353 20 20 ! Current revisions: 21 21 ! ------------------ 22 ! REAL constants provided with KIND-attribute 22 23 ! 23 24 ! Former revisions: … … 206 207 ! 207 208 !-- Send heat flux at bottom surface to the ocean 208 total_2d_a = 0.0 209 total_2d = 0.0 209 total_2d_a = 0.0_wp 210 total_2d = 0.0_wp 210 211 total_2d(nys:nyn,nxl:nxr) = shf(nys:nyn,nxl:nxr) 211 212 … … 216 217 !-- Send humidity flux at bottom surface to the ocean 217 218 IF ( humidity ) THEN 218 total_2d_a = 0.0 219 total_2d = 0.0 219 total_2d_a = 0.0_wp 220 total_2d = 0.0_wp 220 221 total_2d(nys:nyn,nxl:nxr) = qsws(nys:nyn,nxl:nxr) 221 222 … … 236 237 ! 237 238 !-- Send momentum flux (u) at bottom surface to the ocean 238 total_2d_a = 0.0 239 total_2d = 0.0 239 total_2d_a = 0.0_wp 240 total_2d = 0.0_wp 240 241 total_2d(nys:nyn,nxl:nxr) = usws(nys:nyn,nxl:nxr) 241 242 CALL MPI_REDUCE( total_2d, total_2d_a, ngp_a, MPI_REAL, MPI_SUM, 0, & … … 244 245 ! 245 246 !-- Send momentum flux (v) at bottom surface to the ocean 246 total_2d_a = 0.0 247 total_2d = 0.0 247 total_2d_a = 0.0_wp 248 total_2d = 0.0_wp 248 249 total_2d(nys:nyn,nxl:nxr) = vsws(nys:nyn,nxl:nxr) 249 250 CALL MPI_REDUCE( total_2d, total_2d_a, ngp_a, MPI_REAL, MPI_SUM, 0, & … … 338 339 ! 339 340 !-- Send surface temperature to atmosphere 340 total_2d_o = 0.0 341 total_2d = 0.0 341 total_2d_o = 0.0_wp 342 total_2d = 0.0_wp 342 343 total_2d(nys:nyn,nxl:nxr) = pt(nzt,nys:nyn,nxl:nxr) 343 344 … … 367 368 ! 368 369 !-- Send u to atmosphere 369 total_2d_o = 0.0 370 total_2d = 0.0 370 total_2d_o = 0.0_wp 371 total_2d = 0.0_wp 371 372 total_2d(nys:nyn,nxl:nxr) = u(nzt,nys:nyn,nxl:nxr) 372 373 CALL MPI_REDUCE( total_2d, total_2d_o, ngp_o, MPI_REAL, MPI_SUM, 0, & … … 375 376 ! 376 377 !-- Send v to atmosphere 377 total_2d_o = 0.0 378 total_2d = 0.0 378 total_2d_o = 0.0_wp 379 total_2d = 0.0_wp 379 380 total_2d(nys:nyn,nxl:nxr) = v(nzt,nys:nyn,nxl:nxr) 380 381 CALL MPI_REDUCE( total_2d, total_2d_o, ngp_o, MPI_REAL, MPI_SUM, 0, & … … 396 397 !-- following Steinhorn (1991), JPO 21, pp. 1681-1683: 397 398 !-- S'w' = -S * evaporation / ( rho_water * ( 1 - S ) ) 398 saswst = -1.0 * sa(nzt,:,:) * qswst_remote / &399 ( rho(nzt,:,:) * ( 1.0 - sa(nzt,:,:) ) )399 saswst = -1.0_wp * sa(nzt,:,:) * qswst_remote / & 400 ( rho(nzt,:,:) * ( 1.0_wp - sa(nzt,:,:) ) ) 400 401 ENDIF 401 402 … … 474 475 dny2 = 2 * ( dny / 2 ) 475 476 476 total_2d_a = 0.0 477 total_2d_a = 0.0_wp 477 478 ! 478 479 !-- Interpolation from ocean-grid-layer to atmosphere-grid-layer
Note: See TracChangeset
for help on using the changeset viewer.