Ignore:
Timestamp:
Mar 20, 2014 8:40:49 AM (7 years ago)
Author:
raasch
Message:

ONLY-attribute added to USE-statements,
kind-parameters added to all INTEGER and REAL declaration statements,
kinds are defined in new module kinds,
old module precision_kind is removed,
revision history before 2012 removed,
comment fields (!:) to be used for variable explanations added to all variable declaration statements

File:
1 edited

Legend:

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

    r1310 r1320  
    2020! Current revisions:
    2121! ------------------
    22 !
     22! ONLY-attribute added to USE-statements,
     23! kind-parameters added to all INTEGER and REAL declaration statements,
     24! kinds are defined in new module kinds,
     25! old module precision_kind is removed,
     26! revision history before 2012 removed,
     27! comment fields (!:) to be used for variable explanations added to
     28! all variable declaration statements
    2329!
    2430! Former revisions:
     
    5157! special treatment in case of mirror velocity boundary condition removed
    5258!
    53 ! 707 2011-03-29 11:39:40Z raasch
    54 ! bc_lr/ns replaced by bc_lr/ns_cyc
    55 !
    56 ! 667 2010-12-23 12:06:00Z suehring/gryschka
    57 ! Exchange of terminate_coupled between ocean and atmosphere via PE0
    58 ! Minimum grid spacing dxyz2_min(k) is now calculated using dzw instead of dzu
    59 !
    60 ! 622 2010-12-10 08:08:13Z raasch
    61 ! optional barriers included in order to speed up collective operations
    62 !
    63 ! 343 2009-06-24 12:59:09Z maronga
    64 ! Additional timestep criterion in case of simulations with plant canopy
    65 ! Output of messages replaced by message handling routine.
    66 !
    67 ! 222 2009-01-12 16:04:16Z letzel
    68 ! Implementation of a MPI-1 Coupling: replaced myid with target_id
    69 ! Bugfix for nonparallel execution
    70 !
    71 ! 108 2007-08-24 15:10:38Z letzel
    72 ! modifications to terminate coupled runs
    73 !
    74 ! RCS Log replace by Id keyword, revision history cleaned up
    75 !
    76 ! Revision 1.21  2006/02/23 12:59:44  raasch
    77 ! nt_anz renamed current_timestep_number
    78 !
    7959! Revision 1.1  1997/08/11 06:26:19  raasch
    8060! Initial revision
     
    8666!------------------------------------------------------------------------------!
    8767
    88     USE arrays_3d
    89     USE cloud_parameters
    90     USE control_parameters
    91     USE cpulog
    92     USE grid_variables
    93     USE indices
     68    USE arrays_3d,                                                             &
     69        ONLY:  cdc, dzu, dzw, kh, km, lad_u, lad_v, lad_w, u, v, w
     70
     71    USE cloud_parameters,                                                      &
     72        ONLY:  dt_precipitation
     73
     74    USE control_parameters,                                                    &
     75        ONLY:  cfl_factor, coupling_mode, dt_3d, dt_fixed, dt_max,             &
     76               galilei_transformation, old_dt, plant_canopy, message_string,   &
     77               stop_dt, terminate_coupled, terminate_coupled_remote,           &
     78               timestep_reason, u_gtrans, use_ug_for_galilei_tr, v_gtrans
     79
     80    USE cpulog,                                                                &
     81        ONLY:  cpu_log, log_point
     82
     83    USE grid_variables,                                                        &
     84        ONLY:  dx, dx2, dy, dy2
     85
     86    USE indices,                                                               &
     87        ONLY:  nxl, nxlg, nxr, nxrg, nyn, nyng, nys, nysg, nzb, nzt
     88
    9489    USE interfaces
     90
     91    USE kinds
     92
    9593    USE pegrid
    96     USE statistics
     94
     95    USE statistics,                                                            &
     96        ONLY:  flow_statistics_called, hom, u_max, u_max_ijk, v_max, v_max_ijk,&
     97               w_max, w_max_ijk
    9798
    9899    IMPLICIT NONE
    99100
    100     INTEGER ::  i, j, k
    101 
    102     REAL ::  div, dt_diff, dt_diff_l, dt_plant_canopy, dt_plant_canopy_l,     &
    103              dt_plant_canopy_u, dt_plant_canopy_v, dt_plant_canopy_w,         &
    104              dt_u, dt_u_l, dt_v, dt_v_l, dt_w, dt_w_l, u_gtrans_l, u_max_l,   &
    105              u_min_l, value, v_gtrans_l, v_max_l, v_min_l, w_max_l, w_min_l
    106 
    107     REAL, DIMENSION(2)         ::  uv_gtrans, uv_gtrans_l
    108     REAL, DIMENSION(3)         ::  reduce, reduce_l
    109     REAL, DIMENSION(nzb+1:nzt) ::  dxyz2_min
     101    INTEGER(iwp) ::  i !:
     102    INTEGER(iwp) ::  j !:
     103    INTEGER(iwp) ::  k !:
     104
     105    REAL(wp) ::  div               !:
     106    REAL(wp) ::  dt_diff           !:
     107    REAL(wp) ::  dt_diff_l         !:
     108    REAL(wp) ::  dt_plant_canopy   !:
     109    REAL(wp) ::  dt_plant_canopy_l !:
     110    REAL(wp) ::  dt_plant_canopy_u !:
     111    REAL(wp) ::  dt_plant_canopy_v !:
     112    REAL(wp) ::  dt_plant_canopy_w !:
     113    REAL(wp) ::  dt_u              !:
     114    REAL(wp) ::  dt_u_l            !:
     115    REAL(wp) ::  dt_v              !:
     116    REAL(wp) ::  dt_v_l            !:
     117    REAL(wp) ::  dt_w              !:
     118    REAL(wp) ::  dt_w_l            !:
     119    REAL(wp) ::  u_gtrans_l        !:
     120    REAL(wp) ::  u_max_l           !:
     121    REAL(wp) ::  u_min_l           !:
     122    REAL(wp) ::  value             !:
     123    REAL(wp) ::  v_gtrans_l        !:
     124    REAL(wp) ::  v_max_l           !:
     125    REAL(wp) ::  v_min_l           !:
     126    REAL(wp) ::  w_max_l           !:
     127    REAL(wp) ::  w_min_l           !:
     128 
     129    REAL(wp), DIMENSION(2)         ::  uv_gtrans   !:
     130    REAL(wp), DIMENSION(2)         ::  uv_gtrans_l !:
     131    REAL(wp), DIMENSION(3)         ::  reduce      !:
     132    REAL(wp), DIMENSION(3)         ::  reduce_l    !:
     133    REAL(wp), DIMENSION(nzb+1:nzt) ::  dxyz2_min   !: 
    110134
    111135
     
    222246    ENDIF
    223247#else
    224     CALL global_min_max( nzb, nzt+1, nysg, nyng, nxlg, nxrg, u, 'abs', 0.0, &
     248    CALL global_min_max( nzb, nzt+1, nysg, nyng, nxlg, nxrg, u, 'abs', 0.0_wp, &
    225249                         u_max, u_max_ijk )
    226     CALL global_min_max( nzb, nzt+1, nysg, nyng, nxlg, nxrg, v, 'abs', 0.0, &
     250    CALL global_min_max( nzb, nzt+1, nysg, nyng, nxlg, nxrg, v, 'abs', 0.0_wp, &
    227251                         v_max, v_max_ijk )
    228     CALL global_min_max( nzb, nzt+1, nysg, nyng, nxlg, nxrg, w, 'abs', 0.0, &
     252    CALL global_min_max( nzb, nzt+1, nysg, nyng, nxlg, nxrg, w, 'abs', 0.0_wp, &
    229253                         w_max, w_max_ijk )
    230254#endif
Note: See TracChangeset for help on using the changeset viewer.