Changeset 3497 for branches/2012/dev_r3452_NOCL02_Smagorinsky/NEMOGCM/NEMO/OPA_SRC/LDF/ldftra.F90

Timestamp:

2012-10-12T18:14:36+02:00 (12 years ago)

Author:

hliu

Message:

Upload code files for Smagorinsky viscosity/diffusivity work, documentation to be converted to latex format and uploaded soon. See the ticket (coming in minutes)

File:

• branches/2012/dev_r3452_NOCL02_Smagorinsky/NEMOGCM/NEMO/OPA_SRC/LDF/ldftra.F90 (modified) (2 diffs)

branches/2012/dev_r3452_NOCL02_Smagorinsky/NEMOGCM/NEMO/OPA_SRC/LDF/ldftra.F90

@@ -66 +66 @@
       NAMELIST/namtra_ldf/ ln_traldf_lap  , ln_traldf_bilap,                  &
          &                 ln_traldf_level, ln_traldf_hor  , ln_traldf_iso,   &
-         &                 ln_traldf_grif , ln_traldf_gdia,                   &
-         &                 ln_triad_iso   , ln_botmix_grif,                   &
+         &                 ln_traldf_grif , ln_traldf_gdia ,                  &
+         &                 ln_triad_iso   , ln_botmix_grif ,                  &
          &                 rn_aht_0       , rn_ahtb_0      , rn_aeiv_0,       &
-         &                 rn_slpmax
+         &                 rn_slpmax      , rn_chsmag      ,    rn_smsh,      &
+         &                 rn_aht_m
       !!----------------------------------------------------------------------
 
@@ -153 +154 @@
       IF(lwp)WRITE(numout,*) '      constant eddy diffusivity coef.   ahtu = ahtv = ahtw = aht0 = ', aht0
       IF( lk_traldf_eiv ) THEN
-         IF(lwp)WRITE(numout,*)
          IF(lwp)WRITE(numout,*) '      constant eddy induced velocity coef.   aeiu = aeiv = aeiw = aeiv0 = ', aeiv0
+      
       ENDIF
 #endif
+
+#if defined key_traldf_smag && ! defined key_traldf_c3d
+        CALL ctl_stop( 'key_traldf_smag can only be used with key_traldf_c3d' )
+#endif
+#if defined key_traldf_smag
+        IF(lwp) WRITE(numout,*)' SMAGORINSKY DIFFUSION'
+        IF(lwp .AND. rn_smsh < 1)  WRITE(numout,*)' only  shear is used '
+        IF(lwp.and.ln_traldf_bilap) CALL ctl_stop(' SMAGORINSKY + BILAPLACIAN - UNSTABLE OR NON_CONSERVATIVE' )
+#endif
+
       !
    END SUBROUTINE ldf_tra_init