model/src/reset_nlfs_vars.F

C $Header: /u/gcmpack/MITgcm/model/src/reset_nlfs_vars.F,v 1.1 2014/04/29 21:07:39 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 .AND. select_rStar.GE.1 ) 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	C $Header: /u/gcmpack/MITgcm/model/src/reset_nlfs_vars.F,v 1.1 2014/04/29 21:07:39 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 .AND. select_rStar.GE.1 ) 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