@@ -210,7 +210,6 @@ subroutine clm_cpl_HYDRO()
210
210
write (6 ,* ) " end of drive ndhms"
211
211
#endif
212
212
213
- return
214
213
end subroutine clm_cpl_HYDRO
215
214
216
215
subroutine Toclm3d (v1d_out ,nn ,kk1 ,z1_in ,ix ,jx ,kk2 ,z2 ,v2_in )
@@ -235,7 +234,6 @@ subroutine Toclm3d (v1d_out,nn,kk1,z1_in,ix,jx,kk2,z2,v2_in)
235
234
call Toclm2d(v1(:,:,k),ix,jx,v1d_out(:,k),nn)
236
235
end do
237
236
238
- return
239
237
end subroutine Toclm3d
240
238
241
239
subroutine Toclm2d (v1 ,ix ,jx ,v1d_out ,nn )
@@ -266,7 +264,6 @@ subroutine Toclm2d(v1,ix,jx,v1d_out,nn)
266
264
write (6 ,* ) " after scatter_1d_r"
267
265
268
266
#endif
269
- return
270
267
end subroutine Toclm2d
271
268
272
269
subroutine TO1d (vg2d ,v1d ,nx ,ny )
@@ -280,7 +277,6 @@ subroutine TO1d(vg2d,v1d,nx,ny)
280
277
v1d(n) = vg2d(i,j)
281
278
enddo
282
279
enddo
283
- return
284
280
end subroutine TO1d
285
281
286
282
subroutine clm2ND3d (v1d_in ,nn ,kk1 ,z1_in ,ix ,jx ,kk2 ,z2 ,vout )
@@ -302,8 +298,6 @@ subroutine clm2ND3d (v1d_in,nn,kk1,z1_in,ix,jx,kk2,z2,vout)
302
298
end do
303
299
call Interp3D (z1(1 :kk1),v1,kk1,z2(1 :kk2),vout,ix,jx,kk2)
304
300
305
-
306
- return
307
301
end subroutine clm2ND3d
308
302
309
303
subroutine clm2ND2d (v1d_in ,nn ,vout ,ix ,jx )
@@ -346,7 +340,6 @@ subroutine clm2ND2d (v1d_in,nn,vout,ix,jx)
346
340
! domain decomposition
347
341
call decompose_data_real(v2d_g,vout)
348
342
349
- return
350
343
end subroutine clm2ND2d
351
344
352
345
subroutine Interp3D (z1 ,v1 ,kk1 ,z ,vout ,ix ,jx ,kk )
@@ -366,8 +359,8 @@ subroutine Interp3D (z1,v1,kk1,z,vout,ix,jx,kk)
366
359
do j = 1 , jx
367
360
do i = 1 , ix
368
361
do k = 1 , kk
369
- ! call interpLayer(abs (Z1),v1(i,j,1 :kk1),kk1,abs ( Z(k) ),vout(i,j,k))
370
- call interpLayer(Z1(1 :kk1),v1(i,j,1 :kk1),kk1,Z(k),vout(i,j,k))
362
+ ! call interpLayer(abs (Z1),v1(i,j,1 :kk1),kk1,abs ( Z(k) ),vout(i,j,k))
363
+ call interpLayer(Z1(1 :kk1),v1(i,j,1 :kk1),kk1,Z(k),vout(i,j,k))
371
364
end do
372
365
end do
373
366
end do
@@ -408,7 +401,7 @@ subroutine interpLayer(inZ,inV,inK,outZ,outV)
408
401
if (inZ(k1) .eq. inZ(k2)) then
409
402
write (6 ,* ) " FATAL ERROR: inZ(k1)=inZ(k2) " , inZ(k1),inZ(k2)
410
403
!stop 99
411
- stop (" FATAL ERROR: In module_clm_HYDRO.F interpLayer()- inZ(k1)=inZ(k2)" )
404
+ stop (" FATAL ERROR: In module_clm_HYDRO.F interpLayer()- inZ(k1)=inZ(k2)" )
412
405
end if
413
406
w1 = (outZ- inZ(k1))/ (inZ(k2)- inZ(k1))
414
407
w2 = (inZ(k2)- outZ)/ (inZ(k2)- inZ(k1))
@@ -692,7 +685,6 @@ subroutine output_nc(array,idim,jdim, var_name, file_name)
692
685
iret = nf_inq_varid(ncid,var_name,varid)
693
686
iret = nf_put_var_real(ncid,varid,array)
694
687
iret= nf_close(ncid)
695
- return
696
688
end subroutine output_nc
697
689
698
690
subroutine g2c_2d (lbc , ubc , lbl , ubl , lbg , ubg , num2d , carr , garr , &
@@ -995,7 +987,7 @@ subroutine g2c_2d_tmp(lbc, ubc, lbl, ubl, lbg, ubg, num2d, carr, garr, &
995
987
!w_yw = scale_c2l(c) * scale_l2g(l) * wtgcell(c)
996
988
w_yw = scale_c2l(c) * scale_l2g(l)
997
989
if (w_yw .ne. 0 ) then
998
- if (abs (garr(g,j)) .gt. minv .and. abs (garr(g,j)) .lt. maxv) then
990
+ if (abs (garr(g,j)) .gt. minv .and. abs (garr(g,j)) .lt. maxv) then
999
991
! carr(c,j) = garr(g,j) / w_yw * sumwt(g)
1000
992
yw_r = garr(g,j) / w_yw
1001
993
if (yw_r .gt. 0 ) then
0 commit comments