/[MITgcm]/MITgcm_contrib/dgoldberg/streamice/eta_gl_streamice.F
ViewVC logotype

Diff of /MITgcm_contrib/dgoldberg/streamice/eta_gl_streamice.F

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

revision 1.1 by dgoldberg, Tue Feb 19 19:13:30 2013 UTC revision 1.2 by dgoldberg, Wed Aug 27 19:29:12 2014 UTC
# Line 1  Line 1 
1    C $Header$
2    C $Name$
3    
4  #include "CPP_OPTIONS.h"  #include "CPP_OPTIONS.h"
5    
6        _RL FUNCTION ETA_GL_STREAMICE (X, M1, M2, Y0, W)        _RL FUNCTION ETA_GL_STREAMICE (X, M1, M2, Y0, W)
# Line 15  Line 18 
18         _RL X, M1, M2, Y0, W         _RL X, M1, M2, Y0, W
19         _RL TMP1, PI         _RL TMP1, PI
20    
21         IF (X<-1.0*W/2.0) THEN         IF (X.le.-0.5*W) THEN
22           ETA_GL_STREAMICE = Y0 + M1 * X           ETA_GL_STREAMICE = Y0 + M1 * X
23         ELSEIF(X>W/2.0) THEN         ELSEIF(X.gt.0.5*W) THEN
24           ETA_GL_STREAMICE = Y0 + M2 * X           ETA_GL_STREAMICE = Y0 + M2 * X
25         ELSE         ELSE
26           PI = 3.14159265358979323844D0           PI = 3.14159265358979323844D0
27           TMP1 = W/PI * COS(PI*X/W)           TMP1 = W/PI * COS(PI*X/W)
28           ETA_GL_STREAMICE = Y0 +           ETA_GL_STREAMICE = Y0 +
29       &     M1/2 * (X-W/2+TMP1) +       &     M1/2 * (X-0.5*W+TMP1) +
30       &     M2/2 * (X+W/2-TMP1)       &     M2/2 * (X+0.5*W-TMP1)
31         ENDIF         ENDIF
32    
33         RETURN         RETURN
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.2
  ViewVC Help
Powered by ViewVC 1.1.22