pkg/ecco/cost_gencost_customize.F

C $Header: /u/gcmpack/MITgcm/pkg/ecco/cost_gencost_customize.F,v 1.1 2010/02/06 02:43:03 heimbach Exp $
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 ==

      jtlo = mybylo(mythid)
      jthi = mybyhi(mythid)
      itlo = mybxlo(mythid)
      ithi = mybxhi(mythid)
      jmin = 1
      jmax = sny
      imin = 1
      imax = snx

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

#endif /* ALLOW_GENCOST_CONTRIBUTION */

      end
1	C $Header: /u/gcmpack/MITgcm/pkg/ecco/cost_gencost_customize.F,v 1.1 2010/02/06 02:43:03 heimbach Exp $
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	jtlo = mybylo(mythid)
43	jthi = mybyhi(mythid)
44	itlo = mybxlo(mythid)
45	ithi = mybxhi(mythid)
46	jmin = 1
47	jmax = sny
48	imin = 1
49	imax = snx
50
51	do bj = jtlo,jthi
52	do bi = itlo,ithi
53	do j = jmin,jmax
54	do i = imin,imax
55	c #1:
56	c gencost_modfld(i,j,bi,bi,1) = qnet(i,j,bi,bi)
57	c #2:
58	gencost_modfld(i,j,bi,bi,1) = empmr(i,j,bi,bi)
59	c #3:
60	c gencost_modfld(i,j,bi,bi,3) = fu(i,j,bi,bi)
61	c #4:
62	c gencost_modfld(i,j,bi,bi,4) = fv(i,j,bi,bi)
63	c
64	enddo
65	enddo
66	enddo
67	enddo
68
69	#endif /* ALLOW_GENCOST_CONTRIBUTION */
70
71	end