/[MITgcm]/MITgcm/pkg/exf/exf_mapfields.F
ViewVC logotype

Diff of /MITgcm/pkg/exf/exf_mapfields.F

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

revision 1.20 by jmc, Mon May 14 19:34:57 2007 UTC revision 1.22 by jmc, Mon Oct 1 13:36:05 2007 UTC
# Line 94  c     == end of interface == Line 94  c     == end of interface ==
94       &                      + act4*max1*max2*max3       &                      + act4*max1*max2*max3
95  #endif /* ALLOW_AUTODIFF_TAMC */  #endif /* ALLOW_AUTODIFF_TAMC */
96    
97    c     Heat flux.
98            do j = jmin,jmax            do j = jmin,jmax
99              do i = imin,imax              do i = imin,imax
 c            Heat flux.  
100               qnet(i,j,bi,bj) = exf_outscal_hflux*hflux(i,j,bi,bj)               qnet(i,j,bi,bj) = exf_outscal_hflux*hflux(i,j,bi,bj)
101               if ( hfluxfile .EQ. ' ' )              enddo
102       &            qnet(i,j,bi,bj) = qnet(i,j,bi,bj) -            enddo
103              if ( hfluxfile .EQ. ' ' ) then
104               do j = jmin,jmax
105                do i = imin,imax
106                      qnet(i,j,bi,bj) = qnet(i,j,bi,bj) -
107       &            exf_outscal_hflux * ( hflux_exfremo_intercept +       &            exf_outscal_hflux * ( hflux_exfremo_intercept +
108       &            hflux_exfremo_slope*(mytime-starttime) )       &            hflux_exfremo_slope*(mytime-starttime) )
109              enddo              enddo
110            enddo             enddo
111              endif
112    
113    c     Salt flux.
114            do j = jmin,jmax            do j = jmin,jmax
115              do i = imin,imax              do i = imin,imax
116  c            Salt flux.               EmPmR(i,j,bi,bj)= exf_outscal_sflux*sflux(i,j,bi,bj)
117               empmr(i,j,bi,bj)= exf_outscal_sflux*sflux(i,j,bi,bj)       &                                          *rhoConstFresh
118               if ( sfluxfile .EQ. ' ' )              enddo
119       &            empmr(i,j,bi,bj) = empmr(i,j,bi,bj) -            enddo
120              if ( sfluxfile .EQ. ' ' ) then
121               do j = jmin,jmax
122                do i = imin,imax
123                     EmPmR(i,j,bi,bj) = EmPmR(i,j,bi,bj) - rhoConstFresh*
124       &            exf_outscal_sflux * ( sflux_exfremo_intercept +       &            exf_outscal_sflux * ( sflux_exfremo_intercept +
125       &            sflux_exfremo_slope*(mytime-starttime) )       &            sflux_exfremo_slope*(mytime-starttime) )
126              enddo              enddo
127            enddo             enddo
128              endif
129    
130  #ifdef ALLOW_AUTODIFF_TAMC  #ifdef ALLOW_AUTODIFF_TAMC
131  CADJ STORE ustress(:,:,bi,bj) = comlev1_bibj, key=ikey, byte=isbyte  CADJ STORE ustress(:,:,bi,bj) = comlev1_bibj, key=ikey, byte=isbyte
Legend:
Removed from v.1.20  
changed lines
  Added in v.1.22
  ViewVC Help
Powered by ViewVC 1.1.22