Changeset 5102 for trunk/NEMOGCM/NEMO/TOP_SRC
- Timestamp:
- 2015-02-23T10:42:04+01:00 (9 years ago)
- Location:
- trunk/NEMOGCM/NEMO/TOP_SRC/TRP
- Files:
-
- 2 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMOGCM/NEMO/TOP_SRC/TRP/trcdmp.F90
r4990 r5102 24 24 USE trdtra 25 25 USE trd_oce 26 USE iom 26 27 27 28 IMPLICIT NONE … … 302 303 !!---------------------------------------------------------------------- 303 304 ! 305 INTEGER :: imask !local file handle 306 304 307 IF( nn_timing == 1 ) CALL timing_start('trc_dmp_init') 305 308 ! 306 SELECT CASE ( nn_hdmp_tr )307 CASE ( -1 ) ; IF(lwp) WRITE(numout,*) ' tracer damping in the Med & Red seas only'308 CASE ( 1:90 ) ; IF(lwp) WRITE(numout,*) ' tracer damping poleward of', nn_hdmp_tr, ' degrees'309 CASE DEFAULT310 WRITE(ctmp1,*) ' bad flag value for nn_hdmp_tr = ', nn_hdmp_tr311 CALL ctl_stop(ctmp1)312 END SELECT313 309 314 310 IF( lzoom ) nn_zdmp_tr = 0 ! restoring to climatology at closed north or south boundaries … … 325 321 & CALL ctl_stop( 'passive trace damping need key_tradmp to compute damping coef.' ) 326 322 ! 327 ! ! Damping coefficients initialization328 IF( lzoom ) THEN ; CALL dtacof_zoom( restotr )329 ELSE ; CALL dtacof( nn_hdmp_tr, rn_surf_tr, rn_bot_tr, rn_dep_tr, &330 & nn_file_tr, 'TRC' , restotr)331 ENDIF323 ! ! Read damping coefficients from file 324 !Read in mask from file 325 CALL iom_open ( cn_resto_tr, imask) 326 CALL iom_get ( imask, jpdom_autoglo, 'resto', restotr) 327 CALL iom_close( imask ) 332 328 ! 333 329 IF( nn_timing == 1 ) CALL timing_stop('trc_dmp_init') -
trunk/NEMOGCM/NEMO/TOP_SRC/TRP/trcnam_trp.F90
r4624 r5102 51 51 ! !!: ** newtonian damping namelist (nam_trcdmp) ** 52 52 ! !!* Namelist namtrc_dmp : passive tracer newtonian damping * 53 INTEGER , PUBLIC :: nn_hdmp_tr ! = 0/-1/'latitude' for damping over passive tracer54 53 INTEGER , PUBLIC :: nn_zdmp_tr ! = 0/1/2 flag for damping in the mixed layer 55 REAL(wp), PUBLIC :: rn_surf_tr ! surface time scale for internal damping [days] 56 REAL(wp), PUBLIC :: rn_bot_tr ! bottom time scale for internal damping [days] 57 REAL(wp), PUBLIC :: rn_dep_tr ! depth of transition between rn_surf and rn_bot [meters] 58 INTEGER , PUBLIC :: nn_file_tr ! = 1 create a damping.coeff NetCDF file 54 CHARACTER(LEN=200) , PUBLIC :: cn_resto_tr !File containing restoration coefficient 59 55 60 56 !!---------------------------------------------------------------------- … … 82 78 NAMELIST/namtrc_zdf/ ln_trczdf_exp , nn_trczdf_exp 83 79 NAMELIST/namtrc_rad/ ln_trcrad 84 NAMELIST/namtrc_dmp/ nn_hdmp_tr, nn_zdmp_tr, rn_surf_tr, & 85 & rn_bot_tr , rn_dep_tr , nn_file_tr 80 NAMELIST/namtrc_dmp/ nn_zdmp_tr , cn_resto_tr 86 81 !!---------------------------------------------------------------------- 87 82 … … 184 179 WRITE(numout,*) '~~~~~~~' 185 180 WRITE(numout,*) ' Namelist namtrc_dmp : set damping parameter' 186 WRITE(numout,*) ' tracer damping option nn_hdmp_tr = ', nn_hdmp_tr187 181 WRITE(numout,*) ' mixed layer damping option nn_zdmp_tr = ', nn_zdmp_tr, '(zoom: forced to 0)' 188 WRITE(numout,*) ' surface time scale (days) rn_surf_tr = ', rn_surf_tr 189 WRITE(numout,*) ' bottom time scale (days) rn_bot_tr = ', rn_bot_tr 190 WRITE(numout,*) ' depth of transition (meters) rn_dep_tr = ', rn_dep_tr 191 WRITE(numout,*) ' create a damping.coeff file nn_file_tr = ', nn_file_tr 182 WRITE(numout,*) ' Restoration coeff file cn_resto_tr = ', cn_resto_tr 192 183 ENDIF 193 184 !
Note: See TracChangeset
for help on using the changeset viewer.