/[MITgcm]/MITgcm/pkg/shap_filt/shap_filt_readparms.F
ViewVC logotype

Diff of /MITgcm/pkg/shap_filt/shap_filt_readparms.F

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

revision 1.4 by adcroft, Fri Aug 17 19:09:23 2001 UTC revision 1.5 by jmc, Tue Dec 11 14:31:59 2001 UTC
# Line 22  C     === Routine arguments === Line 22  C     === Routine arguments ===
22    
23  #ifdef ALLOW_SHAP_FILT  #ifdef ALLOW_SHAP_FILT
24    
25        NAMELIST /SHAP_PARM01/ Shap_funct,        NAMELIST /SHAP_PARM01/
26       &          nShapT,  nShapTrPhys,       & Shap_funct, shap_filt_uvStar, shap_filt_TrStagg,
27       &          nShapUV, nShapUVPhys,       &    nShapT,  nShapTrPhys, Shap_Trtau, Shap_TrLength,
28       &          Shap_Trtau, Shap_TrLength,       &    nShapUV, nShapUVPhys, Shap_uvtau, Shap_uvLength
      &          Shap_uvtau, Shap_uvLength  
29    
30  C     === Local variables ===  C     === Local variables ===
31  C     msgBuf      - Informational/error meesage buffer  C     msgBuf      - Informational/error meesage buffer
# Line 51  c     SHAPIsOn=.TRUE. Line 50  c     SHAPIsOn=.TRUE.
50    
51  C--   Default flags and values for Shapiro Filter  C--   Default flags and values for Shapiro Filter
52        Shap_funct = 2        Shap_funct = 2
53          shap_filt_uvStar =.TRUE.
54          shap_filt_TrStagg=.TRUE.
55        nShapT = 0        nShapT = 0
56        nShapUV = 0        nShapUV = 0
57        nShapTrPhys = 0        nShapTrPhys = 0
# Line 78  C--   Close the open data file Line 79  C--   Close the open data file
79  C--   Everyone else must wait for the parameters to be loaded  C--   Everyone else must wait for the parameters to be loaded
80        _BARRIER        _BARRIER
81    
82    C--   Check the parameters :
83    
84           IF ( .NOT.shap_filt_uvStar ) THEN
85    
86    C- Notes: applying the filter at the end of the time step (after SOLVE_FOR_P)
87    C    affects the barotropic flow divergence ; this might not be consistent
88    C    with some option of the code.
89    
90            IF ( rigidLid ) THEN
91             WRITE(msgBuf,'(2A)') 'SHAP_FILT with rigidLid ',
92         &                         'needs shap_filt_uvStar=.true.'
93             CALL PRINT_ERROR( msgBuf , 1)
94             STOP 'ABNORMAL END: S/R SHAP_FILT_READPARMS'
95            ELSEIF ( .NOT.exactConserv ) THEN
96             WRITE(msgBuf,'(2A)') 'S/R SHAP_FILT_READPARMS: WARNING <<< ',
97         &    'applying Filter after SOLVE_FOR_P (shap_filt_uvStar=FALSE)'
98             CALL PRINT_ERROR( msgBuf , 1)
99             WRITE(msgBuf,'(2A)') 'S/R SHAP_FILT_READPARMS: WARNING <<< ',
100         &    'requires to recompute Eta after ==> turn on exactConserv '
101             CALL PRINT_ERROR( msgBuf , 1)
102            ENDIF
103    
104           ENDIF                    
105    
106  #endif /* ALLOW_SHAP_FILT */  #endif /* ALLOW_SHAP_FILT */
107        RETURN        RETURN
108        END        END
Legend:
Removed from v.1.4  
changed lines
  Added in v.1.5
  ViewVC Help
Powered by ViewVC 1.1.22