/[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.8 by heimbach, Wed Jan 12 20:33:13 2005 UTC revision 1.9 by heimbach, Thu Mar 2 02:53:23 2006 UTC
# Line 7  C $Name$ Line 7  C $Name$
7        subroutine ctrl_get_gen(        subroutine ctrl_get_gen(
8       I          xx_gen_file, xx_genstartdate, xx_genperiod,       I          xx_gen_file, xx_genstartdate, xx_genperiod,
9       I          genmask, genfld, xx_gen0, xx_gen1, xx_gen_dummy,       I          genmask, genfld, xx_gen0, xx_gen1, xx_gen_dummy,
10         I          xx_gen_remo_intercept, xx_gen_remo_slope,
11       I          mytime, myiter, mythid       I          mytime, myiter, mythid
12       &                     )       &                     )
13    
# Line 47  c     == routine arguments == Line 48  c     == routine arguments ==
48        _RL     xx_gen0(1-olx:snx+olx,1-oly:sny+oly,nsx,nsy)        _RL     xx_gen0(1-olx:snx+olx,1-oly:sny+oly,nsx,nsy)
49        _RL     xx_gen1(1-olx:snx+olx,1-oly:sny+oly,nsx,nsy)        _RL     xx_gen1(1-olx:snx+olx,1-oly:sny+oly,nsx,nsy)
50        _RL     xx_gen_dummy        _RL     xx_gen_dummy
51          _RL     xx_gen_remo_intercept
52          _RL     xx_gen_remo_slope
53    
54        _RL     mytime        _RL     mytime
55        integer myiter        integer myiter
# Line 167  c--       Calculate mask for tracer cell Line 170  c--       Calculate mask for tracer cell
170                genfld(i,j,bi,bj) = genfld (i,j,bi,bj)                genfld(i,j,bi,bj) = genfld (i,j,bi,bj)
171       &              + gensign*genfac            *xx_gen0(i,j,bi,bj)       &              + gensign*genfac            *xx_gen0(i,j,bi,bj)
172       &              + gensign*(1. _d 0 - genfac)*xx_gen1(i,j,bi,bj)       &              + gensign*(1. _d 0 - genfac)*xx_gen1(i,j,bi,bj)
173                genfld(i,j,bi,bj) = genfld(i,j,bi,bj)*genmask(i,j,k,bi,bj)                genfld(i,j,bi,bj) =
174         &             genmask(i,j,k,bi,bj)*( genfld (i,j,bi,bj) -
175         &             ( xx_gen_remo_intercept +
176         &               xx_gen_remo_slope*(mytime-starttime) ) )
177              enddo              enddo
178            enddo            enddo
179          enddo          enddo
Legend:
Removed from v.1.8  
changed lines
  Added in v.1.9
  ViewVC Help
Powered by ViewVC 1.1.22