Changeset 1537 for trunk/NEMO/OPA_SRC/TRA/tranpc.F90
- Timestamp:
- 2009-07-24T16:22:23+02:00 (15 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMO/OPA_SRC/TRA/tranpc.F90
r1146 r1537 4 4 !! Ocean active tracers: non penetrative convection scheme 5 5 !!============================================================================== 6 !! History : 1.0 ! 90-09 (G. Madec) Original code 7 !! ! 91-11 (G. Madec) 8 !! ! 92-06 (M. Imbard) periodic conditions on t and s 9 !! ! 93-03 (M. Guyon) symetrical conditions 10 !! ! 96-01 (G. Madec) statement function for e3 11 !! suppression of common work arrays 12 !! 8.5 ! 02-06 (G. Madec) free form F90 13 !! 9.0 ! 08-06 (G. Madec) applied on ta, sa and called before tranxt in step.F90 14 !!---------------------------------------------------------------------- 15 16 !!---------------------------------------------------------------------- 17 !! tra_npc : apply the non penetrative convection scheme 18 !! tra_npc_init : initialization and control of the scheme 6 !! History : 1.0 ! 1990-09 (G. Madec) Original code 7 !! ! 1996-01 (G. Madec) statement function for e3 8 !! NEMO 1.0 ! 2002-06 (G. Madec) free form F90 9 !! 3.0 ! 2008-06 (G. Madec) applied on ta, sa and called before tranxt in step.F90 10 !!---------------------------------------------------------------------- 11 12 !!---------------------------------------------------------------------- 13 !! tra_npc : apply the non penetrative convection scheme 19 14 !!---------------------------------------------------------------------- 20 15 USE oce ! ocean dynamics and active tracers 21 16 USE dom_oce ! ocean space and time domain 17 USE zdf_oce ! ocean vertical physics 22 18 USE trdmod ! ocean active tracer trends 23 19 USE trdmod_oce ! ocean variables trends … … 31 27 PUBLIC tra_npc ! routine called by step.F90 32 28 33 !!* Namelist namnpc: non penetrative convection algorithm34 INTEGER :: nnpc1 = 1 ! nnpc1 non penetrative convective scheme frequency35 INTEGER :: nnpc2 = 15 ! nnpc2 non penetrative convective scheme print frequency36 37 29 !! * Substitutions 38 30 # include "domzgr_substitute.h90" 39 31 !!---------------------------------------------------------------------- 40 !! OPA 9.0 , LOCEAN-IPSL (2005)32 !! NEMO/OPA 3.2 , LOCEAN-IPSL (2009) 41 33 !! $Id$ 42 34 !! Software governed by the CeCILL licence (modipsl/doc/NEMO_CeCILL.txt) … … 78 70 !!---------------------------------------------------------------------- 79 71 80 IF( kt == nit000 ) CALL tra_npc_init ! Initialisation 81 82 IF( MOD( kt, nnpc1 ) == 0 ) THEN 72 IF( MOD( kt, nn_npc ) == 0 ) THEN 83 73 84 74 inpcc = 0 … … 214 204 ! 2. non penetrative convective scheme statistics 215 205 ! ----------------------------------------------- 216 IF( nn pc2 /= 0 .AND. MOD( kt, nnpc2) == 0 ) THEN206 IF( nn_npcp /= 0 .AND. MOD( kt, nn_npcp ) == 0 ) THEN 217 207 IF(lwp) WRITE(numout,*)' kt=',kt, ' number of statically instable', & 218 208 & ' water column : ',inpcc, ' number of iteration : ',inpci … … 223 213 END SUBROUTINE tra_npc 224 214 225 226 SUBROUTINE tra_npc_init227 !!----------------------------------------------------------------------228 !! *** ROUTINE tra_npc_init ***229 !!230 !! ** Purpose : initializations of the non-penetrative adjustment scheme231 !!----------------------------------------------------------------------232 NAMELIST/namnpc/ nnpc1, nnpc2233 !234 REWIND( numnam ) ! Namelist namzdf : vertical diffusion235 READ ( numnam, namnpc )236 !237 IF(lwp) THEN ! Namelist print238 WRITE(numout,*)239 WRITE(numout,*) 'tra_npc_init : Non Penetrative Convection (npc) scheme'240 WRITE(numout,*) '~~~~~~~~~~~~'241 WRITE(numout,*) ' Namelist namnpc : set npc scheme parameters'242 WRITE(numout,*) ' npc scheme frequency nnpc1 = ', nnpc1243 WRITE(numout,*) ' npc scheme print frequency nnpc2 = ', nnpc2244 ENDIF245 !246 IF ( nnpc1 == 0 ) THEN ! Parameter controls247 IF(lwp) WRITE(numout,cform_war)248 IF(lwp) WRITE(numout,*) ' nnpc1 = ', nnpc1, ' is forced to 1'249 nnpc1 = 1250 nwarn = nwarn + 1251 ENDIF252 !253 END SUBROUTINE tra_npc_init254 255 215 !!====================================================================== 256 216 END MODULE tranpc
Note: See TracChangeset
for help on using the changeset viewer.