Changeset 430

Timestamp:

06/14/16 23:37:59 (8 years ago)

Author:

ymipsl

Message:

physiq dcmip2016 : Add column integration diagnostique for chlorines

YM

File:

• codes/icosagcm/trunk/src/physics_dcmip2016.f90 (modified) (3 diffs)

codes/icosagcm/trunk/src/physics_dcmip2016.f90

@@ -12 +12 @@
   REAL(rstd),ALLOCATABLE :: precl_packed(:)
 
+  TYPE(t_field),POINTER  :: f_Q1(:)
+  TYPE(t_field),POINTER  :: f_Q2(:)
+  TYPE(t_field),POINTER  :: f_PS(:)
+  TYPE(t_field),POINTER  :: f_rhodz(:)
+  TYPE(t_field),POINTER  :: f_Q1_col_int(:)
+  TYPE(t_field),POINTER  :: f_Q2_col_int(:)
   PUBLIC :: init_physics, full_physics, write_physics
 
@@ -74 +80 @@
     ALLOCATE(precl_packed(ngrid))
     CALL allocate_field(f_precl, field_t,type_real)
+    CALL allocate_field(f_Q1, field_t,type_real,llm)
+    CALL allocate_field(f_Q2, field_t,type_real,llm)
+    CALL allocate_field(f_PS, field_t,type_real)
+    CALL allocate_field(f_rhodz, field_t,type_real,llm)
+    CALL allocate_field(f_Q1_col_int, field_t,type_real)
+    CALL allocate_field(f_Q2_col_int, field_t,type_real)
 
     PRINT *, 'init_physics_new', SIZE(physics_inout%Ai)
@@ -90 +102 @@
     USE output_field_mod
     USE physics_interface_mod
+    use disvert_mod
+    REAL(rstd), POINTER :: Q1(:,:)
+    REAL(rstd), POINTER :: Q2(:,:)
+    REAL(rstd), POINTER :: PS(:)
+    REAL(rstd), POINTER :: rhodz(:,:)
+    REAL(rstd), POINTER :: Q1_col_int(:)
+    REAL(rstd), POINTER :: Q2_col_int(:)
+    
+    
     CALL unpack_field(f_precl, precl_packed)
     CALL output_field("precl",f_precl)
-       
+    
+    CALL unpack_field(f_Q1,physics_inout%q(:,:,4))
+    CALL unpack_field(f_Q2,physics_inout%q(:,:,5))
+    CALL unpack_field(f_ps,physics_inout%p(:,1))
+    
+    DO ind=1,ndomain
+      IF (.NOT. assigned_domain(ind)) CYCLE
+        Q1=f_Q1(ind)
+        Q2=f_Q2(ind)
+        Q1_col_int=f_Q1_col_int(ind)
+        Q2_col_int=f_Q2_col_int(ind)
+        PS=f_PS(ind)
+        rhodz=f_rhodz(ind)
+        DO l=1,llm
+          rhodz(:,l)= ( ap(l)-ap(l+1) + (bp(l)-bp(l+1))*ps(:) )/g   
+        ENDDO
+        Q1_col_int=SUM(Q1*rhodz,2)/SUM(rhodz,2)
+        Q2_col_int=SUM(Q2*rhodz,2)/SUM(rhodz,2)
+    ENDDO
+     
+    CALL output_field("Q1_col_int",f_Q1_col_int)
+    CALL output_field("Q2_col_int",f_Q2_col_int)
+    
   END SUBROUTINE write_physics