Changeset 1739 for trunk

Timestamp:

2009-11-19T14:24:00+01:00 (14 years ago)

Author:

rblod

Message:

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

Location:

trunk/NEMO

Files:

• OPA_SRC/BDY/bdyvol.F90 (modified) (1 diff)
• OPA_SRC/DOM/phycst.F90 (modified) (1 diff)
• OPA_SRC/DYN/dynspg_flt.F90 (modified) (2 diffs)
• OPA_SRC/DYN/dynspg_ts.F90 (modified) (1 diff)
• OPA_SRC/DYN/sshwzv.F90 (modified) (1 diff)
• OPA_SRC/OBC/obcvol.F90 (modified) (1 diff)
• OPA_SRC/TRA/trasbc.F90 (modified) (2 diffs)
• TOP_SRC/TRP/trcsbc.F90 (modified) (2 diffs)

trunk/NEMO/OPA_SRC/BDY/bdyvol.F90

@@ -85 +85 @@
       ! -----------------------------------------------------------------------
       z_cflxemp = 0.e0
-      zraur = 1.e0 / rauw
+      zraur = 1.e0 / rau0
       z_cflxemp = SUM ( emp(:,:) * bdytmask(:,:) * e1t(:,:) * e2t(:,:) * zraur )
       IF( lk_mpp )   CALL mpp_sum( z_cflxemp )   ! sum over the global domain

trunk/NEMO/OPA_SRC/DOM/phycst.F90

@@ -49 +49 @@
       rt0_ice  = 273.05_wp  ,  &  !: melting point of ice   (Kelvin)
 #endif
-      rau0     = 1020._wp   ,  &  !: volumic mass of reference (kg/m3)
+      rau0     = 1035._wp   ,  &  !: volumic mass of reference (kg/m3)
       rauw     = 1000._wp   ,  &  !: volumic mass of pure water (kg/m3)
       rcp      =    4.e+3_wp,  &  !: ocean specific heat

trunk/NEMO/OPA_SRC/DYN/dynspg_flt.F90

@@ -109 +109 @@
       !!                                   
       INTEGER  ::   ji, jj, jk           ! dummy loop indices
-      REAL(wp) ::   z2dt, z2dtg, zraur   ! temporary scalars
+      REAL(wp) ::   z2dt, z2dtg          ! temporary scalars
       REAL(wp) ::   zgcb, zbtd   !   -          -
       REAL(wp) ::   ztdgu, ztdgv         !   -          -
@@ -135 +135 @@
       IF( neuler == 0 .AND. kt == nit000+1 )   CALL sol_mat( kt )
       z2dtg  = grav * z2dt
-      zraur  = 1. / rauw
 
       ! Evaluate the masked next velocity (effect of the additional force not included)

trunk/NEMO/OPA_SRC/DYN/dynspg_ts.F90

@@ -148 +148 @@
       z1_8 = 0.5 * 0.25                                     ! coefficient for vorticity estimates
       z1_4 = 0.5 * 0.5
-      zraur  = 1. / rauw                                    ! 1 / volumic mass of pure water
+      zraur  = 1. / rau0                                    ! 1 / volumic mass
       !
       zhdiv(:,:) = 0.e0                                     ! barotropic divergence

trunk/NEMO/OPA_SRC/DYN/sshwzv.F90

@@ -139 +139 @@
       IF( neuler == 0 .AND. kt == nit000 )   z2dt =rdt
 
-      zraur = 1. / rauw
+      zraur = 1. / rau0
 
       !                                           !------------------------------!

trunk/NEMO/OPA_SRC/OBC/obcvol.F90

@@ -97 +97 @@
       ! ---------------------------------------------------------------------------
 
-      zCflxemp = SUM ( emp(:,:)*obctmsk(:,:)* e1t(:,:) * e2t(:,:)  / rauw )
+      zCflxemp = SUM ( emp(:,:)*obctmsk(:,:)* e1t(:,:) * e2t(:,:)  / rau0 )
 
       IF( lk_mpp )   CALL mpp_sum( zCflxemp )   ! sum over the global domain

trunk/NEMO/OPA_SRC/TRA/trasbc.F90

@@ -80 +80 @@
       !!         where emp, the surface freshwater budget (evaporation minus
       !!         precipitation minus runoff) given in kg/m2/s is divided
-      !!         by 1000 kg/m3 (density of plain water) to obtain m/s.    
+      !!         by 1035 kg/m3 (density of ocena water) to obtain m/s.    
       !!         Note: even though Fwe does not appear explicitly for 
       !!         temperature in this routine, the heat carried by the water
@@ -113 +113 @@
       ENDIF
 
-      zsrau = 1. / rauw             ! initialization
+      zsrau = 1. / rau0             ! initialization
 #if defined key_zco
       zse3t = 1. / e3t_0(1)

trunk/NEMO/TOP_SRC/TRP/trcsbc.F90

@@ -54 +54 @@
       !!         where emp, the surface freshwater budget (evaporation minus
       !!         precipitation minus runoff) given in kg/m2/s is divided
-      !!         by 1000 kg/m3 (density of plain water) to obtain m/s.
+      !!         by 1035 kg/m3 (density of ocean water) to obtain m/s.
       !!
       !! ** Action  : - Update the 1st level of tra with the trend associated
@@ -79 +79 @@
 
       ! 0. initialization
-      zsrau = 1. / rauw
+      zsrau = 1. / rau0
       IF( .NOT. ln_sco )  zse3t = 1. / fse3t(1,1,1)