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

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

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

revision 1.13 by heimbach, Fri Mar 7 23:46:46 2003 UTC revision 1.14 by jmc, Wed Sep 10 22:21:22 2003 UTC
# Line 133  C     abbreviations Line 133  C     abbreviations
133                 t3 = t2*t1                 t3 = t2*t1
134                                
135                 s1  = salt(i,j,k,bi,bj)                 s1  = salt(i,j,k,bi,bj)
136                  IF ( s1 .GT. 0. _d 0 ) THEN
137                 s3o2 = SQRT(s1*s1*s1)                 s3o2 = SQRT(s1*s1*s1)
138                  ELSE
139                   s1   = 0. _d 0
140                   s3o2 = 0. _d 0
141                  ENDIF
142                                
143                 p1  = locPres(i,j)*SItoBar                 p1  = locPres(i,j)*SItoBar
144                 p2 = p1*p1                 p2 = p1*p1
# Line 222  C     of sea water at p Line 227  C     of sea water at p
227                 t1  = theta(i,j,k,bi,bj)                 t1  = theta(i,j,k,bi,bj)
228                 t2  = t1*t1                 t2  = t1*t1
229                 s1  = salt(i,j,k,bi,bj)                 s1  = salt(i,j,k,bi,bj)
230                  IF ( s1 .GT. 0. _d 0 ) THEN
231                 sp5 = SQRT(s1)                 sp5 = SQRT(s1)
232                  ELSE
233                   s1  = 0. _d 0
234                   sp5 = 0. _d 0
235                  ENDIF
236    
237                 p1   = locPres(i,j)*SItodBar                 p1   = locPres(i,j)*SItodBar
238                 p1t1 = p1*t1                 p1t1 = p1*t1
# Line 230  C     of sea water at p Line 240  C     of sea water at p
240                 dnum_dtheta = eosMDJWFnum(1)                 dnum_dtheta = eosMDJWFnum(1)
241       &              + t1*(2.*eosMDJWFnum(2) + 3.*eosMDJWFnum(3)*t1)       &              + t1*(2.*eosMDJWFnum(2) + 3.*eosMDJWFnum(3)*t1)
242       &              + eosMDJWFnum(5)*s1                                       &              + eosMDJWFnum(5)*s1                                
243       &              + p1t1*(2.*eosMDJWFnum(8) + 2.*eosMDJWFnum(11)*p1)           &              + p1t1*(2.*eosMDJWFnum(8) + 2.*eosMDJWFnum(11)*p1)
244                            
245                 dden_dtheta = eosMDJWFden(1)                 dden_dtheta = eosMDJWFden(1)
246       &              + t1*(2.*eosMDJWFden(2)       &              + t1*(2.*eosMDJWFden(2)
# Line 373  C     abbreviations Line 383  C     abbreviations
383                 t3 = t2*t1                 t3 = t2*t1
384                                
385                 s1  = salt(i,j,k,bi,bj)                 s1  = salt(i,j,k,bi,bj)
386                  IF ( s1 .GT. 0. _d 0 ) THEN
387                 s3o2 = 1.5*SQRT(s1)                 s3o2 = 1.5*SQRT(s1)
388                  ELSE
389                   s1   = 0. _d 0
390                   s3o2 = 0. _d 0
391                  ENDIF
392    
393                 p1  = locPres(i,j)*SItoBar                 p1  = locPres(i,j)*SItoBar
394    
# Line 448  C     of sea water at p Line 463  C     of sea water at p
463                 t1  = theta(i,j,k,bi,bj)                 t1  = theta(i,j,k,bi,bj)
464                 t2  = t1*t1                 t2  = t1*t1
465                 s1  = salt(i,j,k,bi,bj)                 s1  = salt(i,j,k,bi,bj)
466                  IF ( s1 .GT. 0. _d 0 ) THEN
467                 sp5 = SQRT(s1)                 sp5 = SQRT(s1)
468                  ELSE
469                   s1  = 0. _d 0
470                   sp5 = 0. _d 0
471                  ENDIF
472                    
473                 p1   = locPres(i,j)*SItodBar                 p1   = locPres(i,j)*SItodBar
474                 p1t1 = p1*t1                 p1t1 = p1*t1
Legend:
Removed from v.1.13  
changed lines
  Added in v.1.14
  ViewVC Help
Powered by ViewVC 1.1.22