Changeset 454 for trunk/NEMO/OPA_SRC/DOM/domstp.F90

Timestamp:

2006-05-10T18:47:31+02:00 (18 years ago)

Author:

opalod

Message:

nemo_v1_update_047:RB: re-organization of coordinate definition, scale factors are now 3d by default, include file for partial steps has been removed

File:

• trunk/NEMO/OPA_SRC/DOM/domstp.F90 (modified) (6 diffs)

trunk/NEMO/OPA_SRC/DOM/domstp.F90

@@ -7 +7 @@
    !!----------------------------------------------------------------------
    !!   dom_stp        : ocean time domain initialization
+   !!----------------------------------------------------------------------
+   !! History :
+   !!        !  90-10  (O. Marti)  Original code
+   !!        !  96-01  (G. Madec)  terrain following coordinates
+   !!   8.5  !  02-08  (G. Madec)  F90: Free form and module
    !!----------------------------------------------------------------------
    !! * Modules used
@@ -24 +29 @@
    !!   OPA 9.0 , LOCEAN-IPSL (2005) 
    !! $Header$ 
-   !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt 
+   !! Software governed by the CeCILL licence (modipsl/doc/NEMO_CeCILL.txt)
    !!----------------------------------------------------------------------
 
@@ -56 +61 @@
       !! References :
       !!      Bryan, K., 1984, J. Phys. Oceanogr., 14, 666-673.
-      !!
-      !! History :
-      !!        !  90-10  (O. Marti)  Original code
-      !!        !  96-01  (G. Madec)  terrain following coordinates
-      !!   8.5  !  02-08  (G. Madec)  F90: Free form and module
       !!----------------------------------------------------------------------
-      !! * Local declarations
       INTEGER ::   jk              ! dummy loop indice
       !!----------------------------------------------------------------------
@@ -77 +76 @@
       atfp1 = 1. - 2. * atfp
 
-
       SELECT CASE ( nacc )
 
@@ -90 +88 @@
             IF(lwp) WRITE(numout,*)'               accelerating the convergence'
             IF(lwp) WRITE(numout,*)'               dynamics time step = ', rdt/3600., ' hours'
-#if defined key_s_coord
-            IF( rdtmin /= rdtmax ) THEN
+            IF( ln_sco .AND. rdtmin /= rdtmax ) THEN
                IF(lwp) WRITE(numout,cform_err)
                IF(lwp) WRITE(numout,*)' depth dependent acceleration of &
@@ -97 +94 @@
                nstop = nstop + 1
             ENDIF
-#endif
-#if defined key_partial_steps
-            IF( rdtmin /= rdtmax ) THEN
-               IF(lwp) WRITE(numout,cform_err)
-               IF(lwp) WRITE(numout,*)' depth dependent acceleration of &
-                                      &convergence not implemented for partial steps case'
-               nstop = nstop + 1
-            ENDIF
-#endif
             IF(lwp) WRITE(numout,*)'         tracers   time step :  dt (hours)  level'
 
             DO jk = 1, jpk
-               IF( fsdept(1,1,jk) <= rdth ) rdttra(jk) = rdtmin
-               IF( fsdept(1,1,jk) >  rdth ) THEN
+               IF( gdept_0(jk) <= rdth ) rdttra(jk) = rdtmin
+               IF( gdept_0(jk) >  rdth ) THEN
                   rdttra(jk) = rdtmin + ( rdtmax - rdtmin )   &
-                                      * ( EXP( ( fsdept(1,1,jk ) - rdth ) / rdth ) - 1. )   &
-                                      / ( EXP( ( fsdept(1,1,jpk) - rdth ) / rdth ) - 1. )
+                                      * ( EXP( ( gdept_0(jk ) - rdth ) / rdth ) - 1. )   &
+                                      / ( EXP( ( gdept_0(jpk) - rdth ) / rdth ) - 1. )
                ENDIF
-               IF(lwp) WRITE(numout,9200) rdttra(jk)/3600., jk
+               IF(lwp) WRITE(numout,"(36x,f5.2,5x,i3)") rdttra(jk)/3600., jk
             END DO  
- 9200       FORMAT(36x,f5.2,'     ',i3)
 
          CASE DEFAULT              ! E R R O R