/[MITgcm]/MITgcm/model/src/config_summary.F
ViewVC logotype

Diff of /MITgcm/model/src/config_summary.F

Parent Directory Parent Directory | Revision Log Revision Log | View Revision Graph Revision Graph | View Patch Patch

revision 1.69 by baylor, Thu Mar 10 03:45:11 2005 UTC revision 1.129 by jmc, Thu May 12 22:32:54 2011 UTC
# Line 27  C     !USES: Line 27  C     !USES:
27  #include "PARAMS.h"  #include "PARAMS.h"
28  #include "EOS.h"  #include "EOS.h"
29  #include "GRID.h"  #include "GRID.h"
 #include "DYNVARS.h"  
30  #ifdef ALLOW_MNC  #ifdef ALLOW_MNC
31  #include "MNC_PARAMS.h"  #include "MNC_PARAMS.h"
32  #endif  #endif
# Line 37  C     myThid ::  Number of this instance Line 36  C     myThid ::  Number of this instance
36        INTEGER myThid        INTEGER myThid
37  CEOP  CEOP
38    
39    C     !FUNCTIONS:
40          INTEGER  ILNBLNK
41          EXTERNAL ILNBLNK
42    
43  C     !LOCAL VARIABLES:  C     !LOCAL VARIABLES:
44  C     msgBuf :: Temp. for building output string.  C     msgBuf :: Temp. for building output string.
45  C     I,J,K  :: Loop counters.  C     rUnits :: vertical coordinate units
46  C     bi,bj  :: Tile loop counters.  C     ioUnit :: Temp. for fortran I/O unit
47  C     xcoord :: Temps. for building lists of values for uni-dimensionally  C     i, k   :: Loop counters.
 C     ycoord :: varying parameters.  
 C     zcoord ::  
48        CHARACTER*(MAX_LEN_MBUF) msgBuf        CHARACTER*(MAX_LEN_MBUF) msgBuf
49        INTEGER                  I,J,K        CHARACTER*2 rUnits
50        INTEGER                  bi, bj        CHARACTER*10 endList
51        _RL                     xcoord(Nx)        INTEGER ioUnit
52        _RL                     ycoord(Ny)        INTEGER i, k
53        _RL                     rcoord(Nr+1)        _RL     bufRL(Nr+1)
54          INTEGER buffI(1)
55        INTEGER coordLine        INTEGER coordLine
56        INTEGER tileLine        INTEGER tileLine
57    
# Line 57  C     zcoord :: Line 59  C     zcoord ::
59        _BARRIER        _BARRIER
60        _BEGIN_MASTER(myThid)        _BEGIN_MASTER(myThid)
61    
62          ioUnit = standardMessageUnit
63          rUnits = ' m'
64          endList = '    ;     '
65          IF ( usingPCoords ) rUnits = 'Pa'
66    
67        WRITE(msgBuf,'(A)')        WRITE(msgBuf,'(A)')
68       &'// ======================================================='       &'// ======================================================='
69        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,        CALL PRINT_MESSAGE( msgBuf, ioUnit, SQUEEZE_RIGHT, myThid )
      &                    SQUEEZE_RIGHT , 1)  
70        WRITE(msgBuf,'(A)') '// Model configuration'        WRITE(msgBuf,'(A)') '// Model configuration'
71        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,        CALL PRINT_MESSAGE( msgBuf, ioUnit, SQUEEZE_RIGHT, myThid )
      &                    SQUEEZE_RIGHT , 1)  
72        WRITE(msgBuf,'(A)')        WRITE(msgBuf,'(A)')
73       &'// ======================================================='       &'// ======================================================='
74        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,        CALL PRINT_MESSAGE( msgBuf, ioUnit, SQUEEZE_RIGHT, myThid )
      &  SQUEEZE_RIGHT , 1)  
75    
76        WRITE(msgBuf,'(A)') '//  '        WRITE(msgBuf,'(A)') '//  '
77        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,        CALL PRINT_MESSAGE( msgBuf, ioUnit, SQUEEZE_RIGHT, myThid )
78       &                    SQUEEZE_RIGHT , 1)        WRITE(msgBuf,'(A)')
       WRITE(msgBuf,'(A)')  
79       & '// "Physical" paramters ( PARM01 in namelist ) '       & '// "Physical" paramters ( PARM01 in namelist ) '
80        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,        CALL PRINT_MESSAGE( msgBuf, ioUnit, SQUEEZE_RIGHT, myThid )
      &                    SQUEEZE_RIGHT , 1)  
81        WRITE(msgBuf,'(A)') '//  '        WRITE(msgBuf,'(A)') '//  '
82        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,        CALL PRINT_MESSAGE( msgBuf, ioUnit, SQUEEZE_RIGHT, myThid )
83       &                    SQUEEZE_RIGHT , 1)        CALL WRITE_0D_C( buoyancyRelation, -1, INDEX_NONE,
84        WRITE(msgBuf,'(A,A40)') 'buoyancyRelation = ', buoyancyRelation       & 'buoyancyRelation =', ' /* Type of relation to get Buoyancy */')
       CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,  
      &                    SQUEEZE_RIGHT , 1)  
