Changeset 1325 for trunk/NEMO/OFF_SRC/opa.F90

Timestamp:

2009-02-20T11:03:39+01:00 (15 years ago)

Author:

cetlod

Message:

rewrite the offline opa.F90 module, see ticket:349

File:

• trunk/NEMO/OFF_SRC/opa.F90 (modified) (6 diffs)

trunk/NEMO/OFF_SRC/opa.F90

@@ -33 +33 @@
    PRIVATE
 
+   !! * Module variables
+   CHARACTER (len=64) ::        &
+      cform_aaa="( /, 'AAAAAAAA', / ) "     ! flag for output listing
+
    !! * Routine accessibility
    PUBLIC opa_model      ! called by model.F90
+   PUBLIC opa_init
    !!----------------------------------------------------------------------
    !!   OPA 9.0 , LOCEAN-IPSL  (2005)
@@ -47 +52 @@
       !!                     ***  ROUTINE opa  ***
       !!
+      !! ** Purpose :   opa solves the primitive equations on an orthogonal
+      !!      curvilinear mesh on the sphere.
+      !!
+      !! ** Method  : - model general initialization
+      !!              - launch the time-stepping (stp routine)
+      !!
+      !! References :
+      !!      Madec, Delecluse,Imbard, and Levy, 1997: reference manual.
+      !!              internal report, IPSL.
+      !!----------------------------------------------------------------------
+      INTEGER ::   istp       ! time step index
+      !!----------------------------------------------------------------------
+
+      CALL opa_init  ! Initializations
+
+      IF( lk_mpp )   CALL mpp_max( nstop )
+
+      ! check that all process are still there... If some process have an error,
+      ! they will never enter in step and other processes will wait until the end of the cpu time!
+      IF( lk_mpp )   CALL mpp_max( nstop )
+
+      istp = nit000
+         !
+      DO WHILE ( istp <= nitend .AND. nstop == 0 )
+         CALL stp( istp )
+         istp = istp + 1
+         IF( lk_mpp )   CALL mpp_max( nstop )
+      END DO
+      !                                     ! ========= !
+      !                                     !  Job end  !
+      !                                     ! ========= !
+
+      IF(lwp) WRITE(numout,cform_aaa)       ! Flag AAAAAAA
+
+      IF( nstop /= 0 .AND. lwp ) THEN                 ! error print
+         WRITE(numout,cform_err)
+         WRITE(numout,*) nstop, ' error have been found'
+      ENDIF
+
+      CALL opa_closefile
+
+      IF( lk_mpp )   CALL mppstop                          ! Close all files (mpp)
+      !
+   END SUBROUTINE opa_model
+
+
+   SUBROUTINE opa_init
+      !!----------------------------------------------------------------------
+      !!                     ***  ROUTINE opa_init ***
+      !!
       !! ** Purpose :   opa solves the primitive equations on an orthogonal 
       !!      curvilinear mesh on the sphere.
       !!
       !! ** Method  : - model general initialization
-      !!              - launch the time-stepping (stp routine)
       !!
       !! References :
@@ -76 +130 @@
       !!----------------------------------------------------------------------
       !! * Local declarations
-      INTEGER ::   istp       ! time step index
       CHARACTER (len=20) ::   namelistname
       CHARACTER (len=28) ::   file_out
-      CHARACTER (len=64) ::        &
-         cform_aaa="( /, 'AAAAAAAA', / ) "     ! flag for output listing
 
       !!----------------------------------------------------------------------
@@ -150 +201 @@
       CALL istate_init                      ! ocean initial state (Dynamics and tracers)
 
-      CALL trc_ini                           ! Passive tracers
+!!!!      CALL trc_ini                           ! Passive tracers
 
       CALL day_init                          ! Calendar
@@ -157 +208 @@
       !                                     ! Ocean physics
       CALL tra_qsr_init                         ! Solar radiation penetration
+
+      CALL trc_ini                           ! Passive tracers
 
 #if ! defined key_off_degrad
@@ -167 +220 @@
       IF(lwp) WRITE(numout,cform_aaa)       ! Flag AAAAAAA
 
-      istp = nit000
-      DO WHILE ( istp <= nitend .AND. nstop == 0 )
-         CALL stp( istp )
-         istp = istp + 1
-      END DO
-      !                                     ! ========= !
-      !                                     !  Job end  !
-      !                                     ! ========= !
-
-      IF(lwp) WRITE(numout,cform_aaa)       ! Flag AAAAAAA
-
-      IF( nstop /= 0 ) THEN                 ! error print
-      IF(lwp) WRITE(numout,cform_err)
-      IF(lwp) WRITE(numout,*) nstop, ' error have been found' 
-      ENDIF
-
-      IF( lk_mpp )   CALL mppstop                          ! Close all files (mpp)
-
-   END SUBROUTINE opa_model
+   END SUBROUTINE opa_init
+
+   SUBROUTINE opa_closefile
+      !!----------------------------------------------------------------------
+      !!                     ***  ROUTINE opa_closefile  ***
+      !!
+      !! ** Purpose :   Close the files
+      !!
+      !! ** Method  :
+      !!
+      !! History :
+      !!   9.0  !  05-01  (O. Le Galloudec)  Original code
+      !!----------------------------------------------------------------------
+      !!----------------------------------------------------------------------
+
+      IF ( lk_mpp ) CALL mppsync
+
+      ! 1. Unit close
+      ! -------------
+
+      CLOSE( numnam )           ! namelist
+      CLOSE( numout )           ! standard model output file
+
+      IF(lwp) CLOSE( numstp )   ! time-step file
+
+   END SUBROUTINE opa_closefile
 
    !!======================================================================