Changeset 5353 for LMDZ6/trunk/libf/phylmd/regr_horiz_time_climoz_m.f90

Timestamp:

Nov 29, 2024, 3:15:38 PM (2 weeks ago)

Author:

yann meurdesoif

Message:

Nvidia compiler has some difficulties to compile correctly some complex array constructor.
This commit decompose it the several phases in order to achieve the compilation.
Please Lionel and David, have a look to this in order to validate.
Probably, in future, when compiler heuristic will be improved, this commit can be reversed.
YM

File:

• LMDZ6/trunk/libf/phylmd/regr_horiz_time_climoz_m.f90 (modified) (4 diffs)

LMDZ6/trunk/libf/phylmd/regr_horiz_time_climoz_m.f90

@@ -139 +139 @@
   REAL, ALLOCATABLE :: test_o3_in(:,:)
   REAL, ALLOCATABLE :: test_o3_out(:)
+  REAL,ALLOCATABLE :: tmp(:)
 
 
@@ -427 +428 @@
        !--- Regrid in longitude
         ALLOCATE(o3_regr_lon(nlon_ou, nlat_in, nlev_in, ie-ib+1, read_climoz))
+        tmp = [boundslon_reg(1,west),boundslon_reg(:,east)]
         CALL regr_conserv(1, o3_in3, xs = lon_in_edge,                             &
-                            xt = [boundslon_reg(1,west),boundslon_reg(:,east)],    &
-                            vt = o3_regr_lon, slope = slopes(1,o3_in3, lon_in_edge))
+                          xt = tmp ,                                               &
+                          vt = o3_regr_lon, slope = slopes(1,o3_in3, lon_in_edge))
         DEALLOCATE(o3_in3)
 
@@ -436 +438 @@
         !--- (inverted indices in "o3_regr_lonlat" because "rlatu" is decreasing)
         ALLOCATE(o3_regr_lonlat(nlon_ou, nlat_ou, nlev_in, 0:13, read_climoz))
-        CALL regr_conserv(2, o3_regr_lon, xs = sinlat_in_edge,                     &
-                        xt = [- 1., SIN(boundslat_reg(nlat_ou-1:1:-1,south)), 1.], &
-                        vt = o3_regr_lonlat(:,nlat_ou:1:- 1,:,ib:ie,:),            &
-                   slope = slopes(2,o3_regr_lon, sinlat_in_edge))
+        tmp = [- 1., SIN(boundslat_reg(nlat_ou-1:1:-1,south)), 1.]   
+        CALL regr_conserv(2, o3_regr_lon, xs = sinlat_in_edge,                       &
+                          xt = tmp,                                                  &
+                          vt = o3_regr_lonlat(:,nlat_ou:1:- 1,:,ib:ie,:),            &
+                          slope = slopes(2,o3_regr_lon, sinlat_in_edge))
         DEALLOCATE(o3_regr_lon)
 
@@ -519 +522 @@
         !--- (inverted indices in "o3_regr_lat" because "rlatu" is decreasing)
         ALLOCATE(o3_regr_lat(nlat_ou, nlev_in, 0:13, read_climoz))
+        tmp = [- 1., SIN(boundslat_reg(nlat_ou-1:1:-1,south)), 1.]
         CALL regr_conserv(1, o3_in2, xs = sinlat_in_edge,                          &
-                        xt = [- 1., SIN(boundslat_reg(nlat_ou-1:1:-1,south)), 1.], &
+                        xt = tmp,                                                  &
                         vt = o3_regr_lat(nlat_ou:1:- 1,:,ib:ie,:),                 &
                      slope = slopes(1,o3_in2, sinlat_in_edge))