/[MITgcm]/MITgcm/optim/optim_readdata.F
ViewVC logotype

Diff of /MITgcm/optim/optim_readdata.F

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

revision 1.2 by heimbach, Fri Nov 15 04:03:25 2002 UTC revision 1.7 by dfer, Tue Jan 15 16:36:32 2008 UTC
# Line 37  cgg   Include ECCO_CPPOPTIONS because th Line 37  cgg   Include ECCO_CPPOPTIONS because th
37  cgg   have headers with options for OBCS masks.  cgg   have headers with options for OBCS masks.
38  #include "ECCO_CPPOPTIONS.h"  #include "ECCO_CPPOPTIONS.h"
39    
 #include "ecco.h"  
40  #include "ctrl.h"  #include "ctrl.h"
41  #include "optim.h"  #include "optim.h"
42  #include "minimization.h"  #include "minimization.h"
# Line 75  c     == local variables == Line 74  c     == local variables ==
74        integer funit        integer funit
75    
76        integer cbuffindex        integer cbuffindex
77        _RL     cbuff( sNx*nSx*nPx*sNy*nSy*nPy )        real*4 cbuff( sNx*nSx*nPx*sNy*nSy*nPy )
78    
79        character*(128) fname        character*(128) fname
80    
81        integer         filei  c      integer         filei
82        integer         filej  c      integer         filej
83        integer         filek  c      integer         filek
84    c      integer         fileiG
85    c      integer         filejG
86    c      integer         filensx
87    c      integer         filensy
88        integer         filenopt        integer         filenopt
       integer         fileig  
       integer         filejg  
       integer         filensx  
       integer         filensy  
