/[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.1 - (show annotations) (download)
Thu Nov 6 22:10:07 2003 UTC (20 years, 6 months ago) by heimbach
Branch: MAIN
CVS Tags: checkpoint58l_post, checkpoint57t_post, checkpoint57o_post, checkpoint52l_pre, checkpoint52e_pre, hrcube4, checkpoint58e_post, checkpoint57v_post, checkpoint52n_post, checkpoint52j_post, checkpoint53d_post, checkpoint58u_post, checkpoint54a_pre, checkpoint57f_post, checkpoint55c_post, checkpoint54e_post, checkpoint52e_post, checkpoint57s_post, checkpoint53c_post, checkpoint57k_post, checkpoint55d_pre, checkpoint57d_post, checkpoint57g_post, checkpoint57b_post, checkpoint57c_pre, checkpoint58r_post, checkpoint55j_post, checkpoint56b_post, checkpoint57i_post, checkpoint57y_post, checkpoint58g_post, checkpoint57x_post, checkpoint52j_pre, checkpoint54a_post, checkpoint57m_post, branch-netcdf, checkpoint52d_pre, checkpoint52l_post, checkpoint55h_post, checkpoint58n_post, checkpoint58x_post, checkpoint52k_post, checkpoint52b_pre, checkpoint57g_pre, checkpoint54b_post, checkpoint53b_pre, checkpoint55b_post, checkpoint58t_post, checkpoint58h_post, checkpoint57e_post, checkpoint58w_post, checkpoint54d_post, checkpoint53f_post, hrcube_1, checkpoint56c_post, checkpoint52m_post, checkpoint57y_pre, checkpoint55, checkpoint53a_post, checkpoint57f_pre, checkpoint57a_post, checkpoint54, checkpoint58q_post, checkpoint54f_post, checkpoint53b_post, checkpoint55g_post, checkpoint58j_post, checkpoint59e, checkpoint59d, checkpoint59g, checkpoint59f, checkpoint59a, checkpoint55f_post, checkpoint59c, checkpoint59b, checkpoint59h, checkpoint57r_post, checkpoint59, checkpoint58, checkpoint57a_pre, checkpoint55i_post, checkpoint57, checkpoint56, checkpoint53, checkpoint52, checkpoint52d_post, eckpoint57e_pre, checkpoint52a_post, checkpoint57h_done, checkpoint58f_post, checkpoint52b_post, checkpoint53g_post, checkpoint52f_post, checkpoint57n_post, checkpoint52c_post, checkpoint58d_post, checkpoint58c_post, checkpoint57w_post, checkpoint57p_post, checkpint57u_post, checkpoint58a_post, checkpoint58i_post, checkpoint57q_post, ecco_c52_e35, hrcube5, checkpoint58o_post, checkpoint57z_post, checkpoint57c_post, checkpoint58y_post, checkpoint55e_post, checkpoint52a_pre, checkpoint58k_post, checkpoint52i_post, checkpoint58v_post, checkpoint55a_post, checkpoint53d_pre, checkpoint54c_post, checkpoint58s_post, checkpoint58p_post, checkpoint57j_post, checkpoint58b_post, checkpoint57h_pre, checkpoint58m_post, checkpoint57l_post, checkpoint52i_pre, checkpoint51u_post, checkpoint52h_pre, checkpoint52f_pre, checkpoint57h_post, hrcube_2, hrcube_3, checkpoint56a_post, checkpoint55d_post
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_readargot.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_Readargot(
7 I irec,
8 I mythid
9 & )
10
11 c ==================================================================
12 c SUBROUTINE cost_Readargot
13 c ==================================================================
14 c
15 c o Read a given record of the
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_Readargot
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 argotrec
52 integer beginargot
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 beginargot = argotstartdate(1)/10000
73 beginrun = modelstartdate(1)/10000
74 if ( beginargot .eq. beginrun ) then
75 argotrec = mod(modelstartdate(1)/100,100) -
76 & mod(argotstartdate(1)/100,100) + irec
77 else
78 argotrec = ( beginrun - beginargot - 1)*nmonthyear +
79 & (nmonthyear - mod(argotstartdate(1)/100,100) +
80 & 1) + mod(modelstartdate(1)/100,100) - 1 + irec
81 endif
82 if (argotrec.gt.0) then
83 call mdsreadfield( argotfile, cost_iprec, cost_yftype, nr,
84 & argotobs,argotrec, mythid )
85 else
86 do bj = jtlo,jthi
87 do bi = itlo,ithi
88 do k = 1,nr
89 do j = jmin,jmax
90 do i = imin,imax
91 argotobs(i,j,k,bi,bj)=spval
92 enddo
93 enddo
94 enddo
95 enddo
96 enddo
97 endif
98 return
99 end
100
  ViewVC Help
Powered by ViewVC 1.1.22