Changeset 99 for trunk/NEMO/OPA_SRC/eosbn2.F90
- Timestamp:
- 2004-06-25T10:38:42+02:00 (20 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMO/OPA_SRC/eosbn2.F90
r49 r99 38 38 !! * Share module variables 39 39 INTEGER , PUBLIC :: & !: nameos : ocean physical parameters 40 neos ,& !: = 0/1/2 type of eq. of state and Brunt-Vaisala frequ.40 neos = 0, & !: = 0/1/2 type of eq. of state and Brunt-Vaisala frequ. 41 41 neos_init = 0 !: control flag for initialization 42 42 43 43 REAL(wp), PUBLIC :: & !: nameos : ocean physical parameters 44 ralpha ,& !: thermal expension coeff. (linear equation of state)45 rbeta 44 ralpha = 2.0e-4, & !: thermal expension coeff. (linear equation of state) 45 rbeta = 7.7e-4 !: saline expension coeff. (linear equation of state) 46 46 47 47 !! * Substitutions … … 207 207 zs = psal(ji,jj,jk) 208 208 ! ... density and potential volumic mass 209 prd(ji,jj,jk) = ( rbeta * zs - ralpha * zt - 1.) * tmask(ji,jj,jk)209 prd(ji,jj,jk) = ( rbeta * zs - ralpha * zt ) * tmask(ji,jj,jk) 210 210 END DO 211 211 END DO … … 392 392 zs = psal(ji,jj,jk) 393 393 ! ... density and potential volumic mass 394 prd (ji,jj,jk) = ( rbeta * zs - ralpha * zt - 1.) * tmask(ji,jj,jk)394 prd (ji,jj,jk) = ( rbeta * zs - ralpha * zt ) * tmask(ji,jj,jk) 395 395 prhop(ji,jj,jk) = ( rau0 * prd(ji,jj,jk) + rau0 ) * tmask(ji,jj,jk) 396 396 END DO … … 570 570 DO ji = 1, fs_jpim1 ! vector opt. 571 571 #endif 572 prd(ji,jj) = ( rbeta * psal(ji,jj) - ralpha * ptem(ji,jj) - 1.) * tmask(ji,jj,1)572 prd(ji,jj) = ( rbeta * psal(ji,jj) - ralpha * ptem(ji,jj) ) * tmask(ji,jj,1) 573 573 END DO 574 574 ! ! =============== … … 790 790 ! namelist nameos : ocean physical parameters 791 791 792 ! Default values793 neos = 0 ! flag of eos used794 ralpha = 2.e-4_wp ! thermal exp. coef. (linear eos)795 rbeta = 0.029_wp ! haline exp. coef. (linear eos)796 797 792 ! Read Namelist nameos : equation of state 798 793 REWIND( numnam )
Note: See TracChangeset
for help on using the changeset viewer.