12 |
,tra & |
,tra & |
13 |
,h,lv,cpn,p,ph,tv,tp,tvp,clw & |
,h,lv,cpn,p,ph,tv,tp,tvp,clw & |
14 |
,sig,w0 ) |
,sig,w0 ) |
15 |
use cvparam3 |
use cv3_param_m |
16 |
implicit none |
implicit none |
17 |
|
|
18 |
|
|
19 |
! inputs: |
! inputs: |
20 |
integer len,ncum,nd,ntra,nloc |
integer, intent(in):: len,ncum,nd,ntra,nloc |
21 |
integer iflag1(len),nk1(len),icb1(len),icbs1(len) |
integer iflag1(len),nk1(len),icb1(len),icbs1(len) |
22 |
real plcl1(len),tnk1(len),qnk1(len),gznk1(len) |
real plcl1(len),tnk1(len),qnk1(len),gznk1(len) |
23 |
real pbase1(len),buoybase1(len) |
real pbase1(len),buoybase1(len) |
74 |
100 continue |
100 continue |
75 |
110 continue |
110 continue |
76 |
|
|
|
! do 121 j=1,ntra |
|
|
! do 111 k=1,nd |
|
|
! nn=0 |
|
|
! do 101 i=1,len |
|
|
! if(iflag1(i).eq.0)then |
|
|
! nn=nn+1 |
|
|
! tra(nn,k,j)=tra1(i,k,j) |
|
|
! endif |
|
|
! 101 continue |
|
|
! 111 continue |
|
|
! 121 continue |
|
|
|
|
77 |
if (nn.ne.ncum) then |
if (nn.ne.ncum) then |
78 |
print*,'strange! nn not equal to ncum: ',nn,ncum |
print*,'strange! nn not equal to ncum: ',nn,ncum |
79 |
stop |
stop |