/[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.11 by jmc, Fri Jul 16 01:20:57 2004 UTC revision 1.16 by heimbach, Thu Mar 2 15:30:11 2006 UTC
# Line 2  c $Header$ Line 2  c $Header$
2    
3  #include "EXF_OPTIONS.h"  #include "EXF_OPTIONS.h"
4    
5        subroutine exf_mapfields( mythid )        subroutine exf_mapfields( mytime, myiter, mythid )
6    
7  c     ==================================================================  c     ==================================================================
8  c     SUBROUTINE exf_mapfields  c     SUBROUTINE exf_mapfields
# Line 47  c     == global variables == Line 47  c     == global variables ==
47    
48  #include "EEPARAMS.h"  #include "EEPARAMS.h"
49  #include "SIZE.h"  #include "SIZE.h"
50    #include "PARAMS.h"
51  #include "FFIELDS.h"  #include "FFIELDS.h"
52  #include "GRID.h"  #include "GRID.h"
53    
54  #include "exf_param.h"  #include "exf_param.h"
55  #include "exf_constants.h"  #include "exf_constants.h"
56  #include "exf_fields.h"  #include "exf_fields.h"
57    #include "exf_clim_param.h"
58  #include "exf_clim_fields.h"  #include "exf_clim_fields.h"
59  #ifdef ALLOW_AUTODIFF_TAMC  #ifdef ALLOW_AUTODIFF_TAMC
60  # include "tamc.h"  # include "tamc.h"
# Line 63  c     == routine arguments == Line 65  c     == routine arguments ==
65  c     mythid - thread number for this instance of the routine.  c     mythid - thread number for this instance of the routine.
66    
67        integer mythid        integer mythid
68          integer myiter
69          _RL     mytime
70    
71  c     == local variables ==  c     == local variables ==
72    
73        integer bi,bj        integer bi,bj
74        integer i,j        integer i,j,k
75        integer jtlo        integer jtlo
76        integer jthi        integer jthi
77        integer itlo        integer itlo
# Line 108  c     == end of interface == Line 112  c     == end of interface ==
112              do i = imin,imax              do i = imin,imax
113  c            Heat flux.  c            Heat flux.
114               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)
115                 if ( hfluxfile .EQ. ' ' )
116         &            qnet(i,j,bi,bj) = qnet(i,j,bi,bj) -
117         &            exf_outscal_hflux * ( hflux_exfremo_intercept +
118         &            hflux_exfremo_slope*(mytime-starttime) )
119              enddo              enddo
120            enddo            enddo
121    
# Line 116  c            Heat flux. Line 124  c            Heat flux.
124              do i = imin,imax              do i = imin,imax
125  c            Salt flux.  c            Salt flux.
126               empmr(i,j,bi,bj)= exf_outscal_sflux*sflux(i,j,bi,bj)               empmr(i,j,bi,bj)= exf_outscal_sflux*sflux(i,j,bi,bj)
127                 if ( sfluxfile .EQ. ' ' )
128         &            empmr(i,j,bi,bj) = empmr(i,j,bi,bj) -
129         &            exf_outscal_sflux * ( sflux_exfremo_intercept +
130         &            sflux_exfremo_slope*(mytime-starttime) )
131              enddo              enddo
132            enddo            enddo
133    
# Line 194  c             Short wave radiative flux. Line 206  c             Short wave radiative flux.
206               qsw(i,j,bi,bj)  = exf_outscal_swflux*swflux(i,j,bi,bj)               qsw(i,j,bi,bj)  = exf_outscal_swflux*swflux(i,j,bi,bj)
207              enddo              enddo
208            enddo            enddo
 C         net heat flux = heat flux (except SW) + SW flux.  
           do j = jmin,jmax  
             do i = imin,imax  
              qnet(i,j,bi,bj) = qnet(i,j,bi,bj) + qsw(i,j,bi,bj)  
             enddo  
           enddo  
209  #endif  #endif
210    
211  #ifdef ALLOW_CLIMSST_RELAXATION  #ifdef ALLOW_CLIMSST_RELAXATION
Legend:
Removed from v.1.11  
changed lines
  Added in v.1.16
  ViewVC Help
Powered by ViewVC 1.1.22