Changeset 1497 for trunk/NEMO
- Timestamp:
- 2009-07-20T16:34:02+02:00 (15 years ago)
- Location:
- trunk/NEMO/OFF_SRC
- Files:
-
- 2 edited
-
eosbn2.F90 (modified) (27 diffs)
-
opa.F90 (modified) (1 diff)
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMO/OFF_SRC/eosbn2.F90
r1326 r1497 5 5 !! - Brunt-Vaisala frequency 6 6 !!============================================================================== 7 !! History : ! 89-03 (O. Marti) Original code 8 !! 6.0 ! 94-07 (G. Madec, M. Imbard) add bn2 9 !! 6.0 ! 94-08 (G. Madec) Add Jackett & McDougall eos 10 !! 7.0 ! 96-01 (G. Madec) statement function for e3 11 !! 8.1 ! 97-07 (G. Madec) introduction of neos, OPA8.1 12 !! 8.1 ! 97-07 (G. Madec) density instead of volumic mass 13 !! ! 99-02 (G. Madec, N. Grima) semi-implicit pressure gradient 14 !! ! 01-09 (M. Ben Jelloul) bugfix onlinear eos 15 !! 8.5 ! 02-10 (G. Madec) add eos_init 16 !! 8.5 ! 02-11 (G. Madec, A. Bozec) partial step, eos_insitu_2d 17 !! 9.0 ! 03-08 (G. Madec) F90, free form 18 !! 9.0 ! 06-08 (G. Madec) add tfreez function 19 !!---------------------------------------------------------------------- 7 20 8 21 !!---------------------------------------------------------------------- … … 13 26 !! eos_insitu_2d : Compute the in situ density for 2d fields 14 27 !! eos_bn2 : Compute the Brunt-Vaisala frequency 28 !! tfreez : Compute the surface freezing temperature 15 29 !! eos_init : set eos parameters (namelist) 16 30 !!---------------------------------------------------------------------- 17 !! * Modules used18 31 USE dom_oce ! ocean space and time domain 19 32 USE phycst ! physical constants … … 33 46 END INTERFACE 34 47 35 !! * Routine accessibility 36 PUBLIC eos ! called by step.F90, inidtr.F90, tranpc.F90 and intgrd.F90 37 PUBLIC bn2 ! called by step.F90 38 39 !! * Share module variables 40 INTEGER , PUBLIC :: & !: nameos : ocean physical parameters 41 neos = 0, & !: = 0/1/2 type of eq. of state and Brunt-Vaisala frequ. 42 neos_init = 0 !: control flag for initialization 43 44 REAL(wp), PUBLIC :: & !: nameos : ocean physical parameters 45 ralpha = 2.0e-4, & !: thermal expension coeff. (linear equation of state) 46 rbeta = 7.7e-4 !: saline expension coeff. (linear equation of state) 48 PUBLIC eos ! called by step, istate, tranpc and zpsgrd modules 49 PUBLIC eos_init ! called by istate module 50 PUBLIC bn2 ! called by step module 51 PUBLIC tfreez ! called by sbcice_... modules 52 53 !!* Namelist (nameos) 54 INTEGER , PUBLIC :: neos = 0 !: = 0/1/2 type of eq. of state and Brunt-Vaisala frequ. 55 REAL(wp), PUBLIC :: ralpha = 2.0e-4 !: thermal expension coeff. (linear equation of state) 56 REAL(wp), PUBLIC :: rbeta = 7.7e-4 !: saline expension coeff. (linear equation of state) 47 57 48 58 !! * Substitutions … … 50 60 # include "vectopt_loop_substitute.h90" 51 61 !!---------------------------------------------------------------------- 52 !! OPA 9.0 , LOCEAN-IPSL (200 5)53 !! $Id$ 54 !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt62 !! OPA 9.0 , LOCEAN-IPSL (2006) 63 !! $Id$ 64 !! Software governed by the CeCILL licence (modipsl/doc/NEMO_CeCILL.txt) 55 65 !!---------------------------------------------------------------------- 56 66 57 67 CONTAINS 58 68 59 SUBROUTINE eos_insitu 69 SUBROUTINE eos_insitu( ptem, psal, prd ) 60 70 !!---------------------------------------------------------------------- 61 71 !! *** ROUTINE eos_insitu *** … … 91 101 !! ** Action : compute prd , the in situ density (no units) 92 102 !! 93 !! References : 94 !! Jackett, D.R., and T.J. McDougall. J. Atmos. Ocean. Tech., 1994 95 !! 96 !! History : 97 !! ! 89-03 (o. Marti) Original code 98 !! ! 94-08 (G. Madec) 99 !! ! 96-01 (G. Madec) statement function for e3 100 !! ! 97-07 (G. Madec) introduction of neos, OPA8.1 101 !! ! 97-07 (G. Madec) density instead of volumic mass 102 !! ! 99-02 (G. Madec, N. Grima) semi-implicit pressure gradient 103 !! ! 01-09 (M. Ben Jelloul) bugfix 104 !!---------------------------------------------------------------------- 105 !! * Arguments 106 REAL(wp), DIMENSION(jpi,jpj,jpk), INTENT( in ) :: & 107 ptem, & ! potential temperature 108 psal ! salinity 109 REAL(wp), DIMENSION(jpi,jpj,jpk), INTENT( out ) :: & 110 prd ! potential density (surface referenced) 111 112 !! * Local declarations 103 !! References : Jackett and McDougall, J. Atmos. Ocean. Tech., 1994 104 !!---------------------------------------------------------------------- 105 REAL(wp), DIMENSION(jpi,jpj,jpk), INTENT(in ) :: ptem ! potential temperature [Celcius] 106 REAL(wp), DIMENSION(jpi,jpj,jpk), INTENT(in ) :: psal ! salinity [psu] 107 REAL(wp), DIMENSION(jpi,jpj,jpk), INTENT( out) :: prd ! in situ density 108 !! 113 109 INTEGER :: ji, jj, jk ! dummy loop indices 114 110 REAL(wp) :: & … … 118 114 zd , zc , zaw, za , & ! " " 119 115 zb1, za1, zkw, zk0 ! " " 120 REAL(wp), DIMENSION(jpi,jpj,jpk) :: & 121 zws ! temporary workspace 122 !!---------------------------------------------------------------------- 123 124 125 ! initialization (in not already done) 126 IF( neos_init == 0 ) CALL eos_init 127 116 REAL(wp), DIMENSION(jpi,jpj,jpk) :: zws ! temporary workspace 117 !!---------------------------------------------------------------------- 128 118 129 119 SELECT CASE ( neos ) 130 120 ! 131 121 CASE ( 0 ) ! Jackett and McDougall (1994) formulation 132 122 ! 133 123 !CDIR NOVERRCHK 134 124 zws(:,:,:) = SQRT( ABS( psal(:,:,:) ) ) 135 136 125 ! ! =============== 137 126 DO jk = 1, jpkm1 ! Horizontal slab … … 180 169 END DO ! End of slab 181 170 ! ! =============== 182 183 171 ! 184 172 CASE ( 1 ) ! Linear formulation function of temperature only 185 173 ! 186 174 ! ! =============== 187 175 DO jk = 1, jpkm1 ! Horizontal slab … … 198 186 END DO ! End of slab 199 187 ! ! =============== 200 201 188 ! 202 189 CASE ( 2 ) ! Linear formulation function of temperature and salinity 203 190 ! 204 191 ! ! =============== 205 192 DO jk = 1, jpkm1 ! Horizontal slab … … 216 203 END DO ! End of slab 217 204 ! ! =============== 218 205 ! 219 206 CASE DEFAULT 220 207 ! 221 208 WRITE(ctmp1,*) ' bad flag value for neos = ', neos 222 209 CALL ctl_stop( ctmp1 ) 223 210 ! 224 211 END SELECT 225 226 IF(ln_ctl) THEN 227 CALL prt_ctl(tab3d_1=prd, clinfo1=' eos : ', ovlap=1, kdim=jpk) 228 ENDIF 229 212 ! 213 IF(ln_ctl) CALL prt_ctl(tab3d_1=prd, clinfo1=' eos : ', ovlap=1, kdim=jpk) 214 ! 230 215 END SUBROUTINE eos_insitu 231 216 232 217 233 SUBROUTINE eos_insitu_pot ( ptem, psal, prd, prhop)218 SUBROUTINE eos_insitu_pot( ptem, psal, prd, prhop ) 234 219 !!---------------------------------------------------------------------- 235 220 !! *** ROUTINE eos_insitu_pot *** … … 274 259 !! - prhop, the potential volumic mass (Kg/m3) 275 260 !! 276 !! References : 277 !! Jackett, D.R., and T.J. McDougall. J. Atmos. Ocean. Tech., 1994 278 !! Brown, J. A. and K. A. Campana. Mon. Weather Rev., 1978 279 !! 280 !! History : 281 !! 4.0 ! 89-03 (O. Marti) 282 !! ! 94-08 (G. Madec) 283 !! ! 96-01 (G. Madec) statement function for e3 284 !! ! 97-07 (G. Madec) introduction of neos, OPA8.1 285 !! ! 97-07 (G. Madec) density instead of volumic mass 286 !! ! 99-02 (G. Madec, N. Grima) semi-implicit pressure gradient 287 !! ! 01-09 (M. Ben Jelloul) bugfix 288 !! 9.0 ! 03-08 (G. Madec) F90, free form 289 !!---------------------------------------------------------------------- 290 !! * Arguments 291 REAL(wp), DIMENSION(jpi,jpj,jpk), INTENT( in ) :: & 292 ptem, & ! potential temperature 293 psal ! salinity 294 REAL(wp), DIMENSION(jpi,jpj,jpk), INTENT( out ) :: & 295 prd, & ! potential density (surface referenced) 296 prhop ! potential density (surface referenced) 297 298 !! * Local declarations 261 !! References : Jackett and McDougall, J. Atmos. Ocean. Tech., 1994 262 !! Brown and Campana, Mon. Weather Rev., 1978 263 !!---------------------------------------------------------------------- 264 REAL(wp), DIMENSION(jpi,jpj,jpk), INTENT(in ) :: ptem ! potential temperature [Celcius] 265 REAL(wp), DIMENSION(jpi,jpj,jpk), INTENT(in ) :: psal ! salinity [psu] 266 REAL(wp), DIMENSION(jpi,jpj,jpk), INTENT( out) :: prd ! in situ density 267 REAL(wp), DIMENSION(jpi,jpj,jpk), INTENT( out) :: prhop ! potential density (surface referenced) 268 299 269 INTEGER :: ji, jj, jk ! dummy loop indices 300 270 REAL(wp) :: & ! temporary scalars 301 271 zt, zs, zh, zsr, zr1, zr2, zr3, zr4, zrhop, ze, zbw, & 302 272 zb, zd, zc, zaw, za, zb1, za1, zkw, zk0 303 REAL(wp), DIMENSION(jpi,jpj,jpk) :: zws 304 !!---------------------------------------------------------------------- 305 306 ! initialization (in not already done) 307 IF( neos_init == 0 ) CALL eos_init 308 273 REAL(wp), DIMENSION(jpi,jpj,jpk) :: zws 274 !!---------------------------------------------------------------------- 309 275 310 276 SELECT CASE ( neos ) 311 277 ! 312 278 CASE ( 0 ) ! Jackett and McDougall (1994) formulation 313 279 ! 314 280 !CDIR NOVERRCHK 315 281 zws(:,:,:) = SQRT( ABS( psal(:,:,:) ) ) 316 317 282 ! ! =============== 318 283 DO jk = 1, jpkm1 ! Horizontal slab … … 325 290 zh = fsdept(ji,jj,jk) 326 291 ! square root salinity 327 !!Edmee zsr= SQRT( ABS( zs ) )328 292 zsr= zws(ji,jj,jk) 329 293 ! compute volumic mass pure water at atm pressure … … 365 329 END DO ! End of slab 366 330 ! ! =============== 367 368 331 ! 369 332 CASE ( 1 ) ! Linear formulation function of temperature only 370 333 ! 371 334 ! ! =============== 372 335 DO jk = 1, jpkm1 ! Horizontal slab … … 384 347 END DO ! End of slab 385 348 ! ! =============== 386 387 349 ! 388 350 CASE ( 2 ) ! Linear formulation function of temperature and salinity 389 351 ! 390 352 ! ! =============== 391 353 DO jk = 1, jpkm1 ! Horizontal slab … … 403 365 END DO ! End of slab 404 366 ! ! =============== 405 367 ! 406 368 CASE DEFAULT 407 369 ! 408 370 WRITE(ctmp1,*) ' bad flag value for neos = ', neos 409 371 CALL ctl_stop( ctmp1 ) 410 372 ! 411 373 END SELECT 412 413 IF(ln_ctl) THEN 414 CALL prt_ctl(tab3d_1=prd, clinfo1=' eos-p: ', tab3d_2=prhop, clinfo2=' pot : ', ovlap=1, kdim=jpk) 415 ENDIF 416 374 ! 375 IF(ln_ctl) CALL prt_ctl( tab3d_1=prd, clinfo1=' eos-p: ', tab3d_2=prhop, clinfo2=' pot : ', ovlap=1, kdim=jpk ) 376 ! 417 377 END SUBROUTINE eos_insitu_pot 418 378 419 SUBROUTINE eos_insitu_2d ( ptem, psal, pdep, prd ) 379 380 SUBROUTINE eos_insitu_2d( ptem, psal, pdep, prd ) 420 381 !!---------------------------------------------------------------------- 421 382 !! *** ROUTINE eos_insitu_2d *** … … 451 412 !! ** Action : - prd , the in situ density (no units) 452 413 !! 453 !! References : 454 !! Jackett, D.R., and T.J. McDougall. J. Atmos. Ocean. Tech., 1994 455 !! 456 !! History : 457 !! 8.5 ! 02-11 (G. Madec, A. Bozec) partial step 458 !!---------------------------------------------------------------------- 459 !! * Arguments 460 REAL(wp), DIMENSION(jpi,jpj), INTENT( in ) :: & 461 ptem, & ! potential temperature 462 psal, & ! salinity 463 pdep ! depth 464 REAL(wp), DIMENSION(jpi,jpj), INTENT( out ) :: & 465 prd ! potential density (surface referenced) 466 467 !! * Local declarations 468 INTEGER :: ji, jj ! dummy loop indices 414 !! References : Jackett and McDougall, J. Atmos. Ocean. Tech., 1994 415 !!---------------------------------------------------------------------- 416 REAL(wp), DIMENSION(jpi,jpj), INTENT(in ) :: ptem ! potential temperature [Celcius] 417 REAL(wp), DIMENSION(jpi,jpj), INTENT(in ) :: psal ! salinity [psu] 418 REAL(wp), DIMENSION(jpi,jpj), INTENT(in ) :: pdep ! depth [m] 419 REAL(wp), DIMENSION(jpi,jpj), INTENT( out) :: prd ! in situ density 420 !! 421 INTEGER :: ji, jj ! dummy loop indices 469 422 REAL(wp) :: & ! temporary scalars 470 423 zt, zs, zh, zsr, zr1, zr2, zr3, zr4, zrhop, ze, zbw, & 471 424 zb, zd, zc, zaw, za, zb1, za1, zkw, zk0, & 472 425 zmask 473 REAL(wp), DIMENSION(jpi,jpj) :: zws 474 !!---------------------------------------------------------------------- 475 476 477 ! initialization (in not already done) 478 IF( neos_init == 0 ) CALL eos_init 426 REAL(wp), DIMENSION(jpi,jpj) :: zws 427 !!---------------------------------------------------------------------- 479 428 480 429 prd(:,:) = 0.e0 481 430 482 431 SELECT CASE ( neos ) 483 432 ! 484 433 CASE ( 0 ) ! Jackett and McDougall (1994) formulation 485 434 ! 486 435 !CDIR NOVERRCHK 487 436 DO jj = 1, jpjm1 … … 491 440 END DO 492 441 END DO 493 494 442 ! ! =============== 495 443 DO jj = 1, jpjm1 ! Horizontal slab … … 533 481 ! masked in situ density anomaly 534 482 prd(ji,jj) = ( zrhop / ( 1.0 - zh / ( zk0 - zh * ( za - zh * zb ) ) ) - rau0 ) & 535 / rau0 * zmask 536 END DO 537 ! ! =============== 538 END DO ! End of slab 539 ! ! =============== 540 541 483 & / rau0 * zmask 484 END DO 485 ! ! =============== 486 END DO ! End of slab 487 ! ! =============== 488 ! 542 489 CASE ( 1 ) ! Linear formulation function of temperature only 543 490 ! 544 491 ! ! =============== 545 492 DO jj = 1, jpjm1 ! Horizontal slab … … 551 498 END DO ! End of slab 552 499 ! ! =============== 553 554 500 ! 555 501 CASE ( 2 ) ! Linear formulation function of temperature and salinity 556 502 ! 557 503 ! ! =============== 558 504 DO jj = 1, jpjm1 ! Horizontal slab … … 564 510 END DO ! End of slab 565 511 ! ! =============== 566 512 ! 567 513 CASE DEFAULT 568 514 ! 569 515 WRITE(ctmp1,*) ' bad flag value for neos = ', neos 570 516 CALL ctl_stop( ctmp1 ) 571 517 ! 572 518 END SELECT 573 519 574 IF(ln_ctl) CALL prt_ctl( tab2d_1=prd, clinfo1=' eos2d: ')575 520 IF(ln_ctl) CALL prt_ctl( tab2d_1=prd, clinfo1=' eos2d: ' ) 521 ! 576 522 END SUBROUTINE eos_insitu_2d 577 523 … … 605 551 !! ** Action : - pn2 : the brunt-vaisala frequency 606 552 !! 607 !! References : 608 !! McDougall, T. J., J. Phys. Oceanogr., 17, 1950-1964, 1987. 609 !! 610 !! History : 611 !! 6.0 ! 94-07 (G. Madec, M. Imbard) Original code 612 !! 8.0 ! 97-07 (G. Madec) introduction of statement functions 613 !! 8.5 ! 02-07 (G. Madec) Free form, F90 614 !! 8.5 ! 02-08 (G. Madec) introduction of arguments 615 !!---------------------------------------------------------------------- 616 !! * Arguments 617 REAL(wp), DIMENSION(jpi,jpj,jpk), INTENT( in ) :: & 618 ptem, & ! potential temperature 619 psal ! salinity 620 REAL(wp), DIMENSION(jpi,jpj,jpk), INTENT( out ) :: & 621 pn2 ! Brunt-Vaisala frequency 622 623 !! * Local declarations 553 !! References : McDougall, J. Phys. Oceanogr., 17, 1950-1964, 1987. 554 !!---------------------------------------------------------------------- 555 REAL(wp), DIMENSION(jpi,jpj,jpk), INTENT(in ) :: ptem ! potential temperature [Celcius] 556 REAL(wp), DIMENSION(jpi,jpj,jpk), INTENT(in ) :: psal ! salinity [psu] 557 REAL(wp), DIMENSION(jpi,jpj,jpk), INTENT( out) :: pn2 ! Brunt-Vaisala frequency [s-1] 558 624 559 INTEGER :: ji, jj, jk ! dummy loop indices 625 REAL(wp) :: & 626 zgde3w, zt, zs, zh, & ! temporary scalars 627 zalbet, zbeta ! " " 560 REAL(wp) :: zgde3w, zt, zs, zh, & ! temporary scalars 561 & zalbet, zbeta ! " " 628 562 #if defined key_zdfddm 629 563 REAL(wp) :: zds ! temporary scalars 630 564 #endif 631 565 !!---------------------------------------------------------------------- 632 !! OPA8.5, LODYC-IPSL (2002)633 !!----------------------------------------------------------------------634 566 635 567 ! pn2 : first and last levels … … 642 574 643 575 SELECT CASE ( neos ) 644 576 ! 645 577 CASE ( 0 ) ! Jackett and McDougall (1994) formulation 646 578 ! 647 579 ! ! =============== 648 580 DO jk = 2, jpkm1 ! Horizontal slab … … 694 626 END DO ! End of slab 695 627 ! ! =============== 696 697 628 ! 698 629 CASE ( 1 ) ! Linear formulation function of temperature only 699 630 ! 700 631 ! ! =============== 701 632 DO jk = 2, jpkm1 ! Horizontal slab … … 710 641 END DO ! End of slab 711 642 ! ! =============== 712 713 643 ! 714 644 CASE ( 2 ) ! Linear formulation function of temperature and salinity 715 645 ! 716 646 ! ! =============== 717 647 DO jk = 2, jpkm1 ! Horizontal slab … … 738 668 END DO ! End of slab 739 669 ! ! =============== 740 670 ! 741 671 CASE DEFAULT 742 672 ! 743 673 WRITE(ctmp1,*) ' bad flag value for neos = ', neos 744 674 CALL ctl_stop( ctmp1 ) 745 675 ! 746 676 END SELECT 747 677 748 IF(ln_ctl) THEN 749 CALL prt_ctl(tab3d_1=pn2, clinfo1=' bn2 : ', ovlap=1, kdim=jpk) 678 IF(ln_ctl) CALL prt_ctl(tab3d_1=pn2, clinfo1=' bn2 : ', ovlap=1, kdim=jpk) 750 679 #if defined key_zdfddm 751 CALL prt_ctl(tab3d_1=rrau, clinfo1=' rrau : ', ovlap=1, kdim=jpk)680 IF(ln_ctl) CALL prt_ctl(tab3d_1=rrau, clinfo1=' rrau : ', ovlap=1, kdim=jpk) 752 681 #endif 753 ENDIF 754 682 ! 755 683 END SUBROUTINE eos_bn2 756 684 757 685 686 FUNCTION tfreez( psal ) RESULT( ptf ) 687 !!---------------------------------------------------------------------- 688 !! *** ROUTINE eos_init *** 689 !! 690 !! ** Purpose : Compute the sea surface freezing temperature [Celcius] 691 !! 692 !! ** Method : UNESCO freezing point at the surface (pressure = 0???) 693 !! freezing point [Celcius]=(-.0575+1.710523e-3*sqrt(abs(s))-2.154996e-4*s)*s-7.53e-4*p 694 !! checkvalue: tf= -2.588567 Celsius for s=40.0psu, p=500. decibars 695 !! 696 !! Reference : UNESCO tech. papers in the marine science no. 28. 1978 697 !!---------------------------------------------------------------------- 698 REAL(wp), DIMENSION(jpi,jpj), INTENT(in ) :: psal ! salinity [psu] 699 REAL(wp), DIMENSION(jpi,jpj) :: ptf ! freezing temperature [Celcius] 700 !!---------------------------------------------------------------------- 701 ptf(:,:) = ( - 0.0575 + 1.710523e-3 * SQRT( psal(:,:) ) & 702 & - 2.154996e-4 * psal(:,:) ) * psal(:,:) 703 END FUNCTION tfreez 704 705 758 706 SUBROUTINE eos_init 759 707 !!---------------------------------------------------------------------- … … 762 710 !! ** Purpose : initializations for the equation of state 763 711 !! 764 !! ** Method : Read the namelist nameos 765 !! 766 !! ** Action : blahblah.... 767 !! 768 !! History : 769 !! 8.5 ! 02-10 (G. Madec) Original code 712 !! ** Method : Read the namelist nameos and control the parameters 770 713 !!---------------------------------------------------------------------- 771 714 NAMELIST/nameos/ neos, ralpha, rbeta 772 !!---------------------------------------------------------------------- 773 !! OPA 8.5, LODYC-IPSL (2002) 774 !!---------------------------------------------------------------------- 775 776 ! set the initialization flag to 1 777 neos_init = 1 ! indicate that the initialization has been done 778 779 ! namelist nameos : ocean physical parameters 780 781 ! Read Namelist nameos : equation of state 782 REWIND( numnam ) 715 716 REWIND( numnam ) ! Read Namelist nameos : equation of state 783 717 READ ( numnam, nameos ) 784 718 … … 789 723 WRITE(numout,*) '~~~~~~~~' 790 724 WRITE(numout,*) ' Namelist nameos : set eos parameters' 791 WRITE(numout,*)792 725 WRITE(numout,*) ' flag for eq. of state and N^2 neos = ', neos 793 726 WRITE(numout,*) ' thermal exp. coef. (linear) ralpha = ', ralpha 794 727 WRITE(numout,*) ' saline exp. coef. (linear) rbeta = ', rbeta 795 WRITE(numout,*)796 728 ENDIF 797 729 … … 799 731 800 732 CASE ( 0 ) ! Jackett and McDougall (1994) formulation 801 733 IF(lwp) WRITE(numout,*) 802 734 IF(lwp) WRITE(numout,*) ' use of Jackett & McDougall (1994) equation of state and' 803 735 IF(lwp) WRITE(numout,*) ' McDougall (1987) Brunt-Vaisala frequency' 804 736 ! 805 737 CASE ( 1 ) ! Linear formulation function of temperature only 806 738 IF(lwp) WRITE(numout,*) 807 739 IF(lwp) WRITE(numout,*) ' use of linear eos rho(T) = rau0 * ( 1.0285 - ralpha * T )' 808 740 IF( lk_zdfddm ) CALL ctl_stop( ' double diffusive mixing parameterization requires', & 809 741 & ' that T and S are used as state variables' ) 810 742 ! 811 743 CASE ( 2 ) ! Linear formulation function of temperature and salinity 812 744 IF(lwp) WRITE(numout,*) 813 745 IF(lwp) WRITE(numout,*) ' use of linear eos rho(T,S) = rau0 * ( rbeta * S - ralpha * T )' 814 815 CASE DEFAULT 816 746 ! 747 CASE DEFAULT ! E R R O R in neos 817 748 WRITE(ctmp1,*) ' bad flag value for neos = ', neos 818 749 CALL ctl_stop( ctmp1 ) 819 820 750 END SELECT 821 751 -
trunk/NEMO/OFF_SRC/opa.F90
r1457 r1497 218 218 CALL phy_cst ! Physical constants 219 219 220 CALL eos_init ! Equation of state 221 220 222 CALL dom_cfg ! Domain configuration 221 223
Note: See TracChangeset
for help on using the changeset viewer.
