Changeset 6403 for trunk/NEMOGCM/NEMO/TOP_SRC
- Timestamp:
- 2016-03-25T17:24:35+01:00 (8 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMOGCM/NEMO/TOP_SRC/TRP/trcldf.F90
r6140 r6403 40 40 REAL(wp), PUBLIC :: rn_ahtrc_0 !: laplacian diffusivity coefficient for passive tracer [m2/s] 41 41 REAL(wp), PUBLIC :: rn_bhtrc_0 !: bilaplacian - -- - - [m4/s] 42 REAL(wp), PUBLIC :: rn_fact_lap !: Enhanced zonal diffusivity coefficent in the equatorial domain 42 43 ! 43 44 ! !!: ** lateral mixing namelist (nam_trcldf) ** … … 64 65 INTEGER, INTENT( in ) :: kt ! ocean time-step index 65 66 ! 66 INTEGER :: jn 67 INTEGER :: ji, jj, jk, jn 68 REAL(wp) :: zdep 67 69 CHARACTER (len=22) :: charout 68 70 REAL(wp), POINTER, DIMENSION(:,:,:) :: zahu, zahv … … 76 78 ztrtrd(:,:,:,:) = tra(:,:,:,:) 77 79 ENDIF 78 ! 79 ! !* set the lateral diffusivity coef. for passive tracer 80 ! !* set the lateral diffusivity coef. for passive tracer 80 81 CALL wrk_alloc( jpi,jpj,jpk, zahu, zahv ) 81 zahu(:,:,:) = rldf * ahtu(:,:,:) 82 zahu(:,:,:) = rldf * ahtu(:,:,:) 82 83 zahv(:,:,:) = rldf * ahtv(:,:,:) 83 84 ! !* Enhanced zonal diffusivity coefficent in the equatorial domain 85 DO jk= 1, jpk 86 DO jj = 1, jpj 87 DO ji = 1, jpi 88 IF( gdept_n(ji,jj,jk) > 200. .AND. gphit(ji,jj) < 5. .AND. gphit(ji,jj) > -5. ) THEN 89 zdep = MAX( gdept_n(ji,jj,jk) - 1000., 0. ) / 1000. 90 zahu(ji,jj,jk) = zahu(ji,jj,jk) * MAX( 1., rn_fact_lap * EXP( -zdep ) ) 91 ENDIF 92 END DO 93 END DO 94 END DO 95 ! 84 96 SELECT CASE ( nldf ) !* compute lateral mixing trend and add it to the general trend 85 97 ! … … 136 148 NAMELIST/namtrc_ldf/ ln_trcldf_lap, ln_trcldf_blp, & 137 149 & ln_trcldf_lev, ln_trcldf_hor, ln_trcldf_iso, ln_trcldf_triad, & 138 & rn_ahtrc_0 , rn_bhtrc_0 150 & rn_ahtrc_0 , rn_bhtrc_0, rn_fact_lap 139 151 !!---------------------------------------------------------------------- 140 152 ! … … 164 176 WRITE(numout,*) ' laplacian rn_ahtrc_0 = ', rn_ahtrc_0 165 177 WRITE(numout,*) ' bilaplacian rn_bhtrc_0 = ', rn_bhtrc_0 178 WRITE(numout,*) ' enhanced zonal diffusivity rn_fact_lap = ', rn_fact_lap 179 166 180 ENDIF 167 181 !
Note: See TracChangeset
for help on using the changeset viewer.