/[MITgcm]/MITgcm/pkg/generic_advdiff/gad_dst3fl_adv_x.F
ViewVC logotype

Diff of /MITgcm/pkg/generic_advdiff/gad_dst3fl_adv_x.F

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

revision 1.15 by mlosch, Thu Oct 13 15:10:32 2011 UTC revision 1.16 by jmc, Fri Apr 4 20:29:08 2014 UTC
# Line 42  C- jmc: an alternative would be to compu Line 42  C- jmc: an alternative would be to compu
42  C       (if Rj*Rjm < 0 => psiP*Rj = 0 , elsef Rj > 0 ... , else  ... )  C       (if Rj*Rjm < 0 => psiP*Rj = 0 , elsef Rj > 0 ... , else  ... )
43  C       with no need to compute thetaM (might be easier to differentiate)  C       with no need to compute thetaM (might be easier to differentiate)
44    
45        DO j=1-Oly,sNy+Oly        DO j=1-OLy,sNy+OLy
46         uT(1-Olx,j)=0. _d 0         uT(1-OLx,j)=0. _d 0
47         uT(2-Olx,j)=0. _d 0         uT(2-OLx,j)=0. _d 0
48         uT(sNx+Olx,j)=0. _d 0         uT(sNx+OLx,j)=0. _d 0
49        ENDDO        ENDDO
50        DO j=1-Oly,sNy+Oly        DO j=1-OLy,sNy+OLy
51         DO i=1-Olx+2,sNx+Olx-1         DO i=1-OLx+2,sNx+OLx-1
52  #if (defined ALLOW_AUTODIFF_TAMC && defined TARGET_NEC_SX)  #if (defined ALLOW_AUTODIFF && defined TARGET_NEC_SX)
53  C     These lines make TAF create vectorizable code  C     These lines make TAF create vectorizable code
54          thetaP = 0. _d 0          thetaP = 0. _d 0
55          thetaM = 0. _d 0          thetaM = 0. _d 0
Legend:
Removed from v.1.15  
changed lines
  Added in v.1.16
  ViewVC Help
Powered by ViewVC 1.1.22