rp_frc_adjust_mod Module Reference

Functions/Subroutines
subroutine, public	rp_frc_adjust (ng, tile, linp)
subroutine	rp_frc_adjust_tile (ng, tile, lbi, ubi, lbj, ubj, imins, imaxs, jmins, jmaxs, tl_ustr, tl_vstr, tl_sustr, tl_svstr, tl_tflux, tl_stflx, linp)

Function/Subroutine Documentation

rp_frc_adjust()

subroutine, public rp_frc_adjust_mod::rp_frc_adjust	(	integer, intent(in)	ng,
		integer, intent(in)	tile,
		integer, intent(in)	linp )

Definition at line 33 of file rp_frc_adjust.F.

!***********************************************************************
!
      USE mod_param
      USE mod_forces
!
!  Imported variable declarations.
!
      integer, intent(in) :: ng, tile, Linp
!
!  Local variable declarations.
!
      character (len=*), parameter :: MyFile =                          &
     &  __FILE__
!
# include "tile.h"
!
# ifdef PROFILE
      CALL wclock_on (ng, irpm, 7, __line__, myfile)
# endif
      CALL rp_frc_adjust_tile (ng, tile,                                &
     &                         lbi, ubi, lbj, ubj,                      &
     &                         imins, imaxs, jmins, jmaxs,              &
# ifdef ADJUST_WSTRESS
     &                         forces(ng) % tl_ustr,                    &
     &                         forces(ng) % tl_vstr,                    &
     &                         forces(ng) % tl_sustr,                   &
     &                         forces(ng) % tl_svstr,                   &
# endif
# if defined ADJUST_STFLUX && defined SOLVE3D
     &                         forces(ng) % tl_tflux,                   &
     &                         forces(ng) % tl_stflx,                   &
# endif
     &                         linp)
# ifdef PROFILE
      CALL wclock_off (ng, irpm, 7, __line__, myfile)
# endif
!
      RETURN

References mod_forces::forces, mod_param::irpm, rp_frc_adjust_tile(), wclock_off(), and wclock_on().

Referenced by rp_main3d().

Here is the call graph for this function:

Here is the caller graph for this function:

rp_frc_adjust_tile()

subroutine rp_frc_adjust_mod::rp_frc_adjust_tile ( integer, intent(in) ng, integer, intent(in) tile, integer, intent(in) lbi, integer, intent(in) ubi, integer, intent(in) lbj, integer, intent(in) ubj, integer, intent(in) imins, integer, intent(in) imaxs, integer, intent(in) jmins, integer, intent(in) jmaxs, real(r8), dimension(lbi:,lbj:,:,:), intent(in) tl_ustr, real(r8), dimension(lbi:,lbj:,:,:), intent(in) tl_vstr, real(r8), dimension(lbi:,lbj:), intent(inout) tl_sustr, real(r8), dimension(lbi:,lbj:), intent(inout) tl_svstr, real(r8), dimension(lbi:,lbj:,:,:,:), intent(inout) tl_tflux, real(r8), dimension(lbi:,lbj:,:), intent(inout) tl_stflx, integer, intent(in) linp )

private

Definition at line 75 of file rp_frc_adjust.F.

!***********************************************************************
!
      USE mod_param
      USE mod_scalars
!
!  Imported variable declarations.
!
      integer, intent(in) :: ng, tile
      integer, intent(in) :: LBi, UBi, LBj, UBj
      integer, intent(in) :: IminS, ImaxS, JminS, JmaxS
      integer, intent(in) :: Linp
!
# ifdef ASSUMED_SHAPE
#  ifdef ADJUST_WSTRESS
      real(r8), intent(in) :: tl_ustr(LBi:,LBj:,:,:)
      real(r8), intent(in) :: tl_vstr(LBi:,LBj:,:,:)
#  endif
#  if defined ADJUST_STFLUX && defined SOLVE3D
      real(r8), intent(inout) :: tl_tflux(LBi:,LBj:,:,:,:)
#  endif
#  ifdef ADJUST_WSTRESS
      real(r8), intent(inout) :: tl_sustr(LBi:,LBj:)
      real(r8), intent(inout) :: tl_svstr(LBi:,LBj:)
