pkg/ecco/cost_gencost_customize.F

C $Header: /u/gcmpack/MITgcm/pkg/ecco/cost_gencost_customize.F,v 1.2 2010/03/13 23:12:48 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,bj,1) = qnet(i,j,bi,bj)
c #2:
          gencost_modfld(i,j,bi,bj,1) = empmr(i,j,bi,bj)
c #3:
c          gencost_modfld(i,j,bi,bj,3) = fu(i,j,bi,bj)
c #4:
c          gencost_modfld(i,j,bi,bj,4) = fv(i,j,bi,bj)
c
         enddo
        enddo
       enddo
      enddo

#endif /* ALLOW_GENCOST_CONTRIBUTION */

      end
1	gforget	1.3	C $Header: /u/gcmpack/MITgcm/pkg/ecco/cost_gencost_customize.F,v 1.2 2010/03/13 23:12:48 heimbach Exp $
2	heimbach	1.1	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	heimbach	1.2	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	heimbach	1.1	do bj = jtlo,jthi
52			do bi = itlo,ithi
53			do j = jmin,jmax
54			do i = imin,imax
55			c #1:
56	gforget	1.3	c gencost_modfld(i,j,bi,bj,1) = qnet(i,j,bi,bj)
57	heimbach	1.1	c #2:
58	gforget	1.3	gencost_modfld(i,j,bi,bj,1) = empmr(i,j,bi,bj)
59	heimbach	1.1	c #3:
60	gforget	1.3	c gencost_modfld(i,j,bi,bj,3) = fu(i,j,bi,bj)
61	heimbach	1.1	c #4:
62	gforget	1.3	c gencost_modfld(i,j,bi,bj,4) = fv(i,j,bi,bj)
63	heimbach	1.1	c
64			enddo
65			enddo
66			enddo
67			enddo
68
69			#endif /* ALLOW_GENCOST_CONTRIBUTION */
70
71			end