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

Diff of /MITgcm/pkg/ecco/cost_hyd.F

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

revision 1.3 by heimbach, Fri Apr 29 10:31:56 2005 UTC revision 1.4 by heimbach, Wed Aug 31 19:00:52 2005 UTC
# Line 3  C $Header$ Line 3  C $Header$
3  #include "COST_CPPOPTIONS.h"  #include "COST_CPPOPTIONS.h"
4    
5    
6        subroutine cost_Hyd(        subroutine cost_hyd( myiter, mytime, mythid )
      I                     myiter,  
      I                     mytime,  
      I                     mythid  
      &                   )  
7    
8  c     ==================================================================  c     ==================================================================
9  c     SUBROUTINE cost_Hyd  c     SUBROUTINE cost_hyd
10  c     ==================================================================  c     ==================================================================
11  c  c
12  c     o Evaluate cost function contributions of temperature, salt, and  c     o Evaluate cost function contributions of temperature, salt, and
# Line 24  c              - Restructured the code i Line 20  c              - Restructured the code i
20  c                for the MITgcmUV.  c                for the MITgcmUV.
21  c  c
22  c     ==================================================================  c     ==================================================================
23  c     SUBROUTINE cost_Hyd  c     SUBROUTINE cost_hyd
24  c     ==================================================================  c     ==================================================================
25    
26        implicit none        implicit none
# Line 33  c     == global variables == Line 29  c     == global variables ==
29  #include "EEPARAMS.h"  #include "EEPARAMS.h"
30  #include "SIZE.h"  #include "SIZE.h"
31  #include "PARAMS.h"  #include "PARAMS.h"
32    #include "GRID.h"
33    
34    #include "cal.h"
35    #include "ctrl.h"
36    #include "ctrl_dummy.h"
37    #include "ecco_cost.h"
38    
39  c     == routine arguments ==  c     == routine arguments ==
40    
# Line 43  c     == routine arguments == Line 45  c     == routine arguments ==
45  c     == local variables ==  c     == local variables ==
46        character*(max_len_mbuf) msgbuf        character*(max_len_mbuf) msgbuf
47    
48          integer nnzbar
49          integer nnzobs
50          _RL spminloc
51          _RL spmaxloc
52          _RL spzeroloc
53          _RL localperiod
54    
55  c     == end of interface ==  c     == end of interface ==
56    
57    
# Line 50  c     == end of interface == Line 59  c     == end of interface ==
59        write(msgbuf,'(a)') 'ph-cost call cost_sst'        write(msgbuf,'(a)') 'ph-cost call cost_sst'
60        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
61       &     SQUEEZE_RIGHT , mythid)       &     SQUEEZE_RIGHT , mythid)
62        call cost_sst( myiter, mytime, mythid )        nnzbar = nr
63          nnzobs = 1
64          spminloc = -1.8
65          spmaxloc = 40.
66          spzeroloc = 0.
67          localperiod = 0.
68          call cost_generic(
69         &     nnzbar, tbarfile, tbar, xx_tbar_mean_dummy,
70         &     nnzobs, sstdatfile, sstdat,
71         &     nmonsrec, modelstartdate, localperiod,
72         &     maskC, wsst,
73         &     spminloc, spmaxloc, spzeroloc,
74         &     objf_sst, num_sst,
75         &     myiter, mytime, mythid )
76    cph      call cost_sst   ( myiter, mytime, mythid )
77  #endif  #endif
78    
79  #ifdef ALLOW_ARGO_THETA_COST_CONTRIBUTION  #ifdef ALLOW_ARGO_THETA_COST_CONTRIBUTION
# Line 78  c     == end of interface == Line 101  c     == end of interface ==
101        write(msgbuf,'(a)') 'ph-cost call cost_tmi'        write(msgbuf,'(a)') 'ph-cost call cost_tmi'
102        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
103       &     SQUEEZE_RIGHT , mythid)       &     SQUEEZE_RIGHT , mythid)
104        call cost_tmi( myiter, mytime, mythid )        nnzbar = nr
105          nnzobs = 1
106          spminloc = -1.8
107          spmaxloc = 40.
108          spzeroloc = 0.
109          localperiod = 0.
110          call cost_generic(
111         &     nnzbar, tbarfile, tbar, xx_tbar_mean_dummy,
112         &     nnzobs, sstdatfile, sstdat,
113         &     nmonsrec, modelstartdate, localperiod,
114         &     maskC, wsst,
115         &     spminloc, spmaxloc, spzeroloc,
116         &     objf_tmi, num_tmi,
117         &     myiter, mytime, mythid )
118    cph      call cost_tmi( myiter, mytime, mythid )
119  #endif  #endif
120    
121  #ifdef ALLOW_SSS_COST_CONTRIBUTION  #ifdef ALLOW_SSS_COST_CONTRIBUTION
122        write(msgbuf,'(a)') 'ph-cost call cost_sss'        write(msgbuf,'(a)') 'ph-cost call cost_sss'
123        call print_message( msgbuf, standardmessageunit,        call print_message( msgbuf, standardmessageunit,
124       &     SQUEEZE_RIGHT , mythid)       &     SQUEEZE_RIGHT , mythid)
125        call cost_sss   ( myiter, mytime, mythid )        nnzbar = nr
126          nnzobs = 1
127          spminloc = 20.
128          spmaxloc = 40.
129          spzeroloc = 0.
130          localperiod = 0.
131          call cost_generic(
132         &     nnzbar, sbarfile, sbar, xx_sbar_mean_dummy,
133         &     nnzobs, sssdatfile, sssdat,
134         &     nmonsrec, modelstartdate, localperiod,
135         &     maskC, wsss,
136         &     spminloc, spmaxloc, spzeroloc,
137         &     objf_sss, num_sss,
138         &     myiter, mytime, mythid )
139    cph      call cost_sss   ( myiter, mytime, mythid )
140  #endif  #endif
141    
142  #ifdef ALLOW_CTDS_COST_CONTRIBUTION  #ifdef ALLOW_CTDS_COST_CONTRIBUTION
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.4
  ViewVC Help
Powered by ViewVC 1.1.22