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

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

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

revision 1.20 by jmc, Sat Sep 6 17:42:27 2008 UTC revision 1.21 by mlosch, Tue Jul 19 12:53:24 2011 UTC
# Line 69  C     i,k    :: Loop counters Line 69  C     i,k    :: Loop counters
69        DO k = 0,12        DO k = 0,12
70           eosMDJWFden(k) = 0. _d 0           eosMDJWFden(k) = 0. _d 0
71        ENDDO        ENDDO
72          DO k = 1,48
73             teos(k)        = 0. _d 0
74          ENDDO
75    
76        IF ( equationOfState .EQ. 'LINEAR' ) THEN        IF ( equationOfState .EQ. 'LINEAR' ) THEN
77           IF ( tAlpha .EQ. UNSET_RL ) tAlpha = 2.  _d -4           IF ( tAlpha .EQ. UNSET_RL ) tAlpha = 2.  _d -4
# Line 249  C     5. secant bulk modulus K of sea wa Line 252  C     5. secant bulk modulus K of sea wa
252           eosMDJWFden(11) = -3.03175128 _d -16           eosMDJWFden(11) = -3.03175128 _d -16
253           eosMDJWFden(12) = -1.27934137 _d -17           eosMDJWFden(12) = -1.27934137 _d -17
254    
255          ELSEIF ( equationOfState .EQ. 'TEOS10' ) THEN
256    
257           teos(01) =  9.998420897506056 _d +02
258           teos(02) =  2.839940833161907 _d  00
259           teos(03) = -3.147759265588511 _d -02
260           teos(04) =  1.181805545074306 _d -03
261           teos(05) = -6.698001071123802 _d  00
262           teos(06) = -2.986498947203215 _d -02
263           teos(07) =  2.327859407479162 _d -04
264           teos(08) = -3.988822378968490 _d -02
265           teos(09) =  5.095422573880500 _d -04
266           teos(10) = -1.426984671633621 _d -05
267           teos(11) =  1.645039373682922 _d -07
268           teos(12) = -2.233269627352527 _d -02
269           teos(13) = -3.436090079851880 _d -04
270           teos(14) =  3.726050720345733 _d -06
271           teos(15) = -1.806789763745328 _d -04
272           teos(16) =  6.876837219536232 _d -07
273           teos(17) = -3.087032500374211 _d -07
274           teos(18) = -1.988366587925593 _d -08
275           teos(19) = -1.061519070296458 _d -11
276           teos(20) =  1.550932729220080 _d -10
277           teos(21) =  1.000000000000000 _d  00
278           teos(22) =  2.775927747785646 _d -03
279           teos(23) = -2.349607444135925 _d -05
280           teos(24) =  1.119513357486743 _d -06
281           teos(25) =  6.743689325042773 _d -10
282           teos(26) = -7.521448093615448 _d -03
283           teos(27) = -2.764306979894411 _d -05
284           teos(28) =  1.262937315098546 _d -07
285           teos(29) =  9.527875081696435 _d -10
286           teos(30) = -1.811147201949891 _d -11
287           teos(31) = -3.303308871386421 _d -05
288           teos(32) =  3.801564588876298 _d -07
289           teos(33) = -7.672876869259043 _d -09
290           teos(34) = -4.634182341116144 _d -11
291           teos(35) =  2.681097235569143 _d -12
292           teos(36) =  5.419326551148740 _d -06
293           teos(37) = -2.742185394906099 _d -05
294           teos(38) = -3.212746477974189 _d -07
295           teos(39) =  3.191413910561627 _d -09
296           teos(40) = -1.931012931541776 _d -12
297           teos(41) = -1.105097577149576 _d -07
298           teos(42) =  6.211426728363857 _d -10
299           teos(43) = -1.119011592875110 _d -10
300           teos(44) = -1.941660213148725 _d -11
301           teos(45) = -1.864826425365600 _d -14
302           teos(46) =  1.119522344879478 _d -14
303           teos(47) = -1.200507748551599 _d -15
304           teos(48) =  6.057902487546866 _d -17
305    
306        ELSEIF( equationOfState .EQ. 'IDEALG' ) THEN        ELSEIF( equationOfState .EQ. 'IDEALG' ) THEN
307    
308        ELSE        ELSE
# Line 311  C     i,j,k Line 365  C     i,j,k
365        _RL bulkMod(1-Olx:sNx+Olx,1-Oly:sNy+Oly)        _RL bulkMod(1-Olx:sNx+Olx,1-Oly:sNy+Oly)
366    
367        INTEGER ncheck, kcheck        INTEGER ncheck, kcheck
368        PARAMETER ( ncheck = 13 )        PARAMETER ( ncheck = 20 )
369        _RL tLoc(ncheck), ptLoc(ncheck), sLoc(ncheck), pLoc(ncheck)        _RL tLoc(ncheck), ptLoc(ncheck), sLoc(ncheck), pLoc(ncheck)
370        _RL rLoc(ncheck), bLoc(ncheck)        _RL rLoc(ncheck), bLoc(ncheck)
371        _RS mskSave        _RS mskSave
# Line 326  C     i,j,k Line 380  C     i,j,k
380       &      5.               _d 0, 25.               _d 0,       &      5.               _d 0, 25.               _d 0,
381       &      5.               _d 0, 25.               _d 0,       &      5.               _d 0, 25.               _d 0,
382       &      5.               _d 0, 25.               _d 0,       &      5.               _d 0, 25.               _d 0,
383       &      5.               _d 0, 25.               _d 0/,       &      5.               _d 0, 25.               _d 0,
384         &      5.               _d 0, 25.               _d 0,
385         &      5.               _d 0, 25.               _d 0,
386         &      5.               _d 0, 25.               _d 0,
387         &      5.               _d 0/
388       &     ptLoc       &     ptLoc
389       &     /3.               _d 0, 20.               _d 0,       &     /3.               _d 0, 20.               _d 0,
390       &     25.               _d 0, 20.               _d 0,       &     25.               _d 0, 20.               _d 0,
# Line 334  C     i,j,k Line 392  C     i,j,k
392       &      5.               _d 0, 25.               _d 0,       &      5.               _d 0, 25.               _d 0,
393       &      5.               _d 0, 25.               _d 0,       &      5.               _d 0, 25.               _d 0,
394       &      4.03692566635316 _d 0, 22.84661726775120 _d 0,       &      4.03692566635316 _d 0, 22.84661726775120 _d 0,
395       &      3.62720389416752 _d 0, 22.62420229124846 _d 0/       &      3.62720389416752 _d 0, 22.62420229124846 _d 0,
396         &     19.5130           _d 0,  3.5588           _d 0,
397         &      1.8157           _d 0,  1.3434           _d 0,
398         &      1.1583           _d 0,  1.0518           _d 0,
399         &      1.0073           _d 0/
400       &     sLoc       &     sLoc
401       &     /35.5 _d 0, 35. _d 0,       &     /35.5 _d 0, 35. _d 0,
402       &      35.0 _d 0, 20. _d 0,       &      35.0 _d 0, 20. _d 0,
# Line 342  C     i,j,k Line 404  C     i,j,k
404       &       0.  _d 0,  0. _d 0,       &       0.  _d 0,  0. _d 0,
405       &      35.  _d 0, 35. _d 0,       &      35.  _d 0, 35. _d 0,
406       &       0.  _d 0,  0. _d 0,       &       0.  _d 0,  0. _d 0,
407       &      35.  _d 0, 35. _d 0/       &      35.  _d 0, 35. _d 0,
408         &      34.7392 _d 0, 34.4652 _d 0,
409         &      34.7738 _d 0, 34.8435 _d 0,
410         &      34.8637 _d 0, 34.8739 _d 0, 34.8776 _d 0/
411       &     pLoc       &     pLoc
412       &     /300. _d 5,  200. _d 5,       &     /300. _d 5,  200. _d 5,
413       &      200. _d 5,  100. _d 5,       &      200. _d 5,  100. _d 5,
# Line 350  C     i,j,k Line 415  C     i,j,k
415       &        0. _d 0,    0. _d 0,       &        0. _d 0,    0. _d 0,
416       &        0. _d 0,    0. _d 0,       &        0. _d 0,    0. _d 0,
417       &     1000. _d 5, 1000. _d 5,       &     1000. _d 5, 1000. _d 5,
418       &     1000. _d 5, 1000. _d 5/       &     1000. _d 5, 1000. _d 5,
419         &        0. _d 0, 1010. _d 4,
420         &     2025. _d 4, 3045. _d 4,
421         &     4069. _d 4, 5098. _d 4,
422         &     6131. _d 4/
423        DATA rLoc        DATA rLoc
424       &     /1041.83267  _d 0, 1033.213387 _d 0,       &     /1041.83267  _d 0, 1033.213387 _d 0,
425       &      1031.654229 _d 0, 1017.726743 _d 0,       &      1031.654229 _d 0, 1017.726743 _d 0,
# Line 358  C     i,j,k Line 427  C     i,j,k
427       &       999.96675  _d 0,  997.04796  _d 0,       &       999.96675  _d 0,  997.04796  _d 0,
428       &      1027.67547  _d 0, 1023.34306  _d 0,       &      1027.67547  _d 0, 1023.34306  _d 0,
429       &      1044.12802  _d 0, 1037.90204  _d 0,       &      1044.12802  _d 0, 1037.90204  _d 0,
430       &      1069.48914  _d 0, 1062.53817  _d 0/       &      1069.48914  _d 0, 1062.53817  _d 0,
431         &      1024.571477039354932 _d 0, 1031.937207908966911 _d 0,
432         &      1037.002326523349893 _d 0, 1041.668901630106348 _d 0,
433         &      1046.179673478751511 _d 0, 1050.591667083023594 _d 0,
434         &      1054.901835905235885 _d 0/
435       &     bLoc       &     bLoc
436       &     /   -1.00000 _d 0,    -1.00000 _d 0,       &     /   -1.00000 _d 0,    -1.00000 _d 0,
437       &         -1.00000 _d 0,    -1.00000 _d 0,       &         -1.00000 _d 0,    -1.00000 _d 0,
# Line 366  C     i,j,k Line 439  C     i,j,k
439       &      20337.80375 _d 0, 22100.72106 _d 0,       &      20337.80375 _d 0, 22100.72106 _d 0,
440       &      22185.93358 _d 0, 23726.34949 _d 0,       &      22185.93358 _d 0, 23726.34949 _d 0,
441       &      23643.52599 _d 0, 25405.09717 _d 0,       &      23643.52599 _d 0, 25405.09717 _d 0,
442       &      25577.49819 _d 0, 27108.94504 _d 0/       &      25577.49819 _d 0, 27108.94504 _d 0,
443         &         -1.00000 _d 0,    -1.00000 _d 0,
444         &         -1.00000 _d 0,    -1.00000 _d 0,
445         &         -1.00000 _d 0,    -1.00000 _d 0,
446         &         -1.00000 _d 0/
447    
448        bi   = 1        bi   = 1
449        bj   = 1        bj   = 1
# Line 414  C-    + set rC accordingly Line 490  C-    + set rC accordingly
490       O           rhoLoc,       O           rhoLoc,
491       I           k, bi, bj, myThid )       I           k, bi, bj, myThid )
492    
493            IF ( equationOfState.EQ.'MDJWF' ) THEN            IF ( equationOfState.EQ.'MDJWF' .OR.
494         &         equationOfState.EQ.'TEOS10' ) THEN
495              bulkMod(i,j) = -1. _d 0              bulkMod(i,j) = -1. _d 0
496            ELSE            ELSE
497              CALL FIND_BULKMOD(              CALL FIND_BULKMOD(
# Line 439  C-    + set rC accordingly Line 516  C-    + set rC accordingly
516       &          'EOS_CHECK: check values for eosType=UNESCO:'       &          'EOS_CHECK: check values for eosType=UNESCO:'
517             CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,             CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
518       &          SQUEEZE_RIGHT , myThid )       &          SQUEEZE_RIGHT , myThid )
519              ELSEIF ( kcheck .EQ. 14 ) THEN
520               WRITE(msgBuf,'(A)')
521         &          'EOS_CHECK: check values for eosType=TEOS10:'
522               CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,
523         &          SQUEEZE_RIGHT , myThid )
524            ENDIF            ENDIF
525            WRITE(msgBuf,'(2(A,F4.1),A,F5.0,2A,F10.5,1X,F11.5)')            WRITE(msgBuf,'(2(A,F4.1),A,F5.0,2A,F10.5,1X,F11.5)')
526       &         'rho(', sFld(i,j), ' PSU,',       &         'rho(', sFld(i,j), ' PSU,',
Legend:
Removed from v.1.20  
changed lines
  Added in v.1.21
  ViewVC Help
Powered by ViewVC 1.1.22