Changeset 8088
 Timestamp:
 20170530T09:47:54+02:00 (6 years ago)
 File:

 1 edited
Legend:
 Unmodified
 Added
 Removed

branches/2015/dev_r5003_MERCATOR6_CRS/NEMOGCM/NEMO/TOP_SRC/TRP/trcldf_crs.F90
r8087 r8088 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 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 21 22 USE ldfslp ! ??? 22 23 USE traldf_bilapg ! lateral mixing (tra_ldf_bilapg routine) … … 58 59 INTEGER, INTENT( in ) :: kt ! ocean timestep index 59 60 !! 60 INTEGER :: jn 61 INTEGER :: ji,jj,jk,jn 62 REAL(wp) :: zdep 61 63 CHARACTER (len=22) :: charout 62 64 REAL(wp), POINTER, DIMENSION(:,:,:,:) :: ztrtrd … … 68 70 69 71 rldf = rldf_rat 72 ! 73 r_fact_lap(:,:,:) = 1. 74 DO jk= 1, jpk 75 DO jj = 1, jpj 76 DO ji = 1, jpi 77 IF( fsdept(ji,jj,jk) > 200. .AND. gphit(ji,jj) < 5. .AND. gphit(ji,jj) > 5. ) THEN 78 zdep = MAX( fsdept(ji,jj,jk)  1000., 0. ) / 1000. 79 r_fact_lap(ji,jj,jk) = MAX( 1., rn_fact_lap * EXP( zdep ) ) 80 ENDIF 81 END DO 82 END DO 83 END DO 70 84 71 85 IF( l_trdtrc ) THEN
Note: See TracChangeset
for help on using the changeset viewer.