85        CALL WRITE_0D_L( fluidIsAir,   INDEX_NONE,        CALL WRITE_0D_L( fluidIsAir,   INDEX_NONE,
86       & 'fluidIsAir  =', '  /* fluid major constituent is Air */')       & 'fluidIsAir   =', '  /* fluid major constituent is Air */')
87        CALL WRITE_0D_L( fluidIsWater, INDEX_NONE,        CALL WRITE_0D_L( fluidIsWater, INDEX_NONE,
88       & 'fluidIsWater=', '  /* fuild major constituent is Water */')       & 'fluidIsWater =', '  /* fluid major constituent is Water */')
89        CALL WRITE_0D_L( usingPCoords, INDEX_NONE,        CALL WRITE_0D_L( usingPCoords, INDEX_NONE,
90       & 'usingPCoords =', '  /* use p (or p*) vertical coordinate */')       & 'usingPCoords =', '  /* use p (or p*) vertical coordinate */')
91        CALL WRITE_0D_L( usingZCoords, INDEX_NONE,        CALL WRITE_0D_L( usingZCoords, INDEX_NONE,
92       & 'usingZCoords =', '  /* use z (or z*) vertical coordinate */')       & 'usingZCoords =', '  /* use z (or z*) vertical coordinate */')
93        CALL WRITE_1D_R8( tRef, Nr, INDEX_K,'tRef =',        CALL WRITE_1D_RL( tRef, Nr, INDEX_K, 'tRef =',
94       &'   /* Reference temperature profile ( oC or oK ) */')       & '   /* Reference temperature profile ( oC or K ) */')
95        CALL WRITE_1D_R8( sRef, Nr, INDEX_K,'sRef =',        CALL WRITE_1D_RL( sRef, Nr, INDEX_K, 'sRef =',
96       &'   /* Reference salinity profile ( ppt ) */')       & '   /* Reference salinity profile ( psu ) */')
97        CALL WRITE_0D_R8( viscAh, INDEX_NONE,'viscAh =',        CALL WRITE_0D_RL( viscAh, INDEX_NONE, 'viscAh  =  ',
98       &'   /* Lateral eddy viscosity ( m^2/s ) */')       & ' /* Lateral eddy viscosity ( m^2/s ) */')
99        IF ( viscAhD.NE.viscAh )        IF ( viscAhD.NE.viscAh )
100       & CALL WRITE_0D_R8( viscAhD, INDEX_NONE,'viscAhD =',       & CALL WRITE_0D_RL( viscAhD, INDEX_NONE, 'viscAhD =  ',
101       & '  /* Lateral eddy viscosity (Divergence)( m^2/s ) */')       & ' /* Lateral eddy viscosity (Divergence)( m^2/s ) */')
102        IF ( viscAhZ.NE.viscAh )        IF ( viscAhZ.NE.viscAh )
103       & CALL WRITE_0D_R8( viscAhZ, INDEX_NONE,'viscAhZ =',       & CALL WRITE_0D_RL( viscAhZ, INDEX_NONE, 'viscAhZ =  ',
104       & '  /* Lateral eddy viscosity (Vorticity) ( m^2/s ) */')       & ' /* Lateral eddy viscosity (Vorticity) ( m^2/s ) */')
105        CALL WRITE_0D_R8( viscAhMax, INDEX_NONE,'viscAhMax =',        CALL WRITE_0D_RL( viscAhMax, INDEX_NONE, 'viscAhMax =',
106       &'   /* Maximum lateral eddy viscosity ( m^2/s ) */')       & ' /* Maximum lateral eddy viscosity ( m^2/s ) */')
107        CALL WRITE_0D_R8( viscAhGrid, INDEX_NONE,'viscAhGrid =',        CALL WRITE_0D_RL( viscAhGrid, INDEX_NONE, 'viscAhGrid =',
108       &'   /* Grid dependent lateral eddy viscosity ( non-dim. ) */')       & ' /* Grid dependent lateral eddy viscosity ( non-dim. ) */')
109        CALL WRITE_0D_L( useFullLeith, INDEX_NONE,        CALL WRITE_0D_L( useFullLeith, INDEX_NONE, 'useFullLeith =',
110       &'useFullLeith =',       & ' /* Use Full Form of Leith Viscosity on/off flag*/')
111       &'   /* Use Full Form of Leith Viscosity on/off flag*/')        CALL WRITE_0D_L( useStrainTensionVisc, INDEX_NONE,
112        CALL WRITE_0D_R8( viscC2leith, INDEX_NONE,'viscC2leith =',       & 'useStrainTensionVisc=',
113       &'   /* Leith harmonic viscosity factor ( on grad(vort), non-dim. ) */')       & ' /* Use StrainTension Form of Viscous Operator flag*/')
114        CALL WRITE_0D_R8( viscC2leithD, INDEX_NONE,'viscC2leithD =',        CALL WRITE_0D_L( useAreaViscLength, INDEX_NONE,
115       &'   /* Leith harmonic viscosity factor ( on grad(div), non-dim. ) */')       & 'useAreaViscLength =',
116        CALL WRITE_0D_R8( viscC2smag, INDEX_NONE,'viscC2smag =',       & ' /* Use area for visc length instead of geom. mean*/')
117       &'   /* Smagorinsky harmonic viscosity factor (non-dim.) */')        CALL WRITE_0D_RL( viscC2leith, INDEX_NONE, 'viscC2leith =',
118        CALL WRITE_0D_R8( viscA4, INDEX_NONE,'viscA4 =',       & ' /* Leith harmonic visc. factor (on grad(vort),non-dim.) */')
119       &'   /* Lateral biharmonic viscosity ( m^4/s ) */')        CALL WRITE_0D_RL( viscC2leithD, INDEX_NONE, 'viscC2leithD =',
120        IF ( viscA4D.NE.viscA4 )       & ' /* Leith harmonic viscosity factor (on grad(div),non-dim.)*/')
121       & CALL WRITE_0D_R8( viscA4D, INDEX_NONE,'viscA4D =',        CALL WRITE_0D_RL( viscC2smag, INDEX_NONE, 'viscC2smag =',
122       & '  /* Lateral biharmonic viscosity (Divergence)( m^4/s ) */')       & ' /* Smagorinsky harmonic viscosity factor (non-dim.) */')
123        IF ( viscA4Z.NE.viscA4 )        CALL WRITE_0D_RL( viscA4, INDEX_NONE, 'viscA4  =  ',
124       & CALL WRITE_0D_R8( viscA4Z, INDEX_NONE,'viscA4Z =',       & ' /* Lateral biharmonic viscosity ( m^4/s ) */')
125       & '  /* Lateral biharmonic viscosity (Vorticity) ( m^4/s ) */')        IF ( viscA4D.NE.viscA4 )
126        CALL WRITE_0D_R8( viscA4Max, INDEX_NONE,'viscA4Max =',       & CALL WRITE_0D_RL( viscA4D, INDEX_NONE, 'viscA4D =  ',
127       &'   /* Maximum biharmonic viscosity ( m^2/s ) */')       & ' /* Lateral biharmonic viscosity (Divergence)( m^4/s ) */')
128        CALL WRITE_0D_R8( viscA4Grid, INDEX_NONE,'viscA4Grid =',        IF ( viscA4Z.NE.viscA4 )
129       &'   /* Grid dependent biharmonic viscosity ( non-dim. ) */')       & CALL WRITE_0D_RL( viscA4Z, INDEX_NONE, 'viscA4Z =  ',
130        CALL WRITE_0D_R8( viscC4leith, INDEX_NONE,'viscC4leith =',       & ' /* Lateral biharmonic viscosity (Vorticity) ( m^4/s ) */')
131       &'   /* Leith biharmonic viscosity factor ( on grad(vort), non-dom. ) */')        CALL WRITE_0D_RL( viscA4Max, INDEX_NONE, 'viscA4Max =',
132        CALL WRITE_0D_R8( viscC4leithD, INDEX_NONE,'viscC4leithD =',       & ' /* Maximum biharmonic viscosity ( m^2/s ) */')
133       &'   /* Leith biharmonic viscosity factor ( on grad(div), non-dom. ) */')        CALL WRITE_0D_RL( viscA4Grid, INDEX_NONE, 'viscA4Grid =',
134         & ' /* Grid dependent biharmonic viscosity ( non-dim. ) */')
135          CALL WRITE_0D_RL( viscC4leith, INDEX_NONE,'viscC4leith =',
136         & ' /* Leith biharm viscosity factor (on grad(vort), non-dim.)*/')
137          CALL WRITE_0D_RL( viscC4leithD, INDEX_NONE,'viscC4leithD =',
138         & ' /* Leith biharm viscosity factor (on grad(div), non-dim.) */')
139          CALL WRITE_0D_RL( viscC4Smag, INDEX_NONE,'viscC4Smag =',
140         & ' /* Smagorinsky biharm viscosity factor (non-dim) */')
141        CALL WRITE_0D_L( no_slip_sides, INDEX_NONE,        CALL WRITE_0D_L( no_slip_sides, INDEX_NONE,
142       & 'no_slip_sides =', '  /* Viscous BCs: No-slip sides */')       & 'no_slip_sides =', '  /* Viscous BCs: No-slip sides */')
143        CALL WRITE_0D_R8( viscAr,  INDEX_NONE,'viscAr =',        CALL WRITE_0D_RL( sideDragFactor, INDEX_NONE, 'sideDragFactor =',
144       &'   /* Vertical eddy viscosity ( units of r^2/s ) */')       & ' /* side-drag scaling factor (non-dim) */')
145          CALL WRITE_1D_RL( viscArNr, Nr, INDEX_K, 'viscArNr =',
146         &  ' /* vertical profile of vertical viscosity ('
147         &  //rUnits//'^2/s )*/')
148        CALL WRITE_0D_L( no_slip_bottom, INDEX_NONE,        CALL WRITE_0D_L( no_slip_bottom, INDEX_NONE,
149       & 'no_slip_bottom =', '  /* Viscous BCs: No-slip bottom */')       & 'no_slip_bottom =', '  /* Viscous BCs: No-slip bottom */')
150        CALL WRITE_0D_R8( diffKhT, INDEX_NONE,'diffKhT =',        CALL WRITE_0D_RL( bottomDragLinear, INDEX_NONE,
151         & 'bottomDragLinear =',
152         & ' /* linear bottom-drag coefficient ( m/s ) */')
153          CALL WRITE_0D_RL( bottomDragQuadratic, INDEX_NONE,
154         & 'bottomDragQuadratic =',
155         & ' /* quadratic bottom-drag coefficient (-) */')
156          CALL WRITE_0D_RL( diffKhT, INDEX_NONE,'diffKhT =',
157       &'   /* Laplacian diffusion of heat laterally ( m^2/s ) */')       &'   /* Laplacian diffusion of heat laterally ( m^2/s ) */')
158        CALL WRITE_0D_R8( diffK4T, INDEX_NONE,'diffK4T =',        CALL WRITE_0D_RL( diffK4T, INDEX_NONE,'diffK4T =',
159       &'   /* Bihaarmonic diffusion of heat laterally ( m^4/s ) */')       &'   /* Biharmonic diffusion of heat laterally ( m^4/s ) */')
160        CALL WRITE_0D_R8( diffKhS, INDEX_NONE,'diffKhS =',        CALL WRITE_0D_RL( diffKhS, INDEX_NONE,'diffKhS =',
161       &'   /* Laplacian diffusion of salt laterally ( m^2/s ) */')       &'   /* Laplacian diffusion of salt laterally ( m^2/s ) */')
162        CALL WRITE_0D_R8( diffK4S, INDEX_NONE,'diffK4S =',        CALL WRITE_0D_RL( diffK4S, INDEX_NONE,'diffK4S =',
163       &'   /* Bihaarmonic diffusion of salt laterally ( m^4/s ) */')       &'   /* Biharmonic diffusion of salt laterally ( m^4/s ) */')
164        CALL WRITE_1D_R8( diffKrNrT, Nr, INDEX_K,'diffKrNrT =',        CALL WRITE_1D_RL( diffKrNrT, Nr, INDEX_K, 'diffKrNrT =',
165       & ' /* vertical profile of vertical diffusion of Temp ( m^2/s )*/')       &  ' /* vertical profile of vertical diffusion of Temp ('
166        CALL WRITE_1D_R8( diffKrNrS, Nr, INDEX_K,'diffKrNrS =',       &  //rUnits//'^2/s )*/')
167       & ' /* vertical profile of vertical diffusion of Salt ( m^2/s )*/')        CALL WRITE_1D_RL( diffKrNrS, Nr, INDEX_K, 'diffKrNrS =',
168        CALL WRITE_0D_R8( diffKrBL79surf, INDEX_NONE,'diffKrBL79surf =',       &  ' /* vertical profile of vertical diffusion of Salt ('
169       &'   /* Surface diffusion for Bryan and Lewis 1979 ( m^2/s ) */')       &  //rUnits//'^2/s )*/')
170        CALL WRITE_0D_R8( diffKrBL79deep, INDEX_NONE,'diffKrBL79deep =',        CALL WRITE_0D_RL( diffKrBL79surf, INDEX_NONE,'diffKrBL79surf =',
171       &'   /* Deep diffusion for Bryan and Lewis 1979 ( m^2/s ) */')       &  ' /* Surface diffusion for Bryan and Lewis 79 ( m^2/s ) */')
172        CALL WRITE_0D_R8( diffKrBL79scl, INDEX_NONE,'diffKrBL79scl =',        CALL WRITE_0D_RL( diffKrBL79deep, INDEX_NONE,'diffKrBL79deep =',
173       &'   /* Depth scale for Bryan and Lewis 1979 ( m ) */')       &  ' /* Deep diffusion for Bryan and Lewis 1979 ( m^2/s ) */')
174        CALL WRITE_0D_R8( diffKrBL79Ho, INDEX_NONE,'diffKrBL79Ho =',        CALL WRITE_0D_RL( diffKrBL79scl, INDEX_NONE,'diffKrBL79scl =',
175       &'   /* Turning depth for Bryan and Lewis 1979 ( m ) */')       &  ' /* Depth scale for Bryan and Lewis 1979 ( m ) */')
176        WRITE(msgBuf,'(2A)') ' Equation of State : eosType = ', eosType        CALL WRITE_0D_RL( diffKrBL79Ho, INDEX_NONE,'diffKrBL79Ho =',
177        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,       &  ' /* Turning depth for Bryan and Lewis 1979 ( m ) */')
178       &                    SQUEEZE_RIGHT , 1)        CALL WRITE_0D_RL( ivdc_kappa, INDEX_NONE,'ivdc_kappa =',
179        CALL WRITE_0D_R8( tAlpha,  INDEX_NONE,'tAlpha =',       &  ' /* Implicit Vertical Diffusivity for Convection ('
180       &'   /* Linear EOS thermal expansion coefficient ( 1/degree ) */')       &  //rUnits//'^2/s) */')
181        CALL WRITE_0D_R8( sBeta,   INDEX_NONE,'sBeta =',        CALL WRITE_0D_RL( hMixCriteria, INDEX_NONE,'hMixCriteria=',
182       &'   /* Linear EOS haline contraction coefficient ( 1/ppt ) */')       & '  /* Criteria for mixed-layer diagnostic */')
183          CALL WRITE_0D_RL( dRhoSmall, INDEX_NONE,'dRhoSmall =',
184         & '  /* Parameter for mixed-layer diagnostic */')
185          CALL WRITE_0D_RL( hMixSmooth, INDEX_NONE,'hMixSmooth=',
186         & '  /* Smoothing parameter for mixed-layer diagnostic */')
187          CALL WRITE_0D_C( eosType, 0, INDEX_NONE, 'eosType =',
188         & '  /* Type of Equation of State */')
189          CALL WRITE_0D_RL( tAlpha,  INDEX_NONE,'tAlpha =',
190         &'   /* Linear EOS thermal expansion coefficient ( 1/oC ) */')
191          CALL WRITE_0D_RL( sBeta,   INDEX_NONE,'sBeta =',
192         &'   /* Linear EOS haline contraction coefficient ( 1/psu ) */')
193        IF ( eosType .EQ. 'POLY3' ) THEN        IF ( eosType .EQ. 'POLY3' ) THEN
194          WRITE(msgBuf,'(A)')          WRITE(msgBuf,'(A)')
195       &   '// Polynomial EQS parameters ( from POLY3.COEFFS ) '       &   '// Polynomial EQS parameters ( from POLY3.COEFFS ) '
196          DO K = 1, Nr          DO k = 1, Nr
197           WRITE(msgBuf,'(I3,13F8.3)')           WRITE(msgBuf,'(I3,13F8.3)')
198       &   K,eosRefT(K),eosRefS(K),eosSig0(K), (eosC(I,K),I=1,9)       &   k,eosRefT(k),eosRefS(k),eosSig0(k), (eosC(i,k),i=1,9)
199           CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,           CALL PRINT_MESSAGE( msgBuf, ioUnit, SQUEEZE_RIGHT, myThid )
      &                    SQUEEZE_RIGHT , 1)  
200          ENDDO          ENDDO
201        ENDIF        ENDIF
202        IF ( fluidIsAir ) THEN        IF ( fluidIsAir ) THEN
203         CALL WRITE_0D_R8( atm_Rd, INDEX_NONE, 'atm_Rd =',         CALL WRITE_0D_RL( atm_Rd, INDEX_NONE, 'atm_Rd =',
204       & '  /* gas constant for dry air ( J/kg/K ) */')       & '  /* gas constant for dry air ( J/kg/K ) */')
205        CALL WRITE_0D_R8( atm_Cp, INDEX_NONE, 'atm_Cp =',        CALL WRITE_0D_RL( atm_Cp, INDEX_NONE, 'atm_Cp =',
206       & '  /* specific heat (Cp) of dry air ( J/kg/K ) */')       & '  /* specific heat (Cp) of dry air ( J/kg/K ) */')
207        CALL WRITE_0D_R8( atm_kappa, INDEX_NONE, 'atm_kappa =',        CALL WRITE_0D_RL( atm_kappa, INDEX_NONE, 'atm_kappa =',
208       & '  /* kappa (=Rd/Cp ) of dry air */')       & '  /* kappa (=Rd/Cp ) of dry air */')
209         CALL WRITE_0D_R8( atm_Rq, INDEX_NONE, 'atm_Rq =',         CALL WRITE_0D_RL( atm_Rq, INDEX_NONE, 'atm_Rq =',
210       &  ' /* water vap. specific vol. anomaly relative to dry air */')       &  ' /* water vap. specific vol. anomaly relative to dry air */')
211        CALL WRITE_0D_R8( atm_Po, INDEX_NONE, 'atm_Po =',        CALL WRITE_0D_RL( atm_Po, INDEX_NONE, 'atm_Po =',
212       & '  /* standard reference pressure ( Pa ) */')       & '  /* standard reference pressure ( Pa ) */')
213        CALL WRITE_0D_I( integr_GeoPot, INDEX_NONE, 'integr_GeoPot =',        CALL WRITE_0D_I( integr_GeoPot, INDEX_NONE, 'integr_GeoPot =',
214       & '  /* select how the geopotential is integrated */')       & '  /* select how the geopotential is integrated */')
215        CALL WRITE_0D_I( selectFindRoSurf, INDEX_NONE,        CALL WRITE_0D_I( selectFindRoSurf, INDEX_NONE,
216       & 'selectFindRoSurf=',       & 'selectFindRoSurf=',
217       & '  /* select how Surf.Ref. pressure is defined */')       & '  /* select how Surf.Ref. pressure is defined */')
218        ENDIF        ENDIF
219        CALL WRITE_0D_R8( rhonil,  INDEX_NONE,'rhonil =',        CALL WRITE_0D_RL( rhonil,  INDEX_NONE,'rhonil =',
220       &'   /* Reference density ( kg/m^3 ) */')       &'   /* Reference density ( kg/m^3 ) */')
221        CALL WRITE_0D_R8( rhoConst, INDEX_NONE,'rhoConst =',        CALL WRITE_0D_RL( rhoConst, INDEX_NONE,'rhoConst =',
222       &'   /* Reference density ( kg/m^3 ) */')       &'   /* Reference density ( kg/m^3 ) */')
223        CALL WRITE_0D_R8( rhoConstFresh, INDEX_NONE,'rhoConstFresh =',        CALL WRITE_1D_RL( rhoFacC, Nr,   INDEX_K, 'rhoFacC = ',
224         &  ' /* normalized Reference density @ cell-Center (-) */')
225          CALL WRITE_1D_RL( rhoFacF, Nr+1, INDEX_K, 'rhoFacF = ',
226         &  ' /* normalized Reference density @ W-Interface (-) */')
227          CALL WRITE_0D_RL( rhoConstFresh, INDEX_NONE,'rhoConstFresh =',
228       &'   /* Reference density ( kg/m^3 ) */')       &'   /* Reference density ( kg/m^3 ) */')
229        CALL WRITE_0D_R8( gravity, INDEX_NONE,'gravity =',        CALL WRITE_0D_RL( gravity, INDEX_NONE,'gravity =',
230       &'   /* Gravitational acceleration ( m/s^2 ) */')       &'   /* Gravitational acceleration ( m/s^2 ) */')
231        CALL WRITE_0D_R8( gBaro,   INDEX_NONE,'gBaro =',        CALL WRITE_0D_RL( gBaro,   INDEX_NONE,'gBaro =',
232       &'   /* Barotropic gravity ( m/s^2 ) */')       &'   /* Barotropic gravity ( m/s^2 ) */')
233        CALL WRITE_0D_R8(rotationPeriod,INDEX_NONE,'rotationPeriod =',        CALL WRITE_0D_RL(rotationPeriod,INDEX_NONE,'rotationPeriod =',
234       &'   /* Rotation Period ( s ) */')       &'   /* Rotation Period ( s ) */')
235        CALL WRITE_0D_R8( omega,   INDEX_NONE,'omega =',        CALL WRITE_0D_RL( omega,   INDEX_NONE,'omega =',
236       &'   /* Angular velocity ( rad/s ) */')       &'   /* Angular velocity ( rad/s ) */')
237        CALL WRITE_0D_R8( f0,      INDEX_NONE,'f0 =',        CALL WRITE_0D_RL( f0,      INDEX_NONE,'f0 =',
238       &'   /* Reference coriolis parameter ( 1/s ) */')       &'   /* Reference coriolis parameter ( 1/s ) */')
239        CALL WRITE_0D_R8( beta,    INDEX_NONE,'beta =',        CALL WRITE_0D_RL( beta,    INDEX_NONE,'beta =',
240       &'   /* Beta ( 1/(m.s) ) */')       &'   /* Beta ( 1/(m.s) ) */')
241          CALL WRITE_0D_RL( fPrime,  INDEX_NONE,'fPrime =',
242        CALL WRITE_0D_R8( freeSurfFac, INDEX_NONE,'freeSurfFac =',       &'   /* Second coriolis parameter ( 1/s ) */')
243       &'   /* Implicit free surface factor */')        CALL WRITE_0D_L( rigidLid, INDEX_NONE, 'rigidLid =',
244         &'   /* Rigid lid on/off flag */')
245        CALL WRITE_0D_L( implicitFreeSurface, INDEX_NONE,        CALL WRITE_0D_L( implicitFreeSurface, INDEX_NONE,
246       &                 'implicitFreeSurface =',       &                 'implicitFreeSurface =',
247       &'   /* Implicit free surface on/off flag */')       &'   /* Implicit free surface on/off flag */')
248        CALL WRITE_0D_L( rigidLid, INDEX_NONE,        CALL WRITE_0D_RL( freeSurfFac, INDEX_NONE,'freeSurfFac =',
249       &                 'rigidLid =',       &'   /* Implicit free surface factor */')
250       &'   /* Rigid lid on/off flag */')        CALL WRITE_0D_RL( implicSurfPress, INDEX_NONE,
251        CALL WRITE_0D_R8( implicSurfPress, INDEX_NONE,       & 'implicSurfPress =',
252       &'implicSurfPress =',       & '  /* Surface Pressure implicit factor (0-1)*/')
253       &'   /* Surface Pressure implicit factor (0-1)*/')        CALL WRITE_0D_RL( implicDiv2Dflow, INDEX_NONE,
254        CALL WRITE_0D_R8( implicDiv2Dflow, INDEX_NONE,       & 'implicDiv2Dflow =',
255       &'implicDiv2Dflow =',       & '  /* Barot. Flow Div. implicit factor (0-1)*/')
      &'   /* Barot. Flow Div. implicit factor (0-1)*/')  
256        CALL WRITE_0D_L( exactConserv, INDEX_NONE,        CALL WRITE_0D_L( exactConserv, INDEX_NONE,
257       &'exactConserv =',       & 'exactConserv =',
258       &'   /* Exact Volume Conservation on/off flag*/')       & '  /* Exact Volume Conservation on/off flag*/')
259          CALL WRITE_0D_L( linFSConserveTr, INDEX_NONE,
260         & 'linFSConserveTr =',
261         &  ' /* Tracer correction for Lin Free Surface on/off flag*/')
262        CALL WRITE_0D_L( uniformLin_PhiSurf, INDEX_NONE,        CALL WRITE_0D_L( uniformLin_PhiSurf, INDEX_NONE,
263       &'uniformLin_PhiSurf =',       & 'uniformLin_PhiSurf =',
264       &'   /* use uniform Bo_surf on/off flag*/')       &  ' /* use uniform Bo_surf on/off flag*/')
265        CALL WRITE_0D_I( nonlinFreeSurf, INDEX_NONE,        CALL WRITE_0D_RL( hFacMin, INDEX_NONE, 'hFacMin = ',
266       &'nonlinFreeSurf =',       & '  /* minimum partial cell factor (hFac) */')
267       &'   /* Non-linear Free Surf. options (-1,0,1,2,3)*/')        CALL WRITE_0D_RL( hFacMin, INDEX_NONE, 'hFacMinDr =',
268         &  ' /* minimum partial cell thickness ('//rUnits//') */')
269          WRITE(msgBuf,'(2A)') 'nonlinFreeSurf =',
270         &  ' /* Non-linear Free Surf. options (-1,0,1,2,3)*/'
271          CALL PRINT_MESSAGE( msgBuf, ioUnit, SQUEEZE_RIGHT, myThid )
272          buffI(1) = nonlinFreeSurf
273          CALL PRINT_LIST_I( buffI, 1, 1, INDEX_NONE,
274         &                   .FALSE., .TRUE., ioUnit )
275        WRITE(msgBuf,'(2A)') '     -1,0= Off ; 1,2,3= On,',        WRITE(msgBuf,'(2A)') '     -1,0= Off ; 1,2,3= On,',
276       &  ' 2=+rescale gU,gV, 3=+update cg2d solv.'       &  ' 2=+rescale gU,gV, 3=+update cg2d solv.'
277        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,        CALL PRINT_MESSAGE( msgBuf, ioUnit, SQUEEZE_RIGHT, myThid )
278       &                    SQUEEZE_RIGHT , 1)        CALL PRINT_MESSAGE(endList, ioUnit, SQUEEZE_RIGHT, myThid )
279        CALL WRITE_0D_R8( hFacInf, INDEX_NONE,        CALL WRITE_0D_RL( hFacInf, INDEX_NONE, 'hFacInf = ',
280       &'hFacInf =',       & '  /* lower threshold for hFac (nonlinFreeSurf only)*/')
281       &'   /* lower threshold for hFac (nonlinFreeSurf only)*/')        CALL WRITE_0D_RL( hFacSup, INDEX_NONE, 'hFacSup = ',
282        CALL WRITE_0D_R8( hFacSup, INDEX_NONE,       & '  /* upper threshold for hFac (nonlinFreeSurf only)*/')
      &'hFacSup =',  
      &'   /* upper threshold for hFac (nonlinFreeSurf only)*/')  
283        CALL WRITE_0D_I( select_rStar, INDEX_NONE,        CALL WRITE_0D_I( select_rStar, INDEX_NONE,
284       &'select_rStar =',       & 'select_rStar =',
285       &'   /* r* Coordinate options (not yet implemented)*/')       &  ' /* r* Vertical coord. options (=0 r coord.; >0 uses r*)*/')
286          CALL WRITE_0D_I( selectAddFluid, INDEX_NONE,
287         & 'selectAddFluid =',
288         &  ' /* option for mass source/sink of fluid (=0: off) */')
289        CALL WRITE_0D_L( useRealFreshWaterFlux, INDEX_NONE,        CALL WRITE_0D_L( useRealFreshWaterFlux, INDEX_NONE,
290       &'useRealFreshWaterFlux =',       & 'useRealFreshWaterFlux =',
291       &'   /* Real Fresh Water Flux on/off flag*/')       &  ' /* Real Fresh Water Flux on/off flag*/')
292        IF (useRealFreshWaterFlux .AND. nonlinFreeSurf.GT.0) THEN        CALL WRITE_0D_RL( temp_EvPrRn, INDEX_NONE,
293        CALL WRITE_0D_R8( temp_EvPrRn, INDEX_NONE,       & 'temp_EvPrRn =',
294       &'temp_EvPrRn =',       &  ' /* Temp. of Evap/Prec/R (UNSET=use local T)(oC)*/')
295       &' /* Temp. of Evap/Prec/R (UNSET=use local T)(oC)*/')        CALL WRITE_0D_RL( salt_EvPrRn, INDEX_NONE,
296        CALL WRITE_0D_R8( salt_EvPrRn, INDEX_NONE,       &  'salt_EvPrRn =',
297       &'salt_EvPrRn =',       &  ' /* Salin. of Evap/Prec/R (UNSET=use local S)(psu)*/')
298       &' /* Salin. of Evap/Prec/R (UNSET=use local S)(ppt)*/')        CALL WRITE_0D_RL( temp_addMass, INDEX_NONE,
299        ELSE       & 'temp_addMass =',
300        CALL WRITE_0D_R8( convertFW2Salt, INDEX_NONE,       &  ' /* Temp. of addMass array (UNSET=use local T)(oC)*/')
301       &'convertFW2Salt =',        CALL WRITE_0D_RL( salt_addMass, INDEX_NONE,
302       &' /* convert F.W. Flux to Salt Flux (-1=use local S)(ppt)*/')       & 'salt_addMass =',
303         &  ' /* Salin. of addMass array (UNSET=use local S)(psu)*/')
304          IF ( .NOT.useRealFreshWaterFlux .OR. selectAddFluid.EQ.-1
305         &                                .OR. nonlinFreeSurf.LE.0 ) THEN
306          CALL WRITE_0D_RL( convertFW2Salt, INDEX_NONE,
307         & 'convertFW2Salt =',
308         & ' /* convert F.W. Flux to Salt Flux (-1=use local S)(psu)*/')
309        ENDIF        ENDIF
310    
311          CALL WRITE_0D_L( use3Dsolver, INDEX_NONE,
312         & 'use3Dsolver =', ' /* use 3-D pressure solver on/off flag */')
313        CALL WRITE_0D_L( nonHydrostatic, INDEX_NONE,        CALL WRITE_0D_L( nonHydrostatic, INDEX_NONE,
314       & 'nonHydrostatic =', '  /* Non-Hydrostatic on/off flag */')       & 'nonHydrostatic =', '  /* Non-Hydrostatic on/off flag */')
315          CALL WRITE_0D_RL( nh_Am2, INDEX_NONE, 'nh_Am2 =',
316         & ' /* Non-Hydrostatic terms scaling factor */')
317          CALL WRITE_0D_RL( implicitNHPress, INDEX_NONE,
318         & 'implicitNHPress =',
319         & ' /* Non-Hyd Pressure implicit factor (0-1)*/')
320          CALL WRITE_0D_I( selectNHfreeSurf, INDEX_NONE,
321         & 'selectNHfreeSurf =',
322         & ' /* Non-Hyd (free-)Surface option */')
323          CALL WRITE_0D_L( quasiHydrostatic, INDEX_NONE,
324         & 'quasiHydrostatic =', ' /* Quasi-Hydrostatic on/off flag */')
325        CALL WRITE_0D_L( momStepping,  INDEX_NONE,        CALL WRITE_0D_L( momStepping,  INDEX_NONE,
326       & 'momStepping =', '  /* Momentum equation on/off flag */')       & 'momStepping =', '  /* Momentum equation on/off flag */')
327          CALL WRITE_0D_L( vectorInvariantMomentum,  INDEX_NONE,
328         & 'vectorInvariantMomentum=',
329         & ' /* Vector-Invariant Momentum on/off */')
330        CALL WRITE_0D_L( momAdvection, INDEX_NONE,        CALL WRITE_0D_L( momAdvection, INDEX_NONE,
331       & 'momAdvection =', '  /* Momentum advection on/off flag */')       & 'momAdvection =', '  /* Momentum advection on/off flag */')
332        CALL WRITE_0D_L( momViscosity, INDEX_NONE,        CALL WRITE_0D_L( momViscosity, INDEX_NONE,
333       & 'momViscosity =', '  /* Momentum viscosity on/off flag */')       & 'momViscosity =', '  /* Momentum viscosity on/off flag */')
334        CALL WRITE_0D_L( momImplVertAdv, INDEX_NONE, 'momImplVertAdv =',        CALL WRITE_0D_L( momImplVertAdv, INDEX_NONE, 'momImplVertAdv=',
335       &                '/* Momentum implicit vert. advection on/off*/')       &                ' /* Momentum implicit vert. advection on/off*/')
336        CALL WRITE_0D_L( implicitViscosity, INDEX_NONE,        CALL WRITE_0D_L( implicitViscosity, INDEX_NONE,
337       & 'implicitViscosity =', ' /* Implicit viscosity on/off flag */')       & 'implicitViscosity =', ' /* Implicit viscosity on/off flag */')
338          CALL WRITE_0D_L( metricTerms,  INDEX_NONE, 'metricTerms =',
339         &                '  /* metric-Terms on/off flag */')
340          CALL WRITE_0D_L( useNHMTerms,  INDEX_NONE, 'useNHMTerms =',
341         &              ' /* Non-Hydrostatic Metric-Terms on/off */')
342    c------------
343          WRITE(msgBuf,'(2A)')
344         & 'selectCoriMap =', ' /* Coriolis Map options (0,1,2,3)*/'
345          CALL PRINT_MESSAGE( msgBuf, ioUnit, SQUEEZE_RIGHT, myThid )
346          buffI(1) = selectCoriMap
347          CALL PRINT_LIST_I( buffI, 1, 1, INDEX_NONE,
348         &                   .FALSE., .TRUE., ioUnit )
349          WRITE(msgBuf,'(2A)') '    0= f-Plane ; 1= Beta-Plane ;',
350         &  ' 2= Spherical ; 3= read from file'
351          CALL PRINT_MESSAGE( msgBuf, ioUnit, SQUEEZE_RIGHT, myThid )
352          CALL PRINT_MESSAGE(endList, ioUnit, SQUEEZE_RIGHT, myThid )
353    c------------
354          CALL WRITE_0D_L( use3dCoriolis,  INDEX_NONE,
355         & 'use3dCoriolis =', ' /* 3-D Coriolis on/off flag */')
356        CALL WRITE_0D_L( useCoriolis,  INDEX_NONE,        CALL WRITE_0D_L( useCoriolis,  INDEX_NONE,
357       & 'useCoriolis =', '  /* Coriolis on/off flag */')       & 'useCoriolis =', '  /* Coriolis on/off flag */')
358        CALL WRITE_0D_L( useCDscheme,  INDEX_NONE,        CALL WRITE_0D_L( useCDscheme,  INDEX_NONE,
359       & 'useCDscheme =', '  /* CD scheme on/off flag */')       & 'useCDscheme =', '  /* CD scheme on/off flag */')
360          CALL WRITE_0D_L( useEnergyConservingCoriolis,  INDEX_NONE,
361         & 'useEnergyConservingCoriolis=',
362         & ' /* Flx-Form Coriolis scheme flag */')
363        CALL WRITE_0D_L( useJamartWetPoints,  INDEX_NONE,        CALL WRITE_0D_L( useJamartWetPoints,  INDEX_NONE,
364       & 'useJamartWetPoints=',' /* Coriolis WetPoints method flag */')       & 'useJamartWetPoints=',' /* Coriolis WetPoints method flag */')
365        CALL WRITE_0D_L( useJamartMomAdv,  INDEX_NONE,        CALL WRITE_0D_L( useJamartMomAdv,  INDEX_NONE,
366       & 'useJamartMomAdv=',' /* V.I. Non-linear terms Jamart flag */')       & 'useJamartMomAdv=',' /* V.I Non-linear terms Jamart flag */')
       CALL WRITE_0D_L( SadournyCoriolis,  INDEX_NONE,  
      & 'SadournyCoriolis=',' /* Sadourny Coriolis discr. flag */')  
       CALL WRITE_0D_L( upwindVorticity,  INDEX_NONE,  
      & 'upwindVorticity=',' /* Upwind bias vorticity flag */')  
367        CALL WRITE_0D_L( useAbsVorticity,  INDEX_NONE,        CALL WRITE_0D_L( useAbsVorticity,  INDEX_NONE,
368       & 'useAbsVorticity=',' /* Work with f+zeta in Coriolis */')       & 'useAbsVorticity=',' /* V.I Works with f+zeta in Coriolis */')
369          WRITE(msgBuf,'(2A)') 'selectVortScheme=',
370         &              ' /* V.I Scheme selector for Vorticity-Term */'
371          CALL PRINT_MESSAGE( msgBuf, ioUnit, SQUEEZE_RIGHT, myThid )
372          buffI(1) = selectVortScheme
373          CALL PRINT_LIST_I( buffI, 1, 1, INDEX_NONE,
374         &                   .FALSE., .TRUE., ioUnit )
375          WRITE(msgBuf,'(2A)') '   = 0 : enstrophy (Shallow-Water Eq.)',
376         &                  ' conserving scheme by Sadourny, JAS 75'
377          CALL PRINT_MESSAGE( msgBuf, ioUnit, SQUEEZE_RIGHT, myThid )
378          WRITE(msgBuf,'(2A)') '   = 1 : same as 0 with modified hFac'
379          CALL PRINT_MESSAGE( msgBuf, ioUnit, SQUEEZE_RIGHT, myThid )
380          WRITE(msgBuf,'(2A)') '   = 2 : energy conserving scheme',
381         &         ' (used by Sadourny in JAS 75 paper)'
382          CALL PRINT_MESSAGE( msgBuf, ioUnit, SQUEEZE_RIGHT, myThid )
383          WRITE(msgBuf,'(2A)') '   = 3 : energy (general)',
384         &             ' and enstrophy (2D, nonDiv.) conserving scheme'
385          CALL PRINT_MESSAGE( msgBuf, ioUnit, SQUEEZE_RIGHT, myThid )
386          WRITE(msgBuf,'(2A)') '         from Sadourny',
387         &                     ' (Burridge & Haseler, ECMWF Rep.4, 1977)'
388          CALL PRINT_MESSAGE( msgBuf, ioUnit, SQUEEZE_RIGHT, myThid )
389    c     WRITE(msgBuf,'(2A)') '   = 4 : energy (general)',
390    c    &             ' and enstrophy (2D, nonDiv.) conserving scheme'
391    c     CALL PRINT_MESSAGE( msgBuf, ioUnit, SQUEEZE_RIGHT, myThid )
392    c     WRITE(msgBuf,'(2A)') '         from Arakawa & Lamb, 77'
393    c     CALL PRINT_MESSAGE( msgBuf, ioUnit, SQUEEZE_RIGHT, myThid )
394          CALL PRINT_MESSAGE(endList, ioUnit, SQUEEZE_RIGHT, myThid )
395          CALL WRITE_0D_L( upwindVorticity,  INDEX_NONE,
396         & 'upwindVorticity=',' /* V.I Upwind bias vorticity flag */')
397        CALL WRITE_0D_L( highOrderVorticity,  INDEX_NONE,        CALL WRITE_0D_L( highOrderVorticity,  INDEX_NONE,
398       & 'highOrderVorticity=',' /* High order interp. of vort. flag */')       & 'highOrderVorticity=',' /* V.I High order vort. advect. flag */')
399          CALL WRITE_0D_L( upwindShear,  INDEX_NONE,
400         & 'upwindShear=',' /* V.I Upwind vertical Shear advection flag */')
401          CALL WRITE_0D_I( selectKEscheme, INDEX_NONE,
402         & 'selectKEscheme=',' /* V.I Kinetic Energy scheme selector */')
403        CALL WRITE_0D_L( momForcing,   INDEX_NONE,        CALL WRITE_0D_L( momForcing,   INDEX_NONE,
404       & 'momForcing =', '  /* Momentum forcing on/off flag */')       & 'momForcing =', '  /* Momentum forcing on/off flag */')
405        CALL WRITE_0D_L( momPressureForcing, INDEX_NONE,        CALL WRITE_0D_L( momPressureForcing, INDEX_NONE,
406       & 'momPressureForcing =',         & 'momPressureForcing =',
407       & '  /* Momentum pressure term on/off flag */')       & '  /* Momentum pressure term on/off flag */')
408          CALL WRITE_0D_L( implicitIntGravWave, INDEX_NONE,
409         &  'implicitIntGravWave=',
410         &  ' /* Implicit Internal Gravity Wave flag */')
411        CALL WRITE_0D_L( staggerTimeStep, INDEX_NONE,        CALL WRITE_0D_L( staggerTimeStep, INDEX_NONE,
412       &                 'staggerTimeStep =',       &                 'staggerTimeStep =',
413       &'   /* Stagger time stepping on/off flag */')       &'   /* Stagger time stepping on/off flag */')
414        CALL WRITE_0D_L( multiDimAdvection, INDEX_NONE,        CALL WRITE_0D_L( multiDimAdvection, INDEX_NONE,
415       & 'multiDimAdvection =',       & 'multiDimAdvection =',
416       &'   /* enable/disable Multi-Dim Advection */')       &'   /* enable/disable Multi-Dim Advection */')
417        CALL WRITE_0D_L( useMultiDimAdvec, INDEX_NONE,        CALL WRITE_0D_L( useMultiDimAdvec, INDEX_NONE,
418       & 'useMultiDimAdvec =',       & 'useMultiDimAdvec =',
419       &'   /* Multi-Dim Advection is/is-not used */')       &'   /* Multi-Dim Advection is/is-not used */')
420        CALL WRITE_0D_L( implicitDiffusion, INDEX_NONE,        CALL WRITE_0D_L( implicitDiffusion, INDEX_NONE,
421       & 'implicitDiffusion =','/* Implicit Diffusion on/off flag */')       & 'implicitDiffusion =',' /* Implicit Diffusion on/off flag */')
422        CALL WRITE_0D_L( tempStepping,  INDEX_NONE,        CALL WRITE_0D_L( tempStepping,  INDEX_NONE,
423       & 'tempStepping =', '  /* Temperature equation on/off flag */')       & 'tempStepping =', '  /* Temperature equation on/off flag */')
424        CALL WRITE_0D_L( tempAdvection,  INDEX_NONE,        CALL WRITE_0D_L( tempAdvection,  INDEX_NONE,
425       & 'tempAdvection=', '  /* Temperature advection on/off flag */')       & 'tempAdvection=', '  /* Temperature advection on/off flag */')
426        CALL WRITE_0D_L( tempImplVertAdv,INDEX_NONE,'tempImplVertAdv =',        CALL WRITE_0D_L( tempImplVertAdv,INDEX_NONE,'tempImplVertAdv =',
427       &                '/* Temp. implicit vert. advection on/off */')       &                ' /* Temp. implicit vert. advection on/off */')
428        CALL WRITE_0D_L( tempForcing,  INDEX_NONE,        CALL WRITE_0D_L( tempForcing,  INDEX_NONE,
429       & 'tempForcing  =', '  /* Temperature forcing on/off flag */')       & 'tempForcing  =', '  /* Temperature forcing on/off flag */')
430          CALL WRITE_0D_L( tempIsActiveTr, INDEX_NONE, 'tempIsActiveTr =',
431         & ' /* Temp. is a dynamically Active Tracer */')
432        CALL WRITE_0D_L( saltStepping,  INDEX_NONE,        CALL WRITE_0D_L( saltStepping,  INDEX_NONE,
433       & 'saltStepping =', '  /* Salinity equation on/off flag */')       & 'saltStepping =', '  /* Salinity equation on/off flag */')
434        CALL WRITE_0D_L( saltAdvection,  INDEX_NONE,        CALL WRITE_0D_L( saltAdvection,  INDEX_NONE,
435       & 'saltAdvection=', '  /* Salinity advection on/off flag */')       & 'saltAdvection=', '  /* Salinity advection on/off flag */')
436        CALL WRITE_0D_L( saltImplVertAdv,INDEX_NONE,'saltImplVertAdv =',        CALL WRITE_0D_L( saltImplVertAdv,INDEX_NONE,'saltImplVertAdv =',
437       &                '/* Sali. implicit vert. advection on/off */')       &                ' /* Sali. implicit vert. advection on/off */')
438        CALL WRITE_0D_L( saltForcing,  INDEX_NONE,        CALL WRITE_0D_L( saltForcing,  INDEX_NONE,
439       & 'saltForcing  =', '  /* Salinity forcing on/off flag */')       & 'saltForcing  =', '  /* Salinity forcing on/off flag */')
440        WRITE(msgBuf,'(A)') '//  '        CALL WRITE_0D_L( saltIsActiveTr, INDEX_NONE, 'saltIsActiveTr =',
441        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,       & ' /* Salt  is a dynamically Active Tracer */')
442       &                    SQUEEZE_RIGHT , 1)        CALL WRITE_0D_I( readBinaryPrec, INDEX_NONE, ' readBinaryPrec =',
443         &  ' /* Precision used for reading binary files */')
444          CALL WRITE_0D_I(writeBinaryPrec, INDEX_NONE, 'writeBinaryPrec =',
445         &  ' /* Precision used for writing binary files */')
446          CALL WRITE_0D_L( globalFiles,  INDEX_NONE,
447         & ' globalFiles =',' /* write "global" (=not per tile) files */')
448          CALL WRITE_0D_L( useSingleCpuIO,  INDEX_NONE,
449         & ' useSingleCpuIO =', ' /* only master MPI process does I/O */')
450          CALL WRITE_0D_L( debugMode,  INDEX_NONE,
451         & ' debugMode  =', '  /* Debug Mode on/off flag */')
452          CALL WRITE_0D_I( debLevA, INDEX_NONE,
453         & '   debLevA  =', '  /* 1rst level of debugging */')
454          CALL WRITE_0D_I( debLevB, INDEX_NONE,
455         & '   debLevB  =', '  /* 2nd  level of debugging */')
456          CALL WRITE_0D_I( debugLevel, INDEX_NONE,
457         & ' debugLevel =', '  /* select debugging level */')
458    
459        WRITE(msgBuf,'(A)')        WRITE(msgBuf,'(A)') '//  '
460          CALL PRINT_MESSAGE( msgBuf, ioUnit, SQUEEZE_RIGHT, myThid )
461          WRITE(msgBuf,'(A)')
462       & '// Elliptic solver(s) paramters ( PARM02 in namelist ) '       & '// Elliptic solver(s) paramters ( PARM02 in namelist ) '
463        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,        CALL PRINT_MESSAGE( msgBuf, ioUnit, SQUEEZE_RIGHT, myThid )
      &                    SQUEEZE_RIGHT , 1)  
