Changeset 403

Timestamp:

06/08/16 02:55:08 (8 years ago)

Author:

dubos

Message:

Fixed output of phis

Location:

codes/icosagcm/trunk/src

Files:

• observable.f90 (modified) (3 diffs)
• timeloop_gcm.f90 (modified) (2 diffs)

codes/icosagcm/trunk/src/observable.f90

@@ -37 +37 @@
   END SUBROUTINE init_observable
   
-  SUBROUTINE write_output_fields_basic(f_ps, f_mass, f_geopot, f_theta_rhodz, f_u, f_W, f_q)
+  SUBROUTINE write_output_fields_basic(f_phis, f_ps, f_mass, f_geopot, f_theta_rhodz, f_u, f_W, f_q)
     USE wind_mod
     USE output_field_mod
@@ -51 +51 @@
     USE omega_mod
     
-    TYPE(t_field),POINTER :: f_ps(:), f_mass(:), f_geopot(:), f_theta_rhodz(:), f_u(:), f_W(:), f_q(:)
+    TYPE(t_field),POINTER :: f_phis(:), f_ps(:), f_mass(:), f_geopot(:), f_theta_rhodz(:), f_u(:), f_W(:), f_q(:)
 !    IF (is_master) PRINT *,'CALL write_output_fields_basic'
     REAL :: scalar(1)
@@ -71 +71 @@
       CALL xios_send_field("mid_ap",mid_ap)
       CALL xios_send_field("mid_bp",mid_bp)
+
+      CALL output_field("phis",f_phis)
         
       first_output=.FALSE.

codes/icosagcm/trunk/src/timeloop_gcm.f90

@@ -150 +150 @@
     CALL transfert_request(f_phis,req_i0) 
     CALL transfert_request(f_phis,req_i1) 
-    CALL writefield("phis",f_phis,once=.TRUE.)
 
     CALL init_message(f_ps,req_i0,req_ps0)
@@ -251 +250 @@
        IF (mod(it,itau_out)==0 ) THEN
           CALL transfert_request(f_u,req_e1_vect)
-          CALL write_output_fields_basic(f_ps, f_mass, f_geopot, f_theta_rhodz, f_u, f_W, f_q)
+          CALL write_output_fields_basic(f_phis, f_ps, f_mass, f_geopot, f_theta_rhodz, f_u, f_W, f_q)
        ENDIF