Changeset 203 for codes/icosagcm/trunk/src/vector.f90

Timestamp:

07/09/14 00:58:30 (10 years ago)

Author:

dubos

Message:

Upgraded JW06 and DCMIP5 to new etat0 interface (tested)

File:

• codes/icosagcm/trunk/src/vector.f90 (modified) (1 diff)

codes/icosagcm/trunk/src/vector.f90

@@ -48 +48 @@
     END SUBROUTINE cross_product2
 
-    FUNCTION arc(lon1,lat1, lon2,lat2)
-      REAL(rstd) :: lon1,lat1, lon2,lat2, arc2, &
-           x1,x2,x3,y1,y2,y3,z1,z2,z3
-      x1=cos(lat1)*cos(lon1)
-      x2=cos(lat2)*cos(lon2)
-      y1=cos(lat1)*sin(lon1)
-      y2=cos(lat2)*sin(lon2)
-      z1=sin(lat1)
-      z2=sin(lat2)
-      x3=y1*z2-y2*z1
-      y3=z1*x2-z2*x1
-      z3=x1*y2-x2*y1
-      arc=ASIN(SQRT(x3*x3+y3*y3+z3*z3))
+    FUNCTION arc(lon,lat, lonc,latc)
+      REAL(rstd) :: lon,lat, lonc,latc, arc
+      arc=ACOS(sin(latc)*sin(lat)+cos(latc)*cos(lat)*cos(lon-lonc))
     END FUNCTION arc