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

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

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

revision 1.10 by mlosch, Mon Apr 12 09:00:34 2010 UTC revision 1.11 by gforget, Mon Oct 25 16:21:58 2010 UTC
# Line 29  C     == global variables == Line 29  C     == global variables ==
29    
30  #ifdef ALLOW_AUTODIFF_TAMC  #ifdef ALLOW_AUTODIFF_TAMC
31  #include "tamc.h"  #include "tamc.h"
32    #include "tamc_keys.h"
33  #endif  #endif
34    
35  C     == routine arguments ==  C     == routine arguments ==
# Line 61  C     Loop over tiles. Line 62  C     Loop over tiles.
62        DO bj = myByLo(myThid),myByHi(myThid)        DO bj = myByLo(myThid),myByHi(myThid)
63         DO bi = myBxLo(myThid),myBxHi(myThid)         DO bi = myBxLo(myThid),myBxHi(myThid)
64    
65    #ifdef ALLOW_AUTODIFF_TAMC
66              act1 = bi - myBxLo(myThid)
67              max1 = myBxHi(myThid) - myBxLo(myThid) + 1
68              act2 = bj - myByLo(myThid)
69              max2 = myByHi(myThid) - myByLo(myThid) + 1
70              act3 = myThid - 1
71              max3 = nTx*nTy
72              act4 = ikey_dynamics - 1
73              ikey = (act1 + 1) + act2*max1
74         &                      + act3*max1*max2
75         &                      + act4*max1*max2*max3
76    #endif /* ALLOW_AUTODIFF_TAMC */
77    
78  C--   Initialise  C--   Initialise
79          DO j = 1,sNy          DO j = 1,sNy
80           DO i = 1,sNx           DO i = 1,sNx
# Line 176  C--   infer wind field from wind-speed & Line 190  C--   infer wind field from wind-speed &
190          ENDDO          ENDDO
191  #endif /* ifndef ALLOW_ATM_WIND */  #endif /* ifndef ALLOW_ATM_WIND */
192    
193    #ifdef ALLOW_AUTODIFF_TAMC
194    CADJ STORE wspeed(:,:,bi,bj) = comlev1_bibj, key=ikey, byte=isbyte
195    #endif
196    
197  C--   set wind-speed lower limit  C--   set wind-speed lower limit
198          DO j = 1,sNy          DO j = 1,sNy
199           DO i = 1,sNx           DO i = 1,sNx
Legend:
Removed from v.1.10  
changed lines
  Added in v.1.11
  ViewVC Help
Powered by ViewVC 1.1.22