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

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

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

revision 1.6 by heimbach, Tue May 4 22:13:08 2004 UTC revision 1.7 by heimbach, Thu Nov 11 06:28:33 2004 UTC
# Line 64  c     == local variables == Line 64  c     == local variables ==
64        integer imin,imax        integer imin,imax
65        integer ilgen        integer ilgen
66    
67          _RL     gensign
68        _RL     genfac        _RL     genfac
69        logical doCtrlUpdate        logical doCtrlUpdate
70        logical genfirst        logical genfirst
# Line 144  cph this flag ported from the SIO code ( Line 145  cph this flag ported from the SIO code (
145        else        else
146           doCtrlUpdate = .TRUE.           doCtrlUpdate = .TRUE.
147        endif        endif
148          if ( xx_gen_file .EQ. xx_tauu_file .OR.
149         &     xx_gen_file .EQ. xx_tauv_file ) then
150             gensign = -1.
151          else
152             gensign = 1.
153          endif
154  c  c
155  cph since the above is ECCO specific, we undo it here:  cph since the above is ECCO specific, we undo it here:
156        doCtrlUpdate = .TRUE.  cph      doCtrlUpdate = .TRUE.
157  c  c
158        if ( doCtrlUpdate ) then        if ( doCtrlUpdate ) then
159  cph)  cph)
# Line 157  c--       Calculate mask for tracer cell Line 164  c--       Calculate mask for tracer cell
164            do j = 1,sny            do j = 1,sny
165              do i = 1,snx              do i = 1,snx
166                genfld(i,j,bi,bj) = genfld (i,j,bi,bj)                genfld(i,j,bi,bj) = genfld (i,j,bi,bj)
167       &                         + genfac            *xx_gen0(i,j,bi,bj)       &              + gensign*genfac            *xx_gen0(i,j,bi,bj)
168       &                         + (1. _d 0 - genfac)*xx_gen1(i,j,bi,bj)       &              + gensign*(1. _d 0 - genfac)*xx_gen1(i,j,bi,bj)
169                genfld(i,j,bi,bj) = genfld(i,j,bi,bj)*genmask(i,j,k,bi,bj)                genfld(i,j,bi,bj) = genfld(i,j,bi,bj)*genmask(i,j,k,bi,bj)
170              enddo              enddo
171            enddo            enddo
Legend:
Removed from v.1.6  
changed lines
  Added in v.1.7
  ViewVC Help
Powered by ViewVC 1.1.22