model/src/reset_nlfs_vars.F

C $Header: /u/gcmpack/MITgcm/model/src/update_r_star.F,v 1.9 2011/05/20 01:09:37 jmc Exp $
C $Name:  $

#include "CPP_OPTIONS.h"

CBOP
C     !ROUTINE: RESET_NLFS_VARS
C     !INTERFACE:
      SUBROUTINE RESET_NLFS_VARS( myTime, myIter, myThid )

C     !DESCRIPTION: \bv
C     *==========================================================*
C     | SUBROUTINE RESET_NLFS_VARS
C     | o Re-set some Non-Linear Free-Surface variables
C     |   in order to facilitate the AD tool task of solving
C     |   dependency rules.
C     *==========================================================*
C     \ev

C     !USES:
      IMPLICIT NONE
C     == Global variables
#include "SIZE.h"
#include "EEPARAMS.h"
#include "PARAMS.h"
#include "GRID.h"
#include "SURFACE.h"

C     !INPUT/OUTPUT PARAMETERS:
C     == Routine arguments ==
C     myTime    :: Current time in simulation
C     myIter    :: Current iteration number in simulation
C     myThid    :: my Thread Id number
      _RL myTime
      INTEGER myIter
      INTEGER myThid

C     !LOCAL VARIABLES:
#ifdef NONLIN_FRSURF
C     Local variables
C     i,j,k,bi,bj :: loop counter
      INTEGER i,j,bi,bj
CEOP

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

      DO bj=myByLo(myThid), myByHi(myThid)
       DO bi=myBxLo(myThid), myBxHi(myThid)

#ifndef DISABLE_RSTAR_CODE
        IF ( fluidIsAir ) THEN
         DO j=1-OLy,sNy+OLy
          DO i=1-OLx,sNx+OLx
           pStarFacK(i,j,bi,bj) = rStarFacC(i,j,bi,bj)**atm_kappa
          ENDDO
         ENDDO
        ELSE
         DO j=1-OLy,sNy+OLy
          DO i=1-OLx,sNx+OLx
           pStarFacK(i,j,bi,bj) = 1. _d 0
          ENDDO
         ENDDO
        ENDIF
#endif /* DISABLE_RSTAR_CODE */

C- end bi,bj loop
       ENDDO
      ENDDO

#endif /* NONLIN_FRSURF */

      RETURN
      END
1	jmc	1.1	C $Header: /u/gcmpack/MITgcm/model/src/update_r_star.F,v 1.9 2011/05/20 01:09:37 jmc Exp $
2			C $Name: $
3
4			#include "CPP_OPTIONS.h"
5
6			CBOP
7			C !ROUTINE: RESET_NLFS_VARS
8			C !INTERFACE:
9			SUBROUTINE RESET_NLFS_VARS( myTime, myIter, myThid )
10
11			C !DESCRIPTION: \bv
12			C ==========================================================
13			C \| SUBROUTINE RESET_NLFS_VARS
14			C \| o Re-set some Non-Linear Free-Surface variables
15			C \| in order to facilitate the AD tool task of solving
16			C \| dependency rules.
17			C ==========================================================
18			C \ev
19
20			C !USES:
21			IMPLICIT NONE
22			C == Global variables
23			#include "SIZE.h"
24			#include "EEPARAMS.h"
25			#include "PARAMS.h"
26			#include "GRID.h"
27			#include "SURFACE.h"
28
29			C !INPUT/OUTPUT PARAMETERS:
30			C == Routine arguments ==
31			C myTime :: Current time in simulation
32			C myIter :: Current iteration number in simulation
33			C myThid :: my Thread Id number
34			_RL myTime
35			INTEGER myIter
36			INTEGER myThid
37
38			C !LOCAL VARIABLES:
39			#ifdef NONLIN_FRSURF
40			C Local variables
41			C i,j,k,bi,bj :: loop counter
42			INTEGER i,j,bi,bj
43			CEOP
44
45			C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-\|--+----\|
46
47			DO bj=myByLo(myThid), myByHi(myThid)
48			DO bi=myBxLo(myThid), myBxHi(myThid)
49
50			#ifndef DISABLE_RSTAR_CODE
51			IF ( fluidIsAir ) THEN
52			DO j=1-OLy,sNy+OLy
53			DO i=1-OLx,sNx+OLx
54			pStarFacK(i,j,bi,bj) = rStarFacC(i,j,bi,bj)**atm_kappa
55			ENDDO
56			ENDDO
57			ELSE
58			DO j=1-OLy,sNy+OLy
59			DO i=1-OLx,sNx+OLx
60			pStarFacK(i,j,bi,bj) = 1. _d 0
61			ENDDO
62			ENDDO
63			ENDIF
64			#endif /* DISABLE_RSTAR_CODE */
65
66			C- end bi,bj loop
67			ENDDO
68			ENDDO
69
70			#endif /* NONLIN_FRSURF */
71
72			RETURN
73			END