Changeset 97 for codes/icosagcm/trunk/src/caldyn_gcm.f90

Timestamp:

08/04/12 00:13:44 (12 years ago)

Author:

ymipsl

Message:

Put conforming IO for dcmip testcase

YM

File:

• codes/icosagcm/trunk/src/caldyn_gcm.f90 (modified) (4 diffs)

codes/icosagcm/trunk/src/caldyn_gcm.f90

@@ -643 +643 @@
     USE omega_mod
     USE write_field
+    USE vertical_interp_mod
     TYPE(t_field),POINTER :: f_ps(:), f_phis(:), f_u(:), f_theta_rhodz(:), f_dps(:), &
          f_buf_i(:), f_buf_v(:), f_buf_i3(:), f_buf1_i(:), f_buf2_i(:), f_buf_s(:), f_buf_p(:)
+    
+    REAL(rstd) :: out_pression_lev
+    CHARACTER(LEN=255) :: str_pression
+    
+    out_pression_level=0
+    CALL getin("out_pression_level",out_pression_level) 
+    WRITE(str_pression,*) INT(out_pression_level/100)
+    str_pression=ADJUSTL(str_pression)
     
     CALL writefield("ps",f_ps)
@@ -654 +663 @@
     CALL w_omega(f_ps, f_u, f_buf_i)
     CALL writefield('omega', f_buf_i)
+    IF (out_pression_level<=preff) THEN
+      CALL vertical_interp(f_ps,f_buf_i,f_buf_s,out_pression_level)
+      CALL writefield("omega"//TRIM(str_pression),f_buf_s)
+    ENDIF
     
     ! Temperature
     CALL theta_rhodz2temperature(f_ps,f_theta_rhodz,f_buf_i) ;
     CALL writefield("T",f_buf_i)
+    
+    IF (out_pression_level<=preff) THEN
+      CALL vertical_interp(f_ps,f_buf_i,f_buf_s,out_pression_level)
+      CALL writefield("T"//TRIM(str_pression),f_buf_s)
+    ENDIF
     
     ! velocity components
@@ -663 +681 @@
     CALL writefield("ulon",f_buf1_i)
     CALL writefield("ulat",f_buf2_i)
+
+    IF (out_pression_level<=preff) THEN
+      CALL vertical_interp(f_ps,f_buf1_i,f_buf_s,out_pression_level)
+      CALL writefield("ulon"//TRIM(str_pression),f_buf_s)
+      CALL vertical_interp(f_ps,f_buf2_i,f_buf_s,out_pression_level)
+      CALL writefield("vlon"//TRIM(str_pression),f_buf_s)
+    ENDIF
     
     ! geopotential
@@ -670 +695 @@
     CALL writefield("theta",f_buf1_i) ! potential temperature
     CALL writefield("pk",f_buf2_i) ! Exner pressure
+  
   
   END SUBROUTINE write_output_fields