/[MITgcm]/MITgcm/pkg/diagnostics/diagnostics_main_init.F
ViewVC logotype

Diff of /MITgcm/pkg/diagnostics/diagnostics_main_init.F

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

revision 1.26 by jmc, Sun Jun 25 23:06:24 2006 UTC revision 1.30 by jmc, Fri Apr 27 15:46:04 2007 UTC
# Line 396  c    I          diagName, diagCode, diag Line 396  c    I          diagName, diagCode, diag
396       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
397    
398        diagName  = 'UVELPHI '        diagName  = 'UVELPHI '
399        diagTitle = DIAGS_MK_TITLE( 'Zonal Transp of '        diagTitle = DIAGS_MK_TITLE( 'Zonal Mass-Weight Transp of '
400       I                 //pTitle//' Anomaly', myThid )       I                 //pTitle//' Anomaly', myThid )
401        diagUnits = 'm^3/s^3         '        diagUnits = 'm^3/s^3         '
402        WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '        WRITE(diagCode,'(A,I3.3,A)') 'UU   ',diagNum+2,'MR      '
# Line 404  c    I          diagName, diagCode, diag Line 404  c    I          diagName, diagCode, diag
404       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
405    
406        diagName  = 'VVELPHI '        diagName  = 'VVELPHI '
407        diagTitle = DIAGS_MK_TITLE( 'Meridional Transp of '        diagTitle = DIAGS_MK_TITLE( 'Meridional Mass-Weight Transp of '
408       I                 //pTitle//' Anomaly', myThid )       I                 //pTitle//' Anomaly', myThid )
409        diagUnits = 'm^3/s^3         '        diagUnits = 'm^3/s^3         '
410        WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '        WRITE(diagCode,'(A,I3.3,A)') 'VV   ', diagNum ,'MR      '
# Line 498  c     diagTitle = 'Square of ocean botto Line 498  c     diagTitle = 'Square of ocean botto
498       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
499    
500  C--   surface fluxes:  C--   surface fluxes:
501        diagName  = 'TAUX    '        diagName  = 'oceTAUX '
502        diagTitle = 'zonal surface wind stress, >0 increases uVel'        diagTitle = 'zonal surface wind stress, >0 increases uVel'
503        diagUnits = 'N/m^2           '        diagUnits = 'N/m^2           '
504        diagCode  = 'SU      U1      '        diagCode  = 'SU      U1      '
505        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
506       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
507    
508        diagName  = 'TAUY    '        diagName  = 'oceTAUY '
509        diagTitle = 'meridional surf. wind stress, >0 increases vVel'        diagTitle = 'meridional surf. wind stress, >0 increases vVel'
510        diagUnits = 'N/m^2           '        diagUnits = 'N/m^2           '
511        diagCode  = 'SV      U1      '        diagCode  = 'SV      U1      '
512        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
513       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
514    
515        diagName  = 'TFLUX   '        diagName  = 'atmPload'
516        diagTitle = 'net surface heat flux, >0 increases theta       '        diagTitle = 'Atmospheric pressure loading'
517          diagUnits = 'Pa              '
518          diagCode  = 'SM      U1      '
519          CALL DIAGNOSTICS_ADD2LIST( diagNum,
520         I          diagName, diagCode, diagUnits, diagTitle, myThid )
521    
522          diagName  = 'sIceLoad'
523          diagTitle = 'sea-ice loading (in Mass of ice+snow / area unit)'
524          diagUnits = 'kg/m^2          '
525          diagCode  = 'SM      U1      '
526          CALL DIAGNOSTICS_ADD2LIST( diagNum,
527         I          diagName, diagCode, diagUnits, diagTitle, myThid )
528    
529          diagName  = 'oceFWflx'
530          diagTitle = 'net surface Fresh-Water flux into the ocean'
531         &          //' (+=down), >0 decreases salinity'
532          diagUnits = 'kg/m^2/s        '
533          diagCode  = 'SM      U1      '
534          CALL DIAGNOSTICS_ADD2LIST( diagNum,
535         I          diagName, diagCode, diagUnits, diagTitle, myThid )
536    
537          diagName  = 'oceSflux'
538          diagTitle = 'net surface Salt flux into the ocean (+=down),'
539         &          //' >0 increases salinity'
540          diagUnits = 'g/m^2/s         '
541          diagCode  = 'SM      U1      '
542          CALL DIAGNOSTICS_ADD2LIST( diagNum,
543         I          diagName, diagCode, diagUnits, diagTitle, myThid )
544    
545          diagName  = 'oceQnet '
546          diagTitle = 'net surface heat flux into the ocean (+=down),'
547         &          //' >0 increases theta'
548        diagUnits = 'W/m^2           '        diagUnits = 'W/m^2           '
549        diagCode  = 'SM      U1      '        diagCode  = 'SM      U1      '
550        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
551       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
552    
553        diagName  = 'TRELAX  '        diagName  = 'oceQsw  '
554        diagTitle = 'surface temperature relaxation, >0 increases theta'        diagTitle = 'net Short-Wave radiation (+=down),'
555         &          //' >0 increases theta'
556        diagUnits = 'W/m^2           '        diagUnits = 'W/m^2           '
557        diagCode  = 'SM      U1      '        diagCode  = 'SM      U1      '
558        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
559       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
560    
561        diagName  = 'TICE    '        diagName  = 'oceFreez'
562        diagTitle = 'heat from melt/freeze of sea-ice, >0 increases theta'        diagTitle = 'heating from freezing of sea-water (allowFreezing=T)'
563        diagUnits = 'W/m^2           '        diagUnits = 'W/m^2           '
564        diagCode  = 'SM      U1      '        diagCode  = 'SM      U1      '
565        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
566       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
567    
568        diagName  = 'SFLUX   '        diagName  = 'TRELAX  '
569        diagTitle = 'net surface salt flux, >0 increases salt'        diagTitle = 'surface temperature relaxation, >0 increases theta'
570        diagUnits = 'g/m^2/s         '        diagUnits = 'W/m^2           '
571        diagCode  = 'SM      U1      '        diagCode  = 'SM      U1      '
572        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
573       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
# Line 547  C--   surface fluxes: Line 579  C--   surface fluxes:
579        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
580       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
581    
582        diagName  = 'SWFLUX  '        diagName  = 'surForcT'
583        diagTitle = 'net upward SW radiation, >0 increases theta'        diagTitle = 'model surface forcing for Temperature,'
584         &          //' >0 increases theta'
585          diagUnits = 'W/m^2           '
586          diagCode  = 'SM      U1      '
587          CALL DIAGNOSTICS_ADD2LIST( diagNum,
588         I          diagName, diagCode, diagUnits, diagTitle, myThid )
589    
590          diagName  = 'surForcS'
591          diagTitle = 'model surface forcing for Salinity,'
592         &          //' >0 increases salinity'
593          diagUnits = 'g/m^2/s         '
594          diagCode  = 'SM      U1      '
595          CALL DIAGNOSTICS_ADD2LIST( diagNum,
596         I          diagName, diagCode, diagUnits, diagTitle, myThid )
597    
598          diagName  = 'TFLUX   '
599          diagTitle = 'total heat flux (match heat-content variations),'
600         &          //' >0 increases theta'
601        diagUnits = 'W/m^2           '        diagUnits = 'W/m^2           '
602        diagCode  = 'SM      U1      '        diagCode  = 'SM      U1      '
603        CALL DIAGNOSTICS_ADD2LIST( diagNum,        CALL DIAGNOSTICS_ADD2LIST( diagNum,
604       I          diagName, diagCode, diagUnits, diagTitle, myThid )       I          diagName, diagCode, diagUnits, diagTitle, myThid )
605    
606          diagName  = 'SFLUX   '
607          diagTitle = 'total salt flux (match salt-content variations),'
608         &          //' >0 increases salt'
609          diagUnits = 'g/m^2/s         '
610          diagCode  = 'SM      U1      '
611          CALL DIAGNOSTICS_ADD2LIST( diagNum,
612         I          diagName, diagCode, diagUnits, diagTitle, myThid )
613    
614  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|  C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
615        diagName  = 'PRESSURE'  
616  C-    will be changed to an isomorphic name, such as:        diagName  = 'RCENTER '
 c     diagName  = 'RCENTER '  
617  c     diagTitle = 'Cell-Center r-Position (Pressure, Height) (Pa,m)'  c     diagTitle = 'Cell-Center r-Position (Pressure, Height) (Pa,m)'
618        diagTitle = DIAGS_MK_TITLE( 'Cell-Center '//rTitle, myThid )        diagTitle = DIAGS_MK_TITLE( 'Cell-Center '//rTitle, myThid )
619        diagUnits = DIAGS_MK_UNITS( rUnit2c, myThid )        diagUnits = DIAGS_MK_UNITS( rUnit2c, myThid )

Legend:
Removed from v.1.26  
changed lines
  Added in v.1.30

  ViewVC Help
Powered by ViewVC 1.1.22