New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 771 for branches/dev_001_GM/NEMO/TOP_SRC/TRP/trctrp_lec.F90 – NEMO

Ignore:
Timestamp:
2007-12-17T11:51:41+01:00 (16 years ago)
Author:
gm
Message:

dev_001_GM - small error corrections

File:
1 edited

Legend:

Unmodified
Added
Removed

  • branches/dev_001_GM/NEMO/TOP_SRC/TRP/trctrp_lec.F90

    r770 r771  
    11MODULE trctrp_lec 
    2    !!============================================================================== 
     2   !!====================================================================== 
    33   !!                       ***  MODULE  trctrp_lec  *** 
    4    !! Ocean passive tracers:  namelist read options for transport 
    5    !!============================================================================== 
     4   !! TOP :   namelist read options for transport 
     5   !!====================================================================== 
    66   !! History :   1.0  !  2004-03  (C. Ethe)  Original code 
    77   !!             2.0  !  2007-12  (C. Ethe, G. Madec)  revised architecture 
    … …  
    99#if defined key_passivetrc 
    1010   !!---------------------------------------------------------------------- 
    11    !!   'key_passivetrc'                                          TOP model 
     11   !!   'key_passivetrc'                                         TOP models 
    1212   !!---------------------------------------------------------------------- 
    1313   !!   trc_trp_lec  : read the passive tracer namelist for transport 
    … …  
    9696      !!---------------------------------------------------------------------- 
    9797 
    98       ! Read Namelist namtrcadv : tracer advection scheme 
    99       ! ------------------------- 
    100       REWIND( numnat ) 
     98      IF(lwp) WRITE(numout,*) 
     99      IF(lwp) WRITE(numout,*) ' trc_trp_lec: read namelist for tracer transport' 
     100      IF(lwp) WRITE(numout,*) ' ~~~~~~~~~~~' 
     101 
     102      !                                ! tracer advection scheme 
     103      REWIND( numnat )                 !   ! Read Namelist namtrcadv 
    101104      READ  ( numnat, namtrcadv ) 
    102105 
    103       ! Parameter control and print 
    104       ! --------------------------- 
    105       ! Control print 
    106       IF(lwp) THEN 
     106      IF(lwp) THEN                     !   ! Control print 
    107107         WRITE(numout,*) 
    108          WRITE(numout,*) 'choice/control of the tracer advection scheme' 
    109          WRITE(numout,*) '~~~~~~~~~~~' 
    110          WRITE(numout,*) '          Namelist namtrcadv : chose a advection scheme for tracers' 
    111          WRITE(numout,*) 
    112          WRITE(numout,*) '             2nd order advection scheme     ln_trcadv_cen2   = ', ln_trcadv_cen2 
    113          WRITE(numout,*) '             TVD advection scheme           ln_trcadv_tvd    = ', ln_trcadv_tvd 
    114          WRITE(numout,*) '             MUSCL  advection scheme        ln_trcadv_muscl  = ', ln_trcadv_muscl 
    115          WRITE(numout,*) '             MUSCL2 advection scheme        ln_trcadv_muscl2 = ', ln_trcadv_muscl2 
    116          WRITE(numout,*) '             SMOLARKIEWICZ advection scheme ln_trcadv_smolar = ', ln_trcadv_smolar 
     108         WRITE(numout,*) '   Namelist namtrcadv : chose a advection scheme for tracers' 
     109         WRITE(numout,*) '      2nd order advection scheme     ln_trcadv_cen2   = ', ln_trcadv_cen2 
     110         WRITE(numout,*) '      TVD advection scheme           ln_trcadv_tvd    = ', ln_trcadv_tvd 
     111         WRITE(numout,*) '      MUSCL  advection scheme        ln_trcadv_muscl  = ', ln_trcadv_muscl 
     112         WRITE(numout,*) '      MUSCL2 advection scheme        ln_trcadv_muscl2 = ', ln_trcadv_muscl2 
     113         WRITE(numout,*) '      SMOLARKIEWICZ advection scheme ln_trcadv_smolar = ', ln_trcadv_smolar 
    117114      ENDIF 
    118115 
    119116#if  defined key_trcbbl_dif 
    120       ! Read Namelist namtrcbbl : Bottom boundary layer coef 
    121       ! -------------------------------------------------- 
    122       REWIND( numnat ) 
     117      !                                ! Bottom boundary layer  
     118      REWIND( numnat )                 !   ! Read Namelist namtrcbbl 
    123119      READ  ( numnat, namtrcbbl ) 
    124120 
    125       ! Parameter control and print 
    126       ! --------------------------- 
    127       IF(lwp) THEN 
    128          WRITE(numout,*) ' Diffusive Bottom Boundary Layer' 
    129          WRITE(numout,*) '~~~~~~~~' 
    130          WRITE(numout,*) ' bottom boundary layer coef.    atrcbbl = ', atrcbbl 
     121      IF(lwp) THEN                     !   ! Control print 
     122         WRITE(numout,*) 
     123         WRITE(numout,*) '   Namelist namtrcbbl : set Diffusive Bottom Boundary Layer parameters' 
     124         WRITE(numout,*) '      bottom boundary layer coef.    atrcbbl = ', atrcbbl 
    131125# if defined key_trcbbl_adv 
    132          WRITE(numout,*) ' * Advective Bottom Boundary Layer' 
     126         WRITE(numout,*) '   * Advective Bottom Boundary Layer' 
    133127# endif 
    134          WRITE(numout,*) 
    135128      ENDIF 
    136129#endif 
    137130 
    138       !  Define the lateral tracer physics parameters 
    139       ! ============================================= 
    140      
    141       ! Read Namelist namtrcldf : Lateral physics on tracers 
    142       REWIND( numnat ) 
     131      !                                ! Lateral physics on tracers    
     132      REWIND( numnat )                 !   ! Read Namelist namtrcldf 
    143133      READ  ( numnat, namtrcldf ) 
    144134 
    145       IF(lwp) THEN 
     135      IF(lwp) THEN                     !   ! Control print 
    146136         WRITE(numout,*) 
    147          WRITE(numout,*) 'lateral passive tracer physics' 
    148          WRITE(numout,*) '~~~~~~~' 
    149137         WRITE(numout,*) '   Namelist namtrcldf : set lateral mixing parameters (type, direction, coefficients)' 
    150          WRITE(numout,*) '     perform lateral diffusion or not               ln_trcldf_diff  = ', ln_trcldf_diff 
    151          WRITE(numout,*) '     laplacian operator                             ln_trcldf_lap   = ', ln_trcldf_lap 
    152          WRITE(numout,*) '     bilaplacian operator                           ln_trcldf_bilap = ', ln_trcldf_bilap 
    153          WRITE(numout,*) '     iso-level                                      ln_trcldf_level = ', ln_trcldf_level 
    154          WRITE(numout,*) '     horizontal (geopotential)                      ln_trcldf_hor   = ', ln_trcldf_hor 
    155          WRITE(numout,*) '     iso-neutral                                    ln_trcldf_iso   = ', ln_trcldf_iso 
    156          WRITE(numout,*) '     lateral eddy diffusivity                              ahtrc0   = ', ahtrc0 
    157          WRITE(numout,*) '     background hor. diffusivity                            ahtrb0  = ', ahtrb0 
    158          WRITE(numout,*) '     eddy induced velocity coef.                           aeivtr0  = ', aeivtr0 
    159          WRITE(numout,*) '     ratio between passive and active tracer diffusion coef  trcrat = ', trcrat 
     138         WRITE(numout,*) '      perform lateral diffusion or not               ln_trcldf_diff  = ', ln_trcldf_diff 
     139         WRITE(numout,*) '      laplacian operator                             ln_trcldf_lap   = ', ln_trcldf_lap 
     140         WRITE(numout,*) '      bilaplacian operator                           ln_trcldf_bilap = ', ln_trcldf_bilap 
     141         WRITE(numout,*) '      iso-level                                      ln_trcldf_level = ', ln_trcldf_level 
     142         WRITE(numout,*) '      horizontal (geopotential)                      ln_trcldf_hor   = ', ln_trcldf_hor 
     143         WRITE(numout,*) '      iso-neutral                                    ln_trcldf_iso   = ', ln_trcldf_iso 
     144         WRITE(numout,*) '      lateral eddy diffusivity                              ahtrc0   = ', ahtrc0 
     145         WRITE(numout,*) '      background hor. diffusivity                            ahtrb0  = ', ahtrb0 
     146         WRITE(numout,*) '      eddy induced velocity coef.                           aeivtr0  = ', aeivtr0 
     147         WRITE(numout,*) '      ratio between passive and active tracer diffusion coef  trcrat = ', trcrat 
    160148      ENDIF 
    161149 
    162       ! Read namtrczdf namelist : vertical mixing parameters 
    163       ! -------------------- 
    164       REWIND( numnat ) 
     150      !                                ! Vertical mixing 
     151      REWIND( numnat )                 !   ! Read namtrczdf namelist 
    165152      READ  ( numnat, namtrczdf ) 
    166153 
    167       ! Parameter print 
    168       ! --------------- 
    169       IF(lwp) THEN 
     154      IF(lwp) THEN                     !   ! Control print 
    170155         WRITE(numout,*) 
    171          WRITE(numout,*) 'vertical physics' 
    172          WRITE(numout,*) '~~~~~~~~' 
    173          WRITE(numout,*) '          Namelist namtrczdf : set vertical diffusion parameters' 
    174          WRITE(numout,*) '             time splitting / backward scheme ln_trczdf_exp = ', ln_trczdf_exp 
    175          WRITE(numout,*) '             number of time step               n_trczdf_exp = ', n_trczdf_exp 
     156         WRITE(numout,*) '   Namelist namtrczdf : set vertical diffusion parameters' 
     157         WRITE(numout,*) '      time splitting / backward scheme ln_trczdf_exp = ', ln_trczdf_exp 
     158         WRITE(numout,*) '      number of time step               n_trczdf_exp = ', n_trczdf_exp 
    176159      ENDIF 
    177160 
    178161# if defined key_trcdmp 
    179       ! Read Namelist namtdp : passive tracres damping term 
    180       ! -------------------- 
    181       REWIND( numnat ) 
     162      !                                ! passive tracres damping term 
     163      REWIND( numnat )                 !   ! Read Namelist namtdp 
    182164      READ  ( numnat, namtrcdmp ) 
    183       IF( lzoom )   nmldmptr = 0      ! restoring to climatology at closed north or south boundaries 
     165      IF( lzoom )   nmldmptr = 0           ! restoring to climatology at closed north or south boundaries 
    184166 
    185       ! Parameter control and print 
    186       ! --------------------------- 
    187       IF(lwp) THEN 
     167      IF(lwp) THEN                     !   ! Control print 
    188168         WRITE(numout,*) 
    189          WRITE(numout,*) 'newtonian damping' 
    190          WRITE(numout,*) '~~~~~~~' 
    191          WRITE(numout,*) '          Namelist namtrcdmp : set damping parameter' 
    192          WRITE(numout,*) '             tracers damping option         ndmptr   = ', ndmptr 
    193          WRITE(numout,*) '             create a damping.coeff file    ndmpftr  = ', ndmpftr 
    194          WRITE(numout,*) '             mixed layer damping option     nmldmptr = ', nmldmptr, '(zoom: forced to 0)' 
    195          WRITE(numout,*) '             surface time scale (days)      sdmptr   = ', sdmptr 
    196          WRITE(numout,*) '             bottom time scale (days)       bdmptr   = ', bdmptr 
    197          WRITE(numout,*) '             depth of transition (meters)   hdmptr   = ', hdmptr 
    198          WRITE(numout,*) 
     169         WRITE(numout,*) '   Namelist namtrcdmp : set damping parameter' 
     170         WRITE(numout,*) '      tracers damping option         ndmptr   = ', ndmptr 
     171         WRITE(numout,*) '      create a damping.coeff file    ndmpftr  = ', ndmpftr 
     172         WRITE(numout,*) '      mixed layer damping option     nmldmptr = ', nmldmptr, '(zoom: forced to 0)' 
     173         WRITE(numout,*) '      surface time scale (days)      sdmptr   = ', sdmptr 
     174         WRITE(numout,*) '      bottom time scale (days)       bdmptr   = ', bdmptr 
     175         WRITE(numout,*) '      depth of transition (meters)   hdmptr   = ', hdmptr 
    199176      ENDIF 
    200177#endif 
Note: See TracChangeset for help on using the changeset viewer.