/[MITgcm]/MITgcm/pkg/thsice/thsice_diagnostics_init.F
ViewVC logotype

Diff of /MITgcm/pkg/thsice/thsice_diagnostics_init.F

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

revision 1.3 by jmc, Sun Feb 25 20:00:30 2007 UTC revision 1.4 by jmc, Wed Apr 4 02:40:42 2007 UTC
# Line 22  C     === Global variables === Line 22  C     === Global variables ===
22  #include "EEPARAMS.h"  #include "EEPARAMS.h"
23  c #include "SIZE.h"  c #include "SIZE.h"
24  c #include "PARAMS.h"  c #include "PARAMS.h"
25    #include "THSICE_SIZE.h"
26  c #include "THSICE_PARAMS.h"  c #include "THSICE_PARAMS.h"
27    
28  C     !INPUT/OUTPUT PARAMETERS:  C     !INPUT/OUTPUT PARAMETERS:
29  C     === Routine arguments ===  C     === Routine arguments ===
30  C     myThid ::  my thread Id number  C     myThid ::  My Thread Id. number
31        INTEGER myThid        INTEGER myThid
32  CEOP  CEOP
33    
# Line 45  C     diagTitle :: local string (80c): d Line 46  C     diagTitle :: local string (80c): d
46        CHARACTER*16  diagUnits        CHARACTER*16  diagUnits
47        CHARACTER*(80) diagTitle        CHARACTER*(80) diagTitle
48        INTEGER       numFract, numThick        INTEGER       numFract, numThick
49          INTEGER       n
50          LOGICAL       flag
51          CHARACTER*16  locName
52          CHARACTER*4   THSICE_DIAG_SUFX, diagSufx
53          EXTERNAL      THSICE_DIAG_SUFX
54    
55  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
56    
# Line 178  c     IF ( useDiagnotics ) THEN Line 184  c     IF ( useDiagnotics ) THEN
184          CALL DIAGNOSTICS_ADD2LIST( diagNum,          CALL DIAGNOSTICS_ADD2LIST( diagNum,
185       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
186    
187           DO n=-5,-9,-1
188            flag = .TRUE.
189            IF     ( n .EQ. GAD_SI_FRAC  ) THEN
190              locName   = 'SeaIce Fraction '
191              diagUnits = 'm^2/s           '
192            ELSEIF ( n .EQ. GAD_SI_HSNOW ) THEN
193              locName   = 'Snow Thickness  '
194              diagUnits = 'm^3/s           '
195            ELSEIF ( n .EQ. GAD_SI_HICE  ) THEN
196              locName   = 'SeaIce Thickness'
197              diagUnits = 'm^3/s           '
198            ELSEIF ( n .EQ. GAD_SI_QICE1 ) THEN
199              locName   = 'Ice enthalpy l=1'
200              diagUnits = 'J/kg.m^3/s      '
201    c         diagUnits = 'W.m^3/kg        '
202            ELSEIF ( n .EQ. GAD_SI_QICE2 ) THEN
203              locName   = 'Ice enthalpy l=2'
204              diagUnits = 'J/kg.m^3/s      '
205            ELSE
206              flag = .FALSE.
207            ENDIF
208            IF ( flag ) THEN
209    C-     Advective flux:
210            diagSufx = THSICE_DIAG_SUFX( n, myThid )
211            diagName  = 'ADVx'//diagSufx
212            diagTitle = 'Zonal      Advective Flux of '//locName
213            WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'M1      '
214            CALL DIAGNOSTICS_ADD2LIST( diagNum,
215         I               diagName, diagCode, diagUnits, diagTitle, myThid )
216            diagName  = 'ADVy'//diagSufx
217            diagTitle = 'Meridional Advective Flux of '//locName
218            WRITE(diagCode,'(A,I3.3,A)') 'VV   ',diagNum , 'M1      '
219            CALL DIAGNOSTICS_ADD2LIST( diagNum,
220         I               diagName, diagCode, diagUnits, diagTitle, myThid )
221            ENDIF
222          ENDDO
223    
224  c     ENDIF  c     ENDIF
225    
226  #endif /* ALLOW_DIAGNOSTICS */  #endif /* ALLOW_DIAGNOSTICS */
227    
228        RETURN        RETURN
229        END        END
230    
231    C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
232    
233    CBOP 0
234    C     !ROUTINE: THSICE_DIAG_SUFX
235    
236    C     !INTERFACE:
237          CHARACTER*4 FUNCTION THSICE_DIAG_SUFX( tracerId, myThid )
238    
239    C     !DESCRIPTION:
240    C     *==========================================================*
241    C     | FUNCTION THSICE_DIAG_SUFX
242    C     | o Return diagnostic suffix (4 character long) for the
243    C     |   "tracerId" tracer (used to build diagnostic names).
244    C     *==========================================================*
245    
246    C     !USES:
247          IMPLICIT NONE
248    #include "EEPARAMS.h"
249    #include "THSICE_SIZE.h"
250    
251    C     !INPUT PARAMETERS:
252    C     tracerId   ::  tracer identifier
253    C     myThid     ::  my thread Id number
254          INTEGER      tracerId
255          INTEGER      myThid
256    CEOP
257    
258    C     !LOCAL VARIABLES:
259    
260    C--   Set diagnostic suffix (4 character long) for the "tracerId" tracer
261          IF    ( tracerId.EQ.GAD_SI_FRAC  ) THEN
262            THSICE_DIAG_SUFX = 'SIfr'
263          ELSEIF( tracerId.EQ.GAD_SI_HSNOW ) THEN
264            THSICE_DIAG_SUFX = 'SIhs'
265          ELSEIF( tracerId.EQ.GAD_SI_HICE  ) THEN
266            THSICE_DIAG_SUFX = 'SIhi'
267          ELSEIF( tracerId.EQ.GAD_SI_QICE1 ) THEN
268            THSICE_DIAG_SUFX = 'SIq1'
269          ELSEIF( tracerId.EQ.GAD_SI_QICE2 ) THEN
270            THSICE_DIAG_SUFX = 'SIq2'
271          ELSE
272            THSICE_DIAG_SUFX = 'iiii'
273          ENDIF
274    
275          RETURN
276          END
Legend:
Removed from v.1.3  
changed lines
  Added in v.1.4
  ViewVC Help
Powered by ViewVC 1.1.22