/[MITgcm]/MITgcm/pkg/cost/cost_initvaria.F
ViewVC logotype

Contents of /MITgcm/pkg/cost/cost_initvaria.F

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


Revision 1.2 - (show annotations) (download)
Mon Jun 23 22:29:05 2003 UTC (21 years ago) by heimbach
Branch: MAIN
CVS Tags: checkpoint51f_pre, branchpoint-genmake2, checkpoint51e_post, checkpoint51b_post, checkpoint51c_post, checkpoint51, checkpoint51b_pre, checkpoint51f_post, checkpoint51d_post, checkpoint51a_post, checkpoint50i_post
Branch point for: branch-genmake2
Changes since 1.1: +87 -0 lines
Preparing next differentiable checkpoint and sync
of MAIN vs. ecco-branch

1 C $Header: /u/gcmpack/MITgcm/pkg/cost/cost_init.F,v 1.4 2002/07/13 03:01:39 heimbach Exp $
2
3 #include "COST_CPPOPTIONS.h"
4
5
6 subroutine cost_initvaria( mythid )
7
8 c ==================================================================
9 c SUBROUTINE cost_initvaria
10 c ==================================================================
11 c
12 c o Initialise the variable cost function part.
13 c
14 c started: Christian Eckert eckert@mit.edu 30-Jun-1999
15 c
16 c changed: Christian Eckert eckert@mit.edu 18-Apr-2000
17 c
18 c - Restructured the code in order to create a package
19 c for the MITgcmUV.
20 c
21 c ==================================================================
22 c SUBROUTINE cost_initvaria
23 c ==================================================================
24
25 implicit none
26
27 c == global variables ==
28
29 #include "EEPARAMS.h"
30 #include "SIZE.h"
31
32 #include "cost.h"
33
34 c == routine arguments ==
35
36 integer mythid
37
38 c == local variables ==
39
40 integer i,j,k
41 integer bi,bj
42 integer itlo,ithi
43 integer jtlo,jthi
44
45 c == external functions ==
46
47 c == end of interface ==
48 jtlo = mybylo(mythid)
49 jthi = mybyhi(mythid)
50 itlo = mybxlo(mythid)
51 ithi = mybxhi(mythid)
52
53 c-- Initialize the tiled cost function contributions.
54 do bj = jtlo,jthi
55 do bi = itlo,ithi
56 objf_hflux(bi,bj) = 0. _d 0
57 objf_sflux(bi,bj) = 0. _d 0
58 objf_tauu(bi,bj) = 0. _d 0
59 objf_tauv(bi,bj) = 0. _d 0
60 objf_temp(bi,bj) = 0. _d 0
61 objf_salt(bi,bj) = 0. _d 0
62 objf_sst(bi,bj) = 0. _d 0
63 objf_h(bi,bj) = 0. _d 0
64 objf_atl(bi,bj) = 0. _d 0
65 objf_ctdt(bi,bj) = 0. _d 0
66 objf_ctds(bi,bj) = 0. _d 0
67 objf_test(bi,bj) = 0. _d 0
68 objf_tracer(bi,bj) = 0. _d 0
69 #ifdef ALLOW_COST_VECTOR
70 do i=1,sNx
71 objf_vector(i,bi,bj) = 0. _d 0
72 end do
73 #endif
74 enddo
75 enddo
76
77 c-- Initialise the "global" parts of the cost function.
78 _BEGIN_MASTER( mythid )
79 objf_hmean = 0. _d 0
80 fc = 0. _d 0
81 _END_MASTER( mythid )
82
83 _BARRIER
84
85 return
86 end
87

  ViewVC Help
Powered by ViewVC 1.1.22