Changeset 2031


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

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

Location:
palm/trunk/SOURCE
Files:
23 edited

Legend:

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

    r2012 r2031  
    2020! Current revisions:
    2121! -----------------
    22 !
     22! renamed variable rho to rho_ocean and rho_av to rho_ocean_av
    2323!
    2424! Former revisions:
     
    318318             ENDDO
    319319
    320           CASE ( 'rho' )
    321              DO  i = nxlg, nxrg
    322                 DO  j = nysg, nyng
    323                    DO  k = nzb, nzt+1
    324                       rho_av(k,j,i) = rho_av(k,j,i) / REAL( average_count_3d, KIND=wp )
     320          CASE ( 'rho_ocean' )
     321             DO  i = nxlg, nxrg
     322                DO  j = nysg, nyng
     323                   DO  k = nzb, nzt+1
     324                      rho_ocean_av(k,j,i) = rho_ocean_av(k,j,i) / REAL( average_count_3d, KIND=wp )
    325325                   ENDDO
    326326                ENDDO
  • palm/trunk/SOURCE/check_parameters.f90

    r2027 r2031  
    2020! Current revisions:
    2121! -----------------
    22 !
     22! renamed variable rho to rho_ocean
    2323!
    2424! Former revisions:
     
    299299! 1179 2013-06-14 05:57:58Z raasch
    300300! checks and settings of buoyancy parameters and switches revised,
    301 ! initial profile for rho added to hom (id=77)
     301! initial profile for rho_ocean added to hom (id=77)
    302302!
    303303! 1174 2013-05-31 10:28:08Z gryschka
     
    25282528             hom(:,2,63,:) = SPREAD( zu, 2, statistic_regions+1 )
    25292529
    2530           CASE ( 'rho' )
     2530          CASE ( 'rho_ocean' )
    25312531             IF (  .NOT.  ocean ) THEN
    25322532                message_string = 'data_output_pr = ' //                        &
     
    30053005             unit = 'kg/kg'
    30063006
    3007           CASE ( 'rho' )
     3007          CASE ( 'rho_ocean' )
    30083008             IF (  .NOT.  ocean )  THEN
    30093009                message_string = 'output of "' // TRIM( var ) // '" requi' //  &
  • palm/trunk/SOURCE/data_output_2d.f90

    r2001 r2031  
    2020! Current revisions:
    2121! -----------------
    22 !
     22! renamed variable rho to rho_ocean and rho_av to rho_ocean_av
    2323!
    2424! Former revisions:
     
    157157    USE arrays_3d,                                                             &
    158158        ONLY:  dzw, e, nr, ol, p, pt, precipitation_amount, precipitation_rate,&
    159                prr,q, qc, ql, ql_c, ql_v, ql_vp, qr, qsws, rho, s, sa, shf,    &
     159               prr,q, qc, ql, ql_c, ql_v, ql_vp, qr, qsws, rho_ocean, s, sa, shf,    &
    160160               ssws, tend, ts, u, us, v, vpt, w, z0, z0h, z0q, zu, zw
    161161       
     
    769769             CASE ( 'rho_xy', 'rho_xz', 'rho_yz' )
    770770                IF ( av == 0 )  THEN
    771                    to_be_resorted => rho
    772                 ELSE
    773                    to_be_resorted => rho_av
     771                   to_be_resorted => rho_ocean
     772                ELSE
     773                   to_be_resorted => rho_ocean_av
    774774                ENDIF
    775775
  • palm/trunk/SOURCE/data_output_3d.f90

    r2012 r2031  
    2020! Current revisions:
    2121! ------------------
    22 !
     22! renamed variable rho to rho_ocean and rho_av to rho_ocean_av
    2323!
    2424! Former revisions:
     
    146146
    147147    USE arrays_3d,                                                             &
    148         ONLY:  e, nr, p, pt, prr, q, qc, ql, ql_c, ql_v, qr, rho, s, sa, tend, &
     148        ONLY:  e, nr, p, pt, prr, q, qc, ql, ql_c, ql_v, qr, rho_ocean, s, sa, tend, &
    149149               u, v, vpt, w
    150150       
     
    531531             ENDIF
    532532
    533           CASE ( 'rho' )
    534              IF ( av == 0 )  THEN
    535                 to_be_resorted => rho
    536              ELSE
    537                 to_be_resorted => rho_av
     533          CASE ( 'rho_ocean' )
     534             IF ( av == 0 )  THEN
     535                to_be_resorted => rho_ocean
     536             ELSE
     537                to_be_resorted => rho_ocean_av
    538538             ENDIF
    539539
  • palm/trunk/SOURCE/data_output_mask.f90

    r2001 r2031  
    2020! Current revisions:
    2121! -----------------
    22 !
     22! renamed variable rho to rho_ocean and rho_av to rho_ocean_av
    2323!
    2424! Former revisions:
     
    102102#if defined( __netcdf )
    103103    USE arrays_3d,                                                             &
    104         ONLY:  e, nr, p, pt, q, qc, ql, ql_c, ql_v, qr, rho, s, sa, tend, u,   &
     104        ONLY:  e, nr, p, pt, q, qc, ql, ql_c, ql_v, qr, rho_ocean, s, sa, tend, u,   &
    105105               v, vpt, w
    106106   
    107107    USE averaging,                                                             &
    108108        ONLY:  e_av, lpt_av, nr_av, p_av, pc_av, pr_av, pt_av, q_av, qc_av,    &
    109                ql_av, ql_c_av, ql_v_av, ql_vp_av, qv_av, qr_av, rho_av, s_av,  &
     109               ql_av, ql_c_av, ql_v_av, ql_vp_av, qv_av, qr_av, rho_ocean_av, s_av,  &
    110110               sa_av, u_av, v_av, vpt_av, w_av
    111111   
     
    436436             ENDIF
    437437
    438           CASE ( 'rho' )
    439              IF ( av == 0 )  THEN
    440                 to_be_resorted => rho
    441              ELSE
    442                 to_be_resorted => rho_av
     438          CASE ( 'rho_ocean' )
     439             IF ( av == 0 )  THEN
     440                to_be_resorted => rho_ocean
     441             ELSE
     442                to_be_resorted => rho_ocean_av
    443443             ENDIF
    444444
  • palm/trunk/SOURCE/diffusivities.f90

    r2001 r2031  
    2020! Current revisions:
    2121! -----------------
    22 !
     22! renamed variable rho to rho_ocean
    2323!
    2424! Former revisions:
     
    160160!
    161161!--             Determine the mixing length
    162                 dvar_dz = atmos_ocean_sign * &  ! inverse effect of pt/rho gradient
     162                dvar_dz = atmos_ocean_sign * &  ! inverse effect of pt/rho_ocean gradient
    163163                          ( var(k+1,j,i) - var(k-1,j,i) ) * dd2zu(k)
    164164                IF ( dvar_dz > 0.0_wp ) THEN
  • palm/trunk/SOURCE/eqn_state_seawater.f90

    r2001 r2031  
    2020! Current revisions:
    2121! -----------------
    22 !
     22! renamed variable rho to rho_ocean
    2323!
    2424! Former revisions:
     
    115115
    116116       USE arrays_3d,                                                          &
    117            ONLY:  hyp, prho, pt_p, rho, sa_p
     117           ONLY:  hyp, prho, pt_p, rho_ocean, sa_p
    118118       USE indices,                                                            &
    119119           ONLY:  nxl, nxr, nyn, nys, nzb_s_inner, nzt
     
    181181!
    182182!--             In-situ density
    183                 rho(k,j,i) = pnom / pden
     183                rho_ocean(k,j,i) = pnom / pden
    184184
    185185             ENDDO
     
    188188             prho(nzt+1,j,i)            = prho(nzt,j,i)
    189189             prho(nzb_s_inner(j,i),j,i) = prho(nzb_s_inner(j,i)+1,j,i)
    190              rho(nzt+1,j,i)             = rho(nzt,j,i)
    191              rho(nzb_s_inner(j,i),j,i)  = rho(nzb_s_inner(j,i)+1,j,i)
     190             rho_ocean(nzt+1,j,i)             = rho_ocean(nzt,j,i)
     191             rho_ocean(nzb_s_inner(j,i),j,i)  = rho_ocean(nzb_s_inner(j,i)+1,j,i)
    192192
    193193          ENDDO
     
    205205
    206206       USE arrays_3d,                                                          &
    207            ONLY:  hyp, prho, pt_p, rho, sa_p
     207           ONLY:  hyp, prho, pt_p, rho_ocean, sa_p
    208208           
    209209       USE indices,                                                            &
     
    255255!
    256256!--       In-situ density
    257           rho(k,j,i) = pnom / pden
     257          rho_ocean(k,j,i) = pnom / pden
    258258
    259259
     
    264264       prho(nzt+1,j,i)            = prho(nzt,j,i)
    265265       prho(nzb_s_inner(j,i),j,i) = prho(nzb_s_inner(j,i)+1,j,i)
    266        rho(nzt+1,j,i)             = rho(nzt,j,i)
    267        rho(nzb_s_inner(j,i),j,i)  = rho(nzb_s_inner(j,i)+1,j,i)
     266       rho_ocean(nzt+1,j,i)             = rho_ocean(nzt,j,i)
     267       rho_ocean(nzb_s_inner(j,i),j,i)  = rho_ocean(nzb_s_inner(j,i)+1,j,i)
    268268
    269269    END SUBROUTINE eqn_state_seawater_ij
  • palm/trunk/SOURCE/flow_statistics.f90

    r2027 r2031  
    2020! Current revisions:
    2121! -----------------
    22 !
     22! renamed variable rho to rho_ocean
    2323!
    2424! Former revisions:
     
    212212    USE arrays_3d,                                                             &
    213213        ONLY:  ddzu, ddzw, e, hyp, km, kh, nr, ol, p, prho, prr, pt, q, qc, ql,&
    214                qr, qs, qsws, qswst, rho, s, sa, ss, ssws, sswst, saswsb,       &
     214               qr, qs, qsws, qswst, rho_ocean, s, sa, ss, ssws, sswst, saswsb,       &
    215215               saswst, shf, td_lsa_lpt, td_lsa_q, td_sub_lpt, td_sub_q,        &
    216216               time_vert, ts, tswst, u, ug, us, usws, uswst, vsws, v, vg, vpt, &
     
    926926                                        rmask(j,i,sr)
    927927                   ENDIF
    928                    sums_l(k,64,tn) = sums_l(k,64,tn) + rho(k,j,i) *            &
     928                   sums_l(k,64,tn) = sums_l(k,64,tn) + rho_ocean(k,j,i) *            &
    929929                                                       rmask(j,i,sr)
    930930                   sums_l(k,71,tn) = sums_l(k,71,tn) + prho(k,j,i) *           &
     
    14001400       hom(:,1,55,sr) = sums(:,55)     ! w*u*u*/dz
    14011401       hom(:,1,56,sr) = sums(:,56)     ! w*p*/dz
    1402        hom(:,1,57,sr) = sums(:,57)     ! ( w"e + w"p"/rho )/dz
     1402       hom(:,1,57,sr) = sums(:,57)     ! ( w"e + w"p"/rho_ocean )/dz
    14031403       hom(:,1,58,sr) = sums(:,58)     ! u"pt"
    14041404       hom(:,1,59,sr) = sums(:,59)     ! u*pt*
     
    14071407       hom(:,1,62,sr) = sums(:,62)     ! v*pt*
    14081408       hom(:,1,63,sr) = sums(:,61) + sums(:,62)    ! vpt_t
    1409        hom(:,1,64,sr) = sums(:,64)     ! rho
     1409       hom(:,1,64,sr) = sums(:,64)     ! rho_ocean
    14101410       hom(:,1,65,sr) = sums(:,65)     ! w"sa"
    14111411       hom(:,1,66,sr) = sums(:,66)     ! w*sa*
    14121412       hom(:,1,67,sr) = sums(:,65) + sums(:,66)    ! wsa
    14131413       hom(:,1,68,sr) = sums(:,68)     ! w*p*
    1414        hom(:,1,69,sr) = sums(:,69)     ! w"e + w"p"/rho
     1414       hom(:,1,69,sr) = sums(:,69)     ! w"e + w"p"/rho_ocean
    14151415       hom(:,1,70,sr) = sums(:,70)     ! q*2
    14161416       hom(:,1,71,sr) = sums(:,71)     ! prho
     
    17061706    USE arrays_3d,                                                             &
    17071707        ONLY:  ddzu, ddzw, e, hyp, km, kh, nr, p, prho, pt, q, qc, ql, qr, qs, &
    1708                qsws, qswst, rho, s, sa, saswsb, saswst, shf, ss, ssws, sswst,  &
     1708               qsws, qswst, rho_ocean, s, sa, saswsb, saswst, shf, ss, ssws, sswst,  &
    17091709               td_lsa_lpt, td_lsa_q, td_sub_lpt, td_sub_q, time_vert, ts,      &
    17101710               tswst, u, ug, us, usws, uswst, vsws, v, vg, vpt, vswst, w,      &
     
    27862786          ENDIF
    27872787
    2788           !$acc parallel loop gang present( rflags_invers, rho, prho, rmask, sums_l ) create( s1, s2 )
     2788          !$acc parallel loop gang present( rflags_invers, rho_ocean, prho, rmask, sums_l ) create( s1, s2 )
    27892789          DO  k = nzb, nzt_diff
    27902790             s1 = 0
     
    27932793             DO  i = nxl, nxr
    27942794                DO  j = nys, nyn
    2795                    s1 = s1 + rho(k,j,i)  * rmask(j,i,sr) * rflags_invers(j,i,k+1)
     2795                   s1 = s1 + rho_ocean(k,j,i)  * rmask(j,i,sr) * rflags_invers(j,i,k+1)
    27962796                   s2 = s2 + prho(k,j,i) * rmask(j,i,sr) * rflags_invers(j,i,k+1)
    27972797                ENDDO
     
    34243424       hom(:,1,55,sr) = sums(:,55)     ! w*u*u*/dz
    34253425       hom(:,1,56,sr) = sums(:,56)     ! w*p*/dz
    3426        hom(:,1,57,sr) = sums(:,57)     ! ( w"e + w"p"/rho )/dz
     3426       hom(:,1,57,sr) = sums(:,57)     ! ( w"e + w"p"/rho_ocean )/dz
    34273427       hom(:,1,58,sr) = sums(:,58)     ! u"pt"
    34283428       hom(:,1,59,sr) = sums(:,59)     ! u*pt*
     
    34313431       hom(:,1,62,sr) = sums(:,62)     ! v*pt*
    34323432       hom(:,1,63,sr) = sums(:,61) + sums(:,62)    ! vpt_t
    3433        hom(:,1,64,sr) = sums(:,64)     ! rho
     3433       hom(:,1,64,sr) = sums(:,64)     ! rho_ocean
    34343434       hom(:,1,65,sr) = sums(:,65)     ! w"sa"
    34353435       hom(:,1,66,sr) = sums(:,66)     ! w*sa*
    34363436       hom(:,1,67,sr) = sums(:,65) + sums(:,66)    ! wsa
    34373437       hom(:,1,68,sr) = sums(:,68)     ! w*p*
    3438        hom(:,1,69,sr) = sums(:,69)     ! w"e + w"p"/rho
     3438       hom(:,1,69,sr) = sums(:,69)     ! w"e + w"p"/rho_ocean
    34393439       hom(:,1,70,sr) = sums(:,70)     ! q*2
    34403440       hom(:,1,71,sr) = sums(:,71)     ! prho
  • palm/trunk/SOURCE/init_3d_model.f90

    r2012 r2031  
    2020! Current revisions:
    2121! ------------------
    22 !
     22! renamed variable rho to rho_ocean
    2323!
    2424! Former revisions:
     
    623623#if defined( __nopointer )
    624624       ALLOCATE( prho(nzb:nzt+1,nysg:nyng,nxlg:nxrg),                          &
    625                  rho(nzb:nzt+1,nysg:nyng,nxlg:nxrg),                           &
     625                 rho_ocean(nzb:nzt+1,nysg:nyng,nxlg:nxrg),                           &
    626626                 sa(nzb:nzt+1,nysg:nyng,nxlg:nxrg),                            &
    627627                 sa_p(nzb:nzt+1,nysg:nyng,nxlg:nxrg),                          &
     
    634634                 sa_3(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
    635635       prho => prho_1
    636        rho  => rho_1  ! routines calc_mean_profile and diffusion_e require
     636       rho_ocean  => rho_1  ! routines calc_mean_profile and diffusion_e require
    637637                      ! density to be apointer
    638638#endif
  • palm/trunk/SOURCE/init_masks.f90

    r2001 r2031  
    2020! Current revisions:
    2121! -----------------
    22 !
     22! renamed variable rho to rho_ocean
    2323!
    2424! Former revisions:
     
    350350                unit = 'kg/kg'
    351351
    352              CASE ( 'rho' )
     352             CASE ( 'rho_ocean' )
    353353                IF ( .NOT. ocean )  THEN
    354354                   WRITE ( message_string, * ) 'output of "', TRIM( var ), &
  • palm/trunk/SOURCE/init_ocean.f90

    r2001 r2031  
    2020! Current revisions:
    2121! -----------------
    22 !
     22! renamed variable rho_init to rho_ocean_init
    2323!
    2424! Former revisions:
     
    9090    REAL(wp)     ::  sa_l !<
    9191
    92     REAL(wp), DIMENSION(nzb:nzt+1) ::  rho_init !<
     92    REAL(wp), DIMENSION(nzb:nzt+1) ::  rho_ocean_init !<
    9393
    9494    ALLOCATE( hyp(nzb:nzt+1) )
     
    105105
    106106    hyp(nzt)      = hyp(nzt+1) + rho_surface * g * 0.5_wp * dzu(nzt+1)
    107     rho_init(nzt) = rho_surface
     107    rho_ocean_init(nzt) = rho_surface
    108108
    109109    DO  k = nzt-1, 1, -1
     
    124124          pt_l = 0.5_wp * ( pt_init(k) + pt_init(k+1) )
    125125
    126           rho_init(k) = eqn_state_seawater_func( hyp(k), pt_l, sa_l )
     126          rho_ocean_init(k) = eqn_state_seawater_func( hyp(k), pt_l, sa_l )
    127127
    128           rho_reference = rho_reference + rho_init(k) * dzu(k+1)
     128          rho_reference = rho_reference + rho_ocean_init(k) * dzu(k+1)
    129129
    130130       ENDDO
     
    133133
    134134       DO  k = nzt-1, 0, -1
    135           hyp(k) = hyp(k+1) + g * 0.5_wp * ( rho_init(k) + rho_init(k+1 ) ) * &
     135          hyp(k) = hyp(k+1) + g * 0.5_wp * ( rho_ocean_init(k) + rho_ocean_init(k+1 ) ) * &
    136136                              dzu(k+1)
    137137       ENDDO
     
    161161!
    162162!-- Store initial density profile
    163     hom(:,1,77,:)  = SPREAD( rho_init(:), 2, statistic_regions+1 )
     163    hom(:,1,77,:)  = SPREAD( rho_ocean_init(:), 2, statistic_regions+1 )
    164164
    165165!
     
    168168       ref_state(:) = prho_reference
    169169    ELSE
    170        ref_state(:) = rho_init(:)
     170       ref_state(:) = rho_ocean_init(:)
    171171    ENDIF
    172172
  • palm/trunk/SOURCE/land_surface_model_mod.f90

    r2001 r2031  
    2020! Current revisions:
    2121! -----------------
    22 !
     22! renamed variable rho to rho_ocean
    2323!
    2424! Former revisions:
     
    257257                residual_moisture = 9999999.9_wp,       & !< NAMELIST m_res
    258258                rho_cp,                                 & !< rho_surface * cp
    259                 rho_lv,                                 & !< rho * l_v
     259                rho_lv,                                 & !< rho_ocean * l_v
    260260                rd_d_rv,                                & !< r_d / r_v
    261261                saturation_moisture = 9999999.9_wp,     & !< NAMELIST m_sat
  • palm/trunk/SOURCE/microphysics_mod.f90

    r2001 r2031  
    2020! Current revisions:
    2121! ------------------
    22 !
     22! renamed variable rho to rho_ocean
    2323!
    2424! Former revisions:
     
    8282!
    8383! 1241 2013-10-30 11:36:58Z heinze
    84 ! hyp and rho have to be calculated at each time step if data from external
     84! hyp and rho_ocean have to be calculated at each time step if data from external
    8585! file LSF_DATA are used
    8686!
  • palm/trunk/SOURCE/modules.f90

    r2012 r2031  
    2020! Current revisions:
    2121! ------------------
    22 !
     22! renamed variable rho to rho_ocean and rho_av to rho_ocean_av
    2323!
    2424! Former revisions:
     
    473473    REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET ::                         &
    474474          e, e_p, nr, nr_p, p, prho, pt, pt_p, q, q_p, qc, ql, ql_c, ql_v,     &
    475           ql_vp, qr, qr_p, rho, s, s_p, sa, sa_p, te_m, tnr_m, tpt_m, tq_m,    &
     475          ql_vp, qr, qr_p, rho_ocean, s, s_p, sa, sa_p, te_m, tnr_m, tpt_m, tq_m,    &
    476476          tqr_m, ts_m, tsa_m, tu_m, tv_m, tw_m, u, u_p, v, v_p, vpt, w, w_p
    477477#else
     
    484484    REAL(wp), DIMENSION(:,:,:), POINTER ::                                     &
    485485          e, e_p, nr, nr_p, prho, pt, pt_p, q, q_p, qc, ql, ql_c, qr, qr_p,    &
    486           rho, s, s_p, sa, sa_p, te_m, tnr_m, tpt_m, tq_m, tqr_m, ts_m,        &
     486          rho_ocean, s, s_p, sa, sa_p, te_m, tnr_m, tpt_m, tq_m, tqr_m, ts_m,        &
    487487          tsa_m, tu_m, tv_m, tw_m, u, u_p, v, v_p, vpt, w, w_p
    488488#endif
     
    519519    REAL(wp), DIMENSION(:,:,:), ALLOCATABLE, TARGET ::                         &
    520520          e_av, lpt_av, nr_av, p_av, pc_av, pr_av, prr_av, pt_av, q_av, qc_av, &
    521           ql_av, ql_c_av, ql_v_av, ql_vp_av, qr_av, qv_av, rho_av, s_av, sa_av,&
     521          ql_av, ql_c_av, ql_v_av, ql_vp_av, qr_av, qv_av, rho_ocean_av, s_av, sa_av,&
    522522          u_av, v_av, vpt_av, w_av
    523523 
  • palm/trunk/SOURCE/netcdf_interface_mod.f90

    r2012 r2031  
    2020! Current revisions:
    2121! ------------------
    22 !
     22! renamed variable rho to rho_ocean
    2323!
    2424! Former revisions:
     
    778778                CASE ( 'e', 'lpt', 'nr', 'p', 'pc', 'pr', 'prr', 'pt', 'q',    &
    779779                       'qc', 'ql', 'ql_c', 'ql_v', 'ql_vp', 'qr', 'qv',        &
    780                        'rho', 's', 'sa', 'vpt'  )
     780                       'rho_ocean', 's', 'sa', 'vpt'  )
    781781
    782782                   grid_x = 'x'
     
    12991299!--             Most variables are defined on the scalar grid
    13001300                CASE ( 'e', 'lpt', 'nr', 'p', 'pc', 'pr', 'prr', 'pt', 'q',    &
    1301                        'qc', 'ql', 'ql_c', 'ql_v', 'ql_vp', 'qr', 'qv', 'rho', &
     1301                       'qc', 'ql', 'ql_c', 'ql_v', 'ql_vp', 'qr', 'qv', 'rho_ocean', &
    13021302                       's', 'sa', 'vpt' )
    13031303
     
    44034403!--             Most variables are defined on the zu levels
    44044404                CASE ( 'e', 'lpt', 'nr', 'p', 'pc', 'pr', 'prr', 'pt', 'q',   &
    4405                        'qc', 'ql', 'ql_c', 'ql_v', 'ql_vp', 'qr', 'qv', 'rho',&
     4405                       'qc', 'ql', 'ql_c', 'ql_v', 'ql_vp', 'qr', 'qv', 'rho_ocean',&
    44064406                       's', 'sa', 'u', 'v', 'vpt' )
    44074407
  • palm/trunk/SOURCE/production_e.f90

    r2001 r2031  
    2020! Current revisions:
    2121! -----------------
    22 !
     22! renamed variable rho to rho_ocean
    2323!
    2424! Former revisions:
     
    136136
    137137       USE arrays_3d,                                                          &
    138            ONLY:  ddzw, dd2zu, kh, km, pt, q, ql, qsws, qswst, rho, shf,       &
     138           ONLY:  ddzw, dd2zu, kh, km, pt, q, ql, qsws, qswst, rho_ocean, shf,       &
    139139                  tend, tswst, u, v, vpt, w
    140140
     
    537537                            tend(k,j,i) = tend(k,j,i) +                     &
    538538                                          kh(k,j,i) * g / rho_reference *   &
    539                                           ( rho(k+1,j,i) - rho(k-1,j,i) ) * &
     539                                          ( rho_ocean(k+1,j,i) - rho_ocean(k-1,j,i) ) * &
    540540                                          dd2zu(k)
    541541                         ENDDO
     
    576576                         DO  k = nzb_s_inner(j,i)+1, nzt
    577577                            tend(k,j,i) = tend(k,j,i) +                     &
    578                                           kh(k,j,i) * g / rho(k,j,i) *      &
    579                                           ( rho(k+1,j,i) - rho(k-1,j,i) ) * &
     578                                          kh(k,j,i) * g / rho_ocean(k,j,i) *      &
     579                                          ( rho_ocean(k+1,j,i) - rho_ocean(k-1,j,i) ) * &
    580580                                          dd2zu(k)
    581581                         ENDDO
     
    742742
    743743       USE arrays_3d,                                                          &
    744            ONLY:  ddzw, dd2zu, kh, km, pt, q, ql, qsws, qswst, rho, shf,       &
     744           ONLY:  ddzw, dd2zu, kh, km, pt, q, ql, qsws, qswst, rho_ocean, shf,       &
    745745                  tend, tswst, u, v, vpt, w
    746746
     
    806806!--    Calculate TKE production by shear
    807807       !$acc kernels present( ddzw, dd2zu, kh, km, nzb_diff_s_inner, nzb_diff_s_outer ) &
    808        !$acc         present( nzb_s_inner, pt, q, ql, qsws, qswst, rho )                &
     808       !$acc         present( nzb_s_inner, pt, q, ql, qsws, qswst, rho_ocean )                &
    809809       !$acc         present( shf, tend, tswst, u, v, vpt, w, wall_e_x, wall_e_y )      &
    810810       !$acc         copyin( u_0, v_0 )
     
    11661166                               tend(k,j,i) = tend(k,j,i) +                     &
    11671167                                             kh(k,j,i) * g / rho_reference *   &
    1168                                              ( rho(k+1,j,i) - rho(k-1,j,i) ) * &
     1168                                             ( rho_ocean(k+1,j,i) - rho_ocean(k-1,j,i) ) * &
    11691169                                             dd2zu(k)
    11701170                            ENDIF
     
    12121212                            IF ( k > nzb_s_inner(j,i) )  THEN
    12131213                               tend(k,j,i) = tend(k,j,i) +                     &
    1214                                              kh(k,j,i) * g / rho(k,j,i) *      &
    1215                                              ( rho(k+1,j,i) - rho(k-1,j,i) ) * &
     1214                                             kh(k,j,i) * g / rho_ocean(k,j,i) *      &
     1215                                             ( rho_ocean(k+1,j,i) - rho_ocean(k-1,j,i) ) * &
    12161216                                             dd2zu(k)
    12171217                            ENDIF
     
    14101410
    14111411       USE arrays_3d,                                                          &
    1412            ONLY:  ddzw, dd2zu, kh, km, pt, q, ql, qsws, qswst, rho, shf,       &
     1412           ONLY:  ddzw, dd2zu, kh, km, pt, q, ql, qsws, qswst, rho_ocean, shf,       &
    14131413                  tend, tswst, u, v, vpt, w
    14141414
     
    17661766                      tend(k,j,i) = tend(k,j,i) +                   &
    17671767                                    kh(k,j,i) * g / rho_reference * &
    1768                                     ( rho(k+1,j,i) - rho(k-1,j,i) ) * dd2zu(k)
     1768                                    ( rho_ocean(k+1,j,i) - rho_ocean(k-1,j,i) ) * dd2zu(k)
    17691769                   ENDDO
    17701770
     
    17971797                   DO  k = nzb_s_inner(j,i)+1, nzt
    17981798                      tend(k,j,i) = tend(k,j,i) +                &
    1799                                     kh(k,j,i) * g / rho(k,j,i) * &
    1800                                     ( rho(k+1,j,i) - rho(k-1,j,i) ) * dd2zu(k)
     1799                                    kh(k,j,i) * g / rho_ocean(k,j,i) * &
     1800                                    ( rho_ocean(k+1,j,i) - rho_ocean(k-1,j,i) ) * dd2zu(k)
    18011801                   ENDDO
    18021802
  • palm/trunk/SOURCE/prognostic_equations.f90

    r2012 r2031  
    2020! Current revisions:
    2121! ------------------
    22 !
     22! renamed variable rho to rho_ocean
    2323!
    2424! Former revisions:
     
    215215               nr, nr_p, nrsws, nrswst, pt, ptdf_x, ptdf_y, pt_init, pt_p,     &
    216216               prho, q, q_init, q_p, qsws, qswst, qr, qr_p, qrsws, qrswst, rdf,&
    217                rdf_sc, ref_state, rho, s, s_init, s_p, sa, sa_init, sa_p,      &
     217               rdf_sc, ref_state, rho_ocean, s, s_init, s_p, sa, sa_init, sa_p,      &
    218218               saswsb, saswst, shf, ssws, sswst, tend,                         &
    219219               te_m, tnr_m, tpt_m, tq_m, tqr_m, ts_m, tsa_m, tswst, tu_m, tv_m,&
     
    562562          IF ( .NOT. neutral )  THEN
    563563             IF ( ocean )  THEN
    564                 CALL buoyancy( i, j, rho, 3 )
     564                CALL buoyancy( i, j, rho_ocean, 3 )
    565565             ELSE
    566566                IF ( .NOT. humidity )  THEN
     
    12671267    IF ( .NOT. neutral )  THEN
    12681268       IF ( ocean )  THEN
    1269           CALL buoyancy( rho, 3 )
     1269          CALL buoyancy( rho_ocean, 3 )
    12701270       ELSE
    12711271          IF ( .NOT. humidity )  THEN
     
    22052205    IF ( .NOT. neutral )  THEN
    22062206       IF ( ocean )  THEN
    2207           CALL buoyancy( rho, 3 )
     2207          CALL buoyancy( rho_ocean, 3 )
    22082208       ELSE
    22092209          IF ( .NOT. humidity )  THEN
  • palm/trunk/SOURCE/read_3d_binary.f90

    r2001 r2031  
    2020! Current revisions:
    2121! -----------------
    22 !
     22! renamed variable rho_av to rho_ocean_av
    2323!
    2424! Former revisions:
     
    712712                   IF ( k == 1 )  READ ( 13 )  seq_random_array
    713713
    714                 CASE ( 'rho_av' )
    715                    IF ( .NOT. ALLOCATED( rho_av ) )  THEN
    716                       ALLOCATE( rho_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
    717                    ENDIF
    718                    IF ( k == 1 )  READ ( 13 )  tmp_3d
    719                    rho_av(:,nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) = &
     714                CASE ( 'rho_ocean_av' )
     715                   IF ( .NOT. ALLOCATED( rho_ocean_av ) )  THEN
     716                      ALLOCATE( rho_ocean_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
     717                   ENDIF
     718                   IF ( k == 1 )  READ ( 13 )  tmp_3d
     719                   rho_ocean_av(:,nysc-nbgp:nync+nbgp,nxlc-nbgp:nxrc+nbgp) = &
    720720                                tmp_3d(:,nysf-nbgp:nynf+nbgp,nxlf-nbgp:nxrf+nbgp)
    721721
  • palm/trunk/SOURCE/sum_up_3d_data.f90

    r2025 r2031  
    2020! Current revisions:
    2121! -----------------
    22 !
     22! renamed variable rho to rho_ocean and rho_av to rho_ocean_av
    2323!
    2424! Former revisions:
     
    130130    USE arrays_3d,                                                             &
    131131        ONLY:  dzw, e, nr, ol, p, pt, precipitation_rate, q, qc, ql, ql_c,     &
    132                ql_v, qr, qsws, rho, s, sa, shf, ssws, ts, u, us, v, vpt, w, z0,&
     132               ql_v, qr, qsws, rho_ocean, s, sa, shf, ssws, ts, u, us, v, vpt, w, z0,&
    133133               z0h, z0q
    134134
     
    136136        ONLY:  e_av, lpt_av, lwp_av, nr_av, ol_av, p_av, pc_av, pr_av, prr_av, &
    137137               precipitation_rate_av, pt_av, q_av, qc_av, ql_av, ql_c_av,      &
    138                ql_v_av, ql_vp_av, qr_av, qsws_av, qv_av, rho_av, s_av, sa_av,  &
     138               ql_v_av, ql_vp_av, qr_av, qsws_av, qv_av, rho_ocean_av, s_av, sa_av,  &
    139139               shf_av, ssws_av, ts_av, u_av, us_av, v_av, vpt_av, w_av, z0_av, &
    140140               z0h_av, z0q_av
     
    323323                qv_av = 0.0_wp
    324324
    325              CASE ( 'rho' )
    326                 IF ( .NOT. ALLOCATED( rho_av ) )  THEN
    327                    ALLOCATE( rho_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
    328                 ENDIF
    329                 rho_av = 0.0_wp
     325             CASE ( 'rho_ocean' )
     326                IF ( .NOT. ALLOCATED( rho_ocean_av ) )  THEN
     327                   ALLOCATE( rho_ocean_av(nzb:nzt+1,nysg:nyng,nxlg:nxrg) )
     328                ENDIF
     329                rho_ocean_av = 0.0_wp
    330330
    331331             CASE ( 's' )
     
    657657             ENDDO
    658658
    659           CASE ( 'rho' )
    660              DO  i = nxlg, nxrg
    661                 DO  j = nysg, nyng
    662                    DO  k = nzb, nzt+1
    663                       rho_av(k,j,i) = rho_av(k,j,i) + rho(k,j,i)
     659          CASE ( 'rho_ocean' )
     660             DO  i = nxlg, nxrg
     661                DO  j = nysg, nyng
     662                   DO  k = nzb, nzt+1
     663                      rho_ocean_av(k,j,i) = rho_ocean_av(k,j,i) + rho_ocean(k,j,i)
    664664                   ENDDO
    665665                ENDDO
  • 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
  • palm/trunk/SOURCE/time_integration.f90

    r2012 r2031  
    2020! Current revisions:
    2121! ------------------
    22 !
     22! renamed variable rho to rho_ocean
    2323!
    2424! Former revisions:
     
    242242    USE arrays_3d,                                                             &
    243243        ONLY:  diss, dzu, e, e_p, nr_p, prho, pt, pt_p, pt_init, q_init, q,    &
    244                ql, ql_c, ql_v, ql_vp, qr_p, q_p, ref_state, rho, s, s_p, sa_p, &
     244               ql, ql_c, ql_v, ql_vp, qr_p, q_p, ref_state, rho_ocean, s, s_p, sa_p, &
    245245               tend, u, u_p, v, vpt, v_p, w, w_p
    246246
     
    482482             ENDIF
    483483             IF ( ocean )  THEN
    484                 CALL calc_mean_profile( rho, 64 )
     484                CALL calc_mean_profile( rho_ocean, 64 )
    485485                ref_state(:)  = hom(:,1,64,0)
    486486             ENDIF
     
    552552                IF ( ocean )  THEN
    553553                   CALL exchange_horiz( sa_p, nbgp )
    554                    CALL exchange_horiz( rho, nbgp )
     554                   CALL exchange_horiz( rho_ocean, nbgp )
    555555                   CALL exchange_horiz( prho, nbgp )
    556556                ENDIF
     
    620620                IF ( ocean )  THEN
    621621                   CALL exchange_horiz( sa_p, nbgp )
    622                    CALL exchange_horiz( rho, nbgp )
     622                   CALL exchange_horiz( rho_ocean, nbgp )
    623623                  CALL exchange_horiz( prho, nbgp )
    624624                ENDIF
     
    710710             IF ( ocean )  THEN
    711711                CALL exchange_horiz( sa_p, nbgp )
    712                 CALL exchange_horiz( rho, nbgp )
     712                CALL exchange_horiz( rho_ocean, nbgp )
    713713                CALL exchange_horiz( prho, nbgp )
    714714             ENDIF
  • palm/trunk/SOURCE/urban_surface_mod.f90

    r2025 r2031  
    2121! Current revisions:
    2222! ------------------
    23 !
     23! renamed variable rho to rho_ocean
    2424!
    2525! Former revisions:
     
    419419    INTEGER(iwp), PARAMETER                        :: icsurf   = 5         !< Surface skin layer heat capacity (J m−2 K−1 )
    420420    INTEGER(iwp), PARAMETER                        :: ithick   = 6         !< thickness of the surface (wall, roof, land)  ( m )
    421     INTEGER(iwp), PARAMETER                        :: irhoC    = 7         !< volumetric heat capacity rho*C of the material ( J m−3 K−1 )
     421    INTEGER(iwp), PARAMETER                        :: irhoC    = 7         !< volumetric heat capacity rho_ocean*C of the material ( J m−3 K−1 )
    422422    INTEGER(iwp), PARAMETER                        :: ilambdah = 8         !< thermal conductivity λH of the wall (W m−1 K−1 )
    423423    CHARACTER(12), DIMENSION(:), ALLOCATABLE       :: surface_type_names   !< names of wall types (used only for reports)
     
    28862886        IF ( plant_canopy )  THEN
    28872887            pchf_prep(:) = r_d * (hyp(nzub:nzut) / 100000.0_wp)**0.286_wp &
    2888                         / (cp * hyp(nzub:nzut) * dx*dy*dz) !< equals to 1 / (rho * c_p * Vbox * T)
     2888                        / (cp * hyp(nzub:nzut) * dx*dy*dz) !< equals to 1 / (rho_ocean * c_p * Vbox * T)
    28892889        ENDIF
    28902890
     
    38633863            ENDIF
    38643864           
    3865 !--         volumetric heat capacity rho*C of the wall ( J m−3 K−1 )
     3865!--         volumetric heat capacity rho_ocean*C of the wall ( J m−3 K−1 )
    38663866            rho_c_wall(:,l) = surface_params(irhoC, ip)
    38673867           
     
    39303930            pt1  = pt(k,j,i)
    39313931
    3932 !--         calculate rho * cp coefficient at surface layer
     3932!--         calculate rho_ocean * cp coefficient at surface layer
    39333933            rho_cp  = cp * hyp(k) / ( r_d * pt1 * exn(k) )
    39343934
  • palm/trunk/SOURCE/write_3d_binary.f90

    r2001 r2031  
    2020! Current revisions:
    2121! -----------------
    22 !
     22! renamed variable rho_av to rho_ocean_av
    2323!
    2424! Former revisions:
     
    256256    ENDIF   
    257257    IF ( ocean )  THEN
    258        IF ( ALLOCATED( rho_av ) )  THEN
    259           WRITE ( 14 )  'rho_av              ';  WRITE ( 14 )  rho_av
     258       IF ( ALLOCATED( rho_ocean_av ) )  THEN
     259          WRITE ( 14 )  'rho_ocean_av              ';  WRITE ( 14 )  rho_ocean_av
    260260       ENDIF
    261261       WRITE ( 14 )  'sa                  ';  WRITE ( 14 )  sa
Note: See TracChangeset for help on using the changeset viewer.