Changeset 97 for palm/trunk/SOURCE/init_ocean.f90
- Timestamp:
- Jun 21, 2007 8:23:15 AM (17 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
palm/trunk/SOURCE/init_ocean.f90
r96 r97 40 40 hyp(nzt+1) = surface_pressure * 100.0 41 41 42 hyp(nzt) = hyp(nzt+1) + rho_surface * g * 0.5 * dzu(nzt+1)43 rho_ref 42 hyp(nzt) = hyp(nzt+1) + rho_surface * g * 0.5 * dzu(nzt+1) 43 rho_reference = rho_surface * 0.5 * dzu(nzt+1) 44 44 45 45 DO k = nzt-1, 0, -1 … … 48 48 pt_l = 0.5 * ( pt_init(k) + pt_init(k+1) ) 49 49 50 rho_l 50 rho_l = eqn_state_seawater_func( hyp(k+1), pt_l, sa_l ) 51 51 52 hyp(k) = hyp(k+1) + rho_l * g * dzu(k+1)53 rho_ref = rho_ref+ rho_l * dzu(k+1)52 hyp(k) = hyp(k+1) + rho_l * g * dzu(k+1) 53 rho_reference = rho_reference + rho_l * dzu(k+1) 54 54 55 55 ENDDO 56 56 57 rho_ref = rho_ref / ( zw(nzt) - zu(nzb) ) 58 print*, '*** rho_ref = ', rho_ref 57 rho_reference = rho_reference / ( zw(nzt) - zu(nzb) ) 58 59 ! 60 !-- Calculate the reference potential density 61 prho_reference = 0.0 62 DO k = 0, nzt 63 64 sa_l = 0.5 * ( sa_init(k) + sa_init(k+1) ) 65 pt_l = 0.5 * ( pt_init(k) + pt_init(k+1) ) 66 67 prho_reference = prho_reference + dzu(k+1) * & 68 eqn_state_seawater_func( hyp(0), pt_l, sa_l ) 69 70 ENDDO 71 72 prho_reference = prho_reference / ( zu(nzt) - zu(nzb) ) 59 73 60 74
Note: See TracChangeset
for help on using the changeset viewer.