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

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

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

revision 1.2 by dgoldberg, Thu Mar 7 15:10:00 2013 UTC revision 1.3 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_PRIME (X, M1, M2, Y0, W)        _RL FUNCTION ETA_GL_STREAMICE_PRIME (X, M1, M2, Y0, W)
# Line 16  Line 19 
19         _RL TMP1, PI         _RL TMP1, PI
20    
21    
22         IF (X<-1.0*W/2.0) THEN         IF (X.le.-0.5*W) THEN
23           ETA_GL_STREAMICE_PRIME = M1           ETA_GL_STREAMICE_PRIME = M1
24         ELSEIF(X>W/2.0) THEN         ELSEIF(X.gt.0.5*W) THEN
25           ETA_GL_STREAMICE_PRIME = M2           ETA_GL_STREAMICE_PRIME = M2
26         ELSE         ELSE
27           PI = 3.14159265358979323844D0           PI = 3.14159265358979323844D0
28           TMP1 = SIN (PI*X/W)           TMP1 = SIN (PI*X/W)
29           ETA_GL_STREAMICE_PRIME =           ETA_GL_STREAMICE_PRIME =
30       &     M1/2 * (1-TMP1) +       &     0.5*M1 * (1-TMP1) +
31       &     M2/2 * (1+TMP1)       &     0.5*M2 * (1+TMP1)
32         ENDIF         ENDIF
33    
34         RETURN         RETURN
Legend:
Removed from v.1.2  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22