464        WRITE(msgBuf,'(A)') '//  '        WRITE(msgBuf,'(A)') '//  '
465        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,        CALL PRINT_MESSAGE( msgBuf, ioUnit, SQUEEZE_RIGHT, myThid )
      &                    SQUEEZE_RIGHT , 1)  
466        CALL WRITE_0D_I( cg2dMaxIters,   INDEX_NONE,'cg2dMaxIters =',        CALL WRITE_0D_I( cg2dMaxIters,   INDEX_NONE,'cg2dMaxIters =',
467       &'   /* Upper limit on 2d con. grad iterations  */')       &'   /* Upper limit on 2d con. grad iterations  */')
468        CALL WRITE_0D_I( cg2dChkResFreq, INDEX_NONE,'cg2dChkResFreq =',        CALL WRITE_0D_I( cg2dChkResFreq, INDEX_NONE,'cg2dChkResFreq =',
469       &'   /* 2d con. grad convergence test frequency */')       &'   /* 2d con. grad convergence test frequency */')
470        CALL WRITE_0D_R8( cg2dTargetResidual, INDEX_NONE,        CALL WRITE_0D_RL( cg2dTargetResidual, INDEX_NONE,
471       & 'cg2dTargetResidual =',       & 'cg2dTargetResidual =',
472       &'   /* 2d con. grad target residual  */')       &'   /* 2d con. grad target residual  */')
473        CALL WRITE_0D_R8( cg2dTargetResWunit, INDEX_NONE,        CALL WRITE_0D_RL( cg2dTargetResWunit, INDEX_NONE,
474       & 'cg2dTargetResWunit =',       & 'cg2dTargetResWunit =',
475       &'   /* CG2d target residual [W units] */')       &'   /* CG2d target residual [W units] */')
476        CALL WRITE_0D_I( cg2dPreCondFreq, INDEX_NONE,'cg2dPreCondFreq =',        CALL WRITE_0D_I( cg2dPreCondFreq, INDEX_NONE,'cg2dPreCondFreq =',
477       &'   /* Freq. for updating cg2d preconditioner */')       &'   /* Freq. for updating cg2d preconditioner */')
478          CALL WRITE_0D_L( useSRCGSolver, INDEX_NONE,
479         & 'useSRCGSolver =', '  /* use single reduction CG solver(s) */')
480    
481        WRITE(msgBuf,'(A)') '//  '        WRITE(msgBuf,'(A)') '//  '
482        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,        CALL PRINT_MESSAGE( msgBuf, ioUnit, SQUEEZE_RIGHT, myThid )
483       &                    SQUEEZE_RIGHT , 1)        WRITE(msgBuf,'(A)')
       WRITE(msgBuf,'(A)')  
