pkg/exf/exf_check.F

C $Header: /u/gcmpack/MITgcm/pkg/exf/exf_check.F,v 1.9 2008/01/25 11:49:48 dimitri Exp $
C $Name:  $

#include "EXF_OPTIONS.h"

      subroutine exf_check( mythid )

c     ==================================================================
c     SUBROUTINE exf_check
c     ==================================================================
c
      implicit none

c     == global variables ==

#include "EEPARAMS.h"
#include "SIZE.h"
#include "PARAMS.h"
c#include "FFIELDS.h"
c#include "GRID.h"

#include "EXF_PARAM.h"
#include "EXF_CONSTANTS.h"
c#include "EXF_FIELDS.h"
c     == routine arguments ==

c     mythid - thread number for this instance of the routine.

      integer mythid

c     == local variables ==

c     == end of interface ==

c     check for consistency
      if (.NOT. 
     &     (exf_iprec .EQ. 32 .OR. exf_iprec .EQ. 64)
     &     ) then
       stop 'S/R EXF_CHECK: value of exf_iprec not allowed'
      elseif (.NOT. 
     &      (exf_yftype .EQ. 'RS' .OR. 
     &      exf_yftype .EQ. 'RL')
     &      ) then
       stop 'S/R EXF_CHECK: value of exf_yftype not allowed'
      endif

      if (repeatPeriod.lt.0.) then
       stop 'S/R EXF_CHECK: repeatPeriod must be positive'
      endif

      if (useExfYearlyFields.and.repeatPeriod.ne.0.) then
       print*,'Use of usefldyearlyfields AND repeatPeriod',
     &      ' not implemented'
       stop 'ABNORMAL END: S/R EXF_CHECK'
      endif

#if ( defined (ALLOW_BULKFORMULAE) && defined (ALLOW_ATM_WIND) )
      IF ( ustressfile .NE. ' ' .OR. vstressfile .NE. ' ' ) THEN
        STOP
     & 'S/R EXF_CHECK: use u,v_wind components but not wind-stress'
      ENDIF
#endif

#ifndef ALLOW_ATM_WIND
      IF ( uwindfile .NE. ' ' .OR. vwindfile .NE. ' ' ) THEN
        STOP
     & 'S/R EXF_CHECK: read-in wind-stress but not u,v_wind components'
      ENDIF
#endif

#ifdef USE_EXF_INTERPOLATION
      if ( climsst_nlat .GT. MAX_LAT_INC )
     & stop 'stopped in exf_readparms: climsst_nlat > MAX_LAT_INC'
      if ( climsss_nlat .GT. MAX_LAT_INC )
     & stop 'stopped in exf_readparms: climsss_nlat > MAX_LAT_INC'
      if ( usingCartesianGrid ) then
       print*,'USE_EXF_INTERPOLATION assumes latitude/longitude'
       print*,'input and output coordinates.  Trivial to extend to'
       print*,'cartesian coordinates, but has not yet been done.'
       stop
      endif
      IF ( readStressOnAgrid.OR.readStressOnCgrid ) THEN
        STOP 'wind-stress position defined by Long/Lat'
      ENDIF
      if ( usingCurvilinearGrid ) then
c     some restrictions that can be relaxed later on
       if ( ustressfile .NE. ' ' .and. vstressfile .NE. ' ' ) then
        if ( ustressstartdate .ne. vstressstartdate .or.
     &       ustressperiod    .ne. vstressperiod ) then
         print*,'For useCubedSphereExchange, S/R exf_set_uv.F'
         print*,'assumes that the u and v wind stress files'
         print*,'have the same startdate and period.'
         stop
        endif
       endif
       if ( uwindfile .NE. ' ' .and. vwindfile .NE. ' ' ) then
        if ( uwindstartdate .ne. vwindstartdate .or.
     &       uwindperiod    .ne. vwindperiod ) then
         print*,'For useCubedSphereExchange, S/R exf_set_uv.F'
         print*,'assumes that the u and v wind files'
         print*,'have the same startdate and period.'
         stop
        endif
       endif
      endif
#else /* ifndef USE_EXF_INTERPOLATION */
      IF ( ustressfile .NE. ' ' .OR. vstressfile .NE. ' ' ) THEN
        IF ( (readStressOnAgrid.AND.readStressOnCgrid) .OR.
     &   .NOT.(readStressOnAgrid.OR.readStressOnCgrid) ) THEN
          STOP
     & 'S/R EXF_CHECK: Select 1 wind-stress position: A or C-grid'
        ENDIF
      ELSE
        IF ( readStressOnAgrid .OR. readStressOnCgrid ) THEN
          STOP 'S/R EXF_CHECK: wind-stress position irrelevant'
        ENDIF
      ENDIF
#endif /* USE_EXF_INTERPOLATION */

#ifdef ALLOW_CLIMTEMP_RELAXATION
      STOP 'ALLOW_CLIMTEMP_RELAXATION no longer supported. Use pkg/rbcs'
#endif

#ifdef ALLOW_CLIMSALT_RELAXATION
      STOP 'ALLOW_CLIMSALT_RELAXATION no longer supported. Use pkg/rbcs'
#endif

      RETURN
      END
