/[MITgcm]/MITgcm_contrib/submesoscale/code/gmredi_diagnostics_init.F
ViewVC logotype

Diff of /MITgcm_contrib/submesoscale/code/gmredi_diagnostics_init.F

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

revision 1.1 by dimitri, Fri May 30 21:51:24 2008 UTC revision 1.3 by dimitri, Fri May 30 22:24:25 2008 UTC
# Line 170  c     IF ( useDiagnotics ) THEN Line 170  c     IF ( useDiagnotics ) THEN
170          CALL DIAGNOSTICS_ADD2LIST( diagNum,          CALL DIAGNOSTICS_ADD2LIST( diagNum,
171       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
172    
173  #ifdef ALLOW_TAU_EDDY  #ifdef GM_SUBMESO
174    CBFK These are diagnosed whether or not Bolus Advection is used...
175    CBFK They are the quasi-streamfunction, since they are in the
176    CBFK Visbeck -dPsix/dz=u form, not the Fox-Kemper curl(psi)=u form.
177    CBFK They are always included in the GM part above
178    CBFK (i.e., GM diagnostic=GM+Submeso)
179            diagName  = 'SM_PsiX '
180            diagTitle = 'Submeso Bolus transport quasi-streamfunction : X'
181            diagUnits = 'm^2/s           '
182            WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'LR      '
183            CALL DIAGNOSTICS_ADD2LIST( diagNum,
184         I          diagName, diagCode, diagUnits, diagTitle, myThid )
185    
186            diagName  = 'SM_PsiY '
187            diagTitle = 'Submeso Bolus transport quasi-streamfunction : Y'
188            diagUnits = 'm^2/s           '
189            WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'LR      '
190            CALL DIAGNOSTICS_ADD2LIST( diagNum,
191         I          diagName, diagCode, diagUnits, diagTitle, myThid )
192    
193            diagName  = 'SM_ubT  '
194            diagTitle = 'Zon. Submeso Mass-Weight Bolus Transp of Pot Temp'
195            diagUnits = 'degC.m^3/s      '
196            WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '
197            CALL DIAGNOSTICS_ADD2LIST( diagNum,
198         I          diagName, diagCode, diagUnits, diagTitle, myThid )
199    
200            diagName  = 'SM_vbT  '
201            diagTitle = 'Mer. Submeso Mass-Weight Bolus Transp of Pot Temp'
202            diagUnits = 'degC.m^3/s      '
203            WRITE(diagCode,'(A,I3.3,A)') 'VV   ',diagNum,'MR      '
204            CALL DIAGNOSTICS_ADD2LIST( diagNum,
205         I          diagName, diagCode, diagUnits, diagTitle, myThid )
206    
207            diagName  = 'SM_wbT  '
208            diagTitle = 'Rvel Submeso Mass-Weight Bolus Transp of Pot Temp'
209            diagUnits = 'degC.m^3/s      '
210            WRITE(diagCode,'(A,I3.3,A)') 'WM   ',diagNum,'LR      '
211            CALL DIAGNOSTICS_ADD2LIST( diagNum,
212         I          diagName, diagCode, diagUnits, diagTitle, myThid )
213    
214            diagName  = 'SM_KuzTz'
215            diagTitle = 'Zon. Submeso Mass-Weight Kappa Transp of Pot Temp'
216            diagUnits = 'degC.m^3/s      '
217            WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '
218            CALL DIAGNOSTICS_ADD2LIST( diagNum,
219         I          diagName, diagCode, diagUnits, diagTitle, myThid )
220    
221            diagName  = 'SM_KvzTz'
222            diagTitle = 'Mer. Submeso Mass-Weight Kappa Transp of Pot Temp'
223            diagUnits = 'degC.m^3/s      '
224            WRITE(diagCode,'(A,I3.3,A)') 'VV   ',diagNum,'MR      '
225            CALL DIAGNOSTICS_ADD2LIST( diagNum,
226         I          diagName, diagCode, diagUnits, diagTitle, myThid )
227    
228            diagName  = 'SM_KrddT'
229            diagTitle = 'Rvel Submeso Mass-Weight Kappa Transp of Pot Temp'
230            diagUnits = 'degC.m^3/s      '
231            WRITE(diagCode,'(A,I3.3,A)') 'WM   ',diagNum,'LR      '
232            CALL DIAGNOSTICS_ADD2LIST( diagNum,
233         I          diagName, diagCode, diagUnits, diagTitle, myThid )
234    #endif
235    
236    #ifdef ALLOW_EDDYPSI
237          diagName  = 'GMEdTauX'          diagName  = 'GMEdTauX'
238          diagTitle = 'eddy-induced stress X-comp. estimated from Kwx'          diagTitle = 'eddy-induced stress X-comp. estimated from Kwx'
239          diagUnits = 'N/m^2           '          diagUnits = 'N/m^2           '
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.22