model/src/reset_nlfs_vars.F

ccbf39d1a4 Jean*0001 #include "CPP_OPTIONS.h"
                0002
                0003 CBOP
                0004 C     !ROUTINE: RESET_NLFS_VARS
                0005 C     !INTERFACE:
                0006       SUBROUTINE RESET_NLFS_VARS( myTime, myIter, myThid )
                0007
                0008 C     !DESCRIPTION: \bv
                0009 C     *==========================================================*
                0010 C     | SUBROUTINE RESET_NLFS_VARS
                0011 C     | o Re-set some Non-Linear Free-Surface variables
                0012 C     |   in order to facilitate the AD tool task of solving
                0013 C     |   dependency rules.
                0014 C     *==========================================================*
                0015 C     \ev
                0016
                0017 C     !USES:
                0018       IMPLICIT NONE
                0019 C     == Global variables
                0020 #include "SIZE.h"
                0021 #include "EEPARAMS.h"
                0022 #include "PARAMS.h"
                0023 #include "GRID.h"
                0024 #include "SURFACE.h"
                0025
                0026 C     !INPUT/OUTPUT PARAMETERS:
                0027 C     myTime    :: Current time in simulation
                0028 C     myIter    :: Current iteration number in simulation
                0029 C     myThid    :: my Thread Id number
                0030       _RL myTime
                0031       INTEGER myIter
                0032       INTEGER myThid
                0033
                0034 C     !LOCAL VARIABLES:
                0035 #ifdef NONLIN_FRSURF
9f5240b52a Jean*0036 #ifndef DISABLE_RSTAR_CODE
                0037 C     i, j,     :: loop counter
                0038 C     bi, bj    :: tile indices
                0039       INTEGER i, j, bi, bj
ccbf39d1a4 Jean*0040 CEOP
                0041
                0042 C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
                0043
                0044       DO bj=myByLo(myThid), myByHi(myThid)
                0045        DO bi=myBxLo(myThid), myBxHi(myThid)
                0046
04f533b838 Jean*0047         IF ( fluidIsAir .AND. select_rStar.GE.1 ) THEN
ccbf39d1a4 Jean*0048          DO j=1-OLy,sNy+OLy
                0049           DO i=1-OLx,sNx+OLx
                0050            pStarFacK(i,j,bi,bj) = rStarFacC(i,j,bi,bj)**atm_kappa
                0051           ENDDO
                0052          ENDDO
                0053         ELSE
                0054          DO j=1-OLy,sNy+OLy
                0055           DO i=1-OLx,sNx+OLx
                0056            pStarFacK(i,j,bi,bj) = 1. _d 0
                0057           ENDDO
                0058          ENDDO
                0059         ENDIF
                0060
                0061 C- end bi,bj loop
                0062        ENDDO
                0063       ENDDO
                0064
9f5240b52a Jean*0065 #endif /* DISABLE_RSTAR_CODE */
ccbf39d1a4 Jean*0066 #endif /* NONLIN_FRSURF */
                0067
                0068       RETURN
                0069       END