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 |
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 |
|
|
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 |
|
|
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 |
|
|
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 |
|
|