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

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

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

Revision 1.2 - (show annotations) (download)
Tue Oct 9 00:02:50 2007 UTC (16 years, 7 months ago) by jmc
Branch: MAIN
CVS Tags: checkpoint63p, checkpoint63q, checkpoint63l, checkpoint63m, checkpoint63n, checkpoint63o, checkpoint63h, checkpoint63i, checkpoint63j, checkpoint63k, checkpoint63d, checkpoint63e, checkpoint63f, checkpoint63g, checkpoint63a, checkpoint63b, checkpoint63c, checkpoint60, checkpoint61, checkpoint62, checkpoint63, checkpoint59q, checkpoint59p, checkpoint59r, checkpoint59m, checkpoint59l, checkpoint59o, checkpoint59n, checkpoint59i, checkpoint59k, checkpoint59j, checkpoint62c, checkpoint62b, checkpoint62a, checkpoint62g, checkpoint62f, checkpoint62e, checkpoint62d, checkpoint62k, checkpoint62j, checkpoint62i, checkpoint62h, checkpoint62o, checkpoint62n, checkpoint62m, checkpoint62l, checkpoint62s, checkpoint62r, checkpoint62q, checkpoint62p, checkpoint62w, checkpoint62v, checkpoint62u, checkpoint62t, checkpoint62z, checkpoint62y, checkpoint62x, checkpoint61f, checkpoint61g, checkpoint61d, checkpoint61e, checkpoint61b, checkpoint61c, checkpoint61a, checkpoint61n, checkpoint61o, checkpoint61l, checkpoint61m, checkpoint61j, checkpoint61k, checkpoint61h, checkpoint61i, checkpoint61v, checkpoint61w, checkpoint61t, checkpoint61u, checkpoint61r, checkpoint61s, checkpoint61p, checkpoint61q, checkpoint61z, checkpoint61x, checkpoint61y
Changes since 1.1: +7 -6 lines 
add missing cvs $Header:$ or $Name:$
1 C $Header: $
2 C $Name: $
3
4 #include "COST_CPPOPTIONS.h"
5
6
7 subroutine cost_Readargot(
8 I irec,
9 I mythid
10 & )
11
12 c ==================================================================
13 c SUBROUTINE cost_Readargot
14 c ==================================================================
15 c
16 c o Read a given record of the
17 c
18 c started:Armin Koehl akoehl@ucsd.edu
19 c - Restructured the code in order to create a package
20 c for the MITgcmUV.
21 c
22 c ==================================================================
23 c SUBROUTINE cost_Readargot
24 c ==================================================================
25
26 implicit none
27
28 c == global variables ==
29
30 #include "EEPARAMS.h"
31 #include "SIZE.h"
32 #include "PARAMS.h"
33 #include "GRID.h"
34
35 #include "cal.h"
36 #include "ecco_cost.h"
37
38 c == routine arguments ==
39
40 integer irec
41 integer mythid
42
43 c == local variables ==
44
45 integer bi,bj
46 integer i,j,k
47 integer itlo,ithi
48 integer jtlo,jthi
49 integer jmin,jmax
50 integer imin,imax
51 integer nobs
52 integer argotrec
53 integer beginargot
54 integer beginrun
55
56 _RL spval
57 _RL vartile
58
59 c == end of interface ==
60
61 parameter (spval = -9999. )
62 ce --> there is certainly a better place for this.
63
64 jtlo = mybylo(mythid)
65 jthi = mybyhi(mythid)
66 itlo = mybxlo(mythid)
67 ithi = mybxhi(mythid)
68 jmin = 1
69 jmax = sny
70 imin = 1
71 imax = snx
72
73 beginargot = argotstartdate(1)/10000
74 beginrun = modelstartdate(1)/10000
75 if ( beginargot .eq. beginrun ) then
76 argotrec = mod(modelstartdate(1)/100,100) -
77 & mod(argotstartdate(1)/100,100) + irec
78 else
79 argotrec = ( beginrun - beginargot - 1)*nmonthyear +
80 & (nmonthyear - mod(argotstartdate(1)/100,100) +
81 & 1) + mod(modelstartdate(1)/100,100) - 1 + irec
82 endif
83 if (argotrec.gt.0) then
84 call mdsreadfield( argotfile, cost_iprec, cost_yftype, nr,
85 & argotobs,argotrec, mythid )
86 else
87 do bj = jtlo,jthi
88 do bi = itlo,ithi
89 do k = 1,nr
90 do j = jmin,jmax
91 do i = imin,imax
92 argotobs(i,j,k,bi,bj)=spval
93 enddo
94 enddo
95 enddo
96 enddo
97 enddo
98 endif
99 return
100 end
101
  ViewVC Help
Powered by ViewVC 1.1.22