pkg/ecco/cost_gencost_customize.F

C $Header: $
C $Name:  $

#include "COST_CPPOPTIONS.h"

      subroutine cost_gencost_customize( mythid )
      
c     ==================================================================
c     SUBROUTINE cost_gencost_assignperiod
c     ==================================================================

      implicit none

c     == global variables ==

#include "EEPARAMS.h"
#include "SIZE.h"
#include "PARAMS.h"
#include "DYNVARS.h"
#include "FFIELDS.h"

#include "optim.h"
#include "ecco_cost.h"
#include "ctrl_dummy.h"

c     == routine arguments ==

      integer mythid

#ifdef ALLOW_GENCOST_CONTRIBUTION
c     == local variables ==

      integer bi,bj
      integer i,j
      integer itlo,ithi
      integer jtlo,jthi
      integer jmin,jmax
      integer imin,imax

c     == end of interface ==

      do bj = jtlo,jthi
       do bi = itlo,ithi
        do j = jmin,jmax
         do i =  imin,imax
c #1:
          gencost_modfld(i,j,bi,bi,1) = qnet(i,j,bi,bi)
c #2:
          gencost_modfld(i,j,bi,bi,2) = empmr(i,j,bi,bi)
c #3:
          gencost_modfld(i,j,bi,bi,3) = fu(i,j,bi,bi)
c #4:
          gencost_modfld(i,j,bi,bi,4) = fv(i,j,bi,bi)
c
         enddo
        enddo
       enddo
      enddo

#endif /* ALLOW_GENCOST_CONTRIBUTION */

      end
1	heimbach	1.1	C $Header: $
2			C $Name: $
3
4			#include "COST_CPPOPTIONS.h"
5
6			subroutine cost_gencost_customize( mythid )
7
8			c ==================================================================
9			c SUBROUTINE cost_gencost_assignperiod
10			c ==================================================================
11
12			implicit none
13
14			c == global variables ==
15
16			#include "EEPARAMS.h"
17			#include "SIZE.h"
18			#include "PARAMS.h"
19			#include "DYNVARS.h"
20			#include "FFIELDS.h"
21
22			#include "optim.h"
23			#include "ecco_cost.h"
24			#include "ctrl_dummy.h"
25
26			c == routine arguments ==
27
28			integer mythid
29
30			#ifdef ALLOW_GENCOST_CONTRIBUTION
31			c == local variables ==
32
33			integer bi,bj
34			integer i,j
35			integer itlo,ithi
36			integer jtlo,jthi
37			integer jmin,jmax
38			integer imin,imax
39
40			c == end of interface ==
41
42			do bj = jtlo,jthi
43			do bi = itlo,ithi
44			do j = jmin,jmax
45			do i = imin,imax
46			c #1:
47			gencost_modfld(i,j,bi,bi,1) = qnet(i,j,bi,bi)
48			c #2:
49			gencost_modfld(i,j,bi,bi,2) = empmr(i,j,bi,bi)
50			c #3:
51			gencost_modfld(i,j,bi,bi,3) = fu(i,j,bi,bi)
52			c #4:
53			gencost_modfld(i,j,bi,bi,4) = fv(i,j,bi,bi)
54			c
55			enddo
56			enddo
57			enddo
58			enddo
59
60			#endif /* ALLOW_GENCOST_CONTRIBUTION */
61
62			end