484       & '// Time stepping paramters ( PARM03 in namelist ) '       & '// Time stepping paramters ( PARM03 in namelist ) '
485        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,        CALL PRINT_MESSAGE( msgBuf, ioUnit, SQUEEZE_RIGHT, myThid )
      &                    SQUEEZE_RIGHT , 1)  
486        WRITE(msgBuf,'(A)') '//  '        WRITE(msgBuf,'(A)') '//  '
487        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,        CALL PRINT_MESSAGE( msgBuf, ioUnit, SQUEEZE_RIGHT, myThid )
488       &                    SQUEEZE_RIGHT , 1)        CALL WRITE_0D_RL( deltaTmom, INDEX_NONE,'deltaTmom =',
       CALL WRITE_0D_I( nIter0, INDEX_NONE,'nIter0 =',  
      &'   /* Base timestep number  */')  
       CALL WRITE_0D_I( nTimeSteps, INDEX_NONE,'nTimeSteps =',  
      &'   /* Number of timesteps */')  
       CALL WRITE_0D_R8( deltaTmom, INDEX_NONE,'deltatTmom =',  
489       &'   /* Momentum equation timestep ( s ) */')       &'   /* Momentum equation timestep ( s ) */')
490        CALL WRITE_0D_R8( deltaTfreesurf,INDEX_NONE,'deltaTfreesurf =',        CALL WRITE_0D_RL( deltaTfreesurf,INDEX_NONE,'deltaTfreesurf =',
491       &'   /* FreeSurface equation timestep ( s ) */')       &  ' /* FreeSurface equation timestep ( s ) */')
492        CALL WRITE_1D_R8( dTtracerLev, Nr, INDEX_K, 'dTtracerLev =',        CALL WRITE_1D_RL( dTtracerLev, Nr, INDEX_K, 'dTtracerLev =',
493       &'   /* Tracer equation timestep ( s ) */')       & '  /* Tracer equation timestep ( s ) */')
494        CALL WRITE_0D_R8( deltaTClock, INDEX_NONE,'deltatTClock  =',        CALL WRITE_0D_RL( deltaTClock, INDEX_NONE,'deltaTClock  =',
495       &'   /* Model clock timestep ( s ) */')       &'   /* Model clock timestep ( s ) */')
496        CALL WRITE_0D_R8( cAdjFreq, INDEX_NONE,'cAdjFreq =',        CALL WRITE_0D_RL( cAdjFreq, INDEX_NONE,'cAdjFreq =',
497       &'   /* Convective adjustment interval ( s ) */')       &'   /* Convective adjustment interval ( s ) */')
498        CALL WRITE_0D_L( forcing_In_AB,INDEX_NONE,'forcing_In_AB =',        CALL WRITE_0D_I( momForcingOutAB, INDEX_NONE, 'momForcingOutAB =',
499       &'   /* put T,S Forcing in Adams-Bash. stepping */')       & ' /* =1: take Momentum Forcing out of Adams-Bash. stepping */')
500        CALL WRITE_0D_R8( abeps, INDEX_NONE,'abeps =',        CALL WRITE_0D_I( tracForcingOutAB, INDEX_NONE,
501       &'   /* Adams-Bashforth stabilizing weight */')       & 'tracForcingOutAB =',
502         & ' /* =1: take T,S,pTr Forcing out of Adams-Bash. stepping */')
503          CALL WRITE_0D_L( momDissip_In_AB,INDEX_NONE,'momDissip_In_AB =',
504         & ' /* put Dissipation Tendency in Adams-Bash. stepping */')
505          CALL WRITE_0D_L( doAB_onGtGs, INDEX_NONE, 'doAB_onGtGs =',
506         &  ' /* apply AB on Tendencies (rather than on T,S)*/')
507          CALL WRITE_0D_RL( abEps, INDEX_NONE,'abEps =',
508         &'   /* Adams-Bashforth-2 stabilizing weight */')
509    #ifdef ALLOW_ADAMSBASHFORTH_3
510          CALL WRITE_0D_RL( alph_AB, INDEX_NONE,'alph_AB =',
511         &'   /* Adams-Bashforth-3 primary factor */')
512          CALL WRITE_0D_RL( beta_AB, INDEX_NONE,'beta_AB =',
513         &'   /* Adams-Bashforth-3 secondary factor */')
514          CALL WRITE_0D_L( startFromPickupAB2, INDEX_NONE,
515         & 'startFromPickupAB2=',' /* start from AB-2 pickup */')
516    #endif
517        IF (useCDscheme) THEN        IF (useCDscheme) THEN
518        CALL WRITE_0D_R8( tauCD, INDEX_NONE,'tauCD =',        CALL WRITE_0D_RL( tauCD, INDEX_NONE,'tauCD =',
519       &'   /* CD coupling time-scale ( s ) */')       &'   /* CD coupling time-scale ( s ) */')
520        CALL WRITE_0D_R8( rCD, INDEX_NONE,'rCD =',        CALL WRITE_0D_RL( rCD, INDEX_NONE,'rCD =',
521       &'   /* Normalised CD coupling parameter */')       &'   /* Normalised CD coupling parameter */')
522          CALL WRITE_0D_RL( epsAB_CD, INDEX_NONE,'epsAB_CD =',
523         &  ' /* AB-2 stabilizing weight for CD-scheme*/')
524          ENDIF
525          i = ILNBLNK(pickupSuff)
526          IF ( i.GT.0 ) THEN
527            CALL WRITE_0D_C( pickupSuff, 0, INDEX_NONE,
528         & 'pickupSuff =', ' /* Suffix of pickup-file to restart from */')
529        ENDIF        ENDIF
530        CALL WRITE_0D_R8( startTime, INDEX_NONE,'startTime =',        CALL WRITE_0D_L( pickupStrictlyMatch, INDEX_NONE,
531       &'   /* Run start time ( s ). */')       & 'pickupStrictlyMatch=',
532        CALL WRITE_0D_R8( endTime, INDEX_NONE,'endTime =',       & ' /* stop if pickup do not strictly match */')
533       &'   /* Integration ending time ( s ). */')        CALL WRITE_0D_I( nIter0, INDEX_NONE, 'nIter0   =',
534        CALL WRITE_0D_R8( pChkPtFreq, INDEX_NONE,'pChkPtFreq =',       &'   /* Run starting timestep number */')
535       &'   /* Permanent restart/checkpoint file interval ( s ). */')        CALL WRITE_0D_I( nTimeSteps, INDEX_NONE,'nTimeSteps =',
536        CALL WRITE_0D_R8( chkPtFreq, INDEX_NONE,'chkPtFreq =',       &  ' /* Number of timesteps */')
537       &'   /* Rolling restart/checkpoint file interval ( s ). */')        CALL WRITE_0D_I( nEndIter, INDEX_NONE, 'nEndIter =',
538         &'   /* Run ending timestep number */')
539          CALL WRITE_0D_RL( baseTime, INDEX_NONE,'baseTime =',
540         &'   /* Model base time ( s ) */')
541          CALL WRITE_0D_RL( startTime, INDEX_NONE,'startTime =',
542         & '  /* Run start time ( s ) */')
543          CALL WRITE_0D_RL( endTime, INDEX_NONE,'endTime  =',
544         &'   /* Integration ending time ( s ) */')
545          CALL WRITE_0D_RL( pChkPtFreq, INDEX_NONE,'pChkPtFreq =',
546         &  ' /* Permanent restart/pickup file interval ( s ) */')
547          CALL WRITE_0D_RL( chkPtFreq, INDEX_NONE,'chkPtFreq  =',
548         &  ' /* Rolling restart/pickup file interval ( s ) */')
549        CALL WRITE_0D_L(pickup_write_mdsio,INDEX_NONE,        CALL WRITE_0D_L(pickup_write_mdsio,INDEX_NONE,
550       &     'pickup_write_mdsio =', '   /* Model IO flag. */')       &     'pickup_write_mdsio =', '   /* Model IO flag. */')
551        CALL WRITE_0D_L(pickup_read_mdsio,INDEX_NONE,        CALL WRITE_0D_L(pickup_read_mdsio,INDEX_NONE,
# Line 400  C     zcoord :: Line 558  C     zcoord ::
558  #endif  #endif
559        CALL WRITE_0D_L(pickup_write_immed,INDEX_NONE,        CALL WRITE_0D_L(pickup_write_immed,INDEX_NONE,
560       &     'pickup_write_immed =','   /* Model IO flag. */')       &     'pickup_write_immed =','   /* Model IO flag. */')
561        CALL WRITE_0D_R8( dumpFreq, INDEX_NONE,'dumpFreq =',        CALL WRITE_0D_L(writePickupAtEnd,INDEX_NONE,
562         &     'writePickupAtEnd =','   /* Model IO flag. */')
563          CALL WRITE_0D_RL( dumpFreq, INDEX_NONE,'dumpFreq =',
564       &'   /* Model state write out interval ( s ). */')       &'   /* Model state write out interval ( s ). */')
565          CALL WRITE_0D_L(dumpInitAndLast,INDEX_NONE,'dumpInitAndLast=',
566         &  ' /* write out Initial & Last iter. model state */')
567        CALL WRITE_0D_L(snapshot_mdsio,INDEX_NONE,        CALL WRITE_0D_L(snapshot_mdsio,INDEX_NONE,
568       &     'snapshot_mdsio =', '   /* Model IO flag. */')       &     'snapshot_mdsio =', '   /* Model IO flag. */')
569  #ifdef ALLOW_MNC  #ifdef ALLOW_MNC
570        CALL WRITE_0D_L(snapshot_mnc,INDEX_NONE,        CALL WRITE_0D_L(snapshot_mnc,INDEX_NONE,
571       &     'snapshot_mnc =', '   /* Model IO flag. */')       &     'snapshot_mnc =', '   /* Model IO flag. */')
572  #endif  #endif
573        CALL WRITE_0D_R8( monitorFreq, INDEX_NONE,'monitorFreq =',        CALL WRITE_0D_RL( monitorFreq, INDEX_NONE,'monitorFreq =',
574       &'   /* Monitor output interval ( s ). */')       &'   /* Monitor output interval ( s ). */')
575          CALL WRITE_0D_I( monitorSelect, INDEX_NONE, 'monitorSelect =',
576         & ' /* select group of variables to monitor */')
577        CALL WRITE_0D_L(monitor_stdio,INDEX_NONE,        CALL WRITE_0D_L(monitor_stdio,INDEX_NONE,
578       &     'monitor_stdio =', '   /* Model IO flag. */')       &     'monitor_stdio =', '   /* Model IO flag. */')
579  #ifdef ALLOW_MNC  #ifdef ALLOW_MNC
580        CALL WRITE_0D_L(monitor_mnc,INDEX_NONE,        CALL WRITE_0D_L(monitor_mnc,INDEX_NONE,
581       &     'monitor_mnc =', '   /* Model IO flag. */')       &     'monitor_mnc =', '   /* Model IO flag. */')
582  #endif  #endif
583        CALL WRITE_0D_R8( externForcingPeriod, INDEX_NONE,        CALL WRITE_0D_RL( externForcingPeriod, INDEX_NONE,
584       &   'externForcingPeriod =', '   /* forcing period (s) */')       &   'externForcingPeriod =', '   /* forcing period (s) */')
585        CALL WRITE_0D_R8( externForcingCycle, INDEX_NONE,        CALL WRITE_0D_RL( externForcingCycle, INDEX_NONE,
586       &   'externForcingCycle =', '   /* period of the cyle (s). */')       &   'externForcingCycle =', '   /* period of the cyle (s). */')
587        CALL WRITE_0D_R8( tauThetaClimRelax, INDEX_NONE,        CALL WRITE_0D_RL( tauThetaClimRelax, INDEX_NONE,
588       &   'tauThetaClimRelax =', '   /* relaxation time scale (s) */')       &   'tauThetaClimRelax =', '   /* relaxation time scale (s) */')
589        CALL WRITE_0D_R8( tauSaltClimRelax, INDEX_NONE,        CALL WRITE_0D_RL( tauSaltClimRelax, INDEX_NONE,
590       &   'tauSaltClimRelax =',  '   /* relaxation time scale (s) */')       &   'tauSaltClimRelax =',  '   /* relaxation time scale (s) */')
591        CALL WRITE_0D_R8( latBandClimRelax, INDEX_NONE,        CALL WRITE_0D_RL( latBandClimRelax, INDEX_NONE,
592       &   'latBandClimRelax =', '   /* max. Lat. where relaxation */')       &   'latBandClimRelax =', '   /* max. Lat. where relaxation */')
593    
594        WRITE(msgBuf,'(A)') '//  '        WRITE(msgBuf,'(A)') '//  '
595        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,        CALL PRINT_MESSAGE( msgBuf, ioUnit, SQUEEZE_RIGHT, myThid )
596       &                    SQUEEZE_RIGHT , 1)        WRITE(msgBuf,'(A)')
       WRITE(msgBuf,'(A)')  
597       & '// Gridding paramters ( PARM04 in namelist ) '       & '// Gridding paramters ( PARM04 in namelist ) '
598        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,        CALL PRINT_MESSAGE( msgBuf, ioUnit, SQUEEZE_RIGHT, myThid )
      &                    SQUEEZE_RIGHT , 1)  