89        _RL             fileff        _RL             fileff
90    
91  cgg(  cgg(
# Line 97  cgg) Line 96  cgg)
96    
97  c     == end of interface ==  c     == end of interface ==
98    
99          print *, 'pathei-lsopt in optim_readdata'
100    
101  c--   The reference i/o unit.  c--   The reference i/o unit.
102        funit = 20        funit = 20
103    
# Line 129  c--   Read the data. Line 130  c--   Read the data.
130    
131  c--   Generate file name and open the file.  c--   Generate file name and open the file.
132        write(fname(1:128),'(4a,i4.4)')        write(fname(1:128),'(4a,i4.4)')
133       &     dfile,'_',expId(1:10),'.opt', nopt       &     dfile,'_',yctrlid(1:10),'.opt', nopt
134        open( funit, file   = fname,        open( funit, file   = fname,
135       &     status = 'old',       &     status = 'old',
136       &     form   = 'unformatted',       &     form   = 'unformatted',
# Line 139  c--   Generate file name and open the fi Line 140  c--   Generate file name and open the fi
140  c--   Read the header.  c--   Read the header.
141        read( funit ) nvartype        read( funit ) nvartype
142        read( funit ) nvarlength        read( funit ) nvarlength
143        read( funit ) expId        read( funit ) yctrlid
144        read( funit ) filenopt        read( funit ) filenopt
145        read( funit ) fileff        read( funit ) fileff
146        read( funit ) fileiG        read( funit ) fileiG
# Line 147  c--   Read the header. Line 148  c--   Read the header.
148        read( funit ) filensx        read( funit ) filensx
149        read( funit ) filensy        read( funit ) filensy
150    
151  cph(        read( funit ) (nWetcGlobal(k), k=1,nr)
152        print *,'ph-opt 1 ', nvartype, nvarlength, filensx, filensy        read( funit ) (nWetsGlobal(k), k=1,nr)
153  cph)        read( funit ) (nWetwGlobal(k), k=1,nr)
154    #ifdef ALLOW_CTRL_WETV
155        read( funit ) (((nWetcTile(i,j,k), i=1,nsx), j=1,nsy),        read( funit ) (nWetvGlobal(k), k=1,nr)
156       &     k=1,nr)  #endif
       read( funit ) (((nWetsTile(i,j,k), i=1,nsx), j=1,nsy),  
      &     k=1,nr)  
       read( funit ) (((nWetwTile(i,j,k), i=1,nsx), j=1,nsy),  
      &     k=1,nr)  
157    
158  cgg(    Add OBCS Mask information into the header section for optimization.  cgg(    Add OBCS Mask information into the header section for optimization.
159  #ifdef ALLOW_OBCSN_CONTROL  #ifdef ALLOW_OBCSN_CONTROL
160            read(funit) ((((nWetobcsn(i,j,k,iobcs), k=1,nr),        read( funit ) ((nWetobcsnGlo(k,iobcs), k=1,nr),iobcs= 1,nobcs)
      &          iobcs= 1,nobcs), i=1,nsx) , j=1,nsy)  
161  #endif  #endif
162  #ifdef ALLOW_OBCSS_CONTROL  #ifdef ALLOW_OBCSS_CONTROL
163            read(funit) ((((nWetobcss(i,j,k,iobcs), k=1,nr),        read( funit ) ((nWetobcssGlo(k,iobcs), k=1,nr),iobcs= 1,nobcs)
      &          iobcs= 1,nobcs), i=1,nsx) , j=1,nsy)  
164  #endif  #endif
165  #ifdef ALLOW_OBCSW_CONTROL  #ifdef ALLOW_OBCSW_CONTROL
166            read(funit) ((((nWetobcsw(i,j,k,iobcs), k=1,nr),        read( funit ) ((nWetobcswGlo(k,iobcs), k=1,nr),iobcs= 1,nobcs)
      &          iobcs= 1,nobcs), i=1,nsx) , j=1,nsy)  
167  #endif  #endif
168  #ifdef ALLOW_OBCSE_CONTROL  #ifdef ALLOW_OBCSE_CONTROL
169            read(funit) ((((nWetobcse(i,j,k,iobcs), k=1,nr),        read( funit ) ((nWetobcseGlo(k,iobcs), k=1,nr),iobcs= 1,nobcs)
      &          iobcs= 1,nobcs), i=1,nsx) , j=1,nsy)  
170  #endif  #endif
171  cgg)  cgg)
172        read( funit ) (ncvarindex(i), i=1,maxcvars)        read( funit ) (ncvarindex(i), i=1,maxcvars)
# Line 188  cph( Line 181  cph(
181  cph      if (lheaderonly) then  cph      if (lheaderonly) then
182           print *, 'pathei: nvartype ', nvartype           print *, 'pathei: nvartype ', nvartype
183           print *, 'pathei: nvarlength ', nvarlength           print *, 'pathei: nvarlength ', nvarlength
184           print *, 'pathei: expId ', expId           print *, 'pathei: yctrlid ', yctrlid
185           print *, 'pathei: filenopt ', filenopt           print *, 'pathei: filenopt ', filenopt
186           print *, 'pathei: fileff ', fileff           print *, 'pathei: fileff ', fileff
187           print *, 'pathei: fileiG ', fileiG           print *, 'pathei: fileiG ', fileiG
# Line 196  cph      if (lheaderonly) then Line 189  cph      if (lheaderonly) then
189           print *, 'pathei: filensx ', filensx           print *, 'pathei: filensx ', filensx
190           print *, 'pathei: filensy ', filensy           print *, 'pathei: filensy ', filensy
191                    
192           print *, 'pathei: nWetcTile ',           print *, 'pathei: nWetcGlobal ',
193       &        (((nWetcTile(i,j,k), i=1,nsx), j=1,nsy), k=1,nr)       &        (nWetcGlobal(k),  k=1,nr)
194           print *, 'pathei: nWetsTile ',           print *, 'pathei: nWetsGlobal ',
195       &        (((nWetsTile(i,j,k), i=1,nsx), j=1,nsy), k=1,nr)       &        (nWetsGlobal(k),  k=1,nr)
196           print *, 'pathei: nWetwTile ',           print *, 'pathei: nWetwGlobal ',
197       &        (((nWetwTile(i,j,k), i=1,nsx), j=1,nsy), k=1,nr)       &        (nWetwGlobal(k),  k=1,nr)
198             print *, 'pathei: nWetvGlobal ',
199         &        (nWetvGlobal(k),  k=1,nr)
200           print *, 'pathei: ncvarindex ',           print *, 'pathei: ncvarindex ',
201       &        (ncvarindex(i), i=1,maxcvars)       &        (ncvarindex(i), i=1,maxcvars)
202           print *, 'pathei: ncvarrecs ',           print *, 'pathei: ncvarrecs ',
# Line 267  c--   Read the data. Line 262  c--   Read the data.
262                          do k = 1,ncvarnrmax(icvar)                          do k = 1,ncvarnrmax(icvar)
263                             cbuffindex = 0                             cbuffindex = 0
264                             if (ncvargrd(icvar) .eq. 'c') then                             if (ncvargrd(icvar) .eq. 'c') then
265                                cbuffindex = nwetctile(bi,bj,k)                                cbuffindex = nWetcGlobal(k)
266                             else if (ncvargrd(icvar) .eq. 's') then                             else if (ncvargrd(icvar) .eq. 's') then
267                                cbuffindex = nwetstile(bi,bj,k)                                cbuffindex = nWetsGlobal(k)
268                             else if (ncvargrd(icvar) .eq. 'w') then                             else if (ncvargrd(icvar) .eq. 'w') then
269                                cbuffindex = nwetwtile(bi,bj,k)                                cbuffindex = nWetwGlobal(k)
270                               else if (ncvargrd(icvar) .eq. 'v') then
271                                  cbuffindex = nWetvGlobal(k)
272  cgg(   O.B. points have the grid mask "m".  cgg(   O.B. points have the grid mask "m".
273                             else if (ncvargrd(icvar) .eq. 'm') then                             else if (ncvargrd(icvar) .eq. 'm') then
274  cgg    From "icvrec", calculate what iobcs must be.  cgg    From "icvrec", calculate what iobcs must be.
# Line 280  cgg    From "icvrec", calculate what iob Line 277  cgg    From "icvrec", calculate what iob
277                               iobcs= icvrec - igg*nobcs                               iobcs= icvrec - igg*nobcs
278  #ifdef ALLOW_OBCSN_CONTROL  #ifdef ALLOW_OBCSN_CONTROL
279                               if (icvar .eq. 11) then                                                   if (icvar .eq. 11) then                    
280                                 cbuffindex = nwetobcsn(bi,bj,k,iobcs)                                 cbuffindex = nWetobcsnGlo(k,iobcs)
281                               endif                               endif
282  #endif  #endif
283  #ifdef ALLOW_OBCSS_CONTROL  #ifdef ALLOW_OBCSS_CONTROL
284                               if (icvar .eq. 12) then                               if (icvar .eq. 12) then
285                                 cbuffindex = nwetobcss(bi,bj,k,iobcs)                                 cbuffindex = nWetobcssGlo(k,iobcs)
286                               endif                               endif
287  #endif  #endif
288  #ifdef ALLOW_OBCSW_CONTROL  #ifdef ALLOW_OBCSW_CONTROL
289                               if (icvar .eq. 13) then                               if (icvar .eq. 13) then
290                                 cbuffindex = nwetobcsw(bi,bj,k,iobcs)                                 cbuffindex = nWetobcswGlo(k,iobcs)
291                               endif                               endif
292  #endif  #endif
293  #ifdef ALLOW_OBCSE_CONTROL  #ifdef ALLOW_OBCSE_CONTROL
294                               if (icvar .eq. 14) then                               if (icvar .eq. 14) then
295                                 cbuffindex = nwetobcse(bi,bj,k,iobcs)                                 cbuffindex = nWetobcseGlo(k,iobcs)
296                               endif                               endif
297  #endif  #endif
298  cgg)  cgg)
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.7
  ViewVC Help
Powered by ViewVC 1.1.22