Changeset 196 for codes/icosagcm/trunk/src/wind.f90

Timestamp:

07/06/14 09:10:13 (10 years ago)

Author:

dubos

Message:

First draft of generic dynamics-physics interface - works with DCMIP5.1

File:

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

codes/icosagcm/trunk/src/wind.f90

@@ -185 +185 @@
   END SUBROUTINE compute_wind_from_lonlat_compound
  
+  SUBROUTINE compute_wind_centered_from_lonlat_compound(ulon, ulat, u)
+  USE icosa  
+    
+  IMPLICIT NONE
+  REAL(rstd) :: u(iim*jjm,3,llm)
+  REAL(rstd) :: ulon(iim*jjm,llm)
+  REAL(rstd) :: ulat(iim*jjm,llm)
+
+  INTEGER :: i,j,ij,l     
+  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
+          u(ij,:,l)=ulon(ij,l)*elon_i(ij,:)+ ulat(ij,l)*elat_i(ij,:)
+        ENDDO
+      ENDDO
+    ENDDO
+ 
+  END SUBROUTINE compute_wind_centered_from_lonlat_compound
+ 
   SUBROUTINE compute_wind2D_from_lonlat_compound(ulon, ulat, u)
   USE icosa