/[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.5 - (hide annotations) (download)
Sun Dec 16 18:54:49 2001 UTC (22 years, 6 months ago) by jmc
Branch: MAIN
CVS Tags: checkpoint44e_post, checkpoint44f_post, checkpoint43a-release1mods, chkpt44d_post, checkpoint44e_pre, release1-branch_tutorials, chkpt44a_post, chkpt44c_pre, checkpoint44g_post, release1-branch-end, release1_final_v1, checkpoint44b_post, chkpt44a_pre, checkpoint44b_pre, checkpoint44, chkpt44c_post, checkpoint44f_pre, release1-branch_branchpoint
Branch point for: release1_final, release1-branch
Changes since 1.4: +20 -2 lines 
Modification to the GMREDI package :
 change units of tensor-K arrays, scale now like diffusivity
 initialise all common block arrays in S/R gmredi_init
 add option to use different isopycnal(Redi) & GM diffusivity
 add option to use the advective GM form or the skew-flux form (=default)
 bug in non_unity_diagonal part fixed.
1 jmc 1.5 C $Header: /u/gcmpack/models/MITgcmUV/pkg/gmredi/gmredi_check.F,v 1.4 2001/02/04 14:38:49 cnh Exp $
2     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     CALL WRITE_0D_L( GM_ExtraDiag, INDEX_NONE,
47     & 'GM_ExtraDiag =',' /* Tensor Extra Diag (line 1&2) non 0 */')
48     CALL WRITE_0D_R8( GM_isopycK, INDEX_NONE,'GM_isopycK =',
49     & ' /* Background Isopyc. Diffusivity ( m^2/s ) */')
50     CALL WRITE_0D_R8( GM_background_K*GM_skewflx, INDEX_NONE,
51     & ' GM_skewflx*K =',
52     & ' /* Background GM_SkewFlx Diffusivity ( m^2/s ) */')
53     CALL WRITE_0D_R8( GM_background_K*GM_advect, INDEX_NONE,
54     & ' GM_advec*K =',
55     & ' /* Backg. GM-Advec(=Bolus) Diffusivity ( m^2/s ) */')
56     CALL WRITE_0D_R8( GM_Visbeck_alpha, INDEX_NONE,
57     & ' GM_Visbeck_alpha =',' /* Visbeck alpha coeff. ( ) */')
58     WRITE(msgBuf,'(A,A40)')' Tapering/Cliping : ',GM_taper_scheme
59     CALL PRINT_MESSAGE(msgBuf,standardMessageUnit,SQUEEZE_RIGHT,1)
60 heimbach 1.1
61     ENDIF
62    
63     #endif
64 adcroft 1.3 RETURN
65     END
  ViewVC Help
Powered by ViewVC 1.1.22