/[MITgcm]/MITgcm/pkg/sphere/sphere.F
ViewVC logotype

Diff of /MITgcm/pkg/sphere/sphere.F

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

revision 1.1 by heimbach, Thu Feb 17 18:55:10 2005 UTC revision 1.2 by heimbach, Thu May 5 23:47:57 2005 UTC
# Line 313  c     ================================== Line 313  c     ==================================
313        REAL    A, B, E, F, Q, DA, DB, XLAT        REAL    A, B, E, F, Q, DA, DB, XLAT
314    
315  ce      IF (LMAX.GT.LMAXLMT) CALL ABORT('LMAX.GT.LMAXLMT')  ce      IF (LMAX.GT.LMAXLMT) CALL ABORT('LMAX.GT.LMAXLMT')
316        IF (LMAX.GT.LMAXLMT) CALL shcError('LMAX.GT.LMAXLMT',1)  cph      IF (LMAX.GT.LMAXLMT) CALL shcError('LMAX.GT.LMAXLMT',1)
317    
318        DEG2RAD = ACOS(-1.0)/180.0        DEG2RAD = ACOS(-1.0)/180.0
319        XLAT = 0.5*(XLAT1+XLAT2)        XLAT = 0.5*(XLAT1+XLAT2)
# Line 403  C Line 403  C
403        INTEGER JP, I1, I2        INTEGER JP, I1, I2
404    
405  ce      IF(LMAX.GT.LMAXLMT) CALL ABORT('LMAX.GT.LMAXLMT')  ce      IF(LMAX.GT.LMAXLMT) CALL ABORT('LMAX.GT.LMAXLMT')
406        IF(LMAX.GT.LMAXLMT) CALL shcError('LMAX.GT.LMAXLMT',1)  cph      IF(LMAX.GT.LMAXLMT) CALL shcError('LMAX.GT.LMAXLMT',1)
407    
408        IF(XLAT1 .LT. -90.0+1.E-10) THEN        IF(XLAT1 .LT. -90.0+1.E-10) THEN
409          do i = 1,(1+lmax)*(1+lmax)          do i = 1,(1+lmax)*(1+lmax)
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22