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

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

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


Revision 1.1 - (show annotations) (download)
Thu Nov 6 22:10:07 2003 UTC (20 years, 7 months ago) by heimbach
Branch: MAIN
CVS Tags: checkpoint53f_post, checkpoint54a_pre, checkpoint55c_post, checkpoint53b_pre, checkpoint52l_pre, checkpoint52e_pre, hrcube4, hrcube5, checkpoint52j_post, checkpoint52e_post, checkpoint52d_pre, checkpoint53c_post, checkpoint53d_post, checkpoint55d_pre, checkpoint52j_pre, checkpoint54a_post, branch-netcdf, checkpoint52b_pre, checkpoint52n_post, checkpoint54b_post, checkpoint54d_post, checkpoint54e_post, checkpoint55b_post, checkpoint52m_post, checkpoint55, checkpoint53a_post, checkpoint55a_post, checkpoint52l_post, checkpoint52k_post, checkpoint54, checkpoint53b_post, checkpoint53, checkpoint52, checkpoint52d_post, checkpoint52a_post, checkpoint52b_post, checkpoint53g_post, checkpoint52f_post, checkpoint52c_post, ecco_c52_e35, checkpoint54f_post, checkpoint52a_pre, checkpoint53d_pre, checkpoint54c_post, checkpoint52i_post, checkpoint52i_pre, checkpoint51u_post, checkpoint52h_pre, checkpoint52f_pre, hrcube_1, hrcube_2, hrcube_3
Branch point for: netcdf-sm0
o merging from ecco-branch
o pkg/ecco now containes ecco-specific part of cost function
o top level routines the_main_loop, forward_step
  supersede those in model/src/
  previous input data.cost now in data.ecco
  (new namelist ecco_cost_nml)

1 C $Header: /u/gcmpack/MITgcm/pkg/cost/Attic/cost_readtmifields.F,v 1.1.2.1 2003/06/19 15:21:16 heimbach Exp $
2
3 #include "COST_CPPOPTIONS.h"
4
5
6 subroutine cost_ReadTMIFields(
7 I irec,
8 I mythid
9 & )
10
11 c ==================================================================
12 c SUBROUTINE cost_ReadTMIFields
13 c ==================================================================
14 c
15 c o Read a given record of the TMI SST data.
16 c
17 c started: ARmin Koehl akoehl@ucsd.edu
18 c - Restructured the code in order to create a package
19 c for the MITgcmUV.
20 c
21 c ==================================================================
22 c SUBROUTINE cost_ReadTMIFields
23 c ==================================================================
24
25 implicit none
26
27 c == global variables ==
28
29 #include "EEPARAMS.h"
30 #include "SIZE.h"
31 #include "PARAMS.h"
32 #include "GRID.h"
33
34 #include "cal.h"
35 #include "ecco_cost.h"
36
37 c == routine arguments ==
38
39 integer irec
40 integer mythid
41
42 c == local variables ==
43
44 integer bi,bj
45 integer i,j,k
46 integer itlo,ithi
47 integer jtlo,jthi
48 integer jmin,jmax
49 integer imin,imax
50 integer nobs
51 integer tmirec
52 integer begintmi
53 integer beginrun
54
55 _RL spval
56 _RL vartile
57
58 c == end of interface ==
59
60 parameter (spval = -9999. )
61 ce --> there is certainly a better place for this.
62
63 jtlo = mybylo(mythid)
64 jthi = mybyhi(mythid)
65 itlo = mybxlo(mythid)
66 ithi = mybxhi(mythid)
67 jmin = 1
68 jmax = sny
69 imin = 1
70 imax = snx
71
72 begintmi = tmistartdate(1)/10000
73 beginrun = modelstartdate(1)/10000
74 if ( begintmi .eq. beginrun ) then
75 tmirec = mod(modelstartdate(1)/100,100) -
76 & mod(tmistartdate(1)/100,100) + irec
77 else
78 tmirec = ( beginrun - begintmi - 1)*nmonthyear +
79 & (nmonthyear - mod(tmistartdate(1)/100,100) +
80 & 1) + mod(modelstartdate(1)/100,100) - 1 + irec
81 endif
82 if (tmirec.gt.0) then
83 call mdsreadfield( tmidatfile, cost_iprec, cost_yftype, 1, tmidat,
84 & tmirec, mythid )
85 else
86 do bj = jtlo,jthi
87 do bi = itlo,ithi
88 do j = jmin,jmax
89 do i = imin,imax
90 tmidat(i,j,bi,bj)=spval
91 enddo
92 enddo
93 enddo
94 enddo
95 endif
96 return
97 end
98

  ViewVC Help
Powered by ViewVC 1.1.22