Ignore:
Timestamp:
Oct 21, 2016 3:11:58 PM (5 years ago)
Author:
knoop
Message:

Renamed variable rho to rho_ocean, rho_init to rho_ocean_init and rho_av to rho_ocean_av

File:
1 edited

Legend:

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

    r2001 r2031  
    2020! Current revisions:
    2121! ------------------
    22 !
     22! renamed variable rho to rho_ocean
    2323!
    2424! Former revisions:
     
    8080
    8181    USE arrays_3d,                                                             &
    82         ONLY:  pt, shf, qsws, qswst_remote, rho, sa, saswst, total_2d_a,       &
     82        ONLY:  pt, shf, qsws, qswst_remote, rho_ocean, sa, saswst, total_2d_a,       &
    8383               total_2d_o, tswst, u, usws, uswst, v, vsws, vswst
    8484
     
    421421!--       S'w' = -S * evaporation / ( rho_water * ( 1 - S ) )
    422422          saswst = -1.0_wp * sa(nzt,:,:) * 0.001_wp * qswst_remote /  &
    423                     ( rho(nzt,:,:) * ( 1.0_wp - sa(nzt,:,:) * 0.001_wp ) )
     423                    ( rho_ocean(nzt,:,:) * ( 1.0_wp - sa(nzt,:,:) * 0.001_wp ) )
    424424       ENDIF
    425425
     
    428428!--    (constants are the specific heat capacities for air and water)
    429429!--    now tswst is the ocean top heat flux
    430        tswst = tswst / rho(nzt,:,:) * cp / cpw
     430       tswst = tswst / rho_ocean(nzt,:,:) * cp / cpw
    431431
    432432!
    433433!--    Adjust the momentum fluxes with respect to ocean density
    434        uswst = uswst / rho(nzt,:,:)
    435        vswst = vswst / rho(nzt,:,:)
     434       uswst = uswst / rho_ocean(nzt,:,:)
     435       vswst = vswst / rho_ocean(nzt,:,:)
    436436
    437437    ENDIF
Note: See TracChangeset for help on using the changeset viewer.