Changeset 4366 for branches/2013/dev_MERGE_2013/NEMOGCM/NEMO/OPA_SRC/DOM/domhgr.F90

Timestamp:

2014-01-22T15:57:03+01:00 (10 years ago)

Author:

jchanut

Message:

correct vvl initialization, ticket #1227

File:

• branches/2013/dev_MERGE_2013/NEMOGCM/NEMO/OPA_SRC/DOM/domhgr.F90 (modified) (1 diff)

branches/2013/dev_MERGE_2013/NEMOGCM/NEMO/OPA_SRC/DOM/domhgr.F90

@@ -458 +458 @@
       ! --------------
       e1e2t(:,:) = e1t(:,:) * e2t(:,:)
-
+    
+      ! Useful shortcuts (JC: note the duplicated e2e2t array ! Need some cleaning)
+      ! ---------------------------------------------------------------------------
+      e12t    (:,:) = e1t(:,:) * e2t(:,:)
+      e12u    (:,:) = e1u(:,:) * e2u(:,:)
+      e12v    (:,:) = e1v(:,:) * e2v(:,:)
+      e12f    (:,:) = e1f(:,:) * e2f(:,:)
+      r1_e12t (:,:) = 1._wp    / e12t(:,:)
+      r1_e12u (:,:) = 1._wp    / e12u(:,:)
+      r1_e12v (:,:) = 1._wp    / e12v(:,:)
+      r1_e12f (:,:) = 1._wp    / e12f(:,:)
+      re2u_e1u(:,:) = e2u(:,:) / e1u(:,:)
+      re1v_e2v(:,:) = e1v(:,:) / e2v(:,:)
 
       ! Control printing : Grid informations (if not restart)