Changes in XMLIO_SERVER/trunk/src/XMLIO/mod_grid.f90 [17:29]

File:

• XMLIO_SERVER/trunk/src/XMLIO/mod_grid.f90 (modified) (8 diffs)

XMLIO_SERVER/trunk/src/XMLIO/mod_grid.f90

@@ -3 +3 @@
   USE mod_sorted_list
   USE mod_domain
+  USE mod_zoom
+
   IMPLICIT NONE
 
@@ -19 +21 @@
     INTEGER                     :: nj
     LOGICAL                     :: has_dimension
+    TYPE(vector_zoom),POINTER   :: associated_zoom
+    TYPE(zoom),POINTER          :: global_zoom
   END TYPE grid
 
@@ -29 +33 @@
   INCLUDE 'vector_grid_contains.inc'
 
+  SUBROUTINE grid__swap_context(saved_grid_Ids,saved_Ids)
+  IMPLICIT NONE
+    TYPE(vector_grid),POINTER          :: saved_grid_Ids
+    TYPE(sorted_list),POINTER          :: saved_Ids 
+    
+    grid_ids=>saved_grid_ids
+    ids=>saved_ids
+  END SUBROUTINE grid__swap_context
+  
+  
   SUBROUTINE grid__init
   IMPLICIT NONE
-    
-    ALLOCATE(grid_Ids)
-    ALLOCATE(Ids)
     
     CALL vector_grid__new(grid_Ids)
@@ -68 +79 @@
    ALLOCATE(pt_grid%subdomain)
    ALLOCATE(pt_grid%rank_ids)
+   ALLOCATE(pt_grid%associated_zoom)
+   
    CALL domain__new(pt_grid%domain)
    CALL vector_domain__new(pt_grid%subdomain)
    CALL sorted_list__new(pt_grid%rank_ids)
+   CALL vector_zoom__new(pt_grid%associated_zoom)
    
    pt_grid%has_id          = .FALSE.
@@ -83 +97 @@
      CALL sorted_list__Add(Ids,hash(id),Pos)
    ENDIF
+   
+   CALL grid__get_new_zoom(pt_grid,pt_grid%global_zoom,id)
 
  END SUBROUTINE grid__new
@@ -153 +169 @@
     TYPE(grid), POINTER  :: pt_grid
     TYPE(domain),POINTER :: subdomain
+    TYPE(zoom),POINTER :: pt_zoom
     
     REAL,ALLOCATABLE :: lon(:,:)
     REAL,ALLOCATABLE :: lat(:,:)
-    INTEGER :: ib,ie,jb,je,ni,nj,ibegin,jbegin
+    INTEGER :: ib,ie,jb,je,ni,nj,ibegin,jbegin,iend,jend
     INTEGER :: i
     
@@ -196 +213 @@
       
       CALL domain__set(pt_grid%domain,0,ni,nj,ibegin,jbegin,lon,lat)
-
-      
+      iend=ibegin+ni-1
+      jend=jbegin+nj-1
+      
+      
+      pt_grid%global_zoom%ni_glo=pt_grid%ni
+      pt_grid%global_zoom%nj_glo=pt_grid%nj
+      pt_grid%global_zoom%ibegin_glo=1      
+      pt_grid%global_zoom%jbegin_glo=1
+      
+      DO i=1,pt_grid%associated_zoom%size
+        pt_zoom=>pt_grid%associated_zoom%at(i)%pt
+        
+        ib=MAX(pt_zoom%ibegin_glo-ibegin+1,1)
+        ie=MIN(pt_zoom%ibegin_glo+pt_zoom%ni_glo-ibegin,ni)
+        pt_zoom%ni_loc=MAX(ie-ib+1,0)
+        pt_zoom%ibegin_loc=ib
+
+        jb=MAX(pt_zoom%jbegin_glo-jbegin+1,1)
+        je=MIN(pt_zoom%jbegin_glo+pt_zoom%nj_glo-jbegin,nj)
+        pt_zoom%nj_loc=MAX(je-jb+1,0)
+        pt_zoom%jbegin_loc=jb
+      ENDDO
+                 
+            
       DEALLOCATE(lon)
       DEALLOCATE(lat)
@@ -204 +243 @@
       
             
-        
+  SUBROUTINE grid__get_new_zoom(pt_grid,pt_zoom,zoom_id)
+  USE string_function
+  IMPLICIT NONE
+    TYPE(grid), POINTER                      :: pt_grid
+    TYPE(zoom),POINTER                       :: pt_zoom
+    CHARACTER(LEN=*),INTENT(IN),OPTIONAL     :: zoom_id
+    LOGICAL                                  :: success
+    
+     CALL vector_zoom__get_new(pt_grid%associated_zoom,Pt_zoom)
+     CALL zoom__new(Pt_zoom,zoom_id)
+    
+   END SUBROUTINE grid__get_new_zoom
+     
   SUBROUTINE grid__print(pt_grid)
   IMPLICIT NONE