#  endif
#  if defined ADJUST_STFLUX && defined SOLVE3D
      real(r8), intent(inout) :: tl_stflx(LBi:,LBj:,:)
#  endif
# else
#  ifdef ADJUST_WSTRESS
      real(r8), intent(in) :: tl_ustr(LBi:UBi,LBj:UBj,Nfrec(ng),2)
      real(r8), intent(in) :: tl_vstr(LBi:UBi,LBj:UBj,Nfrec(ng),2)
#  endif
#  if defined ADJUST_STFLUX && defined SOLVE3D
      real(r8), intent(in) :: tl_tflux(LBi:UBi,LBj:UBj,                 &
     &                                 Nfrec(ng),2,NT(ng))
#  endif
#  ifdef ADJUST_WSTRESS
      real(r8), intent(inout) :: tl_sustr(LBi:UBi,LBj:UBj)
      real(r8), intent(inout) :: tl_svstr(LBi:UBi,LBj:UBj)
#  endif
#  if defined ADJUST_STFLUX && defined SOLVE3D
      real(r8), intent(inout) :: tl_stflx(LBi:UBi,LBj:UBj,NT(ng))
#  endif
# endif
!
!  Local variable declarations.
!
      integer :: i, it1, it2, j
# ifdef SOLVE3D
      integer :: itrc
# endif
      real(r8) :: fac, fac1, fac2
 
# include "set_bounds.h"
!
!-----------------------------------------------------------------------
!  Adjust representer model surface forcing fields with 4DVar
!  increments.
!-----------------------------------------------------------------------
!
!  Set time records and interpolation factor, if any.
!
      IF (nfrec(ng).eq.1) THEN
        it1=1
        it2=1
        fac1=1.0_r8
        fac2=0.0_r8
      ELSE
# ifdef GENERIC_DSTART
        it1=max(0,(iic(ng)-ntstart(ng))/nsff(ng))+1
# else
        it1=max(0,(iic(ng)-1)/nsff(ng))+1
# endif
        it2=min(it1+1,nfrec(ng))
        fac1=sf_time(it2,ng)-(time(ng)+dt(ng))
        fac2=(time(ng)+dt(ng))-sf_time(it1,ng)
        fac=1.0_r8/(fac1+fac2)
        fac1=fac*fac1
        fac2=fac*fac2
      END IF
 
# ifdef ADJUST_WSTRESS
!
!  Adjust surface wind stress. Interpolate between surface forcing
!  increments, if appropriate.
!
      DO j=jstrr,jendr
        DO i=istr,iendr
           tl_sustr(i,j)=tl_sustr(i,j)+                                 &
     &                   fac1*tl_ustr(i,j,it1,linp)+                    &
     &                   fac2*tl_ustr(i,j,it2,linp)
        END DO
      END DO
      DO j=jstr,jendr
        DO i=istrr,iendr
           tl_svstr(i,j)=tl_svstr(i,j)+                                 &
     &                   fac1*tl_vstr(i,j,it1,linp)+                    &
     &                   fac2*tl_vstr(i,j,it2,linp)
        END DO
      END DO
# endif
# if defined ADJUST_STFLUX && defined SOLVE3D
!
!  Adjust surface tracer fluxes. Interpolate between surface forcing
!  increments, if appropriate.
!
      DO itrc=1,nt(ng)
        IF (lstflux(itrc,ng)) THEN
          DO j=jstrr,jendr
            DO i=istrr,iendr
               tl_stflx(i,j,itrc)=tl_stflx(i,j,itrc)+                   &
     &                            fac1*tl_tflux(i,j,it1,linp,itrc)+     &
     &                            fac2*tl_tflux(i,j,it2,linp,itrc)
            END DO
          END DO
        END IF
      END DO
# endif
!
      RETURN

References mod_scalars::dt, mod_scalars::iic, mod_scalars::lstflux, mod_scalars::nsff, mod_scalars::ntstart, mod_scalars::sf_time, and mod_scalars::time.

Referenced by rp_frc_adjust().

Here is the caller graph for this function: