Changeset 2779 for trunk/NEMOGCM/NEMO/OPA_SRC

Timestamp:

2011-06-07T19:09:43+02:00 (13 years ago)

Author:

acc

Message:

#832. Corrections to before scale factor calculations (fse3u_b and fse3v_b).

Location:

trunk/NEMOGCM/NEMO/OPA_SRC

Files:

• DOM/domvvl.F90 (modified) (1 diff)
• DYN/dynnxt.F90 (modified) (1 diff)

trunk/NEMOGCM/NEMO/OPA_SRC/DOM/domvvl.F90

@@ -141 +141 @@
       ! surface at t-points and inverse surface at (u/v)-points used in surface averaging computations
       ! for ssh and scale factors
-      zs_t  (:,:) =       e1t(:,:) * e2t(:,:)
-      zs_u_1(:,:) = 0.5 / e1u(:,:) * e2u(:,:)
-      zs_v_1(:,:) = 0.5 / e1v(:,:) * e2v(:,:)
+      zs_t  (:,:) =         e1t(:,:) * e2t(:,:)
+      zs_u_1(:,:) = 0.5 / ( e1u(:,:) * e2u(:,:) )
+      zs_v_1(:,:) = 0.5 / ( e1v(:,:) * e2v(:,:) )
 
       DO jj = 1, jpjm1                          ! initialise before and now Sea Surface Height at u-, v-, f-points

trunk/NEMOGCM/NEMO/OPA_SRC/DYN/dynnxt.F90

@@ -250 +250 @@
             ! surface at t-points and inverse surface at (u/v)-points used in surface averaging computations
             zs_t  (:,:) =       e1t(:,:) * e2t(:,:)
-            zs_u_1(:,:) = 0.5 / e1u(:,:) * e2u(:,:)
-            zs_v_1(:,:) = 0.5 / e1v(:,:) * e2v(:,:)
+            zs_u_1(:,:) = 0.5 / ( e1u(:,:) * e2u(:,:) )
+            zs_v_1(:,:) = 0.5 / ( e1v(:,:) * e2v(:,:) )
             !
             IF( ln_dynadv_vec ) THEN