dgoldberg/streamice_oad_files/streamice_vel_phi.F

C $Header: /u/gcmpack/MITgcm/pkg/streamice/streamice_vel_solve.F,v 1.9 2014/09/09 23:09:21 jmc Exp $
C $Name:  $

#include "STREAMICE_OPTIONS.h"
#ifdef ALLOW_AUTODIFF
# include "AUTODIFF_OPTIONS.h"
#endif

C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|

CBOP
      SUBROUTINE STREAMICE_VEL_PHI ( 
     I           myThid, 
     I           maxNLIter, 
     I           maxCGiter, 
     I           cgtol, 
     O           cg_iters)
C     /============================================================\
C     | SUBROUTINE                                                 |
C     | o                                                          |
C     |============================================================|
C     |                                                            |
C     \============================================================/
      IMPLICIT NONE

C     === Global variables ===
#include "SIZE.h"
#include "EEPARAMS.h"
#include "PARAMS.h"
#include "STREAMICE.h"
#include "STREAMICE_CG.h"
!#ifdef ALLOW_PETSC
!#include "finclude/petsc.h"
!#endif

#ifdef ALLOW_AUTODIFF_TAMC
# include "tamc.h"
#endif

C     !INPUT/OUTPUT ARGUMENTS
      INTEGER myThid
      INTEGER maxNLIter
      INTEGER maxCGIter
      _RL cgtol
      INTEGER cg_iters

#ifdef ALLOW_STREAMICE

C     LOCAL VARIABLES

      CHARACTER*(MAX_LEN_MBUF) msgBuf

#ifdef STREAMICE_HYBRID_STRESS
       CALL STREAMICE_VISC_BETA_HYBRID ( myThid )
#else
       CALL STREAMICE_VISC_BETA ( myThid )
#endif

      _EXCH_XY_RL( tau_beta_eff_streamice , myThid )
      _EXCH_XY_RL( visc_streamice , myThid )

       CALL STREAMICE_CG_WRAPPER(
     &       U_streamice,
     &       V_streamice,
     &       taudx_SI,
     &       taudy_SI,
     &       cgtol,
     &       cg_iters,
     &       maxCGIter,
     &       myThid )

#ifdef STREAMICE_HYBRID_STRESS
        CALL STREAMICE_TAUB (myThid)
#endif

#endif
      RETURN
      END
1	dgoldberg	1.1	C $Header: /u/gcmpack/MITgcm/pkg/streamice/streamice_vel_solve.F,v 1.9 2014/09/09 23:09:21 jmc Exp $
2			C $Name: $
3
4			#include "STREAMICE_OPTIONS.h"
5			#ifdef ALLOW_AUTODIFF
6			# include "AUTODIFF_OPTIONS.h"
7			#endif
8
9			C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-\|--+----\|
10
11			CBOP
12			SUBROUTINE STREAMICE_VEL_PHI (
13			I myThid,
14			I maxNLIter,
15			I maxCGiter,
16			I cgtol,
17			O cg_iters)
18			C /============================================================\
19			C \| SUBROUTINE \|
20			C \| o \|
21			C \|============================================================\|
22			C \| \|
23			C \============================================================/
24			IMPLICIT NONE
25
26			C === Global variables ===
27			#include "SIZE.h"
28			#include "EEPARAMS.h"
29			#include "PARAMS.h"
30			#include "STREAMICE.h"
31			#include "STREAMICE_CG.h"
32			!#ifdef ALLOW_PETSC
33			!#include "finclude/petsc.h"
34			!#endif
35
36			#ifdef ALLOW_AUTODIFF_TAMC
37			# include "tamc.h"
38			#endif
39
40			C !INPUT/OUTPUT ARGUMENTS
41			INTEGER myThid
42			INTEGER maxNLIter
43			INTEGER maxCGIter
44			_RL cgtol
45			INTEGER cg_iters
46
47			#ifdef ALLOW_STREAMICE
48
49			C LOCAL VARIABLES
50
51			CHARACTER*(MAX_LEN_MBUF) msgBuf
52
53			#ifdef STREAMICE_HYBRID_STRESS
54			CALL STREAMICE_VISC_BETA_HYBRID ( myThid )
55			#else
56			CALL STREAMICE_VISC_BETA ( myThid )
57			#endif
58
59			_EXCH_XY_RL( tau_beta_eff_streamice , myThid )
60			_EXCH_XY_RL( visc_streamice , myThid )
61
62			CALL STREAMICE_CG_WRAPPER(
63			& U_streamice,
64			& V_streamice,
65			& taudx_SI,
66			& taudy_SI,
67			& cgtol,
68			& cg_iters,
69			& maxCGIter,
70			& myThid )
71
72			#ifdef STREAMICE_HYBRID_STRESS
73			CALL STREAMICE_TAUB (myThid)
74			#endif
75
76			#endif
77			RETURN
78			END