- Timestamp:
- 2017-04-19T12:15:07+02:00 (7 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
branches/2017/dev_r7881_HPC09_ZDF/NEMOGCM/NEMO/OPA_SRC/TRA/trazdf_imp.F90
r7753 r7931 44 44 45 45 !! * Substitutions 46 # include "zdfddm_substitute.h90"47 46 # include "vectopt_loop_substitute.h90" 48 47 !!---------------------------------------------------------------------- … … 66 65 !! It is computed using a backward time scheme (t=after field) 67 66 !! which provide directly the after tracer field. 68 !! If l k_zdfddm=T, use avs for salinity or for passive tracers67 !! If ln_zdfddm=T, use avs for salinity or for passive tracers 69 68 !! Surface and bottom boundary conditions: no diffusive flux on 70 69 !! both tracers (bottom, applied through the masked field avt). … … 102 101 ! Build matrix if temperature or salinity (only in double diffusion case) or first passive tracer 103 102 ! 104 IF( ( cdtype == 'TRA' .AND. ( jn == jp_tem .OR. ( jn == jp_sal .AND. l k_zdfddm ) ) ) .OR. &103 IF( ( cdtype == 'TRA' .AND. ( jn == jp_tem .OR. ( jn == jp_sal .AND. ln_zdfddm ) ) ) .OR. & 105 104 & ( cdtype == 'TRC' .AND. jn == 1 ) ) THEN 106 105 ! 107 106 ! vertical mixing coef.: avt for temperature, avs for salinity and passive tracers 108 IF( cdtype == 'TRA' .AND. jn == jp_tem ) THEN ; zwt(:,:,2:jpk) = avt 109 ELSE ; zwt(:,:,2:jpk) = fsavs(:,:,2:jpk)107 IF( cdtype == 'TRA' .AND. jn == jp_tem ) THEN ; zwt(:,:,2:jpk) = avt(:,:,2:jpk) 108 ELSE ; zwt(:,:,2:jpk) = avs(:,:,2:jpk) 110 109 ENDIF 111 110 zwt(:,:,1) = 0._wp
Note: See TracChangeset
for help on using the changeset viewer.