Changeset 581 for codes/icosagcm/trunk/src/output

Timestamp:

10/13/17 16:00:58 (7 years ago)

Author:

dubos

Message:

trunk : upgrading to devel

File:

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

codes/icosagcm/trunk/src/output/write_etat0.f90

@@ -5 +5 @@
 CONTAINS
 
-  SUBROUTINE write_etat0(it,f_ps,f_phis,f_theta_rhodz,f_u, f_q)
+  SUBROUTINE write_etat0(it,f_ps,f_phis,f_theta_rhodz,f_u, f_q, f_geopot, f_W)
   USE icosa
   USE restart_mod
@@ -20 +20 @@
     TYPE(t_field),POINTER :: f_u(:)
     TYPE(t_field),POINTER :: f_q(:)
+    TYPE(t_field),POINTER, OPTIONAL :: f_geopot(:)
+    TYPE(t_field),POINTER, OPTIONAL :: f_W(:)
   
     TYPE(t_field),POINTER,SAVE :: f_ulon(:)
@@ -41 +43 @@
     CALL transfert_request(f_u,req_e1_vect)
     CALL un2ulonlat(f_u, f_ulon, f_ulat)
-    
-    CALL write_restart(it,f_ps,f_phis,f_theta_rhodz_1d,f_u, f_ulon, f_ulat, f_q)
 
+    IF(hydrostatic) THEN
+       CALL write_restart(it,f_ps,f_phis,f_theta_rhodz_1d,f_u, f_ulon, f_ulat, f_q)
+    ELSE
+       CALL write_restart(it,f_ps,f_phis,f_theta_rhodz_1d,f_u, f_ulon, f_ulat, f_q, f_geopot, f_W)
+    END IF
     CALL deallocate_field(f_ulon)
     CALL deallocate_field(f_ulat)