pkg/ecco/cost_gencost_all.F

C $Header: /u/gcmpack/MITgcm/pkg/ecco/cost_gencost_all.F,v 1.2 2007/10/09 00:02:50 jmc Exp $
C $Name:  $

#include "COST_CPPOPTIONS.h"


      subroutine cost_gencost_all( myiter, mytime, mythid )

c     ==================================================================
c     SUBROUTINE cost_gencost_all
c     ==================================================================
c
c     o Evaluate user defined penalty terms
c
c     ==================================================================
c     SUBROUTINE cost_gencost_all
c     ==================================================================

      implicit none

c     == global variables ==

#include "SIZE.h"
#include "EEPARAMS.h"
#include "PARAMS.h"
#include "GRID.h"

#include "cal.h"
#include "ctrl.h"
#include "ctrl_dummy.h"
#include "ecco_cost.h"

c     == routine arguments ==

      integer mythid
      integer myiter
      _RL     mytime

#ifdef ALLOW_GENCOST_CONTRIBUTION
c     == local variables ==

      integer i,j,k
      character*(max_len_mbuf) msgbuf

      integer nnzbar
      integer nnzobs
      integer nrecloc

c     == external functions ==

c     == end of interface ==

      do k = 1, NGENCOST

         write(msgbuf,'(a,x,i2.2)') 'ph-cost gencost #', k
         call print_message( msgbuf, standardmessageunit,
     &        SQUEEZE_RIGHT , mythid)
         nnzbar = 1
         nnzobs = 1
         call cost_generic(
     &        nnzbar, gencost_barfile(k), 
     &        gencost_barfld(1-Olx,1-Oly,1,1,k), 
     &        xx_genbar_dummy(k),
     &        nnzobs, gencost_datafile(k), mult_gencost(k),
     &        gencost_nrec(k), modelstartdate, gencost_period(k),
     &        gencost_mask(k), gencost_weight(1-Olx,1-Oly,1,1,k),
     &        gencost_spmin(k), gencost_spmax(k), 
     &        gencost_spzero(k),
     &        objf_gencost(1,1,k), num_gencost(1,1,k),
     &        myiter, mytime, mythid )

      end do

#endif /* ALLOW_GENCOST_CONTRIBUTION */

      end
1	C $Header: /u/gcmpack/MITgcm/pkg/ecco/cost_gencost_all.F,v 1.2 2007/10/09 00:02:50 jmc Exp $
2	C $Name: $
3
4	#include "COST_CPPOPTIONS.h"
5
6
7	subroutine cost_gencost_all( myiter, mytime, mythid )
8
9	c ==================================================================
10	c SUBROUTINE cost_gencost_all
11	c ==================================================================
12	c
13	c o Evaluate user defined penalty terms
14	c
15	c ==================================================================
16	c SUBROUTINE cost_gencost_all
17	c ==================================================================
18
19	implicit none
20
21	c == global variables ==
22
23	#include "SIZE.h"
24	#include "EEPARAMS.h"
25	#include "PARAMS.h"
26	#include "GRID.h"
27
28	#include "cal.h"
29	#include "ctrl.h"
30	#include "ctrl_dummy.h"
31	#include "ecco_cost.h"
32
33	c == routine arguments ==
34
35	integer mythid
36	integer myiter
37	_RL mytime
38
39	#ifdef ALLOW_GENCOST_CONTRIBUTION
40	c == local variables ==
41
42	integer i,j,k
43	character*(max_len_mbuf) msgbuf
44
45	integer nnzbar
46	integer nnzobs
47	integer nrecloc
48
49	c == external functions ==
50
51	c == end of interface ==
52
53	do k = 1, NGENCOST
54
55	write(msgbuf,'(a,x,i2.2)') 'ph-cost gencost #', k
56	call print_message( msgbuf, standardmessageunit,
57	& SQUEEZE_RIGHT , mythid)
58	nnzbar = 1
59	nnzobs = 1
60	call cost_generic(
61	& nnzbar, gencost_barfile(k),
62	& gencost_barfld(1-Olx,1-Oly,1,1,k),
63	& xx_genbar_dummy(k),
64	& nnzobs, gencost_datafile(k), mult_gencost(k),
65	& gencost_nrec(k), modelstartdate, gencost_period(k),
66	& gencost_mask(k), gencost_weight(1-Olx,1-Oly,1,1,k),
67	& gencost_spmin(k), gencost_spmax(k),
68	& gencost_spzero(k),
69	& objf_gencost(1,1,k), num_gencost(1,1,k),
70	& myiter, mytime, mythid )
71
72	end do
73
74	#endif /* ALLOW_GENCOST_CONTRIBUTION */
75
76	end