- Timestamp:
- 2017-05-30T09:44:02+02:00 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2015/dev_r5003_MERCATOR6_CRS/NEMOGCM/NEMO/TOP_SRC/TRP/trcldf_crs.F90
r8086 r8087 18 18 USE trc ! ocean passive tracers variables 19 19 USE trcnam_trp ! passive tracers transport namelist variables 20 !USE ldftra_oce ! lateral diffusion coefficient on tracers 21 USE ldftra_oce,ONLY: ln_traldf_grif,rn_aht_0,rn_ahtb_0,lk_traldf_eiv ! lateral diffusion coefficient on tracers 20 USE ldftra_oce ! lateral diffusion coefficient on tracers 22 21 USE ldfslp ! ??? 23 22 USE traldf_bilapg ! lateral mixing (tra_ldf_bilapg routine) … … 59 58 INTEGER, INTENT( in ) :: kt ! ocean time-step index 60 59 !! 61 INTEGER :: ji,jj,jk,jn 62 REAL(wp) :: zdep 60 INTEGER :: jn 63 61 CHARACTER (len=22) :: charout 64 62 REAL(wp), POINTER, DIMENSION(:,:,:,:) :: ztrtrd … … 70 68 71 69 rldf = rldf_rat 72 !73 r_fact_lap(:,:,:) = 1.74 DO jk= 1, jpk75 DO jj = 1, jpj76 DO ji = 1, jpi77 IF( fsdept(ji,jj,jk) > 200. .AND. gphit(ji,jj) < 5. .AND. gphit(ji,jj) > -5. ) THEN78 zdep = MAX( fsdept(ji,jj,jk) - 1000., 0. ) / 1000.79 r_fact_lap(ji,jj,jk) = MAX( 1., rn_fact_lap * EXP( -zdep ) )80 ENDIF81 END DO82 END DO83 END DO84 70 85 71 IF( l_trdtrc ) THEN
Note: See TracChangeset
for help on using the changeset viewer.