Changeset 198 for codes/icosagcm/trunk

Timestamp:

07/07/14 14:54:04 (10 years ago)

Author:

ymipsl

Message:

For intel compiler, bad inlining when using multithreading and high optimization level, causing crash.
Call has been replaced (temporarily) by corresponding code.

YM

File:

• codes/icosagcm/trunk/src/check_conserve.f90 (modified) (2 diffs)

codes/icosagcm/trunk/src/check_conserve.f90

@@ -200 +200 @@
     REAL(rstd)::masse(iim*jjm,llm)  
     REAL(rstd)::rad,radomeg,lat,lon
+    REAL(rstd) ::xyz(3)
     INTEGER :: i,j,ij,l,ij2
 
@@ -255 +256 @@
           ij=(j-1)*iim+i
           IF (domain(ind)%own(i,j)) THEN
-            CALL xyz2lonlat(xyz_i(ij,:),lon,lat) 
+!            CALL xyz2lonlat(xyz_i(ij,:),lon,lat) 
+! very bad inlining with intel compiler (>14 ?)
+             xyz(:)=xyz_i(ij,:)
+             xyz(:)=xyz(:)/sqrt(sum(xyz(:)**2))
+             lat=asin(xyz(3))
+             lon=atan2(xyz(2),xyz(1))
+  
             ang=ang+rad*cos(lat)*masse(ij,l)*(ulon(ij,l)+ radomeg*cos(lat))
           END IF