Changeset 5960


Ignore:
Timestamp:
2015-12-01T11:31:27+01:00 (5 years ago)
Author:
timgraham
Message:

Merged in changes from /branches/2014/dev_r4650_UKMO12_CFL_diags_take2

Location:
branches/2015/dev_MetOffice_merge_2015
Files:
10 edited
1 copied

Legend:

Unmodified
Added
Removed
  • branches/2015/dev_MetOffice_merge_2015/DOC/TexFiles/Chapters/Chap_DIA.tex

    • Property svn:keywords set to Id
    r5515 r5960  
    15221522\\ 
    15231523 
     1524\section{Courant numbers} 
     1525Courant numbers provide a theoretical indication of the model's numerical stability. The advective Courant numbers can be calculated according to 
     1526\begin{equation} 
     1527\label{eq:CFL} 
     1528C_u = |u|\frac{\rdt}{e_{1u}}, \quad C_v = |v|\frac{\rdt}{e_{2v}}, \quad C_w = |w|\frac{\rdt}{e_{3w}} 
     1529\end{equation} 
     1530in the zonal, meridional and vertical directions respectively. The vertical component is included although it is not strictly valid as the vertical velocity is calculated from the continuity equation rather than as a prognostic variable. Physically this represents the rate at which information is propogated across a grid cell. Values greater than 1 indicate that information is propagated across more than one grid cell in a single time step. 
     1531 
     1532The variables can be activated by setting the \np{nn\_diacfl} namelist parameter to 1 in the \ngn{namctl} namelist. The diagnostics will be written out to an ascii file named cfl\_diagnostics.ascii. In this file the maximum value of $C_u$, $C_v$, and $C_w$ are printed at each timestep along with the coordinates of where the maximum value occurs. At the end of the model run the maximum value of $C_u$, $C_v$, and $C_w$ for the whole model run is printed along with the coordinates of each. The maximum values from the run are also copied to the ocean.output file.  
    15241533 
    15251534 
  • branches/2015/dev_MetOffice_merge_2015/DOC/TexFiles/Namelist/namctl

    • Property svn:keywords set to Id
    r4147 r5960  
    1313                           !     (no physical validity of the results) 
    1414   nn_timing   =    0      !  timing by routine activated (=1) creates timing.output file, or not (=0) 
     15   nn_diacfl   =    0      !  Write out cfl diagnostics (=1) in cfl_diagnostics.ascii, or not (=0) 
    1516/ 
  • branches/2015/dev_MetOffice_merge_2015/NEMOGCM/CONFIG/SHARED/namelist_ref

    • Property svn:keywords set to Id
    r5930 r5960  
    10281028                           !     (no physical validity of the results) 
    10291029   nn_timing   =    0      !  timing by routine activated (=1) creates timing.output file, or not (=0) 
     1030   nn_diacfl   =    0      !  Write out cfl diagnostics (=1) in cfl_diagnostics.ascii, or not (=0) 
    10301031/ 
    10311032!----------------------------------------------------------------------- 
  • branches/2015/dev_MetOffice_merge_2015/NEMOGCM/NEMO/OFF_SRC/nemogcm.F90

    r5836 r5960  
    148148      NAMELIST/namctl/ ln_ctl  , nn_print, nn_ictls, nn_ictle,   & 
    149149         &             nn_isplt, nn_jsplt, nn_jctls, nn_jctle,   & 
    150          &             nn_bench, nn_timing 
     150         &             nn_bench, nn_timing, nn_diacfl 
    151151      NAMELIST/namcfg/ cp_cfg, cp_cfz, jp_cfg, jpidta, jpjdta, jpkdta, jpiglo, jpjglo, & 
    152152         &             jpizoom, jpjzoom, jperio, ln_use_jattr 
  • branches/2015/dev_MetOffice_merge_2015/NEMOGCM/NEMO/OOO_SRC/nemogcm.F90

    r5407 r5960  
    127127      NAMELIST/namctl/ ln_ctl, nn_print, nn_ictls, nn_ictle,   & 
    128128         &             nn_isplt, nn_jsplt, nn_jctls, nn_jctle,   & 
    129          &             nn_bench, nn_timing 
     129         &             nn_bench, nn_timing, nn_diacfl 
    130130      NAMELIST/namcfg/ cp_cfg, cp_cfz, jp_cfg, jpidta, jpjdta, jpkdta, jpiglo, jpjglo, & 
    131131         &             jpizoom, jpjzoom, jperio, ln_use_jattr 
  • branches/2015/dev_MetOffice_merge_2015/NEMOGCM/NEMO/OPA_SRC/IOM/in_out_manager.F90

    r5518 r5960  
    9999   LOGICAL ::   ln_ctl       !: run control for debugging 
    100100   INTEGER ::   nn_timing    !: run control for timing 
     101   INTEGER ::   nn_diacfl    !: flag whether to create CFL diagnostics 
    101102   INTEGER ::   nn_print     !: level of print (0 no print) 
    102103   INTEGER ::   nn_ictls     !: Start i indice for the SUM control 
  • branches/2015/dev_MetOffice_merge_2015/NEMOGCM/NEMO/OPA_SRC/nemogcm.F90

    r5836 r5960  
    6666   USE diadct          ! sections transports           (dia_dct_init routine) 
    6767   USE diaobs          ! Observation diagnostics       (dia_obs_init routine) 
     68   USE diacfl          ! CFL diagnostics               (dia_cfl_init routine) 
    6869   USE lib_fortran     ! Fortran utilities (allows no signed zero when 'key_nosignedzero' defined) 
    6970   USE step            ! NEMO time-stepping                 (stp     routine) 
     
    231232      NAMELIST/namctl/ ln_ctl  , nn_print, nn_ictls, nn_ictle,   & 
    232233         &             nn_isplt, nn_jsplt, nn_jctls, nn_jctle,   & 
    233          &             nn_bench, nn_timing 
     234         &             nn_bench, nn_timing, nn_diacfl 
    234235      NAMELIST/namcfg/ cp_cfg, cp_cfz, jp_cfg, jpidta, jpjdta, jpkdta, jpiglo, jpjglo, & 
    235236         &             jpizoom, jpjzoom, jperio, ln_use_jattr 
     
    458459      !                                      ! Diagnostics 
    459460      IF( lk_floats     )   CALL     flo_init   ! drifting Floats 
     461                            CALL dia_cfl_init   ! Initialise CFL diagnostics 
    460462      IF( lk_diaar5     )   CALL dia_ar5_init   ! ar5 diag 
    461463                            CALL dia_ptr_init   ! Poleward TRansports initialization 
  • branches/2015/dev_MetOffice_merge_2015/NEMOGCM/NEMO/OPA_SRC/step.F90

    r5930 r5960  
    234234      ! diagnostics and outputs              
    235235      !<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<<< 
    236       IF( lk_floats  )   CALL flo_stp       ( kstp )  ! drifting Floats 
    237       IF( lk_diahth  )   CALL dia_hth       ( kstp )  ! Thermocline depth (20 degres isotherm depth) 
    238       IF(.NOT.ln_cpl )   CALL dia_fwb       ( kstp )  ! Fresh water budget diagnostics 
    239       IF( lk_diadct  )   CALL dia_dct       ( kstp )  ! Transports 
    240       IF( lk_diaar5  )   CALL dia_ar5       ( kstp )  ! ar5 diag 
    241       IF( lk_diaharm )   CALL dia_harm      ( kstp )  ! Tidal harmonic analysis 
    242                          CALL dia_wri       ( kstp )  ! ocean model: outputs 
     236      IF( lk_floats  )   CALL flo_stp( kstp )         ! drifting Floats 
     237      IF( nn_diacfl == 1 )   CALL dia_cfl( kstp )         ! Courant number diagnostics 
     238      IF( lk_diahth  )   CALL dia_hth( kstp )         ! Thermocline depth (20 degres isotherm depth) 
     239      IF(.NOT.ln_cpl )   CALL dia_fwb( kstp )         ! Fresh water budget diagnostics 
     240      IF( lk_diadct  )   CALL dia_dct( kstp )         ! Transports 
     241      IF( lk_diaar5  )   CALL dia_ar5( kstp )         ! ar5 diag 
     242      IF( lk_diaharm )   CALL dia_harm( kstp )        ! Tidal harmonic analysis 
     243                         CALL dia_wri( kstp )         ! ocean model: outputs 
    243244      ! 
    244245      IF( ln_crs     )   CALL crs_fld       ( kstp )  ! ocean model: online field coarsening & output 
  • branches/2015/dev_MetOffice_merge_2015/NEMOGCM/NEMO/OPA_SRC/step_oce.F90

    r5930 r5960  
    8282   USE diahsb           ! heat, salt and volume budgets    (dia_hsb routine) 
    8383   USE diaharm 
     84   USE diacfl 
    8485   USE flo_oce          ! floats variables 
    8586   USE floats           ! floats computation               (flo_stp routine) 
  • branches/2015/dev_MetOffice_merge_2015/NEMOGCM/NEMO/SAS_SRC/nemogcm.F90

    r5510 r5960  
    184184      NAMELIST/namctl/ ln_ctl  , nn_print, nn_ictls, nn_ictle,   & 
    185185         &             nn_isplt, nn_jsplt, nn_jctls, nn_jctle,   & 
    186          &             nn_bench, nn_timing 
     186         &             nn_bench, nn_timing, nn_diacfl 
    187187      NAMELIST/namcfg/ cp_cfg, cp_cfz, jp_cfg, jpidta, jpjdta, jpkdta, jpiglo, jpjglo, & 
    188188         &             jpizoom, jpjzoom, jperio, ln_use_jattr 
Note: See TracChangeset for help on using the changeset viewer.