1	C $Header: /u/gcmpack/MITgcm/pkg/exf/exf_check.F,v 1.9 2008/01/25 11:49:48 dimitri Exp $
2	C $Name: $
3
4	#include "EXF_OPTIONS.h"
5
6	subroutine exf_check( mythid )
7
8	c ==================================================================
9	c SUBROUTINE exf_check
10	c ==================================================================
11	c
12	implicit none
13
14	c == global variables ==
15
16	#include "EEPARAMS.h"
17	#include "SIZE.h"
18	#include "PARAMS.h"
19	c#include "FFIELDS.h"
20	c#include "GRID.h"
21
22	#include "EXF_PARAM.h"
23	#include "EXF_CONSTANTS.h"
24	c#include "EXF_FIELDS.h"
25	c == routine arguments ==
26
27	c mythid - thread number for this instance of the routine.
28
29	integer mythid
30
31	c == local variables ==
32
33	c == end of interface ==
34
35	c check for consistency
36	if (.NOT.
37	& (exf_iprec .EQ. 32 .OR. exf_iprec .EQ. 64)
38	& ) then
39	stop 'S/R EXF_CHECK: value of exf_iprec not allowed'
40	elseif (.NOT.
41	& (exf_yftype .EQ. 'RS' .OR.
42	& exf_yftype .EQ. 'RL')
43	& ) then
44	stop 'S/R EXF_CHECK: value of exf_yftype not allowed'
45	endif
46
47	if (repeatPeriod.lt.0.) then
48	stop 'S/R EXF_CHECK: repeatPeriod must be positive'
49	endif
50
51	if (useExfYearlyFields.and.repeatPeriod.ne.0.) then
52	print*,'Use of usefldyearlyfields AND repeatPeriod',
53	& ' not implemented'
54	stop 'ABNORMAL END: S/R EXF_CHECK'
55	endif
56
57	#if ( defined (ALLOW_BULKFORMULAE) && defined (ALLOW_ATM_WIND) )
58	IF ( ustressfile .NE. ' ' .OR. vstressfile .NE. ' ' ) THEN
59	STOP
60	& 'S/R EXF_CHECK: use u,v_wind components but not wind-stress'
61	ENDIF
62	#endif
63
64	#ifndef ALLOW_ATM_WIND
65	IF ( uwindfile .NE. ' ' .OR. vwindfile .NE. ' ' ) THEN
66	STOP
67	& 'S/R EXF_CHECK: read-in wind-stress but not u,v_wind components'
68	ENDIF
69	#endif
70
71	#ifdef USE_EXF_INTERPOLATION
72	if ( climsst_nlat .GT. MAX_LAT_INC )
73	& stop 'stopped in exf_readparms: climsst_nlat > MAX_LAT_INC'
74	if ( climsss_nlat .GT. MAX_LAT_INC )
75	& stop 'stopped in exf_readparms: climsss_nlat > MAX_LAT_INC'
76	if ( usingCartesianGrid ) then
77	print*,'USE_EXF_INTERPOLATION assumes latitude/longitude'
78	print*,'input and output coordinates. Trivial to extend to'
79	print*,'cartesian coordinates, but has not yet been done.'
80	stop
81	endif
82	IF ( readStressOnAgrid.OR.readStressOnCgrid ) THEN
83	STOP 'wind-stress position defined by Long/Lat'
84	ENDIF
85	if ( usingCurvilinearGrid ) then
86	c some restrictions that can be relaxed later on
87	if ( ustressfile .NE. ' ' .and. vstressfile .NE. ' ' ) then
88	if ( ustressstartdate .ne. vstressstartdate .or.
89	& ustressperiod .ne. vstressperiod ) then
90	print*,'For useCubedSphereExchange, S/R exf_set_uv.F'
91	print*,'assumes that the u and v wind stress files'
92	print*,'have the same startdate and period.'
93	stop
94	endif
95	endif
96	if ( uwindfile .NE. ' ' .and. vwindfile .NE. ' ' ) then
97	if ( uwindstartdate .ne. vwindstartdate .or.
98	& uwindperiod .ne. vwindperiod ) then
99	print*,'For useCubedSphereExchange, S/R exf_set_uv.F'
100	print*,'assumes that the u and v wind files'
101	print*,'have the same startdate and period.'
102	stop
103	endif
104	endif
105	endif
106	#else /* ifndef USE_EXF_INTERPOLATION */
107	IF ( ustressfile .NE. ' ' .OR. vstressfile .NE. ' ' ) THEN
108	IF ( (readStressOnAgrid.AND.readStressOnCgrid) .OR.
109	& .NOT.(readStressOnAgrid.OR.readStressOnCgrid) ) THEN
110	STOP
111	& 'S/R EXF_CHECK: Select 1 wind-stress position: A or C-grid'
112	ENDIF
113	ELSE
114	IF ( readStressOnAgrid .OR. readStressOnCgrid ) THEN
115	STOP 'S/R EXF_CHECK: wind-stress position irrelevant'
116	ENDIF
117	ENDIF
118	#endif /* USE_EXF_INTERPOLATION */
119
120	#ifdef ALLOW_CLIMTEMP_RELAXATION
121	STOP 'ALLOW_CLIMTEMP_RELAXATION no longer supported. Use pkg/rbcs'
122	#endif
123
124	#ifdef ALLOW_CLIMSALT_RELAXATION
125	STOP 'ALLOW_CLIMSALT_RELAXATION no longer supported. Use pkg/rbcs'
126	#endif
127
128	RETURN
129	END