599        WRITE(msgBuf,'(A)') '//  '        WRITE(msgBuf,'(A)') '//  '
600        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,        CALL PRINT_MESSAGE( msgBuf, ioUnit, SQUEEZE_RIGHT, myThid )
      &                    SQUEEZE_RIGHT , 1)  
601        CALL WRITE_0D_L( usingCartesianGrid, INDEX_NONE,        CALL WRITE_0D_L( usingCartesianGrid, INDEX_NONE,
602       & 'usingCartesianGrid =',       & 'usingCartesianGrid =',
603       &'   /* Cartesian coordinates flag ( True / False ) */')       & ' /* Cartesian coordinates flag ( True/False ) */')
       CALL WRITE_0D_L( usingSphericalPolarGrid, INDEX_NONE,  
      & 'usingSphericalPolarGrid =',  
      &'   /* Spherical coordinates flag ( True / False ) */')  
604        CALL WRITE_0D_L( usingCylindricalGrid, INDEX_NONE,        CALL WRITE_0D_L( usingCylindricalGrid, INDEX_NONE,
605       & 'usingCylindricalGrid =',       & 'usingCylindricalGrid =',
606       &'   /* Spherical coordinates flag ( True / False ) */')       & ' /* Cylindrical coordinates flag ( True/False ) */')
607        CALL WRITE_0D_L( groundAtK1, INDEX_NONE, 'groundAtK1 =',        CALL WRITE_0D_L( usingSphericalPolarGrid, INDEX_NONE,
608       &'   /* Lower Boundary (ground) at the surface(k=1) ( T / F ) */')       & 'usingSphericalPolarGrid =',
609        CALL WRITE_0D_R8( Ro_SeaLevel, INDEX_NONE,'Ro_SeaLevel =',       & ' /* Spherical coordinates flag ( True/False ) */')
610       &'   /* r(1) ( units of r ) */')        CALL WRITE_0D_L( usingCurvilinearGrid, INDEX_NONE,
611        CALL WRITE_0D_R8( rkFac, INDEX_NONE,'rkFac =',       & 'usingCurvilinearGrid =',
612       &'   /* minus Vertical index orientation  */')       & ' /* Curvilinear coordinates flag ( True/False ) */')
613        CALL WRITE_0D_R8( horiVertRatio, INDEX_NONE,'horiVertRatio =',        CALL WRITE_0D_I( selectSigmaCoord, INDEX_NONE,
614       &'   /* Ratio on units : Horiz - Vertical */')       & 'selectSigmaCoord =',
615  c     CALL WRITE_1D_R8( delZ,Nr, INDEX_K,'delZ = ',       & ' /* Hybrid-Sigma Vert. Coordinate option */')
616  c    &'   /* W spacing ( m ) */')        CALL WRITE_0D_RL( Ro_SeaLevel, INDEX_NONE,'Ro_SeaLevel =',
617  c     CALL WRITE_1D_R8( delP,Nr, INDEX_K,'delP = ',       & ' /* r(1) ( units of r == '//rUnits//' ) */')
618  c    &'   /* W spacing ( Pa ) */')        CALL WRITE_0D_RL( rSigmaBnd, INDEX_NONE, 'rSigmaBnd =',
619  c     CALL WRITE_1D_R8( delR,Nr, INDEX_K,'delR = ',       & ' /* r/sigma transition ( units of r == '//rUnits//' ) */')
620  c    &'   /* W spacing ( units of r ) */')        CALL WRITE_0D_RL( rkSign, INDEX_NONE,'rkSign =',
621        CALL WRITE_1D_R8( drC,Nr, INDEX_K,'drC = ',       &'   /* index orientation relative to vertical coordinate */')
622          CALL WRITE_0D_RL( gravitySign, INDEX_NONE,'gravitySign =',
623         &  ' /* gravity orientation relative to vertical coordinate */')
624          IF ( usingZCoords ) THEN
625           CALL WRITE_0D_RL( mass2rUnit, INDEX_NONE,'mass2rUnit =',
626         & ' /* convert mass per unit area [kg/m2] to r-units [m] */')
627           CALL WRITE_0D_RL( rUnit2mass, INDEX_NONE,'rUnit2mass =',
628         & ' /* convert r-units [m] to mass per unit area [kg/m2] */')
629          ENDIF
630          IF ( usingPCoords ) THEN
631           CALL WRITE_0D_RL( mass2rUnit, INDEX_NONE,'mass2rUnit =',
632         & ' /* convert mass per unit area [kg/m2] to r-units [Pa] */')
633           CALL WRITE_0D_RL( rUnit2mass, INDEX_NONE,'rUnit2mass =',
634         & ' /* convert r-units [Pa] to mass per unit area [kg/m2] */')
635          ENDIF
636          CALL WRITE_COPY1D_RS( bufRL, drC, Nr, INDEX_K, 'drC = ',
637       &'   /* C spacing ( units of r ) */')       &'   /* C spacing ( units of r ) */')
638        CALL WRITE_1D_R8( drF,Nr, INDEX_K,'drF = ',        CALL WRITE_COPY1D_RS( bufRL, drF, Nr, INDEX_K, 'drF = ',
639       &'   /* W spacing ( units of r ) */')       &'   /* W spacing ( units of r ) */')
640        CALL WRITE_1D_R8( delX, Nx, INDEX_I,'delX = ',        IF ( selectSigmaCoord.NE.0 ) THEN
641       &'   /* U spacing ( m - cartesian, degrees - spherical ) */')         CALL WRITE_COPY1D_RS( bufRL,dAHybSigF,Nr,INDEX_K,'dAHybSigF =',
642        CALL WRITE_1D_R8( delY, Ny, INDEX_J,'delY = ',       &  ' /* vertical increment of Hybrid-sigma Coeff. (-) */')
643       &'   /* V spacing ( m - cartesian, degrees - spherical ) */')         CALL WRITE_COPY1D_RS( bufRL,dBHybSigF,Nr,INDEX_K,'dBHybSigF =',
644        CALL WRITE_0D_R8( phiMin, INDEX_NONE,'phiMin = ',       &  ' /* vertical increment of Hybrid-sigma Coeff. (-) */')
645       &'   /* South edge (ignored - cartesian, degrees - spherical ) */')        ENDIF
646        CALL WRITE_0D_R8( thetaMin, INDEX_NONE,'thetaMin = ',        IF ( .NOT.usingCurvilinearGrid ) THEN
647       &'   /* West edge ( ignored - cartesian, degrees - spherical ) */')         CALL WRITE_1D_RL( delX, Nx, INDEX_I, 'delX = ',
648        CALL WRITE_0D_R8( rSphere, INDEX_NONE,'rSphere = ',       &  ' /* U spacing ( m - cartesian, degrees - spherical ) */')
649       &'   /* Radius ( ignored - cartesian, m - spherical ) */')         CALL WRITE_1D_RL( delY, Ny, INDEX_J, 'delY = ',
650        DO bi=1,nSx       &  ' /* V spacing ( m - cartesian, degrees - spherical ) */')
651         DO I=1,sNx        ENDIF
652          xcoord((bi-1)*sNx+I) = xC(I,1,bi,1)        CALL WRITE_0D_RL( xgOrigin, INDEX_NONE,'xgOrigin = ',
653         ENDDO       &'/* X-axis origin of West  edge (cartesian: m, lat-lon: deg.) */')
654        ENDDO        CALL WRITE_0D_RL( ygOrigin, INDEX_NONE,'ygOrigin = ',
655        CALL WRITE_1D_R8( xcoord, sNx*nSx, INDEX_I,'xcoord = ',       &'/* Y-axis origin of South edge (cartesian: m, lat-lon: deg.) */')
656       &'   /* P-point X coord ( m - cartesian, degrees - spherical ) */')        CALL WRITE_0D_RL( rSphere, INDEX_NONE,'rSphere = ',
657        DO bj=1,nSy       &  ' /* Radius ( ignored - cartesian, m - spherical ) */')
658         DO J=1,sNy        IF ( usingCurvilinearGrid ) THEN
659          ycoord((bj-1)*sNy+J) = yC(1,J,1,bj)         CALL WRITE_0D_RL( radius_fromHorizGrid, INDEX_NONE,
660         ENDDO       &  'radius_fromHorizGrid = ',
661        ENDDO       &  '/* sphere Radius of input horiz. grid */')
662        CALL WRITE_1D_R8( ycoord, sNy*nSy, INDEX_J,'ycoord = ',        ENDIF
663       &'   /* P-point Y coord ( m - cartesian, degrees - spherical ) */')        CALL WRITE_0D_L(deepAtmosphere,INDEX_NONE, 'deepAtmosphere =',
664        DO K=1,Nr       &  ' /* Deep/Shallow Atmosphere flag (True/False) */')
665         rcoord(K) = rC(K)        coordLine = 1
666        ENDDO        tileLine  = 1
667        CALL WRITE_1D_R8( rcoord, Nr, INDEX_K,'rcoord = ',        CALL WRITE_XY_XLINE_RS( xC, coordLine, tileLine, 'xC',
668       &'   /* P-point R coordinate (  units of r ) */')       I         ': P-point X coord ( deg. or m if cartesian)')
669        DO K=1,Nr+1        CALL WRITE_XY_YLINE_RS( yC, coordLine, tileLine, 'yC',
670         rcoord(K) = rF(K)       I         ': P-point Y coord ( deg. or m if cartesian)')
671        ENDDO        CALL WRITE_COPY1D_RS( bufRL, rC, Nr, INDEX_K, 'rcoord =',
672        CALL WRITE_1D_R8( rcoord, Nr+1, INDEX_K,'rF = ',       &  ' /* P-point R coordinate (  units of r ) */')
673          CALL WRITE_COPY1D_RS( bufRL, rF,Nr+1,INDEX_K, 'rF = ',
674       &'   /* W-Interf. R coordinate (  units of r ) */')       &'   /* W-Interf. R coordinate (  units of r ) */')
675          IF ( selectSigmaCoord.NE.0 ) THEN
676           CALL WRITE_COPY1D_RS(bufRL,aHybSigmF,Nr+1,INDEX_K,'aHybSigmF =',
677         &  ' /* Hybrid-sigma vert. Coord coeff. @ W-Interface (-) */')
678           CALL WRITE_COPY1D_RS(bufRL,bHybSigmF,Nr+1,INDEX_K,'bHybSigmF =',
679         &  ' /* Hybrid-sigma vert. Coord coeff. @ W-Interface (-) */')
680          ENDIF
681          CALL WRITE_1D_RL( deepFacC, Nr,   INDEX_K, 'deepFacC = ',
682         &  ' /* deep-model grid factor @ cell-Center (-) */')
683          CALL WRITE_1D_RL( deepFacF, Nr+1, INDEX_K, 'deepFacF = ',
684         &  ' /* deep-model grid factor @ W-Interface (-) */')
685          CALL WRITE_1D_RL(rVel2wUnit,Nr+1, INDEX_K,'rVel2wUnit =',
686         &  ' /* convert units: rVel -> wSpeed (=1 if z-coord)*/')
687          CALL WRITE_1D_RL(wUnit2rVel,Nr+1, INDEX_K,'wUnit2rVel =',
688         &  ' /* convert units: wSpeed -> rVel (=1 if z-coord)*/')
689          CALL WRITE_1D_RL( dBdrRef,  Nr,   INDEX_K, 'dBdrRef =',
690         & ' /* Vertical grad. of reference buoyancy [(m/s/r)^2] */')
691          CALL WRITE_0D_L( rotateGrid, INDEX_NONE,
692         & 'rotateGrid =',' /* use rotated grid ( True/False ) */')
693          CALL WRITE_0D_RL( phiEuler, INDEX_NONE,'phiEuler =',
694         &' /* Euler angle, rotation about original z-coordinate [rad] */')
695          CALL WRITE_0D_RL( thetaEuler, INDEX_NONE,'thetaEuler =',
696         & ' /* Euler angle, rotation about new x-coordinate [rad] */')
697          CALL WRITE_0D_RL( psiEuler, INDEX_NONE,'psiEuler =',
698         & ' /* Euler angle, rotation about new z-coordinate [rad] */')
699    
700  C     Grid along selected grid lines  C     Grid along selected grid lines
701        coordLine = 1        coordLine = 1
702        tileLine  = 1        tileLine  = 1
703        CALL WRITE_XY_XLINE_RS( dxF, coordLine, tileLine,        CALL WRITE_XY_XLINE_RS( dxF, coordLine, tileLine, 'dxF',
704       I 'dxF','( m - cartesian, degrees - spherical )')       I              '( units: m )' )
705        CALL WRITE_XY_YLINE_RS( dxF, coordLine, tileLine,        CALL WRITE_XY_YLINE_RS( dxF, coordLine, tileLine, 'dxF',
706       I 'dxF','( m - cartesian, degrees - spherical )')       I              '( units: m )' )
707        CALL WRITE_XY_XLINE_RS( dyF, coordLine, tileLine,        CALL WRITE_XY_XLINE_RS( dyF, coordLine, tileLine, 'dyF',
708       I 'dyF','( m - cartesian, degrees - spherical )')       I              '( units: m )' )
709        CALL WRITE_XY_YLINE_RS( dyF, coordLine, tileLine,        CALL WRITE_XY_YLINE_RS( dyF, coordLine, tileLine, 'dyF',
710       I 'dyF','( m - cartesian, degrees - spherical )')       I              '( units: m )' )
711        CALL WRITE_XY_XLINE_RS( dxG, coordLine, tileLine,        CALL WRITE_XY_XLINE_RS( dxG, coordLine, tileLine, 'dxG',
712       I 'dxG','( m - cartesian, degrees - spherical )')       I              '( units: m )' )
713        CALL WRITE_XY_YLINE_RS( dxG, coordLine, tileLine,        CALL WRITE_XY_YLINE_RS( dxG, coordLine, tileLine, 'dxG',
714       I 'dxG','( m - cartesian, degrees - spherical )')       I              '( units: m )' )
715        CALL WRITE_XY_XLINE_RS( dyG, coordLine, tileLine,        CALL WRITE_XY_XLINE_RS( dyG, coordLine, tileLine, 'dyG',
716       I 'dyG','( m - cartesian, degrees - spherical )')       I              '( units: m )' )
717        CALL WRITE_XY_YLINE_RS( dyG, coordLine, tileLine,        CALL WRITE_XY_YLINE_RS( dyG, coordLine, tileLine, 'dyG',
718       I 'dyG','( m - cartesian, degrees - spherical )')       I              '( units: m )' )
719        CALL WRITE_XY_XLINE_RS( dxC, coordLine, tileLine,        CALL WRITE_XY_XLINE_RS( dxC, coordLine, tileLine, 'dxC',
720       I 'dxC','( m - cartesian, degrees - spherical )')       I              '( units: m )' )
721        CALL WRITE_XY_YLINE_RS( dxC, coordLine, tileLine,        CALL WRITE_XY_YLINE_RS( dxC, coordLine, tileLine, 'dxC',
722       I 'dxC','( m - cartesian, degrees - spherical )')       I              '( units: m )' )
723        CALL WRITE_XY_XLINE_RS( dyC, coordLine, tileLine,        CALL WRITE_XY_XLINE_RS( dyC, coordLine, tileLine, 'dyC',
724       I 'dyC','( m - cartesian, degrees - spherical )')       I              '( units: m )' )
725        CALL WRITE_XY_YLINE_RS( dyC, coordLine, tileLine,        CALL WRITE_XY_YLINE_RS( dyC, coordLine, tileLine, 'dyC',
726       I 'dyC','( m - cartesian, degrees - spherical )')       I              '( units: m )' )
727        CALL WRITE_XY_XLINE_RS( dxV, coordLine, tileLine,        CALL WRITE_XY_XLINE_RS( dxV, coordLine, tileLine, 'dxV',
728       I 'dxV','( m - cartesian, degrees - spherical )')       I              '( units: m )' )
729        CALL WRITE_XY_YLINE_RS( dxV, coordLine, tileLine,        CALL WRITE_XY_YLINE_RS( dxV, coordLine, tileLine, 'dxV',
730       I 'dxV','( m - cartesian, degrees - spherical )')       I              '( units: m )' )
731        CALL WRITE_XY_XLINE_RS( dyU, coordLine, tileLine,        CALL WRITE_XY_XLINE_RS( dyU, coordLine, tileLine, 'dyU',
732       I 'dyU','( m - cartesian, degrees - spherical )')       I              '( units: m )' )
733        CALL WRITE_XY_YLINE_RS( dyU, coordLine, tileLine,        CALL WRITE_XY_YLINE_RS( dyU, coordLine, tileLine, 'dyU',
734       I 'dyU','( m - cartesian, degrees - spherical )')       I              '( units: m )' )
735        CALL WRITE_XY_XLINE_RS( rA, coordLine, tileLine,        CALL WRITE_XY_XLINE_RS( rA , coordLine, tileLine, 'rA ',
736       I 'rA','( m - cartesian, degrees - spherical )')       I              '( units: m^2 )' )
737        CALL WRITE_XY_YLINE_RS( rA, coordLine, tileLine,        CALL WRITE_XY_YLINE_RS( rA , coordLine, tileLine, 'rA ',
738       I 'rA','( m - cartesian, degrees - spherical )')       I              '( units: m^2 )' )
739        CALL WRITE_XY_XLINE_RS( rAw, coordLine, tileLine,        CALL WRITE_XY_XLINE_RS( rAw, coordLine, tileLine, 'rAw',
740       I 'rAw','( m - cartesian, degrees - spherical )')       I              '( units: m^2 )' )
741        CALL WRITE_XY_YLINE_RS( rAw, coordLine, tileLine,        CALL WRITE_XY_YLINE_RS( rAw, coordLine, tileLine, 'rAw',
742       I 'rAw','( m - cartesian, degrees - spherical )')       I              '( units: m^2 )' )
743        CALL WRITE_XY_XLINE_RS( rAs, coordLine, tileLine,        CALL WRITE_XY_XLINE_RS( rAs, coordLine, tileLine, 'rAs',
744       I 'rAs','( m - cartesian, degrees - spherical )')       I              '( units: m^2 )' )
745        CALL WRITE_XY_YLINE_RS( rAs, coordLine, tileLine,        CALL WRITE_XY_YLINE_RS( rAs, coordLine, tileLine, 'rAs',
746       I 'rAs','( m - cartesian, degrees - spherical )')       I              '( units: m^2 )' )
747    
748          CALL WRITE_0D_RL( globalArea, INDEX_NONE, 'globalArea =',
749         & ' /* Integrated horizontal Area (m^2) */')
750    
751          i = ILNBLNK(the_run_name)
752          IF ( i.GT.0 ) THEN
753            CALL WRITE_0D_C( the_run_name, i, INDEX_NONE,
754         &    'the_run_name = ', '/* Name of this simulation */' )
755          ENDIF
756    
757          WRITE(msgBuf,'(A)')
758         &'// ======================================================='
759          CALL PRINT_MESSAGE( msgBuf, ioUnit, SQUEEZE_RIGHT, myThid )
760          WRITE(msgBuf,'(A)') '// End of Model config. summary'
761          CALL PRINT_MESSAGE( msgBuf, ioUnit, SQUEEZE_RIGHT, myThid )
762          WRITE(msgBuf,'(A)')
763         &'// ======================================================='
764          CALL PRINT_MESSAGE( msgBuf, ioUnit, SQUEEZE_RIGHT, myThid )
765        WRITE(msgBuf,'(A)') ' '        WRITE(msgBuf,'(A)') ' '
766        CALL PRINT_MESSAGE( msgBuf, standardMessageUnit,        CALL PRINT_MESSAGE( msgBuf, ioUnit, SQUEEZE_RIGHT, myThid )
      &  SQUEEZE_RIGHT , 1)  
767    
768        _END_MASTER(myThid)        _END_MASTER(myThid)
769        _BARRIER        _BARRIER
770    
   
771        RETURN        RETURN
   100 FORMAT(A,  
      &' '  
      &)  
772        END        END
   
Legend:
Removed from v.1.69  
changed lines
  Added in v.1.129
  ViewVC Help
Powered by ViewVC 1.1.22