Ignore:
Timestamp:
Apr 8, 2014 3:21:23 PM (10 years ago)
Author:
heinze
Message:

REAL constants provided with KIND-attribute

File:
1 edited

Legend:

Unmodified
Added
Removed
  • palm/trunk/SOURCE/init_ocean.f90

    r1323 r1353  
    2020! Current revisions:
    2121! -----------------
    22 !
     22! REAL constants provided with KIND-attribute
    2323!
    2424! Former revisions:
     
    8585!
    8686!-- Set water density near the ocean surface
    87     rho_surface = 1027.62
     87    rho_surface = 1027.62_wp
    8888
    8989!
     
    9393    hyp(nzt+1) = surface_pressure * 100.0_wp
    9494
    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)
    9696    rho_init(nzt) = rho_surface
    9797
     
    106106    DO  n = 1, 5
    107107
    108        rho_reference = rho_surface * 0.5 * dzu(nzt+1)
     108       rho_reference = rho_surface * 0.5_wp * dzu(nzt+1)
    109109
    110110       DO  k = nzt-1, 0, -1
    111111
    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) )
    114114
    115115          rho_init(k) = eqn_state_seawater_func( hyp(k), pt_l, sa_l )
     
    122122
    123123       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 ) ) * &
    125125                              dzu(k+1)
    126126       ENDDO
     
    130130!
    131131!-- Calculate the reference potential density
    132     prho_reference = 0.0
     132    prho_reference = 0.0_wp
    133133    DO  k = 0, nzt
    134134
    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) )
    137137
    138138       prho_reference = prho_reference + dzu(k+1) * &
Note: See TracChangeset for help on using the changeset viewer.