New URL for NEMO forge!   http://forge.nemo-ocean.eu

Since March 2022 along with NEMO 4.2 release, the code development moved to a self-hosted GitLab.
This present forge is now archived and remained online for history.
Changeset 8088 – NEMO

Changeset 8088


Ignore:
Timestamp:
2017-05-30T09:47:54+02:00 (7 years ago)
Author:
cbricaud
Message:

compute r_fact_lap_crs for passive tracer diffusion in CRS mode

File:
1 edited

Legend:

Unmodified
Added
Removed
  • branches/2015/dev_r5003_MERCATOR6_CRS/NEMOGCM/NEMO/TOP_SRC/TRP/trcldf_crs.F90

    r8087 r8088  
    1818   USE trc             ! ocean passive tracers variables 
    1919   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 
    2122   USE ldfslp          ! ??? 
    2223   USE traldf_bilapg   ! lateral mixing            (tra_ldf_bilapg routine) 
     
    5859      INTEGER, INTENT( in ) ::   kt   ! ocean time-step index 
    5960      !! 
    60       INTEGER            :: jn 
     61      INTEGER            :: ji,jj,jk,jn 
     62      REAL(wp)           :: zdep 
    6163      CHARACTER (len=22) :: charout 
    6264      REAL(wp), POINTER, DIMENSION(:,:,:,:) ::   ztrtrd 
     
    6870 
    6971      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 
    7084 
    7185      IF( l_trdtrc )  THEN 
Note: See TracChangeset for help on using the changeset viewer.