--- MITgcm/pkg/gmredi/gmredi_check.F	2001/02/02 21:36:29	1.3
+++ MITgcm/pkg/gmredi/gmredi_check.F	2003/01/13 19:02:45	1.8
@@ -1,4 +1,5 @@
-C $Header: /home/ubuntu/mnt/e9_copy/MITgcm/pkg/gmredi/gmredi_check.F,v 1.3 2001/02/02 21:36:29 adcroft Exp $
+C $Header: /home/ubuntu/mnt/e9_copy/MITgcm/pkg/gmredi/gmredi_check.F,v 1.8 2003/01/13 19:02:45 jmc Exp $
+C $Name:  $
 
 #include "CPP_OPTIONS.h"
 
@@ -39,6 +40,30 @@
         STOP 'ABNORMAL END: S/R GMREDI_CHECK'
        ENDIF
 
+C- print out some kee parameters :
+       CALL WRITE_0D_L( GM_AdvForm, INDEX_NONE,
+     &  'GM_AdvForm =', '  /* if FALSE => use SkewFlux Form */')     
+       CALL WRITE_0D_L( GM_AdvSeparate, INDEX_NONE,
+     & 'GM_AdvSeparate =',' /* Calc Bolus & Euler Adv. separately */')     
+       CALL WRITE_0D_L( GM_ExtraDiag, INDEX_NONE,
+     &  'GM_ExtraDiag =','  /* Tensor Extra Diag (line 1&2) non 0 */')
+       CALL WRITE_0D_R8( GM_isopycK, INDEX_NONE,'GM_isopycK =',
+     &  '   /* Background Isopyc. Diffusivity ( m^2/s ) */')       
+       CALL WRITE_0D_R8( GM_background_K*GM_skewflx, INDEX_NONE,
+     &  ' GM_skewflx*K =',
+     &  '   /* Background GM_SkewFlx Diffusivity ( m^2/s ) */')       
+       CALL WRITE_0D_R8( GM_background_K*GM_advect, INDEX_NONE,
+     &  ' GM_advec*K =',
+     &  '   /* Backg. GM-Advec(=Bolus) Diffusivity ( m^2/s ) */')       
+       CALL WRITE_0D_R8( GM_Visbeck_alpha, INDEX_NONE,
+     &  ' GM_Visbeck_alpha =','   /* Visbeck alpha coeff. ( ) */')       
+       WRITE(msgBuf,'(A,A40)')' Tapering/Cliping : ',GM_taper_scheme
+       CALL PRINT_MESSAGE(msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
+       CALL WRITE_0D_R8( GM_Small_Number, INDEX_NONE,
+     &  ' GM_Small_Number =','  /* epsilon used in slope calc */')       
+       CALL WRITE_0D_R8( GM_slopeSqCutoff, INDEX_NONE,
+     &  ' GM_slopeSqCutoff =', '  /* Slope^2 cut-off value */')       
+
       ENDIF
 
 #endif