Skip to content

Commit c07c3ae

Browse files
committed
Set ent_type correctly with land is single column
Set ent_type correctly with land is in single column mode for section that zeros river variables.
1 parent 4630bc6 commit c07c3ae

File tree

1 file changed

+10
-3
lines changed

1 file changed

+10
-3
lines changed

driver-moab/main/prep_lnd_mod.F90

Lines changed: 10 additions & 3 deletions
Original file line numberDiff line numberDiff line change
@@ -22,6 +22,7 @@ module prep_lnd_mod
2222

2323
use seq_comm_mct, only: mbaxid ! iMOAB id for atm migrated mesh to coupler pes
2424
use seq_comm_mct, only: atm_pg_active ! whether the atm uses FV mesh or not ; made true if fv_nphys > 0
25+
use seq_comm_mct, only: mb_scm_land
2526
! use dimensions_mod, only: np ! for atmosphere
2627
use seq_comm_mct, only: seq_comm_getinfo => seq_comm_setptrs
2728
use seq_map_type_mod
@@ -400,9 +401,15 @@ subroutine prep_lnd_init(infodata, atm_c2_lnd, rof_c2_lnd, glc_c2_lnd, iac_c2_ln
400401
write(logunit,*) subname,' cant get size of land mesh'
401402
call shr_sys_abort(subname//' ERROR in getting size of land mesh')
402403
endif
403-
! land is now cell mesh on coupler side
404-
mlsize = nvise(1)
405-
ent_type = 1 ! cell
404+
! land is usually cell on coupler but could be point
405+
if(mb_scm_land) then
406+
mlsize = nvert(1)
407+
ent_type = 0 ! point cloud
408+
else
409+
mlsize = nvise(1)
410+
ent_type = 1 ! cell
411+
endif
412+
406413
! set to 0 all fields that are projected from river
407414
nrflds = mct_aVect_nRattr(r2x_lx(1)) ! these are the numbers of fields in seq_flds_r2x_fields
408415
arrsize = nrflds*mlsize

0 commit comments

Comments
 (0)