Ignore:
Timestamp:
Oct 1, 2018 2:37:10 AM (3 years ago)
Author:
raasch
Message:

modularization of the ocean code

File:
1 edited

Legend:

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

    r3274 r3294  
    2525! -----------------
    2626! $Id$
     27! ocean renamed ocean_mode
     28!
     29! 3274 2018-09-24 15:42:55Z knoop
    2730! Modularization of all bulk cloud physics code components
    2831!
     
    296299                dt_3d, humidity, initializing_actions, land_surface,           &
    297300                large_scale_forcing, large_scale_subsidence, max_pr_user,      &
    298                 message_string, neutral, ocean, passive_scalar, simulated_time,&
    299                 simulated_time_at_begin, use_subsidence_tendencies,            &
    300                 use_surface_fluxes, use_top_fluxes, ws_scheme_mom,            &
    301                 ws_scheme_sca
     301                message_string, neutral, ocean_mode, passive_scalar,           &
     302                simulated_time, simulated_time_at_begin,                       &
     303                use_subsidence_tendencies, use_surface_fluxes, use_top_fluxes, &
     304                ws_scheme_mom, ws_scheme_sca
    302305       
    303306    USE cpulog,                                                                &
     
    432435!--       According to the Neumann bc for the horizontal velocity components,
    433436!--       the corresponding fluxes has to satisfiy the same bc.
    434           IF ( ocean )  THEN
     437          IF ( ocean_mode )  THEN
    435438             sums_us2_ws_l(nzt+1,:) = sums_us2_ws_l(nzt,:)
    436439             sums_vs2_ws_l(nzt+1,:) = sums_vs2_ws_l(nzt,:)
     
    459462             sums_l(:,17,i)                        = sums_wspts_ws_l(:,i)      &
    460463                                           * heatflux_output_conversion  ! w*pt*
    461              IF ( ocean          ) sums_l(:,66,i)  = sums_wssas_ws_l(:,i) ! w*sa*
     464             IF ( ocean_mode     ) sums_l(:,66,i)  = sums_wssas_ws_l(:,i) ! w*sa*
    462465             IF ( humidity       ) sums_l(:,49,i)  = sums_wsqs_ws_l(:,i)       &
    463466                                           * waterflux_output_conversion  ! w*q*
     
    490493!
    491494!--    Horizontally averaged profile of salinity
    492        IF ( ocean )  THEN
     495       IF ( ocean_mode )  THEN
    493496          !$OMP DO
    494497          DO  i = nxl, nxr
     
    563566             sums_l(:,2,0) = sums_l(:,2,0) + sums_l(:,2,i)
    564567             sums_l(:,4,0) = sums_l(:,4,0) + sums_l(:,4,i)
    565              IF ( ocean )  THEN
     568             IF ( ocean_mode )  THEN
    566569                sums_l(:,23,0) = sums_l(:,23,0) + sums_l(:,23,i)
    567570             ENDIF
     
    592595       CALL MPI_ALLREDUCE( sums_l(nzb,4,0), sums(nzb,4), nzt+2-nzb, MPI_REAL,  &
    593596                           MPI_SUM, comm2d, ierr )
    594        IF ( ocean )  THEN
     597       IF ( ocean_mode )  THEN
    595598          IF ( collective_wait )  CALL MPI_BARRIER( comm2d, ierr )
    596599          CALL MPI_ALLREDUCE( sums_l(nzb,23,0), sums(nzb,23), nzt+2-nzb,       &
     
    623626       sums(:,2) = sums_l(:,2,0)
    624627       sums(:,4) = sums_l(:,4,0)
    625        IF ( ocean )  sums(:,23) = sums_l(:,23,0)
     628       IF ( ocean_mode )  sums(:,23) = sums_l(:,23,0)
    626629       IF ( humidity ) THEN
    627630          sums(:,44) = sums_l(:,44,0)
     
    648651!
    649652!--    Salinity
    650        IF ( ocean )  THEN
     653       IF ( ocean_mode )  THEN
    651654          sums(:,23) = sums(:,23) / ngp_2dh_s_inner(:,sr)
    652655          hom(:,1,23,sr) = sums(:,23)             ! sa
     
    930933!
    931934!--             Salinity flux w"sa"
    932                 IF ( ocean )  THEN
     935                IF ( ocean_mode )  THEN
    933936                   sums_l(k,65,tn) = sums_l(k,65,tn)                           &
    934937                                         - 0.5_wp * ( kh(k,j,i) + kh(k+1,j,i) )&
     
    10001003                         sums_l(k+ki,61,tn) = sums_l(k+ki,61,tn) + &
    10011004                                    0.0_wp * rmask(j,i,sr)        ! v"pt"
    1002                          IF ( ocean )  THEN
     1005                         IF ( ocean_mode )  THEN
    10031006                            sums_l(k+ki,65,tn) = sums_l(k+ki,65,tn) + &
    10041007                                       surf_def_h(l)%sasws(m) * rmask(j,i,sr)  ! w"sa"
     
    10531056                   sums_l(nzb,61,tn) = sums_l(nzb,61,tn) + &
    10541057                                    0.0_wp * rmask(j,i,sr)        ! v"pt"
    1055                    IF ( ocean )  THEN
     1058                   IF ( ocean_mode )  THEN
    10561059                      sums_l(nzb,65,tn) = sums_l(nzb,65,tn) + &
    10571060                                       surf_lsm_h%sasws(m) * rmask(j,i,sr)  ! w"sa"
     
    11041107                   sums_l(nzb,61,tn) = sums_l(nzb,61,tn) + &
    11051108                                    0.0_wp * rmask(j,i,sr)        ! v"pt"
    1106                    IF ( ocean )  THEN
     1109                   IF ( ocean_mode )  THEN
    11071110                      sums_l(nzb,65,tn) = sums_l(nzb,65,tn) + &
    11081111                                       surf_usm_h%sasws(m) * rmask(j,i,sr)  ! w"sa"
     
    12761279                                    0.0_wp * rmask(j,i,sr)        ! v"pt"
    12771280
    1278                 IF ( ocean )  THEN
     1281                IF ( ocean_mode )  THEN
    12791282                   sums_l(nzt,65,tn) = sums_l(nzt,65,tn) + &
    12801283                                       surf_def_h(2)%sasws(m) * rmask(j,i,sr)  ! w"sa"
     
    13361339!--             Salinity flux and density (density does not belong to here,
    13371340!--             but so far there is no other suitable place to calculate)
    1338                 IF ( ocean )  THEN
     1341                IF ( ocean_mode )  THEN
    13391342                   IF( .NOT. ws_scheme_sca .OR. sr /= 0 )  THEN
    13401343                      pts = 0.5_wp * ( sa(k,j,i)   - hom(k,1,23,sr) +          &
     
    15621565!
    15631566!--    Density at top follows Neumann condition
    1564        IF ( ocean )  THEN
     1567       IF ( ocean_mode )  THEN
    15651568          sums_l(nzt+1,64,tn) = sums_l(nzt,64,tn)
    15661569          sums_l(nzt+1,71,tn) = sums_l(nzt,71,tn)
     
    20322035          first = .TRUE.
    20332036
    2034           IF ( ocean )  THEN
     2037          IF ( ocean_mode )  THEN
    20352038             DO  k = nzt, nzb+1, -1
    20362039                IF ( first  .AND.  hom(k,1,18,sr) < -1.0E-8_wp )  THEN
     
    20802083          dptdz_threshold = 0.2_wp / 100.0_wp
    20812084
    2082           IF ( ocean )  THEN
     2085          IF ( ocean_mode )  THEN
    20832086             DO  k = nzt+1, nzb+5, -1
    20842087                IF ( dptdz(k) > dptdz_threshold  .AND.                         &
Note: See TracChangeset for help on using the changeset viewer.