/[MITgcm]/MITgcm/pkg/ecco/cost_gencost_customize.F
ViewVC logotype

Contents of /MITgcm/pkg/ecco/cost_gencost_customize.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph


Revision 1.1 - (show annotations) (download)
Sat Feb 6 02:43:03 2010 UTC (15 years, 5 months ago) by heimbach
Branch: MAIN
CVS Tags: checkpoint62c
Preparing usage of generic cost function terms.
Enable with CPP option
#ifdef ALLOW_GENCOST_CONTRIBUTION
First usage is adding air-sea flux constraints when using bulk controls.
---> NOT YET READY FOR PRIME TIME <---

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

  ViewVC Help
Powered by ViewVC 1.1.22