C $Header: /u/gcmpack/MITgcm/pkg/exf/exf_getsurfacefluxes.F,v 1.24 2015/05/26 17:08:45 gforget Exp $
C $Name: $
#include "EXF_OPTIONS.h"
#ifdef ALLOW_CTRL
# include "CTRL_OPTIONS.h"
#endif
#ifdef ALLOW_ECCO
# include "ECCO_OPTIONS.h"
#endif
subroutine EXF_GETSURFACEFLUXES(
I mytime,
I myiter,
I mythid
& )
c ==================================================================
c SUBROUTINE exf_GetSurfaceFluxes
c ==================================================================
c
c o Mid-level routine for enabling the use of flux fields as control
c variables.
c
c started: Christian Eckert eckert@mit.edu 30-Jun-1999
c
c changed: Christian Eckert eckert@mit.edu 14-Jan-2000
c - Restructured the code in order to create a package
c for the MITgcmUV.
c
c Christian Eckert eckert@mit.edu 12-Feb-2000
c - Changed Routine names (package prefix: exf_)
c
c mods for pkg/seaice: menemenlis@jpl.nasa.gov 20-Dec-2002
c
c ==================================================================
c SUBROUTINE exf_GetSurfaceFluxes
c ==================================================================
implicit none
c == global variables ==
#include "EEPARAMS.h"
#include "SIZE.h"
#include "PARAMS.h"
#include "GRID.h"
#include "EXF_FIELDS.h"
#include "EXF_PARAM.h"
#ifdef ALLOW_CTRL
# include "CTRL_SIZE.h"
# include "ctrl.h"
# include "ctrl_dummy.h"
# include "CTRL_GENARR.h"
#endif
#if (defined (ALLOW_ECCO) defined (ECCO_CTRL_DEPRECATED))
# include "ecco_cost.h"
#endif
c == routine arguments ==
_RL mytime
integer myiter
integer mythid
c == global variables ==
#ifdef ALLOW_CTRL
#ifdef ALLOW_ROTATE_UV_CONTROLS
_RL tmpUE(1-olx:snx+olx,1-oly:sny+oly,nsx,nsy)
_RL tmpVN(1-olx:snx+olx,1-oly:sny+oly,nsx,nsy)
_RL tmpUX(1-olx:snx+olx,1-oly:sny+oly,nsx,nsy)
_RL tmpVY(1-olx:snx+olx,1-oly:sny+oly,nsx,nsy)
#endif
#ifdef ALLOW_GENTIM2D_CONTROL
integer iarr
#endif
#if (defined (ALLOW_ROTATE_UV_CONTROLS) defined (ALLOW_GENTIM2D_CONTROL))
integer i,j,bi,bj
#endif
c == end of interface ==
#ifdef ALLOW_ROTATE_UV_CONTROLS
if ( useCTRL ) then
do bj = mybylo(mythid),mybyhi(mythid)
do bi = mybxlo(mythid),mybxhi(mythid)
do j = 1-oly,sny+oly
do i = 1-olx,snx+olx
tmpUE(i,j,bi,bj) = 0. _d 0
tmpVN(i,j,bi,bj) = 0. _d 0
tmpUX(i,j,bi,bj) = 0. _d 0
tmpVY(i,j,bi,bj) = 0. _d 0
enddo
enddo
enddo
enddo
endif
#endif
#if (defined (ALLOW_CTRL)
defined (ALLOW_GENTIM2D_CONTROL))
if ( useCTRL.AND.ctrlUseGen) then
DO bj = myByLo(myThid),myByHi(myThid)
DO bi = myBxLo(myThid),mybxhi(myThid)
DO j = 1,sNy
DO i = 1,sNx
do iarr = 1, maxCtrlTim2D
#ifndef ALLOW_ROTATE_UV_CONTROLS
if (xx_gentim2d_file(iarr)(1:7).EQ.'xx_tauu')
& ustress(i,j,bi,bj)=ustress(i,j,bi,bj)+
& xx_gentim2d(i,j,bi,bj,iarr)
if (xx_gentim2d_file(iarr)(1:7).EQ.'xx_tauv')
& vstress(i,j,bi,bj)=vstress(i,j,bi,bj)+
& xx_gentim2d(i,j,bi,bj,iarr)
#else
if (xx_gentim2d_file(iarr)(1:7).EQ.'xx_tauu')
& tmpUE(i,j,bi,bj)=tmpUE(i,j,bi,bj)
& +xx_gentim2d(i,j,bi,bj,iarr)
if (xx_gentim2d_file(iarr)(1:7).EQ.'xx_tauv')
& tmpVN(i,j,bi,bj)=tmpVN(i,j,bi,bj)
& +xx_gentim2d(i,j,bi,bj,iarr)
#endif
enddo
ENDDO
ENDDO
ENDDO
ENDDO
endif !if (ctrlUseGen) then
#endif
#if (!defined (ALLOW_ECCO) || defined (ECCO_CTRL_DEPRECATED))
#ifdef ALLOW_HFLUX_CONTROL
call CTRL_GET_GEN (
& xx_hflux_file, xx_hfluxstartdate, xx_hfluxperiod,
& maskc, hflux, xx_hflux0, xx_hflux1, xx_hflux_dummy,
& xx_hflux_remo_intercept, xx_hflux_remo_slope,
& whflux, mytime, myiter, mythid )
#endif
#ifdef ALLOW_SFLUX_CONTROL
call CTRL_GET_GEN (
& xx_sflux_file, xx_sfluxstartdate, xx_sfluxperiod,
& maskc, sflux, xx_sflux0, xx_sflux1, xx_sflux_dummy,
& xx_sflux_remo_intercept, xx_sflux_remo_slope,
& wsflux, mytime, myiter, mythid )
#endif
IF ( .NOT.useAtmWind ) THEN
IF ( .NOT.ctrlUseGen ) THEN
#ifndef ALLOW_ROTATE_UV_CONTROLS
#ifdef ALLOW_USTRESS_CONTROL
call CTRL_GET_GEN (
& xx_tauu_file, xx_tauustartdate, xx_tauuperiod,
& maskw, ustress, xx_tauu0, xx_tauu1, xx_tauu_dummy,
& xx_tauu_remo_intercept, xx_tauu_remo_slope,
& wtauu, mytime, myiter, mythid )
#endif
#ifdef ALLOW_VSTRESS_CONTROL
call CTRL_GET_GEN (
& xx_tauv_file, xx_tauvstartdate, xx_tauvperiod,
& masks, vstress, xx_tauv0, xx_tauv1, xx_tauv_dummy,
& xx_tauv_remo_intercept, xx_tauv_remo_slope,
& wtauv, mytime, myiter, mythid )
#endif
#else
#if (defined (ALLOW_USTRESS_CONTROL) defined (ALLOW_VSTRESS_CONTROL))
call CTRL_GET_GEN (
& xx_tauu_file, xx_tauustartdate, xx_tauuperiod,
& maskc, tmpUE, xx_tauu0, xx_tauu1, xx_tauu_dummy,
& xx_tauu_remo_intercept, xx_tauu_remo_slope,
& wtauu, mytime, myiter, mythid )
call CTRL_GET_GEN (
& xx_tauv_file, xx_tauvstartdate, xx_tauvperiod,
& maskc, tmpVN, xx_tauv0, xx_tauv1, xx_tauv_dummy,
& xx_tauv_remo_intercept, xx_tauv_remo_slope,
& wtauv, mytime, myiter, mythid )
#endif /* ALLOW_USTRESS_CONTROL and ALLOW_VSTRESS_CONTROL */
#endif /* ALLOW_ROTATE_UV_CONTROLS */
ENDIF ! (.NOT.ctrlUseGen)
#else
IF ( (useCTRL).AND.(.NOT.useAtmWind) ) THEN
#endif /* undef ALLOW_ECCO || def ECCO_CTRL_DEPRECATED */
#ifdef ALLOW_ROTATE_UV_CONTROLS
_EXCH_XY_RL(tmpUE,myThid)
_EXCH_XY_RL(tmpVN,myThid)
call ROTATE_UV2EN_RL(tmpUX,tmpVY,tmpUE,tmpVN,
& .FALSE.,stressIsOnCgrid,.TRUE.,1,mythid)
IF ( stressIsOnCgrid ) THEN
CALL EXCH_UV_XY_RL( tmpUX, tmpVY, .TRUE., myThid )
ELSE
CALL EXCH_UV_AGRID_3D_RL( tmpUX, tmpVY, .TRUE., 1, myThid)
ENDIF
do bj = mybylo(mythid),mybyhi(mythid)
do bi = mybxlo(mythid),mybxhi(mythid)
do j = 1-oly,sny+oly
do i = 1-olx,snx+olx
ustress(i,j,bi,bj)=ustress(i,j,bi,bj)+tmpUX(i,j,bi,bj)
vstress(i,j,bi,bj)=vstress(i,j,bi,bj)+tmpVY(i,j,bi,bj)
enddo
enddo
enddo
enddo
#endif /* ALLOW_ROTATE_UV_CONTROLS */
ENDIF !( .NOT.useAtmWind )
#endif /* ALLOW_CTRL */
end