65 |
INTEGER kd, kd0, ksgn, kStore |
INTEGER kd, kd0, ksgn, kStore |
66 |
CHARACTER*(MAX_LEN_MBUF) msgBuf |
CHARACTER*(MAX_LEN_MBUF) msgBuf |
67 |
integer i,offset,Lena,newindx,jindx |
integer i,offset,Lena,newindx,jindx |
68 |
_RL undef,getcon |
_RL undef |
69 |
INTEGER iSp, ndId, j,l |
INTEGER iSp, ndId, j,l |
70 |
INTEGER region2fill(0:nRegions) |
INTEGER region2fill(0:nRegions) |
71 |
_RL scaleFact |
_RL scaleFact |
72 |
_RL gridField(sNx*sNy,nlevs), gridFrac(sNx*sNy) |
_RL gridField(sNx*sNy,nlevs), gridFrac(sNx*sNy) |
73 |
|
|
74 |
|
#ifdef ALLOW_FIZHI |
75 |
|
_RL getcon |
76 |
|
EXTERNAL getcon |
77 |
|
#endif |
78 |
|
|
79 |
C Run through list of active diagnostics to make sure |
C Run through list of active diagnostics to make sure |
80 |
C we are trying to fill a valid diagnostic |
C we are trying to fill a valid diagnostic |
81 |
|
|
82 |
|
undef = UNSET_RL |
83 |
ndiagnum = 0 |
ndiagnum = 0 |
84 |
ipointer = 0 |
ipointer = 0 |
85 |
DO n=1,nlists |
DO n=1,nlists |
149 |
DO k = kFirst,kLast |
DO k = kFirst,kLast |
150 |
kd = kd0 + ksgn*k |
kd = kd0 + ksgn*k |
151 |
if( check ) then |
if( check ) then |
152 |
|
#ifdef ALLOW_FIZHI |
153 |
undef = getcon('UNDEF') |
undef = getcon('UNDEF') |
154 |
|
#endif |
155 |
do i= 1,Lena |
do i= 1,Lena |
156 |
jindx = 1 + int((indx(i+offset-1)-1)/sNx) |
jindx = 1 + int((indx(i+offset-1)-1)/sNx) |
157 |
newindx = indx(i+offset-1)+(jindx-1)*2*Olx |
newindx = indx(i+offset-1)+(jindx-1)*2*Olx |
236 |
gridField(i,k)= 0. |
gridField(i,k)= 0. |
237 |
ENDDO |
ENDDO |
238 |
if( check ) then |
if( check ) then |
239 |
|
#ifdef ALLOW_FIZHI |
240 |
undef = getcon('UNDEF') |
undef = getcon('UNDEF') |
241 |
|
#endif |
242 |
do i= 1,Lena |
do i= 1,Lena |
243 |
newindx = indx(i+offset-1) |
newindx = indx(i+offset-1) |
244 |
if(gridField(newindx,k).eq.undef |
if(gridField(newindx,k).eq.undef |