Changeset 1353 for palm/trunk/SOURCE/init_ocean.f90
- Timestamp:
- Apr 8, 2014 3:21:23 PM (10 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/init_ocean.f90
r1323 r1353 20 20 ! Current revisions: 21 21 ! ----------------- 22 ! 22 ! REAL constants provided with KIND-attribute 23 23 ! 24 24 ! Former revisions: … … 85 85 ! 86 86 !-- Set water density near the ocean surface 87 rho_surface = 1027.62 87 rho_surface = 1027.62_wp 88 88 89 89 ! … … 93 93 hyp(nzt+1) = surface_pressure * 100.0_wp 94 94 95 hyp(nzt) = hyp(nzt+1) + rho_surface * g * 0.5 * dzu(nzt+1)95 hyp(nzt) = hyp(nzt+1) + rho_surface * g * 0.5_wp * dzu(nzt+1) 96 96 rho_init(nzt) = rho_surface 97 97 … … 106 106 DO n = 1, 5 107 107 108 rho_reference = rho_surface * 0.5 * dzu(nzt+1)108 rho_reference = rho_surface * 0.5_wp * dzu(nzt+1) 109 109 110 110 DO k = nzt-1, 0, -1 111 111 112 sa_l = 0.5 * ( sa_init(k) + sa_init(k+1) )113 pt_l = 0.5 * ( pt_init(k) + pt_init(k+1) )112 sa_l = 0.5_wp * ( sa_init(k) + sa_init(k+1) ) 113 pt_l = 0.5_wp * ( pt_init(k) + pt_init(k+1) ) 114 114 115 115 rho_init(k) = eqn_state_seawater_func( hyp(k), pt_l, sa_l ) … … 122 122 123 123 DO k = nzt-1, 0, -1 124 hyp(k) = hyp(k+1) + g * 0.5 * ( rho_init(k) + rho_init(k+1 ) ) * &124 hyp(k) = hyp(k+1) + g * 0.5_wp * ( rho_init(k) + rho_init(k+1 ) ) * & 125 125 dzu(k+1) 126 126 ENDDO … … 130 130 ! 131 131 !-- Calculate the reference potential density 132 prho_reference = 0.0 132 prho_reference = 0.0_wp 133 133 DO k = 0, nzt 134 134 135 sa_l = 0.5 * ( sa_init(k) + sa_init(k+1) )136 pt_l = 0.5 * ( pt_init(k) + pt_init(k+1) )135 sa_l = 0.5_wp * ( sa_init(k) + sa_init(k+1) ) 136 pt_l = 0.5_wp * ( pt_init(k) + pt_init(k+1) ) 137 137 138 138 prho_reference = prho_reference + dzu(k+1) * &
Note: See TracChangeset
for help on using the changeset viewer.