C $Header: /u/gcmpack/MITgcm/pkg/ctrl/ctrl_map_ini.F,v 1.45 2015/02/18 12:31:10 heimbach Exp $
C $Name: $
#include "CTRL_OPTIONS.h"
#ifdef ALLOW_ECCO
# include "ECCO_OPTIONS.h"
#endif
CBOP
C !ROUTINE: ctrl_map_ini
C !INTERFACE:
subroutine CTRL_MAP_INI( mythid )
C !DESCRIPTION: \bv
c *=================================================================
c | SUBROUTINE ctrl_map_ini
c | Add the temperature, salinity, and diffusivity parts of the
c | control vector to the model state and update the tile halos.
c | The control vector is defined in the header file "ctrl.h".
c *=================================================================
C \ev
C !USES:
implicit none
#ifdef ECCO_CTRL_DEPRECATED
c == global variables ==
#include "SIZE.h"
#include "EEPARAMS.h"
#include "PARAMS.h"
#include "GRID.h"
#include "DYNVARS.h"
#include "FFIELDS.h"
#ifdef ALLOW_CTRL
# include "CTRL_SIZE.h"
# include "ctrl.h"
# include "CTRL_FIELDS.h"
# include "CTRL_GENARR.h"
# include "ctrl_dummy.h"
# include "optim.h"
# ifdef ALLOW_ECCO
# include "ecco_cost.h"
else
# include "ctrl_weights.h"
# endif
#endif
#ifdef ALLOW_PTRACERS
# include "PTRACERS_SIZE.h"
c#include "PTRACERS_PARAMS.h"
# include "PTRACERS_FIELDS.h"
#endif
#endif /* ECCO_CTRL_DEPRECATED */
C !INPUT/OUTPUT PARAMETERS:
c == routine arguments ==
integer mythid
#ifdef ECCO_CTRL_DEPRECATED
C !LOCAL VARIABLES:
c == local variables ==
integer bi,bj
integer i,j,k
integer itlo,ithi
integer jtlo,jthi
integer jmin,jmax
integer imin,imax
integer il
logical equal
logical doglobalread
logical ladinit
character*( 80) fnamegeneric
_RL fac
_RL tmptest
c == external ==
integer ilnblnk
external
c == end of interface ==
CEOP
jtlo = mybylo(mythid)
jthi = mybyhi(mythid)
itlo = mybxlo(mythid)
ithi = mybxhi(mythid)
jmin = 1
jmax = sny
imin = 1
imax = snx
doglobalread = .false.
ladinit = .false.
equal = .true.
if ( equal ) then
fac = 1. _d 0
else
fac = 0. _d 0
endif
#ifdef ALLOW_THETA0_CONTROL
c-- Temperature field.
il=ilnblnk( xx_theta_file )
write(fnamegeneric(1:80),'(2a,i10.10)')
& xx_theta_file(1:il),'.',optimcycle
call ACTIVE_READ_XYZ( fnamegeneric, tmpfld3d, 1,
& doglobalread, ladinit, optimcycle,
& mythid, xx_theta_dummy )
do bj = jtlo,jthi
do bi = itlo,ithi
do k = 1,nr
do j = jmin,jmax
do i = imin,imax
#if (defined (ALLOW_ECCO) defined (ECCO_CTRL_DEPRECATED))
IF (abs(tmpfld3d(i,j,k,bi,bj)).gt.
$ 2.0/sqrt(wtheta(k,bi,bj)))
$ tmpfld3d(i,j,k,bi,bj)=
$ sign(2.0/sqrt(wtheta(k,bi,bj)),tmpfld3d(i,j,k,bi,bj))
#endif
theta(i,j,k,bi,bj) = theta(i,j,k,bi,bj) +
& fac*tmpfld3d(i,j,k,bi,bj)
#ifndef DISABLE_CTRL_THETA_LIMIT
if(theta(i,j,k,bi,bj).lt.-2.0)
& theta(i,j,k,bi,bj)= -2.0
#endif
enddo
enddo
enddo
enddo
enddo
#endif
#ifdef ALLOW_SALT0_CONTROL
c-- Temperature field.
il=ilnblnk( xx_salt_file )
write(fnamegeneric(1:80),'(2a,i10.10)')
& xx_salt_file(1:il),'.',optimcycle
call ACTIVE_READ_XYZ( fnamegeneric, tmpfld3d, 1,
& doglobalread, ladinit, optimcycle,
& mythid, xx_salt_dummy )
do bj = jtlo,jthi
do bi = itlo,ithi
do k = 1,nr
do j = jmin,jmax
do i = imin,imax
#if (defined (ALLOW_ECCO) defined (ECCO_CTRL_DEPRECATED))
IF (abs(tmpfld3d(i,j,k,bi,bj)).gt.
$ 2.0/sqrt(wsalt(k,bi,bj)))
$ tmpfld3d(i,j,k,bi,bj)=
$ sign(2.0/sqrt(wsalt(k,bi,bj)),tmpfld3d(i,j,k,bi,bj))
#endif
salt(i,j,k,bi,bj) = salt(i,j,k,bi,bj) +
& fac*tmpfld3d(i,j,k,bi,bj)
enddo
enddo
enddo
enddo
enddO
#endif
#ifdef ALLOW_TR10_CONTROL
#ifdef ALLOW_PTRACERS
c-- Temperature field.
il=ilnblnk( xx_tr1_file )
write(fnamegeneric(1:80),'(2a,i10.10)')
& xx_tr1_file(1:il),'.',optimcycle
call ACTIVE_READ_XYZ( fnamegeneric, tmpfld3d, 1,
& doglobalread, ladinit, optimcycle,
& mythid, xx_tr1_dummy )
do bj = jtlo,jthi
do bi = itlo,ithi
do k = 1,nr
do j = jmin,jmax
do i = imin,imax
#ifdef ALLOW_OPENAD
ptracer(i,j,k,bi,bj,1) = ptracer(i,j,k,bi,bj,1) +
& fac*xx_tr1(i,j,k,bi,bj) +
& fac*tmpfld3d(i,j,k,bi,bj)
#else
ptracer(i,j,k,bi,bj,1) = ptracer(i,j,k,bi,bj,1) +
& fac*tmpfld3d(i,j,k,bi,bj)
#endif
enddo
enddo
enddo
enddo
enddo
#endif
#endif
#ifdef ALLOW_SST0_CONTROL
c-- sst0.
il=ilnblnk( xx_sst_file )
write(fnamegeneric(1:80),'(2a,i10.10)')
& xx_sst_file(1:il),'.',optimcycle
call ACTIVE_READ_XY ( fnamegeneric, tmpfld2d, 1,
& doglobalread, ladinit, optimcycle,
& mythid, xx_sst_dummy )
do bj = jtlo,jthi
do bi = itlo,ithi
do j = jmin,jmax
do i = imin,imax
cph sst(i,j,bi,bj) = sst(i,j,bi,bj) + tmpfld2d(i,j,bi,bj)
theta(i,j,1,bi,bj) = theta(i,j,1,bi,bj)
& + tmpfld2d(i,j,bi,bj)
enddo
enddo
enddo
enddo
#endif
#ifdef ALLOW_SSS0_CONTROL
c-- sss0.
il=ilnblnk( xx_sss_file )
write(fnamegeneric(1:80),'(2a,i10.10)')
& xx_sss_file(1:il),'.',optimcycle
call ACTIVE_READ_XY ( fnamegeneric, tmpfld2d, 1,
& doglobalread, ladinit, optimcycle,
& mythid, xx_sss_dummy )
do bj = jtlo,jthi
do bi = itlo,ithi
do j = jmin,jmax
do i = imin,imax
cph sss(i,j,bi,bj) = sss(i,j,bi,bj) + tmpfld2d(i,j,bi,bj)
salt(i,j,1,bi,bj) = salt(i,j,1,bi,bj)
& + tmpfld2d(i,j,bi,bj)
enddo
enddo
enddo
enddo
#endif
#ifdef ALLOW_AUTODIFF
#ifdef ALLOW_DIFFKR_CONTROL
c-- diffkr.
il=ilnblnk( xx_diffkr_file )
write(fnamegeneric(1:80),'(2a,i10.10)')
& xx_diffkr_file(1:il),'.',optimcycle
call ACTIVE_READ_XYZ( fnamegeneric, tmpfld3d, 1,
& doglobalread, ladinit, optimcycle,
& mythid, xx_diffkr_dummy )
do bj = jtlo,jthi
do bi = itlo,ithi
do k = 1,nr
do j = jmin,jmax
do i = imin,imax
diffkr(i,j,k,bi,bj) = diffkr(i,j,k,bi,bj) +
& tmpfld3d(i,j,k,bi,bj)
enddo
enddo
enddo
enddo
enddo
#endif
#endif
#ifdef ALLOW_AUTODIFF
#ifdef ALLOW_KAPGM_CONTROL
c-- kapgm.
il=ilnblnk( xx_kapgm_file )
write(fnamegeneric(1:80),'(2a,i10.10)')
& xx_kapgm_file(1:il),'.',optimcycle
call ACTIVE_READ_XYZ( fnamegeneric, tmpfld3d, 1,
& doglobalread, ladinit, optimcycle,
& mythid, xx_kapgm_dummy )
do bj = jtlo,jthi
do bi = itlo,ithi
do k = 1,nr
do j = jmin,jmax
do i = imin,imax
#ifdef ALLOW_OPENAD
kapGM(i,j,k,bi,bj) = kapGM(i,j,k,bi,bj) +
& xx_kapgm(i,j,k,bi,bj) +
& tmpfld3d(i,j,k,bi,bj)
#else
kapGM(i,j,k,bi,bj) = kapGM(i,j,k,bi,bj) +
& tmpfld3d(i,j,k,bi,bj)
#endif
enddo
enddo
enddo
enddo
enddo
#endif
#endif
#ifdef ALLOW_AUTODIFF
#ifdef ALLOW_KAPREDI_CONTROL
c-- kapredi.
il=ilnblnk( xx_kapredi_file )
write(fnamegeneric(1:80),'(2a,i10.10)')
& xx_kapredi_file(1:il),'.',optimcycle
call ACTIVE_READ_XYZ( fnamegeneric, tmpfld3d, 1,
& doglobalread, ladinit, optimcycle,
& mythid, xx_kapredi_dummy )
do bj = jtlo,jthi
do bi = itlo,ithi
do k = 1,nr
do j = jmin,jmax
do i = imin,imax
kapRedi(i,j,k,bi,bj) = kapRedi(i,j,k,bi,bj) +
& tmpfld3d(i,j,k,bi,bj)
enddo
enddo
enddo
enddo
enddo
#endif
#endif
#ifdef ALLOW_EFLUXY0_CONTROL
c-- y-component EP-flux field.
il=ilnblnk( xx_efluxy_file )
write(fnamegeneric(1:80),'(2a,i10.10)')
& xx_efluxy_file(1:il),'.',optimcycle
call ACTIVE_READ_XYZ( fnamegeneric, tmpfld3d, 1,
& doglobalread, ladinit, optimcycle,
& mythid, xx_efluxy_dummy )
do bj = jtlo,jthi
do bi = itlo,ithi
do k = 1,nr
do j = jmin,jmax
do i = imin,imax
EfluxY(i,j,k,bi,bj) = EfluxY(i,j,k,bi,bj)
& - fac*tmpfld3d(i,j,k,bi,bj)
& *maskS(i,j,k,bi,bj)
cph EfluxY(i,j,k,bi,bj) = EfluxY(i,j,k,bi,bj)
cph & - rSphere*cosFacU(J,bi,bj)
cph & *fac*tmpfld3d(i,j,k,bi,bj)
enddo
enddo
enddo
enddo
enddo
#endif
#ifdef ALLOW_EFLUXP0_CONTROL
c-- p-component EP-flux field.
il=ilnblnk( xx_efluxp_file )
write(fnamegeneric(1:80),'(2a,i10.10)')
& xx_efluxp_file(1:il),'.',optimcycle
call ACTIVE_READ_XYZ( fnamegeneric, tmpfld3d, 1,
& doglobalread, ladinit, optimcycle,
& mythid, xx_efluxp_dummy )
do bj = jtlo,jthi
do bi = itlo,ithi
do k = 1,nr
do j = jmin,jmax
do i = imin,imax
EfluxP(i,j,k,bi,bj) = EfluxP(i,j,k,bi,bj)
& + fCori(i,j,bi,bj)
& *fac*tmpfld3d(i,j,k,bi,bj)
& *hFacV(i,j,k,bi,bj)
cph EfluxP(i,j,k,bi,bj) = EfluxP(i,j,k,bi,bj)
cph & + fCori(i,j,bi,bj)
cph & *rSphere*cosFacU(J,bi,bj)
cph & *fac*tmpfld3d(i,j,k,bi,bj)
enddo
enddo
enddo
enddo
enddo
#endif
#ifdef ALLOW_BOTTOMDRAG_CONTROL_NONGENERIC
c-- bottom drag
il=ilnblnk( xx_bottomdrag_file )
write(fnamegeneric(1:80),'(2a,i10.10)')
& xx_bottomdrag_file(1:il),'.',optimcycle
call ACTIVE_READ_XY ( fnamegeneric, tmpfld2d, 1,
& doglobalread, ladinit, optimcycle,
& mythid, xx_bottomdrag_dummy )
do bj = jtlo,jthi
do bi = itlo,ithi
do j = jmin,jmax
do i = imin,imax
bottomDragFld(i,j,bi,bj) = bottomDragFld(i,j,bi,bj)
& + tmpfld2d(i,j,bi,bj)
enddo
enddo
enddo
enddo
#endif
#ifdef ALLOW_EDDYPSI_CONTROL
c-- zonal eddy streamfunction : eddyPsiX
il=ilnblnk( xx_edtaux_file )
write(fnamegeneric(1:80),'(2a,i10.10)')
& xx_edtaux_file(1:il),'.',optimcycle
call ACTIVE_READ_XYZ( fnamegeneric, tmpfld3d, 1,
& doglobalread, ladinit, optimcycle,
& mythid, xx_edtaux_dummy )
do bj = jtlo,jthi
do bi = itlo,ithi
do k = 1,nr
do j = jmin,jmax
do i = imin,imax
eddyPsiX(i,j,k,bi,bj) = eddyPsiX(i,j,k,bi,bj) +
& tmpfld3d(i,j,k,bi,bj)
enddo
enddo
enddo
enddo
enddo
c-- meridional eddy streamfunction : eddyPsiY
il=ilnblnk( xx_edtauy_file )
write(fnamegeneric(1:80),'(2a,i10.10)')
& xx_edtauy_file(1:il),'.',optimcycle
call ACTIVE_READ_XYZ( fnamegeneric, tmpfld3d, 1,
& doglobalread, ladinit, optimcycle,
& mythid, xx_edtauy_dummy )
do bj = jtlo,jthi
do bi = itlo,ithi
do k = 1,nr
do j = jmin,jmax
do i = imin,imax
eddyPsiY(i,j,k,bi,bj) = eddyPsiY(i,j,k,bi,bj) +
& tmpfld3d(i,j,k,bi,bj)
enddo
enddo
enddo
enddo
enddo
#endif
#ifdef ALLOW_UVEL0_CONTROL
c-- initial zonal velocity
il=ilnblnk( xx_uvel_file )
write(fnamegeneric(1:80),'(2a,i10.10)')
& xx_uvel_file(1:il),'.',optimcycle
call ACTIVE_READ_XYZ( fnamegeneric, tmpfld3d, 1,
& doglobalread, ladinit, optimcycle,
& mythid, xx_uvel_dummy )
do bj = jtlo,jthi
do bi = itlo,ithi
do k = 1,nr
do j = jmin,jmax
do i = imin,imax
#ifdef ALLOW_OPENAD
uVel(i,j,k,bi,bj) = uVel(i,j,k,bi,bj) +
& fac*xx_uvel(i,j,k,bi,bj)
#else
uVel(i,j,k,bi,bj) = uVel(i,j,k,bi,bj) +
& fac*tmpfld3d(i,j,k,bi,bj)
#endif
enddo
enddo
enddo
enddo
enddo
#endif
#ifdef ALLOW_VVEL0_CONTROL
c-- initial merid. velocity
il=ilnblnk( xx_vvel_file )
write(fnamegeneric(1:80),'(2a,i10.10)')
& xx_vvel_file(1:il),'.',optimcycle
call ACTIVE_READ_XYZ( fnamegeneric, tmpfld3d, 1,
& doglobalread, ladinit, optimcycle,
& mythid, xx_vvel_dummy )
do bj = jtlo,jthi
do bi = itlo,ithi
do k = 1,nr
do j = jmin,jmax
do i = imin,imax
#ifdef ALLOW_OPENAD
vVel(i,j,k,bi,bj) = vVel(i,j,k,bi,bj) +
& fac*xx_vvel(i,j,k,bi,bj)
#else
vVel(i,j,k,bi,bj) = vVel(i,j,k,bi,bj) +
& fac*tmpfld3d(i,j,k,bi,bj)
#endif
enddo
enddo
enddo
enddo
enddo
#endif
#ifdef ALLOW_ETAN0_CONTROL
c-- initial Eta.
il=ilnblnk( xx_etan_file )
write(fnamegeneric(1:80),'(2a,i10.10)')
& xx_etan_file(1:il),'.',optimcycle
call ACTIVE_READ_XY ( fnamegeneric, tmpfld2d, 1,
& doglobalread, ladinit, optimcycle,
& mythid, xx_etan_dummy )
do bj = jtlo,jthi
do bi = itlo,ithi
do j = jmin,jmax
do i = imin,imax
#ifdef ALLOW_OPENAD
etaN(i,j,bi,bj) = etaN(i,j,bi,bj) +
& fac*xx_etan(i,j,bi,bj)
#else
etaN(i,j,bi,bj) = etaN(i,j,bi,bj) +
& fac*tmpfld2d(i,j,bi,bj)
etaH(i,j,bi,bj) = etaH(i,j,bi,bj) +
& fac*tmpfld2d(i,j,bi,bj)
#endif
enddo
enddo
enddo
enddo
#endif
#ifdef ALLOW_RELAXSST_CONTROL
c-- SST relaxation coefficient.
il=ilnblnk( xx_relaxsst_file )
write(fnamegeneric(1:80),'(2a,i10.10)')
& xx_relaxsst_file(1:il),'.',optimcycle
call ACTIVE_READ_XY ( fnamegeneric, tmpfld2d, 1,
& doglobalread, ladinit, optimcycle,
& mythid, xx_relaxsst_dummy )
do bj = jtlo,jthi
do bi = itlo,ithi
do j = jmin,jmax
do i = imin,imax
lambdaThetaClimRelax(i,j,bi,bj) =
& lambdaThetaClimRelax(i,j,bi,bj)
& + tmpfld2d(i,j,bi,bj)
enddo
enddo
enddo
enddo
#endif
#ifdef ALLOW_RELAXSSS_CONTROL
c-- SSS relaxation coefficient.
il=ilnblnk( xx_relaxsss_file )
write(fnamegeneric(1:80),'(2a,i10.10)')
& xx_relaxsss_file(1:il),'.',optimcycle
call ACTIVE_READ_XY ( fnamegeneric, tmpfld2d, 1,
& doglobalread, ladinit, optimcycle,
& mythid, xx_relaxsss_dummy )
do bj = jtlo,jthi
do bi = itlo,ithi
do j = jmin,jmax
do i = imin,imax
lambdaSaltClimRelax(i,j,bi,bj) =
& lambdaSaltClimRelax(i,j,bi,bj)
& + tmpfld2d(i,j,bi,bj)
enddo
enddo
enddo
enddo
#endif
c-- Update the tile edges.
#if (defined (ALLOW_THETA0_CONTROL) defined (ALLOW_SST0_CONTROL))
_EXCH_XYZ_RL( theta, mythid )
#endif
#if (defined (ALLOW_SALT0_CONTROL) defined (ALLOW_SSS0_CONTROL))
_EXCH_XYZ_RL( salt, mythid )
#endif
#ifdef ALLOW_TR10_CONTROL
#ifdef ALLOW_PTRACERS
_EXCH_XYZ_RL(pTracer(1-Olx,1-Oly,1,1,1,1),myThid)
#endif
#endif
#ifdef ALLOW_AUTODIFF
# ifdef ALLOW_DIFFKR_CONTROL
_EXCH_XYZ_RL( diffKr, mythid)
# endif
# ifdef ALLOW_KAPGM_CONTROL
_EXCH_XYZ_RL( kapGM, mythid)
# endif
# ifdef ALLOW_KAPREDI_CONTROL
_EXCH_XYZ_RL( kapRedi, mythid)
# endif
#endif
#ifdef ALLOW_EFLUXY0_CONTROL
_EXCH_XYZ_RL( EfluxY, mythid )
#endif
#ifdef ALLOW_EFLUXP0_CONTROL
_EXCH_XYZ_RL( EfluxP, mythid )
#endif
#ifdef ALLOW_BOTTOMDRAG_CONTROL
_EXCH_XY_RL( bottomDragFld, mythid )
#endif
#ifdef ALLOW_EDDYPSI_CONTROL
CALL EXCH_UV_XYZ_RS(eddyPsiX,eddyPsiY,.TRUE.,myThid)
#endif
#ifdef ALLOW_UVEL0_CONTROL
_EXCH_XYZ_RL( uVel, mythid)
#endif
#ifdef ALLOW_VVEL0_CONTROL
_EXCH_XYZ_RL( vVel, mythid)
#endif
#ifdef ALLOW_ETAN0_CONTROL
_EXCH_XY_RL( etaN, mythid )
_EXCH_XY_RL( etaH, mythid )
#endif
#ifdef ALLOW_RELAXSST_CONTROL
_EXCH_XY_RS( lambdaThetaClimRelax, mythid )
#endif
#ifdef ALLOW_RELAXSSS_CONTROL
_EXCH_XY_RS( lambdaThetaClimRelax, mythid )
#endif
#endif /* ECCO_CTRL_DEPRECATED */
RETURN
END