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