C $Header: /u/gcmpack/MITgcm/pkg/ctrl/ctrl_cost_driver.F,v 1.9 2015/10/29 03:43:59 gforget Exp $ C $Name: $ #include "CTRL_OPTIONS.h" subroutine CTRL_COST_DRIVER( myThid ) c ================================================================== c SUBROUTINE ctrl_cost_driver c ================================================================== c c ================================================================== c SUBROUTINE ctrl_cost_driver c ================================================================== implicit none c == global variables == #include "EEPARAMS.h" #include "SIZE.h" #include "PARAMS.h" #include "GRID.h" #ifdef ALLOW_CTRL # include "CTRL_SIZE.h" # include "ctrl.h" # include "ctrl_dummy.h" # include "CTRL_GENARR.h" #endif c == routine arguments == integer myThid c == local variables == #ifdef ALLOW_CTRL #ifdef ALLOW_GENTIM2D_CONTROL integer startrec integer endrec #endif #if (defined (ALLOW_GENARR2D_CONTROL) defined (ALLOW_GENARR3D_CONTROL) defined (ALLOW_GENTIM2D_CONTROL)) integer iarr logical dodimensionalcost integer k2 #endif c == end of interface == c-- Evaluate the individual cost function contributions. #ifdef ALLOW_GENTIM2D_CONTROL do iarr = 1, maxCtrlTim2D dodimensionalcost=.FALSE. do k2 = 1, maxCtrlProc if (xx_gentim2d_preproc(k2,iarr).EQ.'noscaling') then dodimensionalcost=.TRUE. endif enddo if (xx_gentim2d_weight(iarr).NE.' ') then startrec = ncvarrecstart(300+iarr) endrec = ncvarrecsend(300+iarr) do k2 = 1, maxCtrlProc if (xx_gentim2d_preproc(k2,iarr).EQ.'replicate') then if (xx_gentim2d_preproc_i(k2,iarr).NE.0) then endrec=min(endrec,xx_gentim2d_preproc_i(k2,iarr)) endif endif enddo call CTRL_COST_GEN2D ( & startrec, endrec, & xx_gentim2d_file(iarr),xx_gentim2d_dummy(iarr), & xx_gentim2d_period(iarr), & wgentim2d(1-Olx,1-Oly,1,1,iarr), & dodimensionalcost, & num_gentim2d(1,1,iarr), & objf_gentim2d(1,1,iarr), #ifdef ECCO_CTRL_DEPRECATED & zeroRL, num_gentim2dm(1,1,iarr), objf_gentim2dm(1,1,iarr), & objf_gentim2dsmoo(1,1,iarr), zeroRL, zeroRL, #endif /* ECCO_CTRL_DEPRECATED */ & maskC, myThid ) endif enddo #endif #ifdef ALLOW_GENARR2D_CONTROL do iarr = 1, maxCtrlArr2D dodimensionalcost=.FALSE. do k2 = 1, maxCtrlProc if (xx_genarr2d_preproc(k2,iarr).EQ.'noscaling') then dodimensionalcost=.TRUE. endif enddo if (xx_genarr2d_weight(iarr).NE.' ') then call CTRL_COST_GEN2D ( & 1,1, & xx_genarr2d_file(iarr),xx_genarr2d_dummy(iarr), & zeroRL, wgenarr2d(1-Olx,1-Oly,1,1,iarr), & dodimensionalcost, & num_genarr2d(1,1,iarr), objf_genarr2d(1,1,iarr), #ifdef ECCO_CTRL_DEPRECATED & zeroRL, num_zero_mean, objf_zero_mean, & objf_zero_smoo, zeroRL, zeroRL, #endif /* ECCO_CTRL_DEPRECATED */ & maskC, myThid ) endif enddo #endif #ifdef ALLOW_GENARR3D_CONTROL do iarr = 1, maxCtrlArr3D dodimensionalcost=.FALSE. do k2 = 1, maxCtrlProc if (xx_genarr3d_preproc(k2,iarr).EQ.'noscaling') then dodimensionalcost=.TRUE. endif enddo if (xx_genarr3d_weight(iarr).NE.' ') then call CTRL_COST_GEN3D ( & xx_genarr3d_file(iarr),xx_genarr3d_dummy(iarr), & wgenarr3d(1-Olx,1-Oly,1,1,1,iarr), & dodimensionalcost, & num_genarr3d(1,1,iarr), objf_genarr3d(1,1,iarr), & maskC, myThid ) endif enddo #endif #endif return end