Changeset 3895 for dynamico_lmdz/aquaplanet/LMDZ5/libf/phylmd/regr_pr_av_m.F90

Timestamp:

Jan 25, 2016, 10:47:14 AM (9 years ago)

Author:

ymipsl

Message:

Make LMDZ5 be compliant to generate initiale state and compute in openMP mode suing unstructured grid.

YM

File:

• dynamico_lmdz/aquaplanet/LMDZ5/libf/phylmd/regr_pr_av_m.F90 (modified) (4 diffs)

dynamico_lmdz/aquaplanet/LMDZ5/libf/phylmd/regr_pr_av_m.F90

@@ -46 +46 @@
 
     integer, intent(in):: ncid ! NetCDF ID of the file
+
+!ym   character(len=*), intent(in):: name(:) ! of the NetCDF variables
+!ym  for strange reason, ifort doens't detect correctly the array size (return 0) 
+!ym   compilator bug , version dependent ?
+
     character(len=*), intent(in):: name(:) ! of the NetCDF variables
     integer, intent(in):: julien ! jour julien, 1 <= julien <= 360
@@ -66 +71 @@
     integer varid, ncerr ! for NetCDF
 
-    real  v1(nbp_lon, nbp_lat, size(press_in_edg) - 1, size(name))
+    real  v1(nbp_lon, nbp_lat, size(press_in_edg) - 1, size(v3,3))
     ! input fields at day "julien", on the global "dynamics" horizontal grid
     ! First dimension is for longitude.
@@ -74 +79 @@
     ! NetCDF variable "name(l)".
 
-    real v2(klon, size(press_in_edg) - 1, size(name))
+    real v2(klon, size(press_in_edg) - 1, size(v3,3))
     ! fields scattered to the partial "physics" horizontal grid
     ! "v2(i, k, l)" is at longitude "xlon(i)", latitude "xlat(i)",
@@ -85 +90 @@
 
     call assert(size(v3, 1) == klon, size(v3, 2) == klev, "regr_pr_av v3 klon")
-    n_var = assert_eq(size(name), size(v3, 3), "regr_pr_av v3 n_var")
+    PRINT *,'size name',size(name)
+!ym    n_var = assert_eq(size(name), size(v3, 3), "regr_pr_av v3 n_var")
+!ym intel bugs compiler : ifort doesn't detect correcte size of nme (return 0)
+
+    n_var = size(v3, 3)
+
     call assert(shape(paprs) == (/klon, klev+1/), "regr_pr_av paprs")