/[MITgcm]/MITgcm/model/src/thermodynamics.F
ViewVC logotype

Diff of /MITgcm/model/src/thermodynamics.F

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

revision 1.40 by adcroft, Tue May 13 17:42:00 2003 UTC revision 1.41 by heimbach, Mon Jun 23 22:32:02 2003 UTC
# Line 173  C--   dummy statement to end declaration Line 173  C--   dummy statement to end declaration
173        itdkey = 1        itdkey = 1
174  #endif /* ALLOW_AUTODIFF_TAMC */  #endif /* ALLOW_AUTODIFF_TAMC */
175    
 C--   Set up work arrays with valid (i.e. not NaN) values  
 C     These inital values do not alter the numerical results. They  
 C     just ensure that all memory references are to valid floating  
 C     point numbers. This prevents spurious hardware signals due to  
 C     uninitialised but inert locations.  
       DO j=1-OLy,sNy+OLy  
        DO i=1-OLx,sNx+OLx  
         xA(i,j)      = 0. _d 0  
         yA(i,j)      = 0. _d 0  
         uTrans(i,j)  = 0. _d 0  
         vTrans(i,j)  = 0. _d 0  
         rhok   (i,j) = 0. _d 0  
         phiSurfX(i,j) = 0. _d 0  
         phiSurfY(i,j) = 0. _d 0  
        ENDDO  
       ENDDO  
   
   
176  #ifdef ALLOW_AUTODIFF_TAMC  #ifdef ALLOW_AUTODIFF_TAMC
177  C--   HPF directive to help TAMC  C--   HPF directive to help TAMC
178  CHPF$ INDEPENDENT  CHPF$ INDEPENDENT
# Line 221  CHPF$&                  ) Line 203  CHPF$&                  )
203       &                      + act4*max1*max2*max3       &                      + act4*max1*max2*max3
204  #endif /* ALLOW_AUTODIFF_TAMC */  #endif /* ALLOW_AUTODIFF_TAMC */
205    
206  C--     Set up work arrays that need valid initial values  C--   Set up work arrays with valid (i.e. not NaN) values
207    C     These inital values do not alter the numerical results. They
208    C     just ensure that all memory references are to valid floating
209    C     point numbers. This prevents spurious hardware signals due to
210    C     uninitialised but inert locations.
211    
212          DO j=1-OLy,sNy+OLy          DO j=1-OLy,sNy+OLy
213           DO i=1-OLx,sNx+OLx           DO i=1-OLx,sNx+OLx
214              xA(i,j)        = 0. _d 0
215              yA(i,j)        = 0. _d 0
216              uTrans(i,j)    = 0. _d 0
217              vTrans(i,j)    = 0. _d 0
218              rhok   (i,j)   = 0. _d 0
219              rhoKM1 (i,j)   = 0. _d 0
220              phiSurfX(i,j)  = 0. _d 0
221              phiSurfY(i,j)  = 0. _d 0
222            rTrans (i,j)   = 0. _d 0            rTrans (i,j)   = 0. _d 0
223            fVerT  (i,j,1) = 0. _d 0            fVerT  (i,j,1) = 0. _d 0
224            fVerT  (i,j,2) = 0. _d 0            fVerT  (i,j,2) = 0. _d 0
# Line 231  C--     Set up work arrays that need val Line 226  C--     Set up work arrays that need val
226            fVerS  (i,j,2) = 0. _d 0            fVerS  (i,j,2) = 0. _d 0
227            fVerTr1(i,j,1) = 0. _d 0            fVerTr1(i,j,1) = 0. _d 0
228            fVerTr1(i,j,2) = 0. _d 0            fVerTr1(i,j,2) = 0. _d 0
           rhoKM1 (i,j)   = 0. _d 0  
229           ENDDO           ENDDO
230          ENDDO          ENDDO
231    
Legend:
Removed from v.1.40  
changed lines
  Added in v.1.41
  ViewVC Help
Powered by ViewVC 1.1.22