Changeset 170 for trunk/NEMO/OPA_SRC/DYN/dynhpg.F90

Timestamp:

2004-11-05T14:52:59+01:00 (19 years ago)

Author:

opalod

Message:

CT : BUGFIX114 : use scale factor of each grid point i.e. fse3w(ji,jj,jk) instead of only one point fse3(1,1,jk)

File:

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

trunk/NEMO/OPA_SRC/DYN/dynhpg.F90

@@ -229 +229 @@
       ! 1. Surface value
       ! ----------------
-      zcoef1 = zcoef0 * fse3w(1,1,1)
       DO jj = 2, jpjm1
          DO ji = fs_2, fs_jpim1   ! vector opt.
+            zcoef1 = zcoef0 * fse3w(ji,jj,1)
             ! hydrostatic pressure gradient
             zhpi(ji,jj,1) = zcoef1 * ( rhd(ji+1,jj,1) - rhd(ji,jj,1) ) / e1u(ji,jj)
@@ -249 +249 @@
       ! -----------------
       DO jk = 2, jpkm1
-         zcoef1 = zcoef0 * fse3w(1,1,jk)
          DO jj = 2, jpjm1
             DO ji = fs_2, fs_jpim1   ! vector opt.
+               zcoef1 = zcoef0 * fse3w(ji,jj,jk)
                ! hydrostatic pressure gradient
                zhpi(ji,jj,jk) = zhpi(ji,jj,jk-1)   &