1 |
subroutine diagnostics_init_vals (myThid) |
C $Header$ |
2 |
C *********************************************************************** |
C $Name$ |
|
C ***** **** |
|
|
C ***** PURPOSE **** |
|
|
C ***** Initialize Diagnostic Names, Levels, and GPARMS **** |
|
|
C ***** **** |
|
|
C ***** **** |
|
|
C ***** Note: GDIAG is defined as character*16 and can be **** |
|
|
C ***** equivalenced to character*1 parse(16) **** |
|
|
C ***** **** |
|
|
C ***** parse(1) => 'S' for a scalar diagnostic **** |
|
|
C ***** => 'U' for a U-vector component diagnostic **** |
|
|
C ***** => 'V' for a V-vector component diagnostic **** |
|
|
C ***** **** |
|
|
C ***** parse(2) => 'U' for a C-grid U-Point **** |
|
|
C ***** => 'V' for a C-grid V-Point **** |
|
|
C ***** => 'M' for a C-grid Mass Point **** |
|
|
C ***** => 'Z' for a C-grid Corner Point **** |
|
|
C ***** **** |
|
|
C ***** parse(3) => not used at present **** |
|
|
C ***** **** |
|
|
C ***** parse(4) => 'P' if diagnostic is positive definite **** |
|
|
C ***** **** |
|
|
C ***** parse(5) => 'C' for a diagnostic with a counter array **** |
|
|
C ***** => 'D' to disable an array for output **** |
|
|
C ***** **** |
|
|
C ***** parse(6-8) => 3-digit integer corresponding to vector **** |
|
|
C ***** component mate if parse(1) = 'U' or 'V' **** |
|
|
C ***** or to the counter if parse(5)='C' **** |
|
|
C ***** **** |
|
|
C ***** parse(9) => 'U' if defined at Model-level plus 1/2 **** |
|
|
C ***** => 'M' if defined at Model-level middle **** |
|
|
C ***** => 'L' if defined at Model-level minus 1/2 **** |
|
|
C ***** **** |
|
|
C ***** **** |
|
|
C *********************************************************************** |
|
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 \begin{center} |
19 |
|
C \begin{tabular}[h]{|c|c|}\hline |
20 |
|
C \textbf{Positions} & \textbf{Characters} |
21 |
|
C & \textbf{Meanings} \\\hline |
22 |
|
C parse(1) & S & scalar \\ |
23 |
|
C & U & U-vector component \\ |
24 |
|
C & V & V-vector component \\ |
25 |
|
C parse(2) & U & C-grid U-Point \\ |
26 |
|
C & V & C-grid V-Point \\ |
27 |
|
C & M & C-grid Mass Point \\ |
28 |
|
C & Z & C-grid Corner Point \\ |
29 |
|
C parse(3) & -- & --Unused-- \\ |
30 |
|
C parse(4) & P & positive definite \\ |
31 |
|
C parse(5 ) & C & with counter array \\ |
32 |
|
C & D & disable an array for output \\ |
33 |
|
C parse(6--8) & '123' & 3-digit mate number \\ |
34 |
|
C parse(9) & U & model-level plus 1/2 \\ |
35 |
|
C & M & model-level middle \\ |
36 |
|
C & L & model-level minus 1/2 \\ |
37 |
|
C parse(10) & 0 & levels = 0 \\ |
38 |
|
C & 1 & levels = 1 \\ |
39 |
|
C & R & levels = Nr \\ |
40 |
|
C & L & levels = MAX(Nr,NrPhys) \\ |
41 |
|
C & M & levels = MAX(Nr,NrPhys) - 1 \\ |
42 |
|
C \end{tabular} |
43 |
|
C \end{center} |
44 |
|
|
45 |
|
C !USES: |
46 |
implicit none |
implicit none |
|
#include "CPP_OPTIONS.h" |
|
47 |
#include "SIZE.h" |
#include "SIZE.h" |
48 |
|
#include "EEPARAMS.h" |
49 |
|
#include "EESUPPORT.h" |
50 |
|
#include "PARAMS.h" |
51 |
|
|
52 |
#ifdef ALLOW_FIZHI |
#ifdef ALLOW_FIZHI |
53 |
#include "fizhi_SIZE.h" |
#include "fizhi_SIZE.h" |
59 |
#include "diagnostics_SIZE.h" |
#include "diagnostics_SIZE.h" |
60 |
#include "diagnostics.h" |
#include "diagnostics.h" |
61 |
|
|
62 |
|
C !INPUT PARAMETERS: |
63 |
integer myThid |
integer myThid |
64 |
integer n |
CEOP |
65 |
|
|
66 |
|
integer n, nlevs |
67 |
|
|
68 |
|
nlevs = max(Nr,Nrphys) |
69 |
|
|
70 |
do n=1,ndiagt |
do n=1,ndiagt |
71 |
kdiag (n) = 0 |
kdiag(n) = 0 |
72 |
cdiag(n) = ' ' |
cdiag(n) = ' ' |
73 |
gdiag(n) = ' ' |
gdiag(n) = ' ' |
74 |
udiag(n) = ' ' |
udiag(n) = ' ' |
75 |
tdiag(n) = ' ' |
tdiag(n) = ' ' |
76 |
enddo |
enddo |
77 |
|
|
78 |
c Diagnostic Names |
c Diagnostic Names |
357 |
cdiag(278) = 'UTRAC5 ' |
cdiag(278) = 'UTRAC5 ' |
358 |
cdiag(279) = 'VTRAC5 ' |
cdiag(279) = 'VTRAC5 ' |
359 |
cdiag(280) = 'WTRAC5 ' |
cdiag(280) = 'WTRAC5 ' |
360 |
|
cdiag(281) = 'TRAC1 ' |
361 |
|
cdiag(282) = 'TRAC2 ' |
362 |
|
cdiag(283) = 'TRAC3 ' |
363 |
|
cdiag(284) = 'TRAC4 ' |
364 |
|
cdiag(285) = 'TRAC5 ' |
365 |
|
cdiag(286) = 'DICBIOA ' |
366 |
|
cdiag(287) = 'DICCARB ' |
367 |
|
cdiag(288) = 'DICTFLX ' |
368 |
|
cdiag(289) = 'DICOFLX ' |
369 |
|
cdiag(290) = 'DICCFLX ' |
370 |
|
cdiag(291) = 'DICPCO2 ' |
371 |
|
cdiag(292) = 'DICPHAV ' |
372 |
|
cdiag(293) = 'DTCONV ' |
373 |
|
cdiag(294) = 'DQCONV ' |
374 |
|
cdiag(295) = 'RELHUM ' |
375 |
|
cdiag(296) = 'PRECLS ' |
376 |
|
cdiag(297) = 'ENPREC ' |
377 |
|
cdiag(298) = 'VISCA4 ' |
378 |
|
|
379 |
c Diagnostic Descriptions |
c Diagnostic Descriptions |
380 |
c ---------------------------------------------------------------------- |
c ---------------------------------------------------------------------- |
731 |
tdiag(278) = 'Zonal Mass-Weighted Transp of Tracer 5 (m-kg/s-kg)' |
tdiag(278) = 'Zonal Mass-Weighted Transp of Tracer 5 (m-kg/s-kg)' |
732 |
tdiag(279) = 'Merid Mass-Weighted Transp of Tracer 5 (m-kg/s-kg)' |
tdiag(279) = 'Merid Mass-Weighted Transp of Tracer 5 (m-kg/s-kg)' |
733 |
tdiag(280) = 'Vert Mass-Weighted Transp of Tracer 5 (m-kg/s-kg)' |
tdiag(280) = 'Vert Mass-Weighted Transp of Tracer 5 (m-kg/s-kg)' |
734 |
|
tdiag(281) = 'Mass-Weighted Tracer 1 ()' |
735 |
|
tdiag(282) = 'Mass-Weighted Tracer 2 ()' |
736 |
|
tdiag(283) = 'Mass-Weighted Tracer 3 ()' |
737 |
|
tdiag(284) = 'Mass-Weighted Tracer 4 ()' |
738 |
|
tdiag(285) = 'Mass-Weighted Tracer 5 ()' |
739 |
|
tdiag(286) = 'Biological Productivity (mol/m3/s)' |
740 |
|
tdiag(287) = 'Carbonate chg-biol prod and remin (mol eq/m3/s)' |
741 |
|
tdiag(288) = 'Tendency of DIC due to air-sea exch (mol/m3/s)' |
742 |
|
tdiag(289) = 'Tendency of O2 due to air-sea exch (mol/m3/s)' |
743 |
|
tdiag(290) = 'Flux of CO2 - air-sea exch (mol/m2/s)' |
744 |
|
tdiag(291) = 'Partial Pressure of CO2 (atm)' |
745 |
|
tdiag(292) = 'pH (dimensionless)' |
746 |
|
tdiag(293) = 'Temp Change due to Convection (K/sec)' |
747 |
|
tdiag(294) = 'Humidity Change due to Convection (g/kg/sec)' |
748 |
|
tdiag(295) = 'Relative Humidity (fraction)' |
749 |
|
tdiag(296) = 'Large Scale Precipitation (g/m2/sec)' |
750 |
|
tdiag(297) = 'Energy of Precipitation (snow, rain Temp) in (J/g)' |
751 |
|
tdiag(298) = 'Biharmonic Viscosity Coefficient in (m4/s) ' |
752 |
|
|
753 |
c Diagnostic Units |
c Diagnostic Units |
754 |
c ---------------- |
c ---------------- |
1032 |
udiag(278) = 'm-kg/s-kg ' |
udiag(278) = 'm-kg/s-kg ' |
1033 |
udiag(279) = 'm-kg/s-kg ' |
udiag(279) = 'm-kg/s-kg ' |
1034 |
udiag(280) = 'm-kg/s-kg ' |
udiag(280) = 'm-kg/s-kg ' |
1035 |
|
udiag(281) = 'kg/kg ' |
1036 |
|
udiag(282) = 'kg/kg ' |
1037 |
|
udiag(283) = 'kg/kg ' |
1038 |
|
udiag(284) = 'kg/kg ' |
1039 |
|
udiag(285) = 'kg/kg ' |
1040 |
|
udiag(286) = 'mol/m3/sec ' |
1041 |
|
udiag(287) = 'mol eq/m3/sec ' |
1042 |
|
udiag(288) = 'mol/m3/sec ' |
1043 |
|
udiag(289) = 'mol/m3/sec ' |
1044 |
|
udiag(290) = 'mol/m2/sec ' |
1045 |
|
udiag(291) = 'atm ' |
1046 |
|
udiag(292) = 'dimensionless ' |
1047 |
|
udiag(293) = 'k/sec ' |
1048 |
|
udiag(294) = 'g/kg/sec ' |
1049 |
|
udiag(295) = 'Fraction ' |
1050 |
|
udiag(296) = 'g/m2/sec ' |
1051 |
|
udiag(297) = 'J/g ' |
1052 |
|
udiag(298) = 'm4/sec ' |
1053 |
|
|
1054 |
c Diagnostic Parms |
c Diagnostic Parms |
1055 |
c ---------------- |
c ---------------- 123456789 |
1056 |
gdiag( 1) = 'UM 002L ' |
gdiag( 1) = 'UM 002L1 ' |
1057 |
gdiag( 2) = 'VM 001L ' |
gdiag( 2) = 'VM 001L1 ' |
1058 |
gdiag( 3) = 'SM L ' |
gdiag( 3) = 'SM L1 ' |
1059 |
gdiag( 4) = 'SM L ' |
gdiag( 4) = 'SM L1 ' |
1060 |
gdiag( 5) = 'SM L ' |
gdiag( 5) = 'SM L1 ' |
1061 |
gdiag( 6) = 'SM L ' |
gdiag( 6) = 'SM L1 ' |
1062 |
gdiag( 7) = 'SM P L ' |
gdiag( 7) = 'SM P L1 ' |
1063 |
gdiag( 8) = 'SM L ' |
gdiag( 8) = 'SM LL ' |
1064 |
gdiag( 9) = 'SM L ' |
gdiag( 9) = 'SM L1 ' |
1065 |
gdiag(10) = 'SM L ' |
gdiag(10) = 'SM L1 ' |
1066 |
gdiag(11) = 'SM L ' |
gdiag(11) = 'SM LL ' |
1067 |
gdiag(12) = 'SM L ' |
gdiag(12) = 'SM LL ' |
1068 |
gdiag(13) = 'UM 014M ' |
gdiag(13) = 'UM 014ML ' |
1069 |
gdiag(14) = 'VM 013M ' |
gdiag(14) = 'VM 013ML ' |
1070 |
gdiag(15) = 'SM M ' |
gdiag(15) = 'SM ML ' |
1071 |
gdiag(16) = 'SM M ' |
gdiag(16) = 'SM ML ' |
1072 |
gdiag(17) = 'SM M ' |
gdiag(17) = 'SM ML ' |
1073 |
gdiag(18) = 'SM M ' |
gdiag(18) = 'SM ML ' |
1074 |
gdiag(19) = 'SM M ' |
gdiag(19) = 'SM ML ' |
1075 |
gdiag(20) = 'SM P M ' |
gdiag(20) = 'SM P ML ' |
1076 |
gdiag(21) = 'SM P L ' |
gdiag(21) = 'SM P L1 ' |
1077 |
gdiag(22) = 'SM P L ' |
gdiag(22) = 'SM P L1 ' |
1078 |
gdiag(23) = 'UM 024L ' |
gdiag(23) = 'UM 024LL ' |
1079 |
gdiag(24) = 'VM 023L ' |
gdiag(24) = 'VM 023LL ' |
1080 |
gdiag(25) = 'SM L ' |
gdiag(25) = 'SM LL ' |
1081 |
gdiag(26) = 'SM L ' |
gdiag(26) = 'SM LL ' |
1082 |
gdiag(27) = 'SM L ' |
gdiag(27) = 'SM L1 ' |
1083 |
gdiag(28) = 'SM L ' |
gdiag(28) = 'SM L1 ' |
1084 |
gdiag(29) = 'SM L ' |
gdiag(29) = 'SM L1 ' |
1085 |
gdiag(30) = 'SM L ' |
gdiag(30) = 'SM L1 ' |
1086 |
gdiag(31) = 'SM L ' |
gdiag(31) = 'SM L1 ' |
1087 |
gdiag(32) = 'SM L ' |
gdiag(32) = 'SM L1 ' |
1088 |
gdiag(33) = 'SM P L ' |
gdiag(33) = 'SM P L1 ' |
1089 |
gdiag(34) = 'SM P L ' |
gdiag(34) = 'SM P L1 ' |
1090 |
gdiag(35) = 'SM L ' |
gdiag(35) = 'SM L1 ' |
1091 |
gdiag(36) = 'SM L ' |
gdiag(36) = 'SM L1 ' |
1092 |
gdiag(37) = 'SM U ' |
gdiag(37) = 'SM U1 ' |
1093 |
gdiag(38) = 'SM U ' |
gdiag(38) = 'SM U1 ' |
1094 |
gdiag(39) = 'SM L ' |
gdiag(39) = 'SM L1 ' |
1095 |
gdiag(40) = 'SM M ' |
gdiag(40) = 'SM ML ' |
1096 |
gdiag(41) = 'SM M ' |
gdiag(41) = 'SM ML ' |
1097 |
gdiag(42) = 'SM P M ' |
gdiag(42) = 'SM P ML ' |
1098 |
gdiag(43) = 'SM M ' |
gdiag(43) = 'SM ML ' |
1099 |
gdiag(44) = 'SM P M ' |
gdiag(44) = 'SM P ML ' |
1100 |
gdiag(45) = 'SM P M ' |
gdiag(45) = 'SM P ML ' |
1101 |
gdiag(46) = 'SM L ' |
gdiag(46) = 'SM L1 ' |
1102 |
gdiag(47) = 'SM M ' |
gdiag(47) = 'SM ML ' |
1103 |
gdiag(48) = 'SM P U ' |
gdiag(48) = 'SM P U1 ' |
1104 |
gdiag(49) = 'SM C050M ' |
gdiag(49) = 'SM C050ML ' |
1105 |
gdiag(50) = 'SM D M ' |
gdiag(50) = 'SM D ML ' |
1106 |
gdiag(51) = 'SM P M ' |
gdiag(51) = 'SM P M1 ' |
1107 |
gdiag(52) = 'SM L ' |
gdiag(52) = 'SM L1 ' |
1108 |
gdiag(53) = 'SM L ' |
gdiag(53) = 'SM L1 ' |
1109 |
gdiag(54) = 'UU 055M ' |
gdiag(54) = 'UU 055MR ' |
1110 |
gdiag(55) = 'VV 054M ' |
gdiag(55) = 'VV 054MR ' |
1111 |
gdiag(56) = 'SM M ' |
gdiag(56) = 'SM MR ' |
1112 |
gdiag(57) = 'SM M ' |
gdiag(57) = 'SM MR ' |
1113 |
gdiag(58) = 'SM M ' |
gdiag(58) = 'SM ML ' |
1114 |
gdiag(59) = 'UU 060M ' |
gdiag(59) = 'UU 060ML ' |
1115 |
gdiag(60) = 'VV 059M ' |
gdiag(60) = 'VV 059ML ' |
1116 |
gdiag(61) = 'SM M ' |
gdiag(61) = 'SM ML ' |
1117 |
gdiag(62) = 'SM M ' |
gdiag(62) = 'SM ML ' |
1118 |
gdiag(63) = 'SM M ' |
gdiag(63) = 'SM ML ' |
1119 |
gdiag(64) = 'SM L ' |
gdiag(64) = 'SM LL ' |
1120 |
gdiag(65) = 'SM M ' |
gdiag(65) = 'SM ML ' |
1121 |
gdiag(66) = 'SM M ' |
gdiag(66) = 'SM ML ' |
1122 |
gdiag(67) = 'SM L ' |
gdiag(67) = 'SM L1 ' |
1123 |
gdiag(68) = 'SM L ' |
gdiag(68) = 'SM L1 ' |
1124 |
gdiag(69) = 'SM P M ' |
gdiag(69) = 'SM P MM ' |
1125 |
gdiag(70) = 'SM P M ' |
gdiag(70) = 'SM P M1 ' |
1126 |
gdiag(71) = 'SM P M ' |
gdiag(71) = 'SM P ML ' |
1127 |
gdiag(72) = 'SM P U ' |
gdiag(72) = 'SM P U1 ' |
1128 |
gdiag(73) = 'SM P U ' |
gdiag(73) = 'SM P U1 ' |
1129 |
gdiag(74) = 'SM P L ' |
gdiag(74) = 'SM P LL ' |
1130 |
gdiag(75) = 'UU 076M ' |
gdiag(75) = 'UU 076ML ' |
1131 |
gdiag(76) = 'VV 075M ' |
gdiag(76) = 'VV 075ML ' |
1132 |
gdiag(77) = 'SM M ' |
gdiag(77) = 'SM ML ' |
1133 |
gdiag(78) = 'SM P M ' |
gdiag(78) = 'SM P ML ' |
1134 |
gdiag(79) = 'SM M ' |
gdiag(79) = 'SM ML ' |
1135 |
gdiag(80) = 'SM L ' |
gdiag(80) = 'SM L1 ' |
1136 |
gdiag(81) = 'SM P L ' |
gdiag(81) = 'SM P LL ' |
1137 |
gdiag(82) = 'SM P L ' |
gdiag(82) = 'SM P L1 ' |
1138 |
gdiag(83) = 'SM L ' |
gdiag(83) = 'SM L1 ' |
1139 |
gdiag(84) = 'SM L ' |
gdiag(84) = 'SM L1 ' |
1140 |
gdiag(85) = 'SM L ' |
gdiag(85) = 'SM L1 ' |
1141 |
gdiag(86) = 'SM M ' |
gdiag(86) = 'SM ML ' |
1142 |
gdiag(87) = 'SM M ' |
gdiag(87) = 'SM ML ' |
1143 |
gdiag(88) = 'UU 089M ' |
gdiag(88) = 'UU 089ML ' |
1144 |
gdiag(89) = 'VV 088M ' |
gdiag(89) = 'VV 088ML ' |
1145 |
gdiag(90) = 'SM M ' |
gdiag(90) = 'SM ML ' |
1146 |
gdiag(91) = 'SM M ' |
gdiag(91) = 'SM ML ' |
1147 |
gdiag(92) = 'UU 093M ' |
gdiag(92) = 'UU 093ML ' |
1148 |
gdiag(93) = 'VV 092M ' |
gdiag(93) = 'VV 092ML ' |
1149 |
gdiag(94) = 'UM 095M ' |
gdiag(94) = 'UM 095ML ' |
1150 |
gdiag(95) = 'VM 094M ' |
gdiag(95) = 'VM 094ML ' |
1151 |
gdiag(96) = 'UM 097L ' |
gdiag(96) = 'UM 097L1 ' |
1152 |
gdiag(97) = 'VM 096L ' |
gdiag(97) = 'VM 096L1 ' |
1153 |
gdiag(98) = 'UM 099U ' |
gdiag(98) = 'UM 099U1 ' |
1154 |
gdiag(99) = 'VM 098U ' |
gdiag(99) = 'VM 098U1 ' |
1155 |
gdiag(100) = 'SM M ' |
gdiag(100) = 'SM ML ' |
1156 |
gdiag(101) = 'SM L ' |
gdiag(101) = 'SM L1 ' |
1157 |
gdiag(102) = 'UM 103M ' |
gdiag(102) = 'UM 103M1 ' |
1158 |
gdiag(103) = 'VM 102M ' |
gdiag(103) = 'VM 102M1 ' |
1159 |
gdiag(104) = 'UM 105M ' |
gdiag(104) = 'UM 105M1 ' |
1160 |
gdiag(105) = 'VM 104M ' |
gdiag(105) = 'VM 104M1 ' |
1161 |
gdiag(106) = 'SM P M ' |
gdiag(106) = 'SM P M1 ' |
1162 |
gdiag(107) = 'SM P M ' |
gdiag(107) = 'SM P M1 ' |
1163 |
gdiag(108) = 'UM 109M ' |
gdiag(108) = 'UM 109M1 ' |
1164 |
gdiag(109) = 'VM 108M ' |
gdiag(109) = 'VM 108M1 ' |
1165 |
gdiag(110) = 'SM M ' |
gdiag(110) = 'SM M1 ' |
1166 |
gdiag(111) = 'SM P M ' |
gdiag(111) = 'SM P M1 ' |
1167 |
gdiag(112) = 'UM 113M ' |
gdiag(112) = 'UM 113M1 ' |
1168 |
gdiag(113) = 'VM 112M ' |
gdiag(113) = 'VM 112M1 ' |
1169 |
gdiag(114) = 'SM M ' |
gdiag(114) = 'SM M1 ' |
1170 |
gdiag(115) = 'SM P M ' |
gdiag(115) = 'SM P M1 ' |
1171 |
gdiag(116) = 'SM P L ' |
gdiag(116) = 'SM P LL ' |
1172 |
gdiag(117) = 'SM M ' |
gdiag(117) = 'SM ML ' |
1173 |
gdiag(118) = 'SM M ' |
gdiag(118) = 'SM M1 ' |
1174 |
gdiag(119) = 'SM P L ' |
gdiag(119) = 'SM P L1 ' |
1175 |
gdiag(120) = 'UU 121M ' |
gdiag(120) = 'UU 121ML ' |
1176 |
gdiag(121) = 'VV 120M ' |
gdiag(121) = 'VV 120ML ' |
1177 |
gdiag(122) = 'SM M ' |
gdiag(122) = 'SM ML ' |
1178 |
gdiag(123) = 'SM M ' |
gdiag(123) = 'SM ML ' |
1179 |
gdiag(124) = 'SM L ' |
gdiag(124) = 'SM L1 ' |
1180 |
gdiag(125) = 'SM L ' |
gdiag(125) = 'SM L1 ' |
1181 |
gdiag(126) = 'SM L ' |
gdiag(126) = 'SM L1 ' |
1182 |
gdiag(127) = 'SM L ' |
gdiag(127) = 'SM L1 ' |
1183 |
gdiag(128) = 'SM L ' |
gdiag(128) = 'SM L1 ' |
1184 |
gdiag(129) = 'SM L ' |
gdiag(129) = 'SM L1 ' |
1185 |
gdiag(130) = 'SM L ' |
gdiag(130) = 'SM L1 ' |
1186 |
gdiag(131) = 'SM L ' |
gdiag(131) = 'SM L1 ' |
1187 |
gdiag(132) = 'SM M ' |
gdiag(132) = 'SM ML ' |
1188 |
gdiag(133) = 'SM M ' |
gdiag(133) = 'SM ML ' |
1189 |
gdiag(134) = 'SM M ' |
gdiag(134) = 'SM ML ' |
1190 |
gdiag(135) = 'SM M ' |
gdiag(135) = 'SM ML ' |
1191 |
gdiag(136) = 'SM M ' |
gdiag(136) = 'SM ML ' |
1192 |
gdiag(137) = 'SM M ' |
gdiag(137) = 'SM ML ' |
1193 |
gdiag(138) = 'SM M ' |
gdiag(138) = 'SM ML ' |
1194 |
gdiag(139) = 'SM M ' |
gdiag(139) = 'SM ML ' |
1195 |
gdiag(140) = 'SM P M ' |
gdiag(140) = 'SM P M1 ' |
1196 |
gdiag(141) = 'SM P M ' |
gdiag(141) = 'SM P M1 ' |
1197 |
gdiag(142) = 'SM PC150M ' |
gdiag(142) = 'SM PC150M1 ' |
1198 |
gdiag(143) = 'SM PC151M ' |
gdiag(143) = 'SM PC151M1 ' |
1199 |
gdiag(144) = 'SM PC152M ' |
gdiag(144) = 'SM PC152M1 ' |
1200 |
gdiag(145) = 'SM M ' |
gdiag(145) = 'SM ML ' |
1201 |
gdiag(146) = 'SM PC147M ' |
gdiag(146) = 'SM PC147M1 ' |
1202 |
gdiag(147) = 'SM D M ' |
gdiag(147) = 'SM D M1 ' |
1203 |
gdiag(148) = 'SM PC149M ' |
gdiag(148) = 'SM PC149M1 ' |
1204 |
gdiag(149) = 'SM D M ' |
gdiag(149) = 'SM D M1 ' |
1205 |
gdiag(150) = 'SM D M ' |
gdiag(150) = 'SM D M1 ' |
1206 |
gdiag(151) = 'SM D M ' |
gdiag(151) = 'SM D M1 ' |
1207 |
gdiag(152) = 'SM D M ' |
gdiag(152) = 'SM D M1 ' |
1208 |
gdiag(153) = 'SM L ' |
gdiag(153) = 'SM L1 ' |
1209 |
gdiag(154) = 'SM L ' |
gdiag(154) = 'SM L1 ' |
1210 |
gdiag(155) = 'SM L ' |
gdiag(155) = 'SM L1 ' |
1211 |
gdiag(156) = 'SM L ' |
gdiag(156) = 'SM L1 ' |
1212 |
gdiag(157) = 'SM L ' |
gdiag(157) = 'SM L1 ' |
1213 |
gdiag(158) = 'SM L ' |
gdiag(158) = 'SM L1 ' |
1214 |
gdiag(159) = 'SM L ' |
gdiag(159) = 'SM L1 ' |
1215 |
gdiag(160) = 'SM L ' |
gdiag(160) = 'SM L1 ' |
1216 |
gdiag(161) = 'SM L ' |
gdiag(161) = 'SM L1 ' |
1217 |
gdiag(162) = 'SM L ' |
gdiag(162) = 'SM L1 ' |
1218 |
gdiag(163) = 'SM P L ' |
gdiag(163) = 'SM P L1 ' |
1219 |
gdiag(164) = 'SM L ' |
gdiag(164) = 'SM L1 ' |
1220 |
gdiag(165) = 'SM L ' |
gdiag(165) = 'SM L1 ' |
1221 |
gdiag(166) = 'SM L ' |
gdiag(166) = 'SM L1 ' |
1222 |
gdiag(167) = 'SM L ' |
gdiag(167) = 'SM L1 ' |
1223 |
gdiag(168) = 'SM L ' |
gdiag(168) = 'SM L1 ' |
1224 |
gdiag(169) = 'SM L ' |
gdiag(169) = 'SM L1 ' |
1225 |
gdiag(170) = 'SM L ' |
gdiag(170) = 'SM L1 ' |
1226 |
gdiag(171) = 'SM L ' |
gdiag(171) = 'SM L1 ' |
1227 |
gdiag(172) = 'SM P L ' |
gdiag(172) = 'SM P L1 ' |
1228 |
gdiag(173) = 'SM P L ' |
gdiag(173) = 'SM P L1 ' |
1229 |
gdiag(174) = 'SM L ' |
gdiag(174) = 'SM L1 ' |
1230 |
gdiag(175) = 'SM L ' |
gdiag(175) = 'SM L1 ' |
1231 |
gdiag(176) = 'SM L ' |
gdiag(176) = 'SM L1 ' |
1232 |
gdiag(177) = 'SM L ' |
gdiag(177) = 'SM L1 ' |
1233 |
gdiag(178) = 'SM L ' |
gdiag(178) = 'SM L1 ' |
1234 |
gdiag(179) = 'SM M ' |
gdiag(179) = 'SM M1 ' |
1235 |
gdiag(180) = 'SM M ' |
gdiag(180) = 'SM M1 ' |
1236 |
gdiag(181) = 'SM M ' |
gdiag(181) = 'SM M1 ' |
1237 |
gdiag(182) = 'SM P M ' |
gdiag(182) = 'SM P M1 ' |
1238 |
gdiag(183) = 'SM M ' |
gdiag(183) = 'SM M1 ' |
1239 |
gdiag(184) = 'SM PC185M ' |
gdiag(184) = 'SM PC185M1 ' |
1240 |
gdiag(185) = 'SM PD M ' |
gdiag(185) = 'SM PD M1 ' |
1241 |
gdiag(186) = 'SM L ' |
gdiag(186) = 'SM L1 ' |
1242 |
gdiag(187) = 'SM C188L ' |
gdiag(187) = 'SM C188L1 ' |
1243 |
gdiag(188) = 'SM D L ' |
gdiag(188) = 'SM D L1 ' |
1244 |
gdiag(189) = ' ' |
gdiag(189) = ' 0 ' |
1245 |
gdiag(190) = 'SM L ' |
gdiag(190) = 'SM L1 ' |
1246 |
gdiag(191) = 'SM P L ' |
gdiag(191) = 'SM P L1 ' |
1247 |
gdiag(192) = 'SM L ' |
gdiag(192) = 'SM L1 ' |
1248 |
gdiag(193) = 'SM L ' |
gdiag(193) = 'SM L1 ' |
1249 |
gdiag(194) = 'SM L ' |
gdiag(194) = 'SM L1 ' |
1250 |
gdiag(195) = 'SM L ' |
gdiag(195) = 'SM L1 ' |
1251 |
gdiag(196) = 'SM L ' |
gdiag(196) = 'SM L1 ' |
1252 |
gdiag(197) = 'SM L ' |
gdiag(197) = 'SM L1 ' |
1253 |
gdiag(198) = 'SM L ' |
gdiag(198) = 'SM L1 ' |
1254 |
gdiag(199) = 'SM L ' |
gdiag(199) = 'SM L1 ' |
1255 |
gdiag(200) = 'SM L ' |
gdiag(200) = 'SM L1 ' |
1256 |
gdiag(201) = 'SM L ' |
gdiag(201) = 'SM L1 ' |
1257 |
gdiag(202) = 'SM M ' |
gdiag(202) = 'SM ML ' |
1258 |
gdiag(203) = 'SM M ' |
gdiag(203) = 'SM ML ' |
1259 |
gdiag(204) = 'SM C210L ' |
gdiag(204) = 'SM C210L1 ' |
1260 |
gdiag(205) = 'SM L ' |
gdiag(205) = 'SM L1 ' |
1261 |
gdiag(206) = 'SM L ' |
gdiag(206) = 'SM L1 ' |
1262 |
gdiag(207) = 'SM L ' |
gdiag(207) = 'SM L1 ' |
1263 |
gdiag(208) = 'SM L ' |
gdiag(208) = 'SM L1 ' |
1264 |
gdiag(209) = 'SM M ' |
gdiag(209) = 'SM ML ' |
1265 |
gdiag(210) = 'SM D L ' |
gdiag(210) = 'SM D L1 ' |
1266 |
gdiag(211) = 'SM M ' |
gdiag(211) = 'SM ML ' |
1267 |
gdiag(212) = 'SM M ' |
gdiag(212) = 'SM ML ' |
1268 |
gdiag(213) = 'SM P M ' |
gdiag(213) = 'SM P ML ' |
1269 |
gdiag(214) = 'SM P L ' |
gdiag(214) = 'SM P L1 ' |
1270 |
gdiag(215) = 'SM P L ' |
gdiag(215) = 'SM P L1 ' |
1271 |
gdiag(216) = 'SM P L ' |
gdiag(216) = 'SM P L1 ' |
1272 |
gdiag(217) = 'SM P L ' |
gdiag(217) = 'SM P L1 ' |
1273 |
gdiag(218) = 'SM P L ' |
gdiag(218) = 'SM P L1 ' |
1274 |
gdiag(219) = 'SM P L ' |
gdiag(219) = 'SM P L1 ' |
1275 |
gdiag(220) = 'SM P M ' |
gdiag(220) = 'SM P ML ' |
1276 |
gdiag(221) = 'SM M ' |
gdiag(221) = 'SM ML ' |
1277 |
gdiag(222) = 'SM M ' |
gdiag(222) = 'SM ML ' |
1278 |
gdiag(223) = 'SM M ' |
gdiag(223) = 'SM ML ' |
1279 |
gdiag(224) = 'SM M ' |
gdiag(224) = 'SM ML ' |
1280 |
gdiag(225) = 'SM M ' |
gdiag(225) = 'SM ML ' |
1281 |
gdiag(226) = 'SM M ' |
gdiag(226) = 'SM ML ' |
1282 |
gdiag(227) = 'SM M ' |
gdiag(227) = 'SM ML ' |
1283 |
gdiag(228) = 'SM M ' |
gdiag(228) = 'SM ML ' |
1284 |
gdiag(229) = 'UM 230L ' |
gdiag(229) = 'UM 230L1 ' |
1285 |
gdiag(230) = 'VM 229L ' |
gdiag(230) = 'VM 229L1 ' |
1286 |
gdiag(231) = 'SM L ' |
gdiag(231) = 'SM L1 ' |
1287 |
gdiag(232) = 'UM 233L ' |
gdiag(232) = 'UM 233L1 ' |
1288 |
gdiag(233) = 'VM 232L ' |
gdiag(233) = 'VM 232L1 ' |
1289 |
gdiag(234) = 'UM 235L ' |
gdiag(234) = 'UM 235L1 ' |
1290 |
gdiag(235) = 'VM 234L ' |
gdiag(235) = 'VM 234L1 ' |
1291 |
gdiag(236) = 'UM 237M ' |
gdiag(236) = 'UM 237ML ' |
1292 |
gdiag(237) = 'VM 236M ' |
gdiag(237) = 'VM 236ML ' |
1293 |
gdiag(238) = 'SM M ' |
gdiag(238) = 'SM M1 ' |
1294 |
gdiag(239) = 'SM M ' |
gdiag(239) = 'SM M1 ' |
1295 |
gdiag(240) = 'SM M ' |
gdiag(240) = 'SM MR ' |
1296 |
gdiag(241) = 'SM M ' |
gdiag(241) = 'SM MR ' |
1297 |
gdiag(242) = 'UU 243M ' |
gdiag(242) = 'UU 243MR ' |
1298 |
gdiag(243) = 'VV 242M ' |
gdiag(243) = 'VV 242MR ' |
1299 |
gdiag(244) = 'WM L ' |
gdiag(244) = 'WM LR ' |
1300 |
gdiag(245) = 'SM M ' |
gdiag(245) = 'SM MR ' |
1301 |
gdiag(246) = 'SM M ' |
gdiag(246) = 'SM MR ' |
1302 |
gdiag(247) = 'UU 248M ' |
gdiag(247) = 'UU 248MR ' |
1303 |
gdiag(248) = 'VV 247M ' |
gdiag(248) = 'VV 247MR ' |
1304 |
gdiag(249) = 'WM L ' |
gdiag(249) = 'WM LR ' |
1305 |
gdiag(250) = 'UZ M ' |
gdiag(250) = 'UZ 250MR ' |
1306 |
gdiag(251) = 'UU 252M ' |
gdiag(251) = 'UU 252MR ' |
1307 |
gdiag(252) = 'VV 251M ' |
gdiag(252) = 'VV 251MR ' |
1308 |
gdiag(253) = 'WM M ' |
gdiag(253) = 'WM MR ' |
1309 |
gdiag(254) = 'UU 255M ' |
gdiag(254) = 'UU 255MR ' |
1310 |
gdiag(255) = 'VV 254M ' |
gdiag(255) = 'VV 254MR ' |
1311 |
gdiag(256) = 'WM M ' |
gdiag(256) = 'WM MR ' |
1312 |
gdiag(257) = 'UU 258M ' |
gdiag(257) = 'UU 258MR ' |
1313 |
gdiag(258) = 'VV 257M ' |
gdiag(258) = 'VV 257MR ' |
1314 |
gdiag(259) = 'WM M ' |
gdiag(259) = 'WM MR ' |
1315 |
gdiag(260) = 'UU 261M ' |
gdiag(260) = 'UU 261MR ' |
1316 |
gdiag(261) = 'VV 260M ' |
gdiag(261) = 'VV 260MR ' |
1317 |
gdiag(262) = 'WM M ' |
gdiag(262) = 'WM MR ' |
1318 |
gdiag(263) = 'UU 264M ' |
gdiag(263) = 'UU 264MR ' |
1319 |
gdiag(264) = 'VV 263M ' |
gdiag(264) = 'VV 263MR ' |
1320 |
gdiag(265) = 'WM M ' |
gdiag(265) = 'WM MR ' |
1321 |
gdiag(266) = 'UU 267M ' |
gdiag(266) = 'UU 267MR ' |
1322 |
gdiag(267) = 'VV 266M ' |
gdiag(267) = 'VV 266MR ' |
1323 |
gdiag(268) = 'WM M ' |
gdiag(268) = 'WM MR ' |
1324 |
gdiag(269) = 'UU 270M ' |
gdiag(269) = 'UU 270MR ' |
1325 |
gdiag(270) = 'VV 269M ' |
gdiag(270) = 'VV 269MR ' |
1326 |
gdiag(271) = 'WM M ' |
gdiag(271) = 'WM MR ' |
1327 |
gdiag(272) = 'UU 273M ' |
gdiag(272) = 'UU 273MR ' |
1328 |
gdiag(273) = 'VV 272M ' |
gdiag(273) = 'VV 272MR ' |
1329 |
gdiag(274) = 'WM M ' |
gdiag(274) = 'WM MR ' |
1330 |
gdiag(275) = 'UU 276M ' |
gdiag(275) = 'UU 276MR ' |
1331 |
gdiag(276) = 'VV 275M ' |
gdiag(276) = 'VV 275MR ' |
1332 |
gdiag(277) = 'WM M ' |
gdiag(277) = 'WM MR ' |
1333 |
gdiag(278) = 'UU 279M ' |
gdiag(278) = 'UU 279MR ' |
1334 |
gdiag(279) = 'VV 278M ' |
gdiag(279) = 'VV 278MR ' |
1335 |
gdiag(280) = 'WM M ' |
gdiag(280) = 'WM MR ' |
1336 |
|
gdiag(281) = 'SM P MR ' |
1337 |
c Diagnostic Levels |
gdiag(282) = 'SM P MR ' |
1338 |
c ----------------- |
gdiag(283) = 'SM P MR ' |
1339 |
KDIAG( 1) = 1 |
gdiag(284) = 'SM P MR ' |
1340 |
KDIAG( 2) = 1 |
gdiag(285) = 'SM P MR ' |
1341 |
KDIAG( 3) = 1 |
gdiag(286) = 'SM P MR ' |
1342 |
KDIAG( 4) = 1 |
gdiag(287) = 'SM P MR ' |
1343 |
KDIAG( 5) = 1 |
gdiag(288) = 'SM P L1 ' |
1344 |
KDIAG( 6) = 1 |
gdiag(289) = 'SM P L1 ' |
1345 |
KDIAG( 7) = 1 |
gdiag(290) = 'SM P L1 ' |
1346 |
KDIAG( 8) = nrphys |
gdiag(291) = 'SM P M1 ' |
1347 |
KDIAG( 9) = 1 |
gdiag(292) = 'SM P M1 ' |
1348 |
KDIAG(10) = 1 |
gdiag(293) = 'SM MR ' |
1349 |
KDIAG(11) = nrphys |
gdiag(294) = 'SM MR ' |
1350 |
KDIAG(12) = nrphys |
gdiag(295) = 'SM MR ' |
1351 |
KDIAG(13) = nrphys |
gdiag(296) = 'SM M1 ' |
1352 |
KDIAG(14) = nrphys |
gdiag(297) = 'SM M1 ' |
1353 |
KDIAG(15) = nrphys |
gdiag(298) = 'SM MR ' |
1354 |
KDIAG(16) = nrphys |
|
1355 |
KDIAG(17) = nrphys |
c Diagnostic Levels |
1356 |
KDIAG(18) = nrphys |
c ----------------- |
1357 |
KDIAG(19) = nrphys |
do n = 1,ndiagt |
1358 |
KDIAG(20) = nrphys |
IF (gdiag(n)(10:10) .eq. '0') KDIAG(n) = 0 |
1359 |
KDIAG(21) = 1 |
IF (gdiag(n)(10:10) .eq. '1') KDIAG(n) = 1 |
1360 |
KDIAG(22) = 1 |
IF (gdiag(n)(10:10) .eq. 'R') KDIAG(n) = Nr |
1361 |
KDIAG(23) = nrphys |
IF (gdiag(n)(10:10) .eq. 'L') KDIAG(n) = nlevs |
1362 |
KDIAG(24) = nrphys |
IF (gdiag(n)(10:10) .eq. 'M') KDIAG(n) = nlevs - 1 |
1363 |
KDIAG(25) = nrphys |
enddo |
|
KDIAG(26) = nrphys |
|
|
KDIAG(27) = 1 |
|
|
KDIAG(28) = 1 |
|
|
KDIAG(29) = 1 |
|
|
KDIAG(30) = 1 |
|
|
KDIAG(31) = 1 |
|
|
KDIAG(32) = 1 |
|
|
KDIAG(33) = 1 |
|
|
KDIAG(34) = 1 |
|
|
KDIAG(35) = 1 |
|
|
KDIAG(36) = 1 |
|
|
KDIAG(37) = 1 |
|
|
KDIAG(38) = 1 |
|
|
KDIAG(39) = 1 |
|
|
KDIAG(40) = nrphys |
|
|
KDIAG(41) = nrphys |
|
|
KDIAG(42) = nrphys |
|
|
KDIAG(43) = nrphys |
|
|
KDIAG(44) = nrphys |
|
|
KDIAG(45) = nrphys |
|
|
KDIAG(46) = 1 |
|
|
KDIAG(47) = nrphys |
|
|
KDIAG(48) = 1 |
|
|
KDIAG(49) = nrphys |
|
|
KDIAG(50) = nrphys |
|
|
KDIAG(51) = 1 |
|
|
KDIAG(52) = 1 |
|
|
KDIAG(53) = 1 |
|
|
KDIAG(54) = Nr |
|
|
KDIAG(55) = Nr |
|
|
KDIAG(56) = Nr |
|
|
KDIAG(57) = Nr |
|
|
KDIAG(58) = nrphys |
|
|
KDIAG(59) = nrphys |
|
|
KDIAG(60) = nrphys |
|
|
KDIAG(61) = nrphys |
|
|
KDIAG(62) = nrphys |
|
|
KDIAG(63) = nrphys |
|
|
KDIAG(64) = nrphys |
|
|
KDIAG(65) = nrphys |
|
|
KDIAG(66) = nrphys |
|
|
KDIAG(67) = 1 |
|
|
KDIAG(68) = 1 |
|
|
KDIAG(69) = nrphys-1 |
|
|
KDIAG(70) = 1 |
|
|
KDIAG(71) = nrphys |
|
|
KDIAG(72) = 1 |
|
|
KDIAG(73) = 1 |
|
|
KDIAG(74) = nrphys |
|
|
KDIAG(75) = nrphys |
|
|
KDIAG(76) = nrphys |
|
|
KDIAG(77) = nrphys |
|
|
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 |
|
|
KDIAG(238)= 1 |
|
|
KDIAG(239)= 1 |
|
|
KDIAG(240)= Nr |
|
|
KDIAG(241)= Nr |
|
|
KDIAG(242)= Nr |
|
|
KDIAG(243)= Nr |
|
|
KDIAG(244)= Nr |
|
|
KDIAG(245)= Nr |
|
|
KDIAG(246)= Nr |
|
|
KDIAG(247)= Nr |
|
|
KDIAG(248)= Nr |
|
|
KDIAG(249)= Nr |
|
|
KDIAG(250)= Nr |
|
|
KDIAG(251)= Nr |
|
|
KDIAG(252)= Nr |
|
|
KDIAG(253)= Nr |
|
|
KDIAG(254)= Nr |
|
|
KDIAG(255)= Nr |
|
|
KDIAG(256)= Nr |
|
|
KDIAG(257)= Nr |
|
|
KDIAG(258)= Nr |
|
|
KDIAG(259)= Nr |
|
|
KDIAG(260)= Nr |
|
|
KDIAG(261)= Nr |
|
|
KDIAG(262)= Nr |
|
|
KDIAG(263)= Nr |
|
|
KDIAG(264)= Nr |
|
|
KDIAG(265)= Nr |
|
|
KDIAG(266)= Nr |
|
|
KDIAG(267)= Nr |
|
|
KDIAG(268)= Nr |
|
|
KDIAG(269)= Nr |
|
|
KDIAG(270)= Nr |
|
|
KDIAG(271)= Nr |
|
|
KDIAG(272)= Nr |
|
|
KDIAG(273)= Nr |
|
|
KDIAG(274)= Nr |
|
|
KDIAG(275)= Nr |
|
|
KDIAG(276)= Nr |
|
|
KDIAG(277)= Nr |
|
|
KDIAG(278)= Nr |
|
|
KDIAG(279)= Nr |
|
|
KDIAG(280)= Nr |
|
1364 |
|
|
1365 |
do n=1,ndiagt |
do n=1,ndiagt |
1366 |
idiag (n) = 0 |
idiag(n) = 0 |
1367 |
ndiag(n) = 0 |
ndiag(n) = 0 |
1368 |
enddo |
enddo |
1369 |
|
|
1370 |
RETURN |
RETURN |
1371 |
END |
END |
1372 |
|
|
1373 |
|
|
1374 |
|
C---+----1----+----2----+----3----+----4----+----5----+----6----+----7-|--+----| |