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 4739 for branches/2014/dev_r4650_UKMO3_masked_damping/NEMOGCM/NEMO/OPA_SRC/TRA/tradmp.F90 – NEMO

Ignore:
Timestamp:
2014-08-13T10:46:04+02:00 (10 years ago)
Author:
timgraham
Message:

Updated C1D/dyndmp.F90 and trcdmp.F90 to read restoration coefficient from a file.
Modified namelist_top_ref to match new options
Bug fixes to DMP_TOOLS tool and addition of custom.F90 to allow users to make modifications. Also changed to use working precision (wp) throughout.

File:
1 edited

Legend:

Unmodified
Added
Removed

  • branches/2014/dev_r4650_UKMO3_masked_damping/NEMOGCM/NEMO/OPA_SRC/TRA/tradmp.F90

    r4738 r4739  
    3636   USE wrk_nemo       ! Memory allocation 
    3737   USE timing         ! Timing 
     38   USE iom 
    3839 
    3940   IMPLICIT NONE 
    … …  
    4647   LOGICAL , PUBLIC ::   ln_tradmp   !: internal damping flag 
    4748   INTEGER , PUBLIC ::   nn_zdmp     ! = 0/1/2 flag for damping in the mixed layer 
    48    CHARACTER(LEN=200) :: cn_resto      ! name of netcdf file containing restoration coefficient field 
     49   CHARACTER(LEN=200) , PUBLIC :: cn_resto      ! name of netcdf file containing restoration coefficient field 
    4950   LOGICAL , PUBLIC ::   ln_miss      ! check for missing data in T/S data file (slow?) 
    5051   REAL(wp), PUBLIC ::   rn_miss      ! Value of missing data 
     52   ! 
    5153 
    5254 
    … …  
    189191      !!---------------------------------------------------------------------- 
    190192 
    191       NAMELIST/namtra_dmp/ ln_tradmp, rn_dmp, nn_zdmp, cn_resto 
     193      NAMELIST/namtra_dmp/ ln_tradmp, nn_zdmp, cn_resto 
    192194      INTEGER ::  ios         ! Local integer for output status of namelist read 
    193       REAL(wp), POINTER, DIMENSION(:,:,:) ::  dmp_mask   ! 3D mask 
     195      INTEGER :: imask        ! File handle  
    194196      !!---------------------------------------------------------------------- 
    195197 
    … …  
    224226         CASE ( 1 )  ; IF(lwp) WRITE(numout,*) '   no tracer damping in the turbocline' 
    225227         CASE ( 2 )  ; IF(lwp) WRITE(numout,*) '   no tracer damping in the mixed layer' 
     228         END SELECT 
    226229 
    227230         !Initialisation of dtatsd - Would it be better to have dmpdta routine 
Note: See TracChangeset for help on using the changeset viewer.