Changeset 1739 for trunk/NEMO/OPA_SRC/DYN/dynspg_flt.F90

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

File:

• trunk/NEMO/OPA_SRC/DYN/dynspg_flt.F90 (modified) (2 diffs)

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)