pkg/cost/cost_readargot.F

C $Header: /u/gcmpack/development/heimbach/ecco_env/pkg/cost/cost_readsstfields.F,v 1.5 2001/02/12 19:40:54 heimbach Exp $

#include "COST_CPPOPTIONS.h"


      subroutine cost_Readargot(
     I                               irec,
     I                               mythid
     &                             )

c     ==================================================================
c     SUBROUTINE cost_Readargot
c     ==================================================================
c
c     o Read a given record of the 
c
c     started:Armin Koehl akoehl@ucsd.edu
c              - Restructured the code in order to create a package
c                for the MITgcmUV.
c
c     ==================================================================
c     SUBROUTINE cost_Readargot
c     ==================================================================

      implicit none

c     == global variables ==

#include "EEPARAMS.h"
#include "SIZE.h"
#include "PARAMS.h"
#include "GRID.h"

#include "cal.h"
#include "cost.h"

c     == routine arguments ==

      integer irec
      integer mythid 

c     == local variables ==

      integer bi,bj
      integer i,j,k
      integer itlo,ithi
      integer jtlo,jthi
      integer jmin,jmax
      integer imin,imax
      integer nobs
      integer argotrec
      integer beginargot
      integer beginrun

      _RL spval
      _RL vartile

c     == end of interface ==

      parameter (spval = -9999. ) 
ce    --> there is certainly a better place for this.

      jtlo = mybylo(mythid)
      jthi = mybyhi(mythid)
      itlo = mybxlo(mythid)
      ithi = mybxhi(mythid)
      jmin = 1
      jmax = sny
      imin = 1
      imax = snx

      beginargot = argotstartdate(1)/10000
      beginrun = modelstartdate(1)/10000
      if ( beginargot .eq. beginrun ) then
        argotrec = mod(modelstartdate(1)/100,100) -
     &           mod(argotstartdate(1)/100,100) + irec
      else
        argotrec = ( beginrun - beginargot - 1)*nmonthyear +
     &           (nmonthyear - mod(argotstartdate(1)/100,100) +
     &            1) + mod(modelstartdate(1)/100,100) - 1 + irec
      endif
      if (argotrec.gt.0) then
      call mdsreadfield( argotfile, cost_iprec, cost_yftype, nr, 
     &                   argotobs,argotrec, mythid )
      else         
         do bj = jtlo,jthi
            do bi = itlo,ithi
               do k = 1,nr
               do j = jmin,jmax
                  do i = imin,imax
                     argotobs(i,j,k,bi,bj)=spval
                  enddo
               enddo
               enddo
            enddo
         enddo
      endif
      return
      end

1	C $Header: /u/gcmpack/development/heimbach/ecco_env/pkg/cost/cost_readsstfields.F,v 1.5 2001/02/12 19:40:54 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 "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