Changeset 36

Timestamp:

07/27/12 16:10:47 (12 years ago)

Author:

ymipsl

Message:

add compute_wind_perp_from_lonlat_compound functionnality

YM

File:

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

codes/icosagcm/trunk/src/wind.f90

@@ -166 +166 @@
   END SUBROUTINE compute_wind_from_lonlat_compound
  
- 
+  SUBROUTINE compute_wind_perp_from_lonlat_compound(ulon, ulat, up)
+  USE icosa  
+    
+  IMPLICIT NONE
+  REAL(rstd) :: up(3*iim*jjm,llm)
+  REAL(rstd) :: ulon(3*iim*jjm,llm)
+  REAL(rstd) :: ulat(3*iim*jjm,llm)
+  REAL(rstd) :: u(3*iim*jjm,3,llm)
+
+  INTEGER :: i,j,ij,l     
+  
+   CALL compute_wind_from_lonlat_compound(ulon, ulat, u)
+
+    DO l=1,llm
+      DO j=jj_begin-1,jj_end+1
+        DO i=ii_begin-1,ii_end+1
+          ij=(j-1)*iim+i
+          up(ij+u_right,l)=sum(u(ij+u_right,:,l)*ep_e(ij+u_right,:))
+          up(ij+u_lup,l)=sum(u(ij+u_lup,:,l)*ep_e(ij+u_lup,:))
+          up(ij+u_ldown,l)=sum(u(ij+u_ldown,:,l)*ep_e(ij+u_ldown,:))
+        ENDDO
+      ENDDO
+    ENDDO
+ 
+  END SUBROUTINE compute_wind_perp_from_lonlat_compound
+   
  SUBROUTINE compute_wind_centered_lonlat_compound(uc, ulon, ulat)
   USE icosa  
     
   IMPLICIT NONE
-  REAL(rstd) :: uc(3*iim*jjm,3,llm)
-  REAL(rstd) :: ulon(3*iim*jjm,3,llm)
-  REAL(rstd) :: ulat(3*iim*jjm,3,llm)
+  REAL(rstd) :: uc(iim*jjm,3,llm)
+  REAL(rstd) :: ulon(iim*jjm,3,llm)
+  REAL(rstd) :: ulat(iim*jjm,3,llm)
 
   INTEGER :: i,j,ij,l     
@@ -182 +207 @@
         DO i=ii_begin,ii_end
           ij=(j-1)*iim+i
-          ulon(ij,:,l)=sum(uc(ij,:,l)*elon_i(ij,:))*elon_i(ij,:) 
-          ulat(ij,:,l)=sum(uc(ij,:,l)*elat_i(ij,:))**elat_i(ij,:) 
+          ulon(ij,l)=sum(uc(ij,:,l)*elon_i(ij,:))
+          ulat(ij,l)=sum(uc(ij,:,l)*elat_i(ij,:)) 
         ENDDO
       ENDDO