/[MITgcm]/MITgcm/pkg/obcs/obcs_check.F
ViewVC logotype

Contents of /MITgcm/pkg/obcs/obcs_check.F

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


Revision 1.7 - (show annotations) (download)
Mon Oct 10 05:53:48 2005 UTC (18 years, 8 months ago) by mlosch
Branch: MAIN
CVS Tags: checkpint57u_post
Changes since 1.6: +12 -1 lines
o OBCS and PTRACERS: add open boundary support for passive tracers
  - either use homogenous (pseudo) v.Neumann conditions or prescribe
    OB-values from file; this is not different from the way theta and salinity
    are treated
  - however, Orlanski-radiation conditions are not supported, and the model
    will stop if you use pTracers and Orlanski at the same time.
  - beefed up the rountine obcs_external_fields_load: now only those open
    boundary values are overwritten with values from files for which there
    are really files, otherwise the OB-fields remain untouched. This makes
    it possible to use different OBs at different ends of the domain (as
    with EXF)
  - TODO: add support for OB?w and OB?eta, which can currently not be read
    from a file.

1 C $Header: /u/gcmpack/MITgcm/pkg/obcs/obcs_check.F,v 1.6 2004/09/24 12:32:09 mlosch Exp $
2 C $Name: $
3
4 #include "OBCS_OPTIONS.h"
5
6 SUBROUTINE OBCS_CHECK( myThid )
7 C /==========================================================\
8 C | SUBROUTINE OBCS_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 "OBCS.h"
19
20 C === Routine arguments ===
21 C myThid - Number of this instances
22 INTEGER myThid
23
24 #ifdef ALLOW_OBCS
25
26 C === Local variables ===
27 C msgBuf - Informational/error meesage buffer
28 CHARACTER*(MAX_LEN_MBUF) msgBuf
29
30 WRITE(msgBuf,'(A)') 'OBCS_CHECK: #define ALLOW_OBCS'
31 CALL PRINT_MESSAGE(msgBuf,standardMessageUnit,
32 & SQUEEZE_RIGHT,myThid)
33
34 #ifdef ALLOW_ORLANSKI
35 WRITE(msgBuf,'(A)') 'OBCS_CHECK: #define ALLOW_ORLANSKI'
36 CALL PRINT_MESSAGE(msgBuf,standardMessageUnit,
37 & SQUEEZE_RIGHT,myThid)
38 #else
39 IF (useOrlanskiNorth.OR.useOrlanskiSouth.OR.
40 & useOrlanskiEast.OR.useOrlanskiWest) THEN
41 WRITE(msgBuf,'(A)')
42 & 'OBCS_CHECK: ERROR: #undef OBCS_RADIATE_ORLANSKI and'
43 CALL PRINT_ERROR( msgBuf , 1)
44 WRITE(msgBuf,'(A)')
45 & 'OBCS_CHECK: ERROR: one of useOrlanski* logicals is true'
46 CALL PRINT_ERROR( msgBuf , 1)
47 STOP 'ABNORMAL END: S/R OBCS_CHECK'
48 ENDIF
49 #endif /* ALLOW_ORLANSKI */
50
51 IF (useOrlanskiNorth.OR.useOrlanskiSouth.OR.
52 & useOrlanskiEast.OR.useOrlanskiWest) THEN
53 IF (nonlinFreeSurf.GT.0) THEN
54 WRITE(msgBuf,'(A)')
55 & 'OBCS_CHECK: ERROR: useOrlanski* Rad OBC with'
56 CALL PRINT_ERROR( msgBuf , 1)
57 WRITE(msgBuf,'(A)')
58 & 'OBCS_CHECK: ERROR: nonlinFreeSurf not yet implemented'
59 CALL PRINT_ERROR( msgBuf , 1)
60 STOP 'ABNORMAL END: S/R OBCS_CHECK'
61 ENDIF
62 #ifdef ALLOW_PTRACERS
63 IF (usePTracers) THEN
64 WRITE(msgBuf,'(A)')
65 & 'OBCS_CHECK: ERROR: useOrlanski* Rad OBC with'
66 CALL PRINT_ERROR( msgBuf , 1)
67 WRITE(msgBuf,'(A)')
68 & 'OBCS_CHECK: ERROR: pTracers not yet implemented'
69 CALL PRINT_ERROR( msgBuf , 1)
70 STOP 'ABNORMAL END: S/R OBCS_CHECK'
71 ENDIF
72 #endif ALLOW_PTRACERS
73 ENDIF
74
75 #ifndef ALLOW_OBCS_PRESCRIBE
76 IF (useOBCSprescribe) THEN
77 WRITE(msgBuf,'(A)')
78 & 'OBCS_CHECK: ERROR: useOBCSprescribe = .TRUE. for'
79 CALL PRINT_ERROR( msgBuf , 1)
80 WRITE(msgBuf,'(A)')
81 & 'OBCS_CHECK: undef ALLOW_OBCS_PRESCRIBE'
82 STOP 'ABNORMAL END: S/R OBCS_CHECK'
83 ENDIF
84 #endif
85
86 WRITE(msgBuf,'(A)') 'OBCS_CHECK: OK'
87 CALL PRINT_MESSAGE(msgBuf,standardMessageUnit,
88 & SQUEEZE_RIGHT,myThid)
89
90 #endif /* ALLOW_OBCS */
91 return
92 end

  ViewVC Help
Powered by ViewVC 1.1.22