/[MITgcm]/MITgcm/pkg/ctrl/ctrl_hfacc_ini.F
ViewVC logotype

Diff of /MITgcm/pkg/ctrl/ctrl_hfacc_ini.F

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

revision 1.4 by heimbach, Mon May 14 22:02:33 2007 UTC revision 1.5 by jmc, Tue Oct 9 00:00:00 2007 UTC
# Line 1  Line 1 
1  C $Header$  C $Header$
2    C $Name$
3    
4  #include "CTRL_CPPOPTIONS.h"  #include "CTRL_CPPOPTIONS.h"
5    
# Line 9  C     !INTERFACE: Line 10  C     !INTERFACE:
10    
11  C     !DESCRIPTION: \bv  C     !DESCRIPTION: \bv
12  c     *=================================================================  c     *=================================================================
13  c     | SUBROUTINE ctrl_hfacc_ini  c     | SUBROUTINE ctrl_hfacc_ini
14  c     | Add the hFacC part of the control vector to the model state  c     | Add the hFacC part of the control vector to the model state
15  c     | and update the tile halos.  c     | and update the tile halos.
16  c     | The control vector is defined in the header file "ctrl.h".  c     | The control vector is defined in the header file "ctrl.h".
17  c     *=================================================================  c     *=================================================================
18  C     \ev  C     \ev
# Line 79  CEOP Line 80  CEOP
80          fac = 0. _d 0          fac = 0. _d 0
81        endif        endif
82    
83  Cml      write(msgbuf,'(a)')  Cml      write(msgbuf,'(a)')
84  Cml     &     'ctrl_hfacc_ini: Re-initialising hFacC,'  Cml     &     'ctrl_hfacc_ini: Re-initialising hFacC,'
85  Cml      call print_message( msgbuf, standardmessageunit,  Cml      call print_message( msgbuf, standardmessageunit,
86  Cml     &                    SQUEEZE_RIGHT , mythid)  Cml     &                    SQUEEZE_RIGHT , mythid)
87  Cml      write(msgbuf,'(a)')  Cml      write(msgbuf,'(a)')
88  Cml     &     '                adding the control vector.'  Cml     &     '                adding the control vector.'
89  Cml      call print_message( msgbuf, standardmessageunit,  Cml      call print_message( msgbuf, standardmessageunit,
90  Cml     &                    SQUEEZE_RIGHT , mythid)  Cml     &                    SQUEEZE_RIGHT , mythid)
91        write(standardmessageunit,'(21x,a)')        write(standardmessageunit,'(21x,a)')
92       &     'ctrl_hfacc_ini: Re-initialising hFacC,'       &     'ctrl_hfacc_ini: Re-initialising hFacC,'
93        write(standardmessageunit,'(21x,a)')        write(standardmessageunit,'(21x,a)')
94       &     '                adding the control vector.'       &     '                adding the control vector.'
95    
96  C     Re-initialize hFacC, so that TAMC/TAF can see it  C     Re-initialize hFacC, so that TAMC/TAF can see it
# Line 113  C     this will be no longer necessary Line 114  C     this will be no longer necessary
114  Cml      _EXCH_XYZ_R8( hFacC ,myThid )  Cml      _EXCH_XYZ_R8( hFacC ,myThid )
115        _EXCH_XYZ_R4( hFacC ,myThid )        _EXCH_XYZ_R4( hFacC ,myThid )
116    
117  C--    C--
118        il=ilnblnk( xx_hfacc_file )        il=ilnblnk( xx_hfacc_file )
119        write(fnamehfacc(1:80),'(2a,i10.10)')        write(fnamehfacc(1:80),'(2a,i10.10)')
120       &     xx_hfacc_file(1:il),'.',optimcycle       &     xx_hfacc_file(1:il),'.',optimcycle
# Line 143  C-- Line 144  C--
144                 do i = imin,imax                 do i = imin,imax
145                    k = k_lowC(i,j,bi,bj)                    k = k_lowC(i,j,bi,bj)
146  c                  if ( k .gt. 0 ) then  c                  if ( k .gt. 0 ) then
147                       hFacC(i,j,k,bi,bj) = hFacC(i,j,k,bi,bj)                       hFacC(i,j,k,bi,bj) = hFacC(i,j,k,bi,bj)
148       &                                  + fac*tmpfld2d(i,j,bi,bj)       &                                  + fac*tmpfld2d(i,j,bi,bj)
149  c                  end if  c                  end if
150                 enddo                 enddo
Legend:
Removed from v.1.4  
changed lines
  Added in v.1.5
  ViewVC Help
Powered by ViewVC 1.1.22