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

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

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

revision 1.5 by molod, Fri Feb 27 23:46:52 2004 UTC revision 1.19 by edhill, Tue Jul 6 05:22:48 2004 UTC
# Line 1  Line 1 
1        subroutine diagnostics_init_vals (myThid)  C $Header$
2    C $Name$
3    
4    #include "DIAG_OPTIONS.h"
5    
6    C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
7    CBOP 0
8    C     !ROUTINE: DIAGNOSTICS_INIT_VALS
9    
10    C     !INTERFACE:
11          SUBROUTINE DIAGNOSTICS_INIT_VALS( myThid )
12    
13    C     !DESCRIPTION:
14    C     Initialize Diagnostic Names, Levels, and GPARMS.  Note that GDIAG
15    C     is defined as character*16 and can be to character*1 parse(16)
16    C     with the following codes currently defined:
17    
18  C ***********************************************************************  C ***********************************************************************
19  C *****                                                              ****  C *****                                                              ****
20  C *****  PURPOSE                                                     ****  C *****  PURPOSE                                                     ****
21  C *****    Initialize Diagnostic Names, Levels, and GPARMS           ****  C *****    Initialize Diagnostic Names, Levels, and GPARMS           ****
22  C *****                                                              ****  C *****                                                              ****
23  C *****                                                              ****  C *****                                                              ****
24  C *****  Note: GDIAG is define as character*16 and can be            ****  C *****  Note: GDIAG is defined as character*16 and can be           ****
25  C *****              equivalenced to character*1 parse(16)           ****  C *****              equivalenced to character*1 parse(16)           ****
26  C *****                                                              ****  C *****                                                              ****
27  C *****  parse(1)   => 'S' for a scalar diagnostic                   ****  C *****  parse(1)   => 'S' for a scalar diagnostic                   ****
# Line 15  C ***** Line 31  C *****
31  C *****  parse(2)   => 'U' for a C-grid U-Point                      ****  C *****  parse(2)   => 'U' for a C-grid U-Point                      ****
32  C *****             => 'V' for a C-grid V-Point                      ****  C *****             => 'V' for a C-grid V-Point                      ****
33  C *****             => 'M' for a C-grid Mass Point                   ****  C *****             => 'M' for a C-grid Mass Point                   ****
34  C *****             => 'Z' for a C-grid Vorticity Point              ****  C *****             => 'Z' for a C-grid Corner Point                 ****
35  C *****                                                              ****  C *****                                                              ****
36  C *****  parse(3)   => 'G' for a field on the geophysical grid       ****  C *****  parse(3)   => not used at present                           ****
37  C *****                                                              ****  C *****                                                              ****
38  C *****  parse(4)   => 'P' if diagnostic is positive definite        ****  C *****  parse(4)   => 'P' if diagnostic is positive definite        ****
39  C *****                                                              ****  C *****                                                              ****
# Line 26  C *****             => 'D' to disable an Line 42  C *****             => 'D' to disable an
42  C *****                                                              ****  C *****                                                              ****
43  C *****  parse(6-8) => 3-digit integer corresponding to vector       ****  C *****  parse(6-8) => 3-digit integer corresponding to vector       ****
44  C *****                component mate if parse(1) = 'U' or 'V'       ****  C *****                component mate if parse(1) = 'U' or 'V'       ****
45    C *****                or to the counter if parse(5)='C'             ****
46  C *****                                                              ****  C *****                                                              ****
47  C *****  parse(9)   => 'U' if defined at Model-level upper-edge      ****  C *****  parse(9)   => 'U' if defined at Model-level plus 1/2        ****
48  C *****             => 'M' if defined at Model-level middle          ****  C *****             => 'M' if defined at Model-level middle          ****
49  C *****             => 'L' if defined at Model-level lower-edge      ****  C *****             => 'L' if defined at Model-level minus 1/2       ****
50  C *****                                                              ****  C *****                                                              ****
51  C *****                                                              ****  C *****                                                              ****
52  C ***********************************************************************  C ***********************************************************************
53    
54    C     \begin{center}
55    C       \begin{tabular}[h]{|c|c|}\hline
56    C         \textbf{Positions}  &  \textbf{Characters}
57    C         &  \textbf{Meanings} \\\hline
58    C         parse(1)  &  S  &  scalar \\
59    C                   &  U  &  U-vector component  \\
60    C                   &  V  &  V-vector component  \\
61    C         parse(2)  &  U  &  C-grid U-Point  \\
62    C                   &  V  &  C-grid V-Point  \\
63    C                   &  M  &  C-grid Mass Point  \\
64    C                   &  Z  &  C-grid Corner Point  \\
65    C         parse(3)  &  --  &  --Unused--  \\
66    C         parse(4)  &  P  &  positive definite  \\
67    C         parse(5 ) &  C  &  with counter array  \\
68    C                   &  D  &  disable an array for output  \\
69    C         parse(6--8) & '123'  &  3-digit mate number \\
70    C         parse(9)  &  U  &  model-level plus 1/2  \\
71    C                   &  M  &  model-level middle  \\
72    C                   &  L  &  model-level minus 1/2  \\
73    C         parse(10) &  0  &  levels = 0  \\
74    C                   &  1  &  levels = 1  \\
75    C                   &  R  &  levels = Nr  \\
76    C                   &  L  &  levels = MAX(Nr,NrPhys)  \\
77    C                   &  M  &  levels = MAX(Nr,NrPhys) - 1  \\
78    C       \end{tabular}
79    C     \end{center}
80    
81    C     !USES:
82        implicit none        implicit none
 #include "CPP_OPTIONS.h"  
