/[MITgcm]/MITgcm/pkg/gmredi/gmredi_check.F
ViewVC logotype

Annotation of /MITgcm/pkg/gmredi/gmredi_check.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph

Revision 1.6 - (hide annotations) (download)
Wed Mar 6 01:56:27 2002 UTC (22 years, 6 months ago) by jmc
Branch: MAIN
CVS Tags: checkpoint46n_post, checkpoint46l_post, checkpoint46g_pre, checkpoint46f_post, checkpoint46b_post, checkpoint46l_pre, checkpoint46d_pre, checkpoint45d_post, checkpoint46j_pre, checkpoint44h_pre, checkpoint46a_post, checkpoint46j_post, checkpoint46k_post, checkpoint45a_post, checkpoint46e_pre, checkpoint45b_post, checkpoint46b_pre, checkpoint46c_pre, checkpoint46, checkpoint46h_pre, checkpoint46m_post, checkpoint46a_pre, checkpoint45c_post, checkpoint44h_post, checkpoint46g_post, checkpoint46i_post, checkpoint46c_post, checkpoint46e_post, checkpoint45, checkpoint46h_post, checkpoint46d_post
Changes since 1.5: +3 -1 lines 
o GM Advective form:
 Tracers are advected using the residual transport (= Eulerian + GM-bolus).
 parameter GM_AdvSeparate=T return to previous form (i.e. compute separately
  Eulerian and Bolus advection fluxes)
1 jmc 1.6 C $Header: /u/gcmpack/MITgcm/pkg/gmredi/gmredi_check.F,v 1.5 2001/12/16 18:54:49 jmc Exp $
2 jmc 1.5 C $Name: $
3 heimbach 1.1
4     #include "CPP_OPTIONS.h"
5    
6     SUBROUTINE GMREDI_CHECK( myThid )
7     C /==========================================================\
8     C | SUBROUTINE GMREDI_CHECK |
9     C | o Check dependances with other packages |
10     C |==========================================================|
11     C \==========================================================/
12     IMPLICIT NONE
13    
14     C === Global variables ===
15     #include "SIZE.h"
16     #include "EEPARAMS.h"
17     #include "PARAMS.h"
18     #include "GMREDI.h"
19    
20     C === Routine arguments ===
21     C myThid - Number of this instances
22     INTEGER myThid
23    
24     C === Local variables ===
25     C msgBuf - Informational/error meesage buffer
26     CHARACTER*(MAX_LEN_MBUF) msgBuf
27    
28    
29     #ifdef ALLOW_GMREDI
30 adcroft 1.3 IF (GMRediIsOn) THEN
31 heimbach 1.1
32 adcroft 1.3 WRITE(msgBuf,'(A)') 'GMREDI_CHECK: #define GMREDI'
33     CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
34     & SQUEEZE_RIGHT , 1)
35 heimbach 1.1
36     C GM/Redi needs implicit diffusion (will be packaged later)
37     IF (.NOT.implicitDiffusion) THEN
38     WRITE(msgBuf,'(A)') 'GM/Redi needs implicitDiffusion=.true.'
39     CALL PRINT_ERROR( msgBuf , 1)
40 adcroft 1.3 STOP 'ABNORMAL END: S/R GMREDI_CHECK'
41 heimbach 1.1 ENDIF
42 jmc 1.5
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 jmc 1.6 CALL WRITE_0D_L( GM_AdvSeparate, INDEX_NONE,
47     & 'GM_AdvSeparate =',' /* Calc Bolus & Euler Adv. separately */')
48 jmc 1.5 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 heimbach 1.1
63     ENDIF
64    
65     #endif
66 adcroft 1.3 RETURN
67     END
  ViewVC Help
Powered by ViewVC 1.1.22