Changeset 254 for codes/icosagcm/branches/SATURN_DYNAMICO/ICOSAGCM/src/physics_lmdz_generic.f90

Timestamp:

07/23/14 18:14:17 (10 years ago)

Author:

ymipsl

Message:

Adding xios output functionnalities

YM

File:

• codes/icosagcm/branches/SATURN_DYNAMICO/ICOSAGCM/src/physics_lmdz_generic.f90 (modified) (7 diffs)

codes/icosagcm/branches/SATURN_DYNAMICO/ICOSAGCM/src/physics_lmdz_generic.f90

@@ -43 +43 @@
   USE mpipara
   USE disvert_mod
+  USE xios_mod
   
   IMPLICIT NONE
@@ -51 +52 @@
   REAL(rstd),ALLOCATABLE :: lonfi(:)
   REAL(rstd),ALLOCATABLE :: airefi(:)
+  REAL(rstd),ALLOCATABLE :: bounds_latfi(:,:)
+  REAL(rstd),ALLOCATABLE :: bounds_lonfi(:,:)
 
     start_day=0
@@ -94 +97 @@
     ALLOCATE(lonfi(nbp_phys))
     ALLOCATE(airefi(nbp_phys))
+    ALLOCATE(bounds_latfi(nbp_phys,6))
+    ALLOCATE(bounds_lonfi(nbp_phys,6))
     
     pos=0
@@ -104 +109 @@
           pos=pos+1
           CALL xyz2lonlat(xyz_i(ij,:),lonfi(pos),latfi(pos))
+          CALL xyz2lonlat(xyz_v(ij+z_rup,:), bounds_lonfi(pos,1), bounds_latfi(pos,1))
+          CALL xyz2lonlat(xyz_v(ij+z_up,:), bounds_lonfi(pos,2), bounds_latfi(pos,2))
+          CALL xyz2lonlat(xyz_v(ij+z_lup,:), bounds_lonfi(pos,3), bounds_latfi(pos,3))
+          CALL xyz2lonlat(xyz_v(ij+z_ldown,:), bounds_lonfi(pos,4), bounds_latfi(pos,4))
+          CALL xyz2lonlat(xyz_v(ij+z_down,:), bounds_lonfi(pos,5), bounds_latfi(pos,5))
+          CALL xyz2lonlat(xyz_v(ij+z_rdown,:), bounds_lonfi(pos,6), bounds_latfi(pos,6))
           airefi(pos)=Ai(ij)  
         ENDDO
@@ -112 +123 @@
     CALL initialize_unstructured_physics(nbp_phys,llm, comm_icosa, mpi_size,distrib,            &
                                          day_length,start_day,itau_physics*dt,                  &
-                                         latfi,lonfi,airefi,radius,g, gas_constant/mu, cpp,     &
+                                         6,latfi,lonfi,airefi,bounds_lonfi,bounds_latfi,         &
+                                         radius,g, gas_constant/mu, cpp,                        &
                                          preff, ap, bp)                                         
     
-    
+    CALL xios_set_context    
     
 !    CALL init_phys_lmdz(128,97,llm, comm_icosa, mpi_size, distrib)
@@ -135 +147 @@
   USE transfert_mod
   USE mpipara
+  USE xios_mod
   IMPLICIT NONE
     INTEGER,INTENT(IN)    :: it
@@ -305 +318 @@
 
 !$OMP BARRIER
+    CALL xios_set_context    
      
   CONTAINS