Changeset 1930 for trunk

Timestamp:

2010-06-15T09:55:52+02:00 (14 years ago)

Author:

smasson

Message:

supress rauw and change rau0 to 1035 in OFF_SRC, see ticket:606

Location:

trunk/NEMO

Files:

• OFF_SRC/DOM/phycst.F90 (modified) (2 diffs)
• OPA_SRC/DOM/phycst.F90 (modified) (2 diffs)

trunk/NEMO/OFF_SRC/DOM/phycst.F90

@@ -37 +37 @@
       rtt      = 273.16_wp  ,  &  !: triple point of temperature (Kelvin)
       rt0      = 273.15_wp  ,  &  !: freezing point of water (Kelvin)
-      rau0     = 1020._wp   ,  &  !: volumic mass of reference (kg/m3)
-      rauw     = 1000._wp   ,  &  !: density of pure water (kg/m3)
+      rau0     = 1035._wp   ,  &  !: volumic mass of reference (kg/m3)
       rcp      =    4.e+3_wp,  &  !: ocean specific heat
       ro0cpr                      !: = 1. / ( rau0 * rcp )
@@ -127 +126 @@
       ro0cpr = 1. / ( rau0 * rcp )
       IF(lwp) WRITE(numout,*)
-      IF(lwp) WRITE(numout,*) '          volumic mass of pure water         rauw   = ', rauw, ' kg/m^3'
       IF(lwp) WRITE(numout,*) '          volumic mass of reference          rau0   = ', rau0, ' kg/m^3'
       IF(lwp) WRITE(numout,*) '          ocean specific heat                rcp    = ', rcp

trunk/NEMO/OPA_SRC/DOM/phycst.F90

@@ -49 +49 @@
 #endif
       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
       ro0cpr                      !: = 1. / ( rau0 * rcp )
@@ -154 +153 @@
       ro0cpr = 1. / ( rau0 * rcp )
       IF(lwp) WRITE(numout,*)
-      IF(lwp) WRITE(numout,*) '          volumic mass of pure water         rauw   = ', rauw, ' kg/m^3'
       IF(lwp) WRITE(numout,*) '          volumic mass of reference          rau0   = ', rau0, ' kg/m^3'
       IF(lwp) WRITE(numout,*) '          ocean specific heat                rcp    = ', rcp