Changeset 1739


Ignore:
Timestamp:
2009-11-19T14:24:00+01:00 (11 years ago)
Author:
rblod
Message:

Use rau0 instead of rauw for sbc, uchange rau0 to 1035, see ticket #606

Location:
trunk/NEMO
Files:
8 edited

Legend:

Unmodified
Added
Removed
  • trunk/NEMO/OPA_SRC/BDY/bdyvol.F90

    r1521 r1739  
    8585      ! ----------------------------------------------------------------------- 
    8686      z_cflxemp = 0.e0 
    87       zraur = 1.e0 / rauw 
     87      zraur = 1.e0 / rau0 
    8888      z_cflxemp = SUM ( emp(:,:) * bdytmask(:,:) * e1t(:,:) * e2t(:,:) * zraur ) 
    8989      IF( lk_mpp )   CALL mpp_sum( z_cflxemp )   ! sum over the global domain 
  • trunk/NEMO/OPA_SRC/DOM/phycst.F90

    r1732 r1739  
    4949      rt0_ice  = 273.05_wp  ,  &  !: melting point of ice   (Kelvin) 
    5050#endif 
    51       rau0     = 1020._wp   ,  &  !: volumic mass of reference (kg/m3) 
     51      rau0     = 1035._wp   ,  &  !: volumic mass of reference (kg/m3) 
    5252      rauw     = 1000._wp   ,  &  !: volumic mass of pure water (kg/m3) 
    5353      rcp      =    4.e+3_wp,  &  !: ocean specific heat 
  • trunk/NEMO/OPA_SRC/DYN/dynspg_flt.F90

    r1683 r1739  
    109109      !!                                    
    110110      INTEGER  ::   ji, jj, jk           ! dummy loop indices 
    111       REAL(wp) ::   z2dt, z2dtg, zraur   ! temporary scalars 
     111      REAL(wp) ::   z2dt, z2dtg          ! temporary scalars 
    112112      REAL(wp) ::   zgcb, zbtd   !   -          - 
    113113      REAL(wp) ::   ztdgu, ztdgv         !   -          - 
     
    135135      IF( neuler == 0 .AND. kt == nit000+1 )   CALL sol_mat( kt ) 
    136136      z2dtg  = grav * z2dt 
    137       zraur  = 1. / rauw 
    138137 
    139138      ! Evaluate the masked next velocity (effect of the additional force not included) 
  • trunk/NEMO/OPA_SRC/DYN/dynspg_ts.F90

    r1708 r1739  
    148148      z1_8 = 0.5 * 0.25                                     ! coefficient for vorticity estimates 
    149149      z1_4 = 0.5 * 0.5 
    150       zraur  = 1. / rauw                                    ! 1 / volumic mass of pure water 
     150      zraur  = 1. / rau0                                    ! 1 / volumic mass 
    151151      ! 
    152152      zhdiv(:,:) = 0.e0                                     ! barotropic divergence 
  • trunk/NEMO/OPA_SRC/DYN/sshwzv.F90

    r1694 r1739  
    139139      IF( neuler == 0 .AND. kt == nit000 )   z2dt =rdt 
    140140 
    141       zraur = 1. / rauw 
     141      zraur = 1. / rau0 
    142142 
    143143      !                                           !------------------------------! 
  • trunk/NEMO/OPA_SRC/OBC/obcvol.F90

    r1528 r1739  
    9797      ! --------------------------------------------------------------------------- 
    9898 
    99       zCflxemp = SUM ( emp(:,:)*obctmsk(:,:)* e1t(:,:) * e2t(:,:)  / rauw ) 
     99      zCflxemp = SUM ( emp(:,:)*obctmsk(:,:)* e1t(:,:) * e2t(:,:)  / rau0 ) 
    100100 
    101101      IF( lk_mpp )   CALL mpp_sum( zCflxemp )   ! sum over the global domain 
  • trunk/NEMO/OPA_SRC/TRA/trasbc.F90

    r1146 r1739  
    8080      !!         where emp, the surface freshwater budget (evaporation minus 
    8181      !!         precipitation minus runoff) given in kg/m2/s is divided 
    82       !!         by 1000 kg/m3 (density of plain water) to obtain m/s.     
     82      !!         by 1035 kg/m3 (density of ocena water) to obtain m/s.     
    8383      !!         Note: even though Fwe does not appear explicitly for  
    8484      !!         temperature in this routine, the heat carried by the water 
     
    113113      ENDIF 
    114114 
    115       zsrau = 1. / rauw             ! initialization 
     115      zsrau = 1. / rau0             ! initialization 
    116116#if defined key_zco 
    117117      zse3t = 1. / e3t_0(1) 
  • trunk/NEMO/TOP_SRC/TRP/trcsbc.F90

    r1562 r1739  
    5454      !!         where emp, the surface freshwater budget (evaporation minus 
    5555      !!         precipitation minus runoff) given in kg/m2/s is divided 
    56       !!         by 1000 kg/m3 (density of plain water) to obtain m/s. 
     56      !!         by 1035 kg/m3 (density of ocean water) to obtain m/s. 
    5757      !! 
    5858      !! ** Action  : - Update the 1st level of tra with the trend associated 
     
    7979 
    8080      ! 0. initialization 
    81       zsrau = 1. / rauw 
     81      zsrau = 1. / rau0 
    8282      IF( .NOT. ln_sco )  zse3t = 1. / fse3t(1,1,1) 
    8383 
Note: See TracChangeset for help on using the changeset viewer.