Changeset 497 for trunk/NEMO/OFF_SRC/DOM/domstp.F90
- Timestamp:
- 2006-09-12T13:03:53+02:00 (18 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/NEMO/OFF_SRC/DOM/domstp.F90
r343 r497 7 7 !!---------------------------------------------------------------------- 8 8 !! dom_stp : ocean time domain initialization 9 !!---------------------------------------------------------------------- 10 !! History : 11 !! ! 90-10 (O. Marti) Original code 12 !! ! 96-01 (G. Madec) terrain following coordinates 13 !! 8.5 ! 02-08 (G. Madec) F90: Free form and module 9 14 !!---------------------------------------------------------------------- 10 15 !! * Modules used … … 22 27 # include "domzgr_substitute.h90" 23 28 !!---------------------------------------------------------------------- 24 !! OPA 9.0 , LOCEAN-IPSL (2005)25 !! $Header$26 !! This software is governed by the CeCILL licence see modipsl/doc/NEMO_CeCILL.txt29 !! OPA 9.0 , LOCEAN-IPSL (2005) 30 !! $Header$ 31 !! Software governed by the CeCILL licence (modipsl/doc/NEMO_CeCILL.txt) 27 32 !!---------------------------------------------------------------------- 28 33 … … 56 61 !! References : 57 62 !! Bryan, K., 1984, J. Phys. Oceanogr., 14, 666-673. 58 !!59 !! History :60 !! ! 90-10 (O. Marti) Original code61 !! ! 96-01 (G. Madec) terrain following coordinates62 !! 8.5 ! 02-08 (G. Madec) F90: Free form and module63 63 !!---------------------------------------------------------------------- 64 !! * Local declarations65 64 INTEGER :: jk ! dummy loop indice 66 65 !!---------------------------------------------------------------------- … … 77 76 atfp1 = 1. - 2. * atfp 78 77 79 80 78 SELECT CASE ( nacc ) 81 79 … … 90 88 IF(lwp) WRITE(numout,*)' accelerating the convergence' 91 89 IF(lwp) WRITE(numout,*)' dynamics time step = ', rdt/3600., ' hours' 92 #if defined key_s_coord 93 IF( rdtmin /= rdtmax ) THEN 94 IF(lwp) WRITE(numout,cform_err) 95 IF(lwp) WRITE(numout,*)' depth dependent acceleration of & 96 &convergence not implemented in s-coordinates' 97 nstop = nstop + 1 98 ENDIF 99 #endif 100 #if defined key_partial_steps 101 IF( rdtmin /= rdtmax ) THEN 102 IF(lwp) WRITE(numout,cform_err) 103 IF(lwp) WRITE(numout,*)' depth dependent acceleration of & 104 &convergence not implemented for partial steps case' 105 nstop = nstop + 1 106 ENDIF 107 #endif 90 IF( ln_sco .AND. rdtmin /= rdtmax ) & 91 & CALL ctl_stop ( ' depth dependent acceleration of convergence not implemented in s-coordinates' ) 108 92 IF(lwp) WRITE(numout,*)' tracers time step : dt (hours) level' 109 93 110 94 DO jk = 1, jpk 111 IF( fsdept(1,1,jk) <= rdth ) rdttra(jk) = rdtmin112 IF( fsdept(1,1,jk) > rdth ) THEN95 IF( gdept_0(jk) <= rdth ) rdttra(jk) = rdtmin 96 IF( gdept_0(jk) > rdth ) THEN 113 97 rdttra(jk) = rdtmin + ( rdtmax - rdtmin ) & 114 * ( EXP( ( fsdept(1,1,jk ) - rdth ) / rdth ) - 1. ) &115 / ( EXP( ( fsdept(1,1,jpk) - rdth ) / rdth ) - 1. )98 * ( EXP( ( gdept_0(jk ) - rdth ) / rdth ) - 1. ) & 99 / ( EXP( ( gdept_0(jpk) - rdth ) / rdth ) - 1. ) 116 100 ENDIF 117 IF(lwp) WRITE(numout, 9200) rdttra(jk)/3600., jk101 IF(lwp) WRITE(numout,"(36x,f5.2,5x,i3)") rdttra(jk)/3600., jk 118 102 END DO 119 9200 FORMAT(36x,f5.2,' ',i3)120 103 121 104 CASE DEFAULT ! E R R O R 122 105 123 IF(lwp) WRITE(numout,cform_err) 124 IF(lwp) WRITE(numout,*) ' nacc value e r r o r, nacc= ',nacc 125 IF(lwp) WRITE(numout,*) ' we stop' 126 nstop = nstop + 1 106 WRITE(ctmp1,*) ' nacc value e r r o r, nacc= ',nacc 107 CALL ctl_stop( ctmp1 ) 127 108 128 109 END SELECT
Note: See TracChangeset
for help on using the changeset viewer.