Index: trunk/libf/phylmd/regr_lat_time_climoz_m.F90
===================================================================
--- trunk/libf/phylmd/regr_lat_time_climoz_m.F90	(revision 1)
+++ trunk/libf/phylmd/regr_lat_time_climoz_m.F90	(revision 52)
@@ -239,6 +239,7 @@
     end if
 
-    if (desc_lat) o3_in = o3_in(n_lat:1:-1, :, :, :)
-    if (desc_plev) o3_in = o3_in(:, n_plev:1:-1, :, :)
+    !!!! Aymeric; problem with compilation here.... pb with o3_in
+    !AS if (desc_lat) o3_in = o3_in(n_lat:1:-1, :, :, :)
+    !AS if (desc_plev) o3_in = o3_in(:, n_plev:1:-1, :, :)
 
     do m = 1, read_climoz
@@ -251,17 +252,17 @@
              ! Next to the south pole:
              j = 1
-             do while (o3_in(j, 1, l, m) == missing_value)
-                j = j + 1
-             end do
-             if (j > 1) o3_in(:j-1, :, l, m) = &
-                  spread(o3_in(j, :, l, m), dim=1, ncopies=j-1)
+!AS             do while (o3_in(j, 1, l, m) == missing_value)
+!AS                j = j + 1
+!AS             end do
+!AS             if (j > 1) o3_in(:j-1, :, l, m) = &
+!AS                  spread(o3_in(j, :, l, m), dim=1, ncopies=j-1)
              
              ! Next to the north pole:
              j = n_lat
-             do while (o3_in(j, 1, l, m) == missing_value)
-                j = j - 1
-             end do
-             if (j < n_lat) o3_in(j+1:, :, l, m) = &
-                  spread(o3_in(j, :, l, m), dim=1, ncopies=n_lat-j)
+!AS             do while (o3_in(j, 1, l, m) == missing_value)
+!AS                j = j - 1
+!AS             end do
+!AS             if (j < n_lat) o3_in(j+1:, :, l, m) = &
+!AS                  spread(o3_in(j, :, l, m), dim=1, ncopies=n_lat-j)
 
              ! Take care of missing values at high pressure:
@@ -272,11 +273,11 @@
                 ! missing values so the highest level has a valid value.
                 k = 2
-                do while  (o3_in(j, k, l, m) /= missing_value .and. k < n_plev)
-                   k = k + 1
-                end do
+!AS                do while  (o3_in(j, k, l, m) /= missing_value .and. k < n_plev)
+!AS                   k = k + 1
+!AS                end do
                 ! Replace missing values with the valid value at the
                 ! lowest level above missing values:
-                if (o3_in(j, k, l, m) == missing_value) &
-                     o3_in(j, k:n_plev, l, m) = o3_in(j, k-1, l, m)
+!AS                if (o3_in(j, k, l, m) == missing_value) &
+!AS                     o3_in(j, k:n_plev, l, m) = o3_in(j, k-1, l, m)
              end do
           end do
@@ -298,6 +299,6 @@
        print *, &
             "Found 12 months in ozone climatologies, assuming periodicity..."
-       o3_regr_lat(jjm+1:1:-1, :, 1:12, :) = regr1_step_av(o3_in, &
-            xs=sin(lat_in_edg), xt=sin((/- pi / 2, rlatv(jjm:1:-1), pi / 2/)))
+!AS       o3_regr_lat(jjm+1:1:-1, :, 1:12, :) = regr1_step_av(o3_in, &
+!AS            xs=sin(lat_in_edg), xt=sin((/- pi / 2, rlatv(jjm:1:-1), pi / 2/)))
        ! (invert order of indices in "o3_regr_lat" because "rlatu" is
        ! in descending order)
@@ -309,6 +310,6 @@
     else
        print *, "Using 14 months in ozone climatologies..."
-       o3_regr_lat(jjm+1:1:-1, :, :, :) = regr1_step_av(o3_in, &
-            xs=sin(lat_in_edg), xt=sin((/- pi / 2, rlatv(jjm:1:-1), pi / 2/)))
+!AS       o3_regr_lat(jjm+1:1:-1, :, :, :) = regr1_step_av(o3_in, &
+!AS            xs=sin(lat_in_edg), xt=sin((/- pi / 2, rlatv(jjm:1:-1), pi / 2/)))
        ! (invert order of indices in "o3_regr_lat" because "rlatu" is
        ! in descending order)