83  #include "SIZE.h"  #include "SIZE.h"
84    #include "EEPARAMS.h"
85    #include "EESUPPORT.h"
86    #include "PARAMS.h"
87    
88  #ifdef ALLOW_FIZHI  #ifdef ALLOW_FIZHI
89  #include "fizhi_SIZE.h"  #include "fizhi_SIZE.h"
90  #else  #else
91         integer Nrphys         integer Nrphys
92         parameter (Nrphys=1)         parameter (Nrphys=0)
93  #endif  #endif
94    
95  #include "diagnostics_SIZE.h"  #include "diagnostics_SIZE.h"
96  #include "diagnostics.h"  #include "diagnostics.h"
97    
98    C     !INPUT PARAMETERS:
99        integer myThid        integer myThid
100        integer  n  CEOP
101    
102          integer n, nlevs
103    
104          nlevs   = max(Nr,Nrphys)
105    
106        do n=1,ndiagt        do n=1,ndiagt
107        kdiag (n) = 0          kdiag(n) = 0
108        cdiag(n) = '        '          cdiag(n) = '        '
109        gdiag(n) = '                '          gdiag(n) = '                '
110        udiag(n) = '                '          udiag(n) = '                '
111        tdiag(n) = ' '          tdiag(n) = ' '
112        enddo        enddo
113    
114  c Diagnostic Names  c Diagnostic Names
# Line 298  c ---------------- Line 350  c ----------------
350        cdiag(235) = 'GWCVS   '        cdiag(235) = 'GWCVS   '
351        cdiag(236) = 'GWCU    '        cdiag(236) = 'GWCU    '
352        cdiag(237) = 'GWCV    '        cdiag(237) = 'GWCV    '
353          cdiag(238) = 'ETAN    '
354          cdiag(239) = 'ETANSQ  '
355          cdiag(240) = 'THETA   '
356          cdiag(241) = 'SALT    '
357          cdiag(242) = 'UVEL    '
358          cdiag(243) = 'VVEL    '
359          cdiag(244) = 'WVEL    '
360          cdiag(245) = 'THETASQ '
361          cdiag(246) = 'SALTSQ  '
362          cdiag(247) = 'UVELSQ  '
363          cdiag(248) = 'VVELSQ  '
364          cdiag(249) = 'WVELSQ  '
365          cdiag(250) = 'UVELVVEL'
366          cdiag(251) = 'UVELMASS'
367          cdiag(252) = 'VVELMASS'
368          cdiag(253) = 'WVELMASS'
369          cdiag(254) = 'UTHMASS '
370          cdiag(255) = 'VTHMASS '
371          cdiag(256) = 'WTHMASS '
372          cdiag(257) = 'USLTMASS'
373          cdiag(258) = 'VSLTMASS'
374          cdiag(259) = 'WSLTMASS'
375          cdiag(260) = 'UVELTH  '
376          cdiag(261) = 'VVELTH  '
377          cdiag(262) = 'WVELTH  '
378          cdiag(263) = 'UVELSLT '
379          cdiag(264) = 'VVELSLT '
380          cdiag(265) = 'WVELSLT '
381          cdiag(266) = 'UTRAC1  '
382          cdiag(267) = 'VTRAC1  '
383          cdiag(268) = 'WTRAC1  '
384          cdiag(269) = 'UTRAC2  '
385          cdiag(270) = 'VTRAC2  '
386          cdiag(271) = 'WTRAC2  '
387          cdiag(272) = 'UTRAC3  '
388          cdiag(273) = 'VTRAC3  '
389          cdiag(274) = 'WTRAC3  '
390          cdiag(275) = 'UTRAC4  '
391          cdiag(276) = 'VTRAC4  '
392          cdiag(277) = 'WTRAC4  '
393          cdiag(278) = 'UTRAC5  '
394          cdiag(279) = 'VTRAC5  '
395          cdiag(280) = 'WTRAC5  '
396          cdiag(281) = 'TRAC1   '
397          cdiag(282) = 'TRAC2   '
398          cdiag(283) = 'TRAC3   '
399          cdiag(284) = 'TRAC4   '
400          cdiag(285) = 'TRAC5   '
401          cdiag(286) = 'DICBIOA '
402          cdiag(287) = 'DICCARB '
403          cdiag(288) = 'DICTFLX '
404          cdiag(289) = 'DICOFLX '
405          cdiag(290) = 'DICCFLX '
406          cdiag(291) = 'DICPCO2 '
407          cdiag(292) = 'DICPHAV '
408          cdiag(293) = 'DTCONV  '
409          cdiag(294) = 'DQCONV  '
410          cdiag(295) = 'RELHUM  '
411          cdiag(296) = 'PRECLS  '
412          cdiag(297) = 'ENPREC  '
413    
414  c Diagnostic Descriptions  c Diagnostic Descriptions
415  c ----------------------------------------------------------------------  c ----------------------------------------------------------------------
# Line 611  c -------------------------------------- Line 723  c --------------------------------------
723       . 'U-Wind Tendency due to Gravity Wave Drag w/Clouds (m/s/day)'       . 'U-Wind Tendency due to Gravity Wave Drag w/Clouds (m/s/day)'
724        tdiag(237) =        tdiag(237) =
725       . 'V-Wind Tendency due to Gravity Wave Drag w/Clouds (m/s/day)'       . 'V-Wind Tendency due to Gravity Wave Drag w/Clouds (m/s/day)'
726          tdiag(238) = 'Perturbation of Surface (pressure, height) (hPa,m)'
727          tdiag(239) = 'Square of Perturbation of Sfc (hPa^2,m^2)'
728          tdiag(240)  = 'Potential Temperature (K)'
729          tdiag(241)  = 'Salt or Water Vapor mixing ratio (g/kg)'
730          tdiag(242)  = 'Zonal Component of Velocity (m/s)'
731          tdiag(243)  = 'Meridional Component of Velocity (m/s)'
732          tdiag(244)  = 'Vertical Component of Velocity (m/s)'
733          tdiag(245)  = 'Square of Potential Temperature (K^2)'
734          tdiag(246)  = 'Square of Salt or Water Vapor mix rat (g^2/kg^2)'
735          tdiag(247)  = 'Square of Zonal Comp of Velocity (m^2/s^2)'
736          tdiag(248)  = 'Square of Meridional Comp of Velocity (m^2/s^2)'
737          tdiag(249)  = 'Square of Vertical Comp of Velocity (m^2/s^2)'
738          tdiag(250)  = 'Meridional Transport of Zonal Momentum (m^2/s^2)'
739          tdiag(251)  = 'Zonal Mass-Weighted Comp of Velocity (m/s)'
740          tdiag(252)  = 'Meridional Mass-Weighted Comp of Velocity (m/s)'
741          tdiag(253)  = 'Vertical Mass-Weighted Comp of Velocity (m/s)'
742          tdiag(254)  = 'Zonal Mass-Weight Transp of Pot Temp (m-K/s)'
743          tdiag(255)  = 'Meridional Mass-Weight Transp of Pot Temp (m-K/s)'
744          tdiag(256)  = 'Vertical Mass-Weight Transp of Pot Temp (m-K/s)'
745          tdiag(257)  = 'Zonal Mass-Weight Transp of Salt (m-kg/s-kg)'
746          tdiag(258)  = 'Meridional Mass-Weight Transp of Salt (m-kg/s-kg)'
747          tdiag(259)  = 'Vertical Mass-Weight Transp of Salt (m-kg/s-kg)'
748          tdiag(260)  = 'Zonal Transp of Pot Temp (m-K/s)'
749          tdiag(261)  = 'Meridional Transp of Pot Temp (m-K/s)'
750          tdiag(262)  = 'Vertical Transp of Pot Temp (m-K/s)'
751          tdiag(263)  = 'Zonal Transp of Salt (m-kg/s-kg)'
752          tdiag(264)  = 'Meridional Transp of Salt (m-kg/s-kg)'
753          tdiag(265)  = 'Vertical Transp of Salt (m-kg/s-kg)'
754          tdiag(266)  = 'Zonal Mass-Weighted Transp of Tracer 1 (m-kg/s-kg)'
755          tdiag(267)  = 'Merid Mass-Weighted Transp of Tracer 1 (m-kg/s-kg)'
756          tdiag(268)  = 'Vert  Mass-Weighted Transp of Tracer 1 (m-kg/s-kg)'
757          tdiag(269)  = 'Zonal Mass-Weighted Transp of Tracer 2 (m-kg/s-kg)'
758          tdiag(270)  = 'Merid Mass-Weighted Transp of Tracer 2 (m-kg/s-kg)'
759          tdiag(271)  = 'Vert  Mass-Weighted Transp of Tracer 2 (m-kg/s-kg)'
760          tdiag(272)  = 'Zonal Mass-Weighted Transp of Tracer 3 (m-kg/s-kg)'
761          tdiag(273)  = 'Merid Mass-Weighted Transp of Tracer 3 (m-kg/s-kg)'
762          tdiag(274)  = 'Vert  Mass-Weighted Transp of Tracer 3 (m-kg/s-kg)'
763          tdiag(275)  = 'Zonal Mass-Weighted Transp of Tracer 4 (m-kg/s-kg)'
764          tdiag(276)  = 'Merid Mass-Weighted Transp of Tracer 4 (m-kg/s-kg)'
765          tdiag(277)  = 'Vert  Mass-Weighted Transp of Tracer 4 (m-kg/s-kg)'
766          tdiag(278)  = 'Zonal Mass-Weighted Transp of Tracer 5 (m-kg/s-kg)'
767          tdiag(279)  = 'Merid Mass-Weighted Transp of Tracer 5 (m-kg/s-kg)'
768          tdiag(280)  = 'Vert  Mass-Weighted Transp of Tracer 5 (m-kg/s-kg)'
769          tdiag(281)  = 'Mass-Weighted Tracer 1 ()'
770          tdiag(282)  = 'Mass-Weighted Tracer 2 ()'
771          tdiag(283)  = 'Mass-Weighted Tracer 3 ()'
772          tdiag(284)  = 'Mass-Weighted Tracer 4 ()'
773          tdiag(285)  = 'Mass-Weighted Tracer 5 ()'
774          tdiag(286)  = 'Biological Productivity (mol/m3/s)'
775          tdiag(287)  = 'Carbonate chg-biol prod and remin (mol eq/m3/s)'
776          tdiag(288)  = 'Tendency of DIC due to air-sea exch (mol/m3/s)'
777          tdiag(289)  = 'Tendency of O2 due to air-sea exch (mol/m3/s)'
778          tdiag(290)  = 'Flux of CO2 - air-sea exch (mol/m2/s)'
779          tdiag(291)  = 'Partial Pressure of CO2 (atm)'
780          tdiag(292)  = 'pH (dimensionless)'
781          tdiag(293)  = 'Temp Change due to Convection (K/sec)'
782          tdiag(294)  = 'Humidity Change due to Convection (g/kg/sec)'
783          tdiag(295)  = 'Relative Humidity (fraction)'
784          tdiag(296)  = 'Large Scale Precipitation (g/m2/sec)'
785          tdiag(297)  = 'Energy of Precipitation (snow, rain Temp) in (J/g)'
786    
787  c Diagnostic Units  c Diagnostic Units
788  c ----------------  c ----------------
# Line 851  c ---------------- Line 1023  c ----------------
1023        udiag(235) = 'N/m2            '        udiag(235) = 'N/m2            '
1024        udiag(236) = 'm/sec/day       '        udiag(236) = 'm/sec/day       '
1025        udiag(237) = 'm/sec/day       '        udiag(237) = 'm/sec/day       '
1026          udiag(238) = 'm or hPa        '
1027          udiag(239) = 'm2 or hPa2      '
1028          udiag(240) = 'K               '
1029          udiag(241) = 'g/kg            '
1030          udiag(242) = 'm/sec           '
1031          udiag(243) = 'm/sec           '
1032          udiag(244) = 'm/sec           '
1033          udiag(245) = 'K2              '
1034          udiag(246) = 'g2/kg2          '
1035          udiag(247) = 'm2/sec2         '
1036          udiag(248) = 'm2/sec2         '
1037          udiag(249) = 'm2/sec2         '
1038          udiag(250) = 'm2/sec2         '
1039          udiag(251) = 'm/sec           '
1040          udiag(252) = 'm/sec           '
1041          udiag(253) = 'm/sec           '
1042          udiag(254) = 'm-K/sec         '
1043          udiag(255) = 'm-K/sec         '
1044          udiag(256) = 'm-K/sec         '
1045          udiag(257) = 'm-kg/s-kg       '
1046          udiag(258) = 'm-kg/s-kg       '
1047          udiag(259) = 'm-kg/s-kg       '
1048          udiag(260) = 'm-K/sec         '
1049          udiag(261) = 'm-K/sec         '
1050          udiag(262) = 'm-K/sec         '
1051          udiag(263) = 'm-kg/s-kg       '
1052          udiag(264) = 'm-kg/s-kg       '
1053          udiag(265) = 'm-kg/s-kg       '
1054          udiag(266) = 'm-kg/s-kg       '
1055          udiag(267) = 'm-kg/s-kg       '
1056          udiag(268) = 'm-kg/s-kg       '
1057          udiag(269) = 'm-kg/s-kg       '
1058          udiag(270) = 'm-kg/s-kg       '
1059          udiag(271) = 'm-kg/s-kg       '
1060          udiag(272) = 'm-kg/s-kg       '
1061          udiag(273) = 'm-kg/s-kg       '
1062          udiag(274) = 'm-kg/s-kg       '
1063          udiag(275) = 'm-kg/s-kg       '
1064          udiag(276) = 'm-kg/s-kg       '
1065          udiag(277) = 'm-kg/s-kg       '
1066          udiag(278) = 'm-kg/s-kg       '
1067          udiag(279) = 'm-kg/s-kg       '
1068          udiag(280) = 'm-kg/s-kg       '
1069          udiag(281) = 'kg/kg           '
1070          udiag(282) = 'kg/kg           '
1071          udiag(283) = 'kg/kg           '
1072          udiag(284) = 'kg/kg           '
1073          udiag(285) = 'kg/kg           '
1074          udiag(286) = 'mol/m3/sec      '
1075          udiag(287) = 'mol eq/m3/sec   '
1076          udiag(288) = 'mol/m3/sec      '
1077          udiag(289) = 'mol/m3/sec      '
1078          udiag(290) = 'mol/m2/sec      '
1079          udiag(291) = 'atm             '
1080          udiag(292) = 'dimensionless   '
1081          udiag(293) = 'k/sec           '
1082          udiag(294) = 'g/kg/sec        '
1083          udiag(295) = 'Fraction        '
1084          udiag(296) = 'g/m2/sec        '
1085          udiag(297) = 'J/g             '
1086    
1087  c Diagnostic Parms  c Diagnostic Parms
1088  c ----------------  c ----------------  123456789
1089        gdiag( 1)  = 'UMG  002L       '        gdiag( 1)  = 'UM   002L1      '
1090        gdiag( 2)  = 'VMG  001L       '        gdiag( 2)  = 'VM   001L1      '
1091        gdiag( 3)  = 'SMG     L       '        gdiag( 3)  = 'SM      L1      '
1092        gdiag( 4)  = 'SMG     L       '        gdiag( 4)  = 'SM      L1      '
1093        gdiag( 5)  = 'SMG     L       '        gdiag( 5)  = 'SM      L1      '
1094        gdiag( 6)  = 'SMG     L       '        gdiag( 6)  = 'SM      L1      '
1095        gdiag( 7)  = 'SMGP    L       '        gdiag( 7)  = 'SM P    L1      '
1096        gdiag( 8)  = 'SMG     L       '        gdiag( 8)  = 'SM      LL      '
1097        gdiag( 9)  = 'SMG     L       '        gdiag( 9)  = 'SM      L1      '
1098        gdiag(10)  = 'SMG     L       '        gdiag(10)  = 'SM      L1      '
1099        gdiag(11)  = 'SMG     L       '        gdiag(11)  = 'SM      LL      '
1100        gdiag(12)  = 'SMG     L       '        gdiag(12)  = 'SM      LL      '
1101        gdiag(13)  = 'UMG  014M       '        gdiag(13)  = 'UM   014ML      '
1102        gdiag(14)  = 'VMG  013M       '        gdiag(14)  = 'VM   013ML      '
1103        gdiag(15)  = 'SMG     M       '        gdiag(15)  = 'SM      ML      '
1104        gdiag(16)  = 'SMG     M       '        gdiag(16)  = 'SM      ML      '
1105        gdiag(17)  = 'SMG     M       '        gdiag(17)  = 'SM      ML      '
1106        gdiag(18)  = 'SMG     M       '        gdiag(18)  = 'SM      ML      '
1107        gdiag(19)  = 'SM      M       '        gdiag(19)  = 'SM      ML      '
1108        gdiag(20)  = 'SM P    M       '        gdiag(20)  = 'SM P    ML      '
1109        gdiag(21)  = 'SMGP    L       '        gdiag(21)  = 'SM P    L1      '
1110        gdiag(22)  = 'SMGP    L       '        gdiag(22)  = 'SM P    L1      '
1111        gdiag(23)  = 'UMG  024L       '        gdiag(23)  = 'UM   024LL      '
1112        gdiag(24)  = 'VMG  023L       '        gdiag(24)  = 'VM   023LL      '
1113        gdiag(25)  = 'SMG     L       '        gdiag(25)  = 'SM      LL      '
1114        gdiag(26)  = 'SMG     L       '        gdiag(26)  = 'SM      LL      '
1115        gdiag(27)  = 'SMG     L       '        gdiag(27)  = 'SM      L1      '
1116        gdiag(28)  = 'SMG     L       '        gdiag(28)  = 'SM      L1      '
1117        gdiag(29)  = 'SMG     L       '        gdiag(29)  = 'SM      L1      '
1118        gdiag(30)  = 'SMG     L       '        gdiag(30)  = 'SM      L1      '
1119        gdiag(31)  = 'SMG     L       '        gdiag(31)  = 'SM      L1      '
1120        gdiag(32)  = 'SMG     L       '        gdiag(32)  = 'SM      L1      '
1121        gdiag(33)  = 'SMGP    L       '        gdiag(33)  = 'SM P    L1      '
1122        gdiag(34)  = 'SMGP    L       '        gdiag(34)  = 'SM P    L1      '
1123        gdiag(35)  = 'SMG     L       '        gdiag(35)  = 'SM      L1      '
1124        gdiag(36)  = 'SMG     L       '        gdiag(36)  = 'SM      L1      '
1125        gdiag(37)  = 'SMG     U       '        gdiag(37)  = 'SM      U1      '
1126        gdiag(38)  = 'SMG     U       '        gdiag(38)  = 'SM      U1      '
1127        gdiag(39)  = 'SMG     L       '        gdiag(39)  = 'SM      L1      '
1128        gdiag(40)  = 'SM      M       '        gdiag(40)  = 'SM      ML      '
1129        gdiag(41)  = 'SMG     M       '        gdiag(41)  = 'SM      ML      '
1130        gdiag(42)  = 'SMGP    M       '        gdiag(42)  = 'SM P    ML      '
1131        gdiag(43)  = 'SMG     M       '        gdiag(43)  = 'SM      ML      '
1132        gdiag(44)  = 'SMGP    M       '        gdiag(44)  = 'SM P    ML      '
1133        gdiag(45)  = 'SMGP    M       '        gdiag(45)  = 'SM P    ML      '
1134        gdiag(46)  = 'SMG     L       '        gdiag(46)  = 'SM      L1      '
1135        gdiag(47)  = 'SMG     M       '        gdiag(47)  = 'SM      ML      '
1136        gdiag(48)  = 'SMGP    U       '        gdiag(48)  = 'SM P    U1      '
1137        gdiag(49)  = 'SMG C050M       '        gdiag(49)  = 'SM  C050ML      '
1138        gdiag(50)  = 'SMG D   M       '        gdiag(50)  = 'SM  D   ML      '
1139        gdiag(51)  = 'SMGP    M       '        gdiag(51)  = 'SM P    M1      '
1140        gdiag(52)  = 'SMG     L       '        gdiag(52)  = 'SM      L1      '
1141        gdiag(53)  = 'SM      L       '        gdiag(53)  = 'SM      L1      '
1142        gdiag(54)  = 'UU   055M       '        gdiag(54)  = 'UU   055MR      '
1143        gdiag(55)  = 'VV   054M       '        gdiag(55)  = 'VV   054MR      '
1144        gdiag(56)  = 'SM      M       '        gdiag(56)  = 'SM      MR      '
1145        gdiag(57)  = 'SM      M       '        gdiag(57)  = 'SM      MR      '
1146        gdiag(58)  = 'SM      M       '        gdiag(58)  = 'SM      ML      '
1147        gdiag(59)  = 'UU   060M       '        gdiag(59)  = 'UU   060ML      '
1148        gdiag(60)  = 'VV   059M       '        gdiag(60)  = 'VV   059ML      '
1149        gdiag(61)  = 'SM      M       '        gdiag(61)  = 'SM      ML      '
1150        gdiag(62)  = 'SM      M       '        gdiag(62)  = 'SM      ML      '
1151        gdiag(63)  = 'SM      M       '        gdiag(63)  = 'SM      ML      '
1152        gdiag(64)  = 'SM      L       '        gdiag(64)  = 'SM      LL      '
1153        gdiag(65)  = 'SMG     M       '        gdiag(65)  = 'SM      ML      '
1154        gdiag(66)  = 'SMG     M       '        gdiag(66)  = 'SM      ML      '
1155        gdiag(67)  = 'SMG     L       '        gdiag(67)  = 'SM      L1      '
1156        gdiag(68)  = 'SMG     L       '        gdiag(68)  = 'SM      L1      '
1157        gdiag(69)  = 'SMGP    M       '        gdiag(69)  = 'SM P    MM      '
1158        gdiag(70)  = 'SMGP    M       '        gdiag(70)  = 'SM P    M1      '
1159        gdiag(71)  = 'SM P    M       '        gdiag(71)  = 'SM P    ML      '
1160        gdiag(72)  = 'SMGP    U       '        gdiag(72)  = 'SM P    U1      '
1161        gdiag(73)  = 'SMGP    U       '        gdiag(73)  = 'SM P    U1      '
1162        gdiag(74)  = 'SMGP    L       '        gdiag(74)  = 'SM P    LL      '
1163        gdiag(75)  = 'UU   076M       '        gdiag(75)  = 'UU   076ML      '
1164        gdiag(76)  = 'VV   075M       '        gdiag(76)  = 'VV   075ML      '
1165        gdiag(77)  = 'SM      M       '        gdiag(77)  = 'SM      ML      '
1166        gdiag(78)  = 'SM P    M       '        gdiag(78)  = 'SM P    ML      '
1167        gdiag(79)  = 'SM      M       '        gdiag(79)  = 'SM      ML      '
1168        gdiag(80)  = 'SM      L       '        gdiag(80)  = 'SM      L1      '
1169        gdiag(81)  = 'SM P    L       '        gdiag(81)  = 'SM P    LL      '
1170        gdiag(82)  = 'SMGP    L       '        gdiag(82)  = 'SM P    L1      '
1171        gdiag(83)  = 'SM      L       '        gdiag(83)  = 'SM      L1      '
1172        gdiag(84)  = 'SMG     L       '        gdiag(84)  = 'SM      L1      '
1173        gdiag(85)  = 'SMG     L       '        gdiag(85)  = 'SM      L1      '
1174        gdiag(86)  = 'SMG     M       '        gdiag(86)  = 'SM      ML      '
1175        gdiag(87)  = 'SMG     M       '        gdiag(87)  = 'SM      ML      '
1176        gdiag(88)  = 'UU   089M       '        gdiag(88)  = 'UU   089ML      '
1177        gdiag(89)  = 'VV   088M       '        gdiag(89)  = 'VV   088ML      '
1178        gdiag(90)  = 'SM      M       '        gdiag(90)  = 'SM      ML      '
1179        gdiag(91)  = 'SM      M       '        gdiag(91)  = 'SM      ML      '
1180        gdiag(92)  = 'UU   093M       '        gdiag(92)  = 'UU   093ML      '
1181        gdiag(93)  = 'VV   092M       '        gdiag(93)  = 'VV   092ML      '
1182        gdiag(94)  = 'UMG  095M       '        gdiag(94)  = 'UM   095ML      '
1183        gdiag(95)  = 'VMG  094M       '        gdiag(95)  = 'VM   094ML      '
1184        gdiag(96)  = 'UMG  097L       '        gdiag(96)  = 'UM   097L1      '
1185        gdiag(97)  = 'VMG  096L       '        gdiag(97)  = 'VM   096L1      '
1186        gdiag(98)  = 'UMG  099U       '        gdiag(98)  = 'UM   099U1      '
1187        gdiag(99)  = 'VMG  098U       '        gdiag(99)  = 'VM   098U1      '
1188        gdiag(100) = 'SMG     M       '        gdiag(100) = 'SM      ML      '
1189        gdiag(101) = 'SM      L       '        gdiag(101) = 'SM      L1      '
1190        gdiag(102) = 'UM   103M       '        gdiag(102) = 'UM   103M1      '
1191        gdiag(103) = 'VM   102M       '        gdiag(103) = 'VM   102M1      '
1192        gdiag(104) = 'UM   105M       '        gdiag(104) = 'UM   105M1      '
1193        gdiag(105) = 'VM   104M       '        gdiag(105) = 'VM   104M1      '
1194        gdiag(106) = 'SMGP    M       '        gdiag(106) = 'SM P    M1      '
1195        gdiag(107) = 'SM P    M       '        gdiag(107) = 'SM P    M1      '
1196        gdiag(108) = 'UMG  109M       '        gdiag(108) = 'UM   109M1      '
1197        gdiag(109) = 'VMG  108M       '        gdiag(109) = 'VM   108M1      '
1198        gdiag(110) = 'SMG     M       '        gdiag(110) = 'SM      M1      '
1199        gdiag(111) = 'SMGP    M       '        gdiag(111) = 'SM P    M1      '
1200        gdiag(112) = 'UMG  113M       '        gdiag(112) = 'UM   113M1      '
1201        gdiag(113) = 'VMG  112M       '        gdiag(113) = 'VM   112M1      '
1202        gdiag(114) = 'SMG     M       '        gdiag(114) = 'SM      M1      '
1203        gdiag(115) = 'SMGP    M       '        gdiag(115) = 'SM P    M1      '
1204        gdiag(116) = 'SMGP    L       '        gdiag(116) = 'SM P    LL      '
1205        gdiag(117) = 'SM      M       '        gdiag(117) = 'SM      ML      '
1206        gdiag(118) = 'SM      M       '        gdiag(118) = 'SM      M1      '
1207        gdiag(119) = 'SM P    L       '        gdiag(119) = 'SM P    L1      '
1208        gdiag(120) = 'UU   121M       '        gdiag(120) = 'UU   121ML      '
1209        gdiag(121) = 'VV   120M       '        gdiag(121) = 'VV   120ML      '
1210        gdiag(122) = 'SM      M       '        gdiag(122) = 'SM      ML      '
1211        gdiag(123) = 'SM      M       '        gdiag(123) = 'SM      ML      '
1212        gdiag(124) = 'SMG     L       '        gdiag(124) = 'SM      L1      '
1213        gdiag(125) = 'SMG     L       '        gdiag(125) = 'SM      L1      '
1214        gdiag(126) = 'SMG     L       '        gdiag(126) = 'SM      L1      '
1215        gdiag(127) = 'SMG     L       '        gdiag(127) = 'SM      L1      '
1216        gdiag(128) = 'SMG     L       '        gdiag(128) = 'SM      L1      '
1217        gdiag(129) = 'SMG     L       '        gdiag(129) = 'SM      L1      '
1218        gdiag(130) = 'SMG     L       '        gdiag(130) = 'SM      L1      '
1219        gdiag(131) = 'SMG     L       '        gdiag(131) = 'SM      L1      '
1220        gdiag(132) = 'SMG     M       '        gdiag(132) = 'SM      ML      '
1221        gdiag(133) = 'SMG     M       '        gdiag(133) = 'SM      ML      '
1222        gdiag(134) = 'SMG     M       '        gdiag(134) = 'SM      ML      '
1223        gdiag(135) = 'SMG     M       '        gdiag(135) = 'SM      ML      '
1224        gdiag(136) = 'SMG     M       '        gdiag(136) = 'SM      ML      '
1225        gdiag(137) = 'SMG     M       '        gdiag(137) = 'SM      ML      '
1226        gdiag(138) = 'SMG     M       '        gdiag(138) = 'SM      ML      '
1227        gdiag(139) = 'SMG     M       '        gdiag(139) = 'SM      ML      '
1228        gdiag(140) = 'SMGP    M       '        gdiag(140) = 'SM P    M1      '
1229        gdiag(141) = 'SMGP    M       '        gdiag(141) = 'SM P    M1      '
1230        gdiag(142) = 'SMGPC150M       '        gdiag(142) = 'SM PC150M1      '
1231        gdiag(143) = 'SMGPC151M       '        gdiag(143) = 'SM PC151M1      '
1232        gdiag(144) = 'SMGPC152M       '        gdiag(144) = 'SM PC152M1      '
1233        gdiag(145) = 'SMG     M       '        gdiag(145) = 'SM      ML      '
1234        gdiag(146) = 'SMGPC147M       '        gdiag(146) = 'SM PC147M1      '
1235        gdiag(147) = 'SMG D   M       '        gdiag(147) = 'SM  D   M1      '
1236        gdiag(148) = 'SMGPC149M       '        gdiag(148) = 'SM PC149M1      '
1237        gdiag(149) = 'SMG D   M       '        gdiag(149) = 'SM  D   M1      '
1238        gdiag(150) = 'SMG D   M       '        gdiag(150) = 'SM  D   M1      '
1239        gdiag(151) = 'SMG D   M       '        gdiag(151) = 'SM  D   M1      '
1240        gdiag(152) = 'SMG D   M       '        gdiag(152) = 'SM  D   M1      '
1241        gdiag(153) = 'SMG     L       '        gdiag(153) = 'SM      L1      '
1242        gdiag(154) = 'SMG     L       '        gdiag(154) = 'SM      L1      '
1243        gdiag(155) = 'SMG     L       '        gdiag(155) = 'SM      L1      '
1244        gdiag(156) = 'SMG     L       '        gdiag(156) = 'SM      L1      '
1245        gdiag(157) = 'SMG     L       '        gdiag(157) = 'SM      L1      '
1246        gdiag(158) = 'SMG     L       '        gdiag(158) = 'SM      L1      '
1247        gdiag(159) = 'SMG     L       '        gdiag(159) = 'SM      L1      '
1248        gdiag(160) = 'SMG     L       '        gdiag(160) = 'SM      L1      '
1249        gdiag(161) = 'SMG     L       '        gdiag(161) = 'SM      L1      '
1250        gdiag(162) = 'SMG     L       '        gdiag(162) = 'SM      L1      '
1251        gdiag(163) = 'SMGP    L       '        gdiag(163) = 'SM P    L1      '
1252        gdiag(164) = 'SMG     L       '        gdiag(164) = 'SM      L1      '
1253        gdiag(165) = 'SMG     L       '        gdiag(165) = 'SM      L1      '
1254        gdiag(166) = 'SMG     L       '        gdiag(166) = 'SM      L1      '
1255        gdiag(167) = 'SMG     L       '        gdiag(167) = 'SM      L1      '
1256        gdiag(168) = 'SMG     L       '        gdiag(168) = 'SM      L1      '
1257        gdiag(169) = 'SMG     L       '        gdiag(169) = 'SM      L1      '
1258        gdiag(170) = 'SMG     L       '        gdiag(170) = 'SM      L1      '
1259        gdiag(171) = 'SMG     L       '        gdiag(171) = 'SM      L1      '
1260        gdiag(172) = 'SMGP    L       '        gdiag(172) = 'SM P    L1      '
1261        gdiag(173) = 'SMGP    L       '        gdiag(173) = 'SM P    L1      '
1262        gdiag(174) = 'SMG     L       '        gdiag(174) = 'SM      L1      '
1263        gdiag(175) = 'SMG     L       '        gdiag(175) = 'SM      L1      '
1264        gdiag(176) = 'SMG     L       '        gdiag(176) = 'SM      L1      '
1265        gdiag(177) = 'SMG     L       '        gdiag(177) = 'SM      L1      '
1266        gdiag(178) = 'SMG     L       '        gdiag(178) = 'SM      L1      '
1267        gdiag(179) = 'SM      M       '        gdiag(179) = 'SM      M1      '
1268        gdiag(180) = 'SM      M       '        gdiag(180) = 'SM      M1      '
1269        gdiag(181) = 'SM      M       '        gdiag(181) = 'SM      M1      '
1270        gdiag(182) = 'SM P    M       '        gdiag(182) = 'SM P    M1      '
1271        gdiag(183) = 'SM      M       '        gdiag(183) = 'SM      M1      '
1272        gdiag(184) = 'SMGPC185M       '        gdiag(184) = 'SM PC185M1      '
1273        gdiag(185) = 'SMGPD   M       '        gdiag(185) = 'SM PD   M1      '
1274        gdiag(186) = 'SMG     L       '        gdiag(186) = 'SM      L1      '
1275        gdiag(187) = 'SMG C188L       '        gdiag(187) = 'SM  C188L1      '
1276        gdiag(188) = 'SMG D   L       '        gdiag(188) = 'SM  D   L1      '
1277        gdiag(189) = '                '        gdiag(189) = '         0      '
1278        gdiag(190) = 'SMG     L       '        gdiag(190) = 'SM      L1      '
1279        gdiag(191) = 'SM P    L       '        gdiag(191) = 'SM P    L1      '
1280        gdiag(192) = 'SMG     L       '        gdiag(192) = 'SM      L1      '
1281        gdiag(193) = 'SMG     L       '        gdiag(193) = 'SM      L1      '
1282        gdiag(194) = 'SMG     L       '        gdiag(194) = 'SM      L1      '
1283        gdiag(195) = 'SMG     L       '        gdiag(195) = 'SM      L1      '
1284        gdiag(196) = 'SMG     L       '        gdiag(196) = 'SM      L1      '
1285        gdiag(197) = 'SMG     L       '        gdiag(197) = 'SM      L1      '
1286        gdiag(198) = 'SMG     L       '        gdiag(198) = 'SM      L1      '
1287        gdiag(199) = 'SMG     L       '        gdiag(199) = 'SM      L1      '
1288        gdiag(200) = 'SMG     L       '        gdiag(200) = 'SM      L1      '
1289        gdiag(201) = 'SMG     L       '        gdiag(201) = 'SM      L1      '
1290        gdiag(202) = 'SMG     M       '        gdiag(202) = 'SM      ML      '
1291        gdiag(203) = 'SMG     M       '        gdiag(203) = 'SM      ML      '
1292        gdiag(204) = 'SMG C210L       '        gdiag(204) = 'SM  C210L1      '
1293        gdiag(205) = 'SMG     L       '        gdiag(205) = 'SM      L1      '
1294        gdiag(206) = 'SMG     L       '        gdiag(206) = 'SM      L1      '
1295        gdiag(207) = 'SMG     L       '        gdiag(207) = 'SM      L1      '
1296        gdiag(208) = 'SMG     L       '        gdiag(208) = 'SM      L1      '
1297        gdiag(209) = 'SMG     M       '        gdiag(209) = 'SM      ML      '
1298        gdiag(210) = 'SMG D   L       '        gdiag(210) = 'SM  D   L1      '
1299        gdiag(211) = 'SM      M       '        gdiag(211) = 'SM      ML      '
1300        gdiag(212) = 'SM      M       '        gdiag(212) = 'SM      ML      '
1301        gdiag(213) = 'SMGP    M       '        gdiag(213) = 'SM P    ML      '
1302        gdiag(214) = 'SMGP    L       '        gdiag(214) = 'SM P    L1      '
1303        gdiag(215) = 'SMGP    L       '        gdiag(215) = 'SM P    L1      '
1304        gdiag(216) = 'SMGP    L       '        gdiag(216) = 'SM P    L1      '
1305        gdiag(217) = 'SMGP    L       '        gdiag(217) = 'SM P    L1      '
1306        gdiag(218) = 'SMGP    L       '        gdiag(218) = 'SM P    L1      '
1307        gdiag(219) = 'SMGP    L       '        gdiag(219) = 'SM P    L1      '
1308        gdiag(220) = 'SMGP    M       '        gdiag(220) = 'SM P    ML      '
1309        gdiag(221) = 'SMG     M       '        gdiag(221) = 'SM      ML      '
1310        gdiag(222) = 'SMG     M       '        gdiag(222) = 'SM      ML      '
1311        gdiag(223) = 'SMG     M       '        gdiag(223) = 'SM      ML      '
1312        gdiag(224) = 'SMG     M       '        gdiag(224) = 'SM      ML      '
1313        gdiag(225) = 'SMG     M       '        gdiag(225) = 'SM      ML      '
1314        gdiag(226) = 'SMG     M       '        gdiag(226) = 'SM      ML      '
1315        gdiag(227) = 'SMG     M       '        gdiag(227) = 'SM      ML      '
1316        gdiag(228) = 'SMG     M       '        gdiag(228) = 'SM      ML      '
1317        gdiag(229) = 'UMG  230L       '        gdiag(229) = 'UM   230L1      '
1318        gdiag(230) = 'VMG  229L       '        gdiag(230) = 'VM   229L1      '
1319        gdiag(231) = 'SMG     L       '        gdiag(231) = 'SM      L1      '
1320        gdiag(232) = 'UMG  233L       '        gdiag(232) = 'UM   233L1      '
1321        gdiag(233) = 'VMG  232L       '        gdiag(233) = 'VM   232L1      '
1322        gdiag(234) = 'UMG  235L       '        gdiag(234) = 'UM   235L1      '
1323        gdiag(235) = 'VMG  234L       '        gdiag(235) = 'VM   234L1      '
1324        gdiag(236) = 'UMG  237M       '        gdiag(236) = 'UM   237ML      '
1325        gdiag(237) = 'VMG  236M       '        gdiag(237) = 'VM   236ML      '
1326                  gdiag(238) = 'SM      M1      '
1327  c Diagnostic Levels        gdiag(239) = 'SM      M1      '
1328  c -----------------        gdiag(240) = 'SM      MR      '
1329        KDIAG( 1) =  1            gdiag(241) = 'SM      MR      '
1330        KDIAG( 2) =  1            gdiag(242) = 'UU   243MR      '
1331        KDIAG( 3) =  1            gdiag(243) = 'VV   242MR      '
1332        KDIAG( 4) =  1            gdiag(244) = 'WM      LR      '
1333        KDIAG( 5) =  1            gdiag(245) = 'SM      MR      '
1334        KDIAG( 6) =  1            gdiag(246) = 'SM      MR      '
1335        KDIAG( 7) =  1            gdiag(247) = 'UU   248MR      '
1336        KDIAG( 8) = nrphys          gdiag(248) = 'VV   247MR      '
1337        KDIAG( 9) =  1            gdiag(249) = 'WM      LR      '
1338        KDIAG(10) =  1            gdiag(250) = 'UZ   250MR      '
1339        KDIAG(11) = nrphys        gdiag(251) = 'UU   252MR      '
1340        KDIAG(12) = nrphys          gdiag(252) = 'VV   251MR      '
1341        KDIAG(13) = nrphys          gdiag(253) = 'WM      MR      '
1342        KDIAG(14) = nrphys          gdiag(254) = 'UU   255MR      '
1343        KDIAG(15) = nrphys          gdiag(255) = 'VV   254MR      '
1344        KDIAG(16) = nrphys          gdiag(256) = 'WM      MR      '
1345        KDIAG(17) = nrphys          gdiag(257) = 'UU   258MR      '
1346        KDIAG(18) = nrphys          gdiag(258) = 'VV   257MR      '
1347        KDIAG(19) = nrphys          gdiag(259) = 'WM      MR      '
1348        KDIAG(20) = nrphys          gdiag(260) = 'UU   261MR      '
1349        KDIAG(21) =  1            gdiag(261) = 'VV   260MR      '
1350        KDIAG(22) =  1            gdiag(262) = 'WM      MR      '
1351        KDIAG(23) = nrphys          gdiag(263) = 'UU   264MR      '
1352        KDIAG(24) = nrphys          gdiag(264) = 'VV   263MR      '
1353        KDIAG(25) = nrphys          gdiag(265) = 'WM      MR      '
1354        KDIAG(26) = nrphys          gdiag(266) = 'UU   267MR      '
1355        KDIAG(27) =  1            gdiag(267) = 'VV   266MR      '
1356        KDIAG(28) =  1            gdiag(268) = 'WM      MR      '
1357        KDIAG(29) =  1            gdiag(269) = 'UU   270MR      '
1358        KDIAG(30) =  1            gdiag(270) = 'VV   269MR      '
1359        KDIAG(31) =  1            gdiag(271) = 'WM      MR      '
1360        KDIAG(32) =  1            gdiag(272) = 'UU   273MR      '
1361        KDIAG(33) =  1            gdiag(273) = 'VV   272MR      '
1362        KDIAG(34) =  1            gdiag(274) = 'WM      MR      '
1363        KDIAG(35) =  1            gdiag(275) = 'UU   276MR      '
1364        KDIAG(36) =  1            gdiag(276) = 'VV   275MR      '
1365        KDIAG(37) =  1            gdiag(277) = 'WM      MR      '
1366        KDIAG(38) =  1            gdiag(278) = 'UU   279MR      '
1367        KDIAG(39) =  1            gdiag(279) = 'VV   278MR      '
1368        KDIAG(40) = nrphys        gdiag(280) = 'WM      MR      '
1369        KDIAG(41) = nrphys        gdiag(281) = 'SM P    MR      '
1370        KDIAG(42) = nrphys        gdiag(282) = 'SM P    MR      '
1371        KDIAG(43) = nrphys        gdiag(283) = 'SM P    MR      '
1372        KDIAG(44) = nrphys        gdiag(284) = 'SM P    MR      '
1373        KDIAG(45) = nrphys        gdiag(285) = 'SM P    MR      '
1374        KDIAG(46) =  1        gdiag(286) = 'SM P    MR      '
1375        KDIAG(47) = nrphys        gdiag(287) = 'SM P    MR      '
1376        KDIAG(48) =  1            gdiag(288) = 'SM P    L1      '
1377        KDIAG(49) = nrphys        gdiag(289) = 'SM P    L1      '
1378        KDIAG(50) = nrphys        gdiag(290) = 'SM P    L1      '
1379        KDIAG(51) =  1            gdiag(291) = 'SM P    M1      '
1380        KDIAG(52) =  1            gdiag(292) = 'SM P    M1      '
1381        KDIAG(53) =  1            gdiag(293) = 'SM      MR      '
1382        KDIAG(54) = Nr        gdiag(294) = 'SM      MR      '
1383        KDIAG(55) = Nr        gdiag(295) = 'SM      MR      '
1384        KDIAG(56) = Nr        gdiag(296) = 'SM      M1      '
1385        KDIAG(57) = Nr        gdiag(297) = 'SM      M1      '
1386        KDIAG(58) = nrphys  
1387        KDIAG(59) = nrphys  c     Diagnostic Levels
1388        KDIAG(60) = nrphys  c     -----------------
1389        KDIAG(61) = nrphys        do n = 1,297
1390        KDIAG(62) = nrphys          IF (gdiag(n)(10:10) .eq. '0')  KDIAG(n) = 0
1391        KDIAG(63) = nrphys          IF (gdiag(n)(10:10) .eq. '1')  KDIAG(n) = 1
1392        KDIAG(64) = nrphys          IF (gdiag(n)(10:10) .eq. 'R')  KDIAG(n) = Nr
1393        KDIAG(65) = nrphys          IF (gdiag(n)(10:10) .eq. 'L')  KDIAG(n) = nlevs
1394        KDIAG(66) = nrphys          IF (gdiag(n)(10:10) .eq. 'M')  KDIAG(n) = nlevs - 1
1395        KDIAG(67) =  1            enddo
1396        KDIAG(68) =  1      
1397        KDIAG(69) = nrphys-1  #ifdef ALLOW_MNC
1398        KDIAG(70) =  1            IF (useMNC) THEN
1399        KDIAG(71) = nrphys  C       Set MNC CW Variable types
1400        KDIAG(72) =  1              do n = 1,297
1401        KDIAG(73) =  1                CALL DIAG_GET_MNC_GTYPE(cdiag(n), gdiag(n),
1402        KDIAG(74) = nrphys       &         tdiag(n), udiag(n), myThid)
1403        KDIAG(75) = nrphys          enddo
1404        KDIAG(76) = nrphys        ENDIF
1405        KDIAG(77) = nrphys  #endif /*  ALLOW_MNC  */
       KDIAG(78) = nrphys  
       KDIAG(79) = nrphys  
       KDIAG(80) =  1      
       KDIAG(81) = nrphys  
       KDIAG(82) =  1      
       KDIAG(83) =  1      
       KDIAG(84) =  1      
       KDIAG(85) =  1      
       KDIAG(86) = nrphys  
       KDIAG(87) = nrphys  
       KDIAG(88) = nrphys  
       KDIAG(89) = nrphys  
       KDIAG(90) = nrphys  
       KDIAG(91) = nrphys  
       KDIAG(92) = nrphys  
       KDIAG(93) = nrphys  
       KDIAG(94) = nrphys  
       KDIAG(95) = nrphys  
       KDIAG(96) =  1      
       KDIAG(97) =  1      
       KDIAG(98) =  1      
       KDIAG(99) =  1      
       KDIAG(100)= nrphys  
       KDIAG(101)=  1  
       KDIAG(102)=  1      
       KDIAG(103)=  1      
       KDIAG(104)=  1      
       KDIAG(105)=  1      
       KDIAG(106)=  1      
       KDIAG(107)=  1      
       KDIAG(108)=  1      
       KDIAG(109)=  1      
       KDIAG(110)=  1      
       KDIAG(111)=  1      
       KDIAG(112)=  1      
       KDIAG(113)=  1      
       KDIAG(114)=  1      
       KDIAG(115)=  1      
       KDIAG(116)= nrphys  
       KDIAG(117)= nrphys  
       KDIAG(118)=  1      
       KDIAG(119)=  1      
       KDIAG(120)= nrphys  
       KDIAG(121)= nrphys  
       KDIAG(122)= nrphys  
       KDIAG(123)= nrphys  
       KDIAG(124)=  1      
       KDIAG(125)=  1      
       KDIAG(126)=  1      
       KDIAG(127)=  1      
       KDIAG(128)=  1      
       KDIAG(129)=  1      
       KDIAG(130)=  1      
       KDIAG(131)=  1      
       KDIAG(132)= nrphys  
       KDIAG(133)= nrphys  
       KDIAG(134)= nrphys  
       KDIAG(135)= nrphys  
       KDIAG(136)= nrphys  
       KDIAG(137)= nrphys  
       KDIAG(138)= nrphys  
       KDIAG(139)= nrphys  
       KDIAG(140)=  1      
       KDIAG(141)=  1      
       KDIAG(142)=  1      
       KDIAG(143)=  1      
       KDIAG(144)=  1      
       KDIAG(145)= nrphys  
       KDIAG(146)=  1      
       KDIAG(147)=  1      
       KDIAG(148)=  1      
       KDIAG(149)=  1      
       KDIAG(150)=  1      
       KDIAG(151)=  1      
       KDIAG(152)=  1      
       KDIAG(153)=  1      
       KDIAG(154)=  1      
       KDIAG(155)=  1      
       KDIAG(156)=  1      
       KDIAG(157)=  1      
       KDIAG(158)=  1      
       KDIAG(159)=  1      
       KDIAG(160)=  1      
       KDIAG(161)=  1      
       KDIAG(162)=  1      
       KDIAG(163)=  1      
       KDIAG(164)=  1      
       KDIAG(165)=  1      
       KDIAG(166)=  1      
       KDIAG(167)=  1      
       KDIAG(168)=  1      
       KDIAG(169)=  1      
       KDIAG(170)=  1      
       KDIAG(171)=  1      
       KDIAG(172)=  1      
       KDIAG(173)=  1      
       KDIAG(174)=  1      
       KDIAG(175)=  1  
       KDIAG(176)=  1  
       KDIAG(177)=  1      
       KDIAG(178)=  1      
       KDIAG(179)=  1      
       KDIAG(180)=  1      
       KDIAG(181)=  1      
       KDIAG(182)=  1      
       KDIAG(183)=  1      
       KDIAG(184)=  1      
       KDIAG(185)=  1      
       KDIAG(186)=  1      
       KDIAG(187)=  1      
       KDIAG(188)=  1      
       KDIAG(189)=  0      
       KDIAG(190)=  1      
       KDIAG(191)=  1      
       KDIAG(192)=  1      
       KDIAG(193)=  1      
       KDIAG(194)=  1      
       KDIAG(195)=  1      
       KDIAG(196)=  1      
       KDIAG(197)=  1      
       KDIAG(198)=  1      
       KDIAG(199)=  1      
       KDIAG(200)=  1      
       KDIAG(201)=  1      
       KDIAG(202)= nrphys  
       KDIAG(203)= nrphys  
       KDIAG(204)=  1      
       KDIAG(205)=  1      
       KDIAG(206)=  1      
       KDIAG(207)=  1      
       KDIAG(208)=  1      
       KDIAG(209)= nrphys  
       KDIAG(210)=  1  
       KDIAG(211)= nrphys  
       KDIAG(212)= nrphys  
       KDIAG(213)= nrphys  
       KDIAG(214)=  1  
       KDIAG(215)=  1  
       KDIAG(216)=  1  
       KDIAG(217)=  1  
       KDIAG(218)=  1  
       KDIAG(219)=  1  
       KDIAG(220)= nrphys  
       KDIAG(221)= nrphys  
       KDIAG(222)= nrphys  
       KDIAG(223)= nrphys  
       KDIAG(224)= nrphys  
       KDIAG(225)= nrphys  
       KDIAG(226)= nrphys  
       KDIAG(227)= nrphys  
       KDIAG(228)= nrphys  
       KDIAG(229)=  1  
       KDIAG(230)=  1  
       KDIAG(231)=  1  
       KDIAG(232)=  1  
       KDIAG(233)=  1  
       KDIAG(234)=  1  
       KDIAG(235)=  1  
       KDIAG(236)= nrphys  
       KDIAG(237)= nrphys  
1406    
1407        do n=1,ndiagt        do n=1,ndiagt
1408        idiag (n) = 0          idiag(n) = 0
1409        ndiag(n) = 0          ndiag(n) = 0
1410        enddo        enddo
1411    
1412        RETURN        RETURN
1413        END        END
1414    
1415    
1416    C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----|
1417    CBOP 1
1418    C     !ROUTINE: DIAG_GET_MNC_GTYPE
1419    
1420    C     !INTERFACE:
1421          SUBROUTINE DIAG_GET_MNC_GTYPE(
1422         I     name,
1423         I     gstr,
1424         I     desc,
1425         I     units,
1426         I     myThid )
1427    
1428    C     !DESCRIPTION:
1429    C     Do the mapping between GDIAG and MNC GTypes
1430    
1431    C     !USES:
1432          implicit none
1433    
1434    C     !INPUT PARAMETERS:
1435          CHARACTER*(*) name, gstr, desc, units
1436          integer myThid
1437    CEOP
1438    
1439          IF (name(1:1) .eq. ' ') THEN
1440            RETURN
1441          ENDIF
1442    
1443          IF (gstr(10:10) .eq. '0') THEN
1444            CALL MNC_CW_ADD_VNAME(name, '-_-_--__-__t', 3,4, myThid)
1445          ENDIF
1446          IF (gstr(10:10) .eq. '1') THEN
1447            CALL MNC_CW_ADD_VNAME(name, 'Cen_xy_Hn__-__t', 3,4, myThid)
1448          ENDIF
1449          IF (gstr(10:10) .eq. 'R') THEN
1450            CALL MNC_CW_ADD_VNAME(name, 'Cen_xy_Hn__C__t', 3,4, myThid)
1451          ENDIF
1452          IF (gstr(10:10) .eq. 'L') THEN
1453            CALL MNC_CW_ADD_VNAME(name, 'Cen_xy_Hn__Phys__t', 3,4, myThid)
1454          ENDIF
1455          IF (gstr(10:10) .eq. 'M') THEN
1456            CALL MNC_CW_ADD_VNAME(name, 'Cen_xy_Hn__PhysI__t', 3,4, myThid)
1457          ENDIF
1458    
1459          IF (desc(1:1) .ne. ' ') THEN
1460            CALL MNC_CW_ADD_VATTR_TEXT(name,1, 'description', desc, myThid)
1461          ENDIF
1462          IF (units(1:1) .ne. ' ') THEN
1463            CALL MNC_CW_ADD_VATTR_TEXT(name,1, 'units', units, myThid)
1464          ENDIF
1465    
1466          RETURN
1467          END
1468    

Legend:
Removed from v.1.5  
changed lines
  Added in v.1.19

  ViewVC Help
Powered by ViewVC 1.1.22