Ignore:
Timestamp:
2010-06-24T17:52:01+02:00 (10 years ago)
Author:
acc
Message:

ticket #684 step 4: Add in changes between the head of the dev_r1821_mixed_ldfdyn branch and the trunk@1821

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/DEV_r1784_mid_year_merge_2010/NEMO/OPA_SRC/LDF/ldfdyn.F90

    r1601 r1954  
    6767      NAMELIST/namdyn_ldf/ ln_dynldf_lap  , ln_dynldf_bilap,                  & 
    6868         &                 ln_dynldf_level, ln_dynldf_hor  , ln_dynldf_iso,   & 
    69          &                 rn_ahm_0       , rn_ahmb_0 
     69         &                 rn_ahm_0       , rn_ahmb_0      , rn_ahm_0_blp 
    7070      !!---------------------------------------------------------------------- 
    7171 
     
    7878         WRITE(numout,*) '~~~~~~~' 
    7979         WRITE(numout,*) '   Namelist nam_dynldf : set lateral mixing parameters' 
    80          WRITE(numout,*) '      laplacian operator          ln_dynldf_lap   = ', ln_dynldf_lap 
    81          WRITE(numout,*) '      bilaplacian operator        ln_dynldf_bilap = ', ln_dynldf_bilap 
    82          WRITE(numout,*) '      iso-level                   ln_dynldf_level = ', ln_dynldf_level 
    83          WRITE(numout,*) '      horizontal (geopotential)   ln_dynldf_hor   = ', ln_dynldf_hor 
    84          WRITE(numout,*) '      iso-neutral                 ln_dynldf_iso   = ', ln_dynldf_iso 
    85          WRITE(numout,*) '      horizontal eddy viscosity   rn_ahm_0        = ', rn_ahm_0 
    86          WRITE(numout,*) '      background viscosity        rn_ahmb_0       = ', rn_ahmb_0 
    87       ENDIF 
    88  
    89       ahm0  = rn_ahm_0                  ! OLD namelist variables defined from DOCTOR namelist variables 
    90       ahmb0 = rn_ahmb_0 
     80         WRITE(numout,*) '      laplacian operator                      ln_dynldf_lap   = ', ln_dynldf_lap 
     81         WRITE(numout,*) '      bilaplacian operator                    ln_dynldf_bilap = ', ln_dynldf_bilap 
     82         WRITE(numout,*) '      iso-level                               ln_dynldf_level = ', ln_dynldf_level 
     83         WRITE(numout,*) '      horizontal (geopotential)               ln_dynldf_hor   = ', ln_dynldf_hor 
     84         WRITE(numout,*) '      iso-neutral                             ln_dynldf_iso   = ', ln_dynldf_iso 
     85         WRITE(numout,*) '      horizontal laplacian eddy viscosity     rn_ahm_0        = ', rn_ahm_0 
     86         WRITE(numout,*) '      background viscosity                    rn_ahmb_0       = ', rn_ahmb_0 
     87         WRITE(numout,*) '      horizontal bilaplacian eddy viscosity   rn_ahm_0        = ', rn_ahm_0 
     88 
     89      ENDIF 
     90 
     91      ahm0     = rn_ahm_0                  ! OLD namelist variables defined from DOCTOR namelist variables 
     92      ahmb0    = rn_ahmb_0 
     93      ahm0_blp = rn_ahm_0_blp 
    9194 
    9295      ! ... check of lateral diffusive operator on tracers 
     
    117120      IF( ln_dynldf_bilap ) THEN 
    118121         IF(lwp) WRITE(numout,*) '   biharmonic momentum diffusion' 
    119          IF( ahm0 > 0 .AND. .NOT. lk_esopa )   CALL ctl_stop( 'The horizontal viscosity coef. ahm0 must be negative' ) 
     122         IF( ahm0_blp == 0.0 ) ahm0_blp = ahm0       ! Old namelist method: bilap specified with ahm0 
     123         IF( .NOT. ln_dynldf_lap ) ahm0 = ahm0_blp   ! Allow spatially varying coefs, which use ahm0 as input 
     124         IF( ahm0_blp > 0 .AND. .NOT. lk_esopa )   CALL ctl_stop( 'The horizontal viscosity coef. ahm0 must be negative' ) 
    120125      ELSE 
    121126         IF(lwp) WRITE(numout,*) '   harmonic momentum diff. (default)' 
Note: See TracChangeset for help on using the changeset viewer.