1 |
C $Header$ |
C $Header$ |
2 |
|
C $Name$ |
3 |
|
|
4 |
#include "CPP_OPTIONS.h" |
#include "GMREDI_OPTIONS.h" |
5 |
|
|
6 |
SUBROUTINE GMREDI_CHECK( myThid ) |
SUBROUTINE GMREDI_CHECK( myThid ) |
7 |
C /==========================================================\ |
C /==========================================================\ |
40 |
STOP 'ABNORMAL END: S/R GMREDI_CHECK' |
STOP 'ABNORMAL END: S/R GMREDI_CHECK' |
41 |
ENDIF |
ENDIF |
42 |
|
|
43 |
|
C- print out some kee parameters : |
44 |
|
CALL WRITE_0D_L( GM_AdvForm, INDEX_NONE, |
45 |
|
& 'GM_AdvForm =', ' /* if FALSE => use SkewFlux Form */') |
46 |
|
CALL WRITE_0D_L( GM_AdvSeparate, INDEX_NONE, |
47 |
|
& 'GM_AdvSeparate =',' /* Calc Bolus & Euler Adv. separately */') |
48 |
|
CALL WRITE_0D_L( GM_ExtraDiag, INDEX_NONE, |
49 |
|
& 'GM_ExtraDiag =',' /* Tensor Extra Diag (line 1&2) non 0 */') |
50 |
|
CALL WRITE_0D_R8( GM_isopycK, INDEX_NONE,'GM_isopycK =', |
51 |
|
& ' /* Background Isopyc. Diffusivity ( m^2/s ) */') |
52 |
|
CALL WRITE_0D_R8( GM_background_K*GM_skewflx, INDEX_NONE, |
53 |
|
& ' GM_skewflx*K =', |
54 |
|
& ' /* Background GM_SkewFlx Diffusivity ( m^2/s ) */') |
55 |
|
CALL WRITE_0D_R8( GM_background_K*GM_advect, INDEX_NONE, |
56 |
|
& ' GM_advec*K =', |
57 |
|
& ' /* Backg. GM-Advec(=Bolus) Diffusivity ( m^2/s ) */') |
58 |
|
CALL WRITE_0D_R8( GM_Visbeck_alpha, INDEX_NONE, |
59 |
|
& ' GM_Visbeck_alpha =',' /* Visbeck alpha coeff. ( ) */') |
60 |
|
WRITE(msgBuf,'(A,A40)')' Tapering/Cliping : ',GM_taper_scheme |
61 |
|
CALL PRINT_MESSAGE(msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1) |
62 |
|
CALL WRITE_0D_R8( GM_Small_Number, INDEX_NONE, |
63 |
|
& ' GM_Small_Number =',' /* epsilon used in slope calc */') |
64 |
|
CALL WRITE_0D_R8( GM_slopeSqCutoff, INDEX_NONE, |
65 |
|
& ' GM_slopeSqCutoff =', ' /* Slope^2 cut-off value */') |
66 |
|
|
67 |
ENDIF |
ENDIF |
68 |
|
|
69 |
#endif |
#endif |