Changeset 476 for codes/icosagcm/trunk/src/write_etat0.f90

Timestamp:

09/12/16 14:39:01 (8 years ago)

Author:

ymipsl

Message:

Bug fix : read and write restart file was brocken due to new extra dimension to theta_rhodz

YM

File:

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

codes/icosagcm/trunk/src/write_etat0.f90

@@ -10 +10 @@
   USE wind_mod
   USE write_field_mod
+  USE domain_mod
+  USE omp_para
   IMPLICIT NONE
     INTEGER,INTENT(IN)    :: it
@@ -20 +22 @@
     TYPE(t_field),POINTER,SAVE :: f_ulon(:)
     TYPE(t_field),POINTER,SAVE :: f_ulat(:)
+    TYPE(t_field),POINTER,SAVE :: f_theta_rhodz_1d(:)
+    REAL(rstd), POINTER :: theta_rhodz(:,:,:),theta_rhodz_1d(:,:)
+    INTEGER :: ind
     
     
     CALL allocate_field(f_ulon,field_t,type_real,llm,name='ulon')
     CALL allocate_field(f_ulat,field_t,type_real,llm,name='ulat')
+    CALL allocate_field(f_theta_rhodz_1d,field_t,type_real,llm,name='theta_rhodz')
 
+!$OMP BARRIER    
+    DO ind=1, ndomain
+       IF (.NOT. assigned_domain(ind) .AND. is_omp_level_master) CYCLE
+       theta_rhodz=f_theta_rhodz(ind) ; theta_rhodz_1d=f_theta_rhodz_1d(ind)
+       theta_rhodz_1d(:,:)=theta_rhodz(:,:,1)
+    ENDDO
+    
     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,f_u, f_ulon, f_ulat, f_q)
+    CALL write_restart(it,f_ps,f_phis,f_theta_rhodz_1d,f_u, f_ulon, f_ulat, f_q)
 
     CALL deallocate_field(f_ulon)