Changeset 155 for codes/icosagcm/trunk/src/metric.f90

Timestamp:

05/20/13 16:35:23 (11 years ago)

Author:

dubos

Message:

Schmidt transform in metric.f90

File:

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

codes/icosagcm/trunk/src/metric.f90

@@ -52 +52 @@
   INTEGER,PARAMETER :: vdown=5
   INTEGER,PARAMETER :: vrdown=6
-
   
 CONTAINS
-  
+
+  SUBROUTINE remap_schmidt
+    USE spherical_geom_mod
+    USE ioipsl
+    IMPLICIT NONE
+    INTEGER :: nf,i,j
+    REAL(rstd) :: schmidt_factor, schmidt_lon, schmidt_lat
+
+    ! Schmidt transform parameters
+    schmidt_factor = 1.
+    CALL getin('schmidt_factor', schmidt_factor)
+    schmidt_factor =  schmidt_factor**2.
+    
+    schmidt_lon = 0.
+    CALL getin('schmidt_lon', schmidt_lon)
+    schmidt_lon = schmidt_lon * pi/180.
+
+    schmidt_lat = 45.
+    CALL getin('schmidt_lat', schmidt_lat)
+    schmidt_lat = schmidt_lat * pi/180.
+
+    DO nf=1,nb_face
+       DO j=1,jjm_glo
+          DO i=1,iim_glo
+             CALL schmidt_transform(vertex_glo(i,j,nf)%xyz, schmidt_factor, schmidt_lon, schmidt_lat)
+          END DO
+       END DO
+    END DO
+  END SUBROUTINE remap_schmidt
+
   SUBROUTINE allocate_metric
   IMPLICIT NONE
@@ -830 +858 @@
 !    CALL compute_face
     CALL compute_face_projection
+    CALL remap_schmidt
+
     CALL set_index
     CALL set_cell