Changes between Version 2 and Version 3 of DevelopmentActivities/Isotopes

Timestamp:

2012-07-17T10:31:53+02:00 (12 years ago)

Author:

dsolyga

Comment:

--

DevelopmentActivities/Isotopes

@@ -102 +102 @@
 
 
+== Answers from Yann (17/07) : ==
 
+   * Yann prefers the explicit loops except for initialization because : a) the compiler knows the bounds so it can optimize the loops b) It gives more information to the developer. But we can use the vector notation
+     for the first index. For example :
+{{{
+DO k = 1,niso
+     DO j = 1, nvm
+          transpir(:,j,k) = 
+     END DO
+END DO
+}}}
 
+   * Yann advises to set the dimension of the isotopes as the last dimension. The explicit loop over the last dimension will be mandatory because of Fortran optimization.
+     This last example is not recommended : 
+{{{
+    DO ji = 1,kjpindex
+       param(ji,:) = alpha * param(ji,:) + (1. - alpha)
+    END DO
+}}}
+    The same advice if we pass a variable to a subroutine. We can pass an array to a subroutine by indicating the full dimensions as input variables (it is currently done in the model).
+   * Concerning the equations : No specific recommendation.