Changeset 131 for codes/icosagcm/trunk/src/disvert_dcmip3.f90

Timestamp:

02/08/13 12:09:35 (11 years ago)

Author:

ymipsl

Message:

Some operations must be only done by the mpi master task.

YM

File:

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

codes/icosagcm/trunk/src/disvert_dcmip3.f90

@@ -23 +23 @@
   SUBROUTINE disvert(ap,bp,presnivs)
   USE icosa
+  USE mpipara
   IMPLICIT NONE
   REAL(rstd),INTENT(OUT) :: ap(:)
@@ -40 +41 @@
   
   eta_top = (GG/Teq*exp(-N**2*llm*ncar_dz/g)+1-GG/Teq)**(1./kappa)
-  PRINT *,'eta_top ->', eta_top
+  IF (is_mpi_root) PRINT *,'eta_top ->', eta_top
   do l = 1,llm+1
      eta = (GG/Teq*exp(-N**2*(l-1)*ncar_dz/g)+1-GG/Teq)**(1./kappa)
-     PRINT *,'eta ->', eta
+     IF (is_mpi_root) PRINT *,'eta ->', eta
      bp(l) = ((eta - eta_top)/(1 - eta_top))**cindx
      ap(l) = preff * ( eta - bp(l) )
   ENDDO
-  PRINT *,'eta ->', eta
+  IF (is_mpi_root) PRINT *,'eta ->', eta
   bp(1)=1.
   ap(1)=0.
@@ -56 +57 @@
   ENDDO
 
-  PRINT *, 'Vertical placement of model levels according to DCMIP Appendix E.3'
-  PRINT *, 'Parameters : ncar_dz=', ncar_dz, '  ncar_p0=',ncar_p0, '  ncar_disvert_c=',cindx
-  PRINT *, 'Isothermal amtosphere with ncar_T0=',ncar_T0 
+  IF (is_mpi_root) PRINT *, 'Vertical placement of model levels according to DCMIP Appendix E.3'
+  IF (is_mpi_root) PRINT *, 'Parameters : ncar_dz=', ncar_dz, '  ncar_p0=',ncar_p0, '  ncar_disvert_c=',cindx
+  IF (is_mpi_root) PRINT *, 'Isothermal amtosphere with ncar_T0=',ncar_T0 
 
 END SUBROUTINE disvert