/[MITgcm]/MITgcm/pkg/admtlm/admtlm_dsvd2model.F
ViewVC logotype

Diff of /MITgcm/pkg/admtlm/admtlm_dsvd2model.F

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

revision 1.2 by heimbach, Tue Nov 1 04:09:46 2005 UTC revision 1.3 by heimbach, Thu Apr 27 12:49:02 2006 UTC
# Line 5  C $Name$ Line 5  C $Name$
5  #include "PACKAGES_CONFIG.h"  #include "PACKAGES_CONFIG.h"
6  #include "CTRL_CPPOPTIONS.h"  #include "CTRL_CPPOPTIONS.h"
7    
8        subroutine admtlm_dsvd2model( first, mythid )        subroutine admtlm_dsvd2model(
9         &     first, postprocev, mythid )
10    
11  c     ==================================================================  c     ==================================================================
12  c     SUBROUTINE ctrl_unpack  c     SUBROUTINE ctrl_unpack
# Line 61  c     == global variables == Line 62  c     == global variables ==
62  c     == routine arguments ==  c     == routine arguments ==
63    
64        logical first        logical first
65          logical postprocev
66        integer mythid        integer mythid
67    
68  #ifndef EXCLUDE_CTRL_PACK  #ifndef EXCLUDE_CTRL_PACK
# Line 112  c--   Initialise Line 114  c--   Initialise
114        nbuffGlobal = 0        nbuffGlobal = 0
115    
116  cph-new(  cph-new(
117        yadprefix = 'g_'        if ( postprocev ) then
118             yadprefix = 'ev'
119          else
120             yadprefix = 'g_'
121          endif
122        nveccount = 0        nveccount = 0
123  cph-new)  cph-new)
124    
# Line 179  c     >>> Write gradient vector <<< Line 185  c     >>> Write gradient vector <<<
185    
186  #ifdef ALLOW_ADMTLM  #ifdef ALLOW_ADMTLM
187    
188         write(cfile(1:128),'(a)') ' '         if (postprocev) then
189         write(cfile,'(a,i4.4)')  cph do a dummy read of initialized EV fields
190       &      'admtlm_vector.it', optimcycle  cph they will be overwritten by array phtmpadmtlm
191              write(cfile(1:128),'(a)') ' '
192              write(cfile,'(a,i4.4)')
193         &         'admtlm_eigen', optimcycle
194           else
195              write(cfile(1:128),'(a)') ' '
196              write(cfile,'(a,i4.4)')
197         &         'admtlm_vector.it', optimcycle
198           endif
199         print *, 'ph-pack: unpacking ', cfile         print *, 'ph-pack: unpacking ', cfile
200         open( cunit, file   = cfile,  cph       open( cunit, file   = cfile,
201       &      status = 'old',  cph     &      status = 'old',
202       &      form   = 'unformatted',  cph     &      form   = 'unformatted',
203       &      access  = 'sequential'   )  cph     &      access  = 'sequential'   )
204    
205  #else /* ndef ALLOW_ADMTLM */  #else /* ndef ALLOW_ADMTLM */
206    
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22