subroutine user_init(comp, importState, exportState, clock, rc)
type(ESMF_GridComp) :: comp
type(ESMF_State) :: importState, exportState
type(ESMF_Clock) :: clock
integer, intent(out) :: rc
! Local variables
type(ESMF_VM) :: vm
integer :: de_id, npets
integer, allocatable :: indexList(:)
type(ESMF_DistGrid) :: distgrid
type(ESMF_LocStream) :: locs
type(ESMF_Field) :: field(3)
type(ESMF_FieldBundle) :: fieldbundle
! Initialize return code
rc = ESMF_SUCCESS
! Query component for VM and create a layout with the right breakdown
call ESMF_GridCompGet(comp, vm=vm, rc=rc)
if(rc/=ESMF_SUCCESS) return
call ESMF_VMGet(vm, localPet=de_id, petCount=npets, rc=rc)
if(rc/=ESMF_SUCCESS) return
print *, de_id, "User Comp 2 Init starting"
! Check for correct number of PETs
if ( npets .ne. 2 ) then
call ESMF_LogSetError(ESMF_RC_ARG_BAD,&
msg="This component must run on 2 PETs.",&
line=__LINE__, file=__FILE__, rcToReturn=rc)
return
endif
if ( de_id .eq. 0 ) then
allocate(indexList(6))
allocate(tmp_data(6))
allocate(hum_data(6))
allocate(prs_data(6))
indexList = (/1,2,3,4,5,6/)
tmp_data(1:6) = 0
hum_data(1:6) = 0
prs_data(1:6) = 0
elseif ( de_id .eq. 1 ) then
allocate(indexList(10))
allocate(tmp_data(10))
allocate(hum_data(10))
allocate(prs_data(10))
indexList = (/7,8,9,10,11,12,13,14,15,16/)
tmp_data(1:10) = 0
hum_data(1:10) = 0
prs_data(1:10) = 0
endif
! Add "temperature" "humidity" "pressure" fields to the export state.
distgrid = ESMF_DistGridCreate(arbSeqIndexList=indexList, rc=rc)
if (rc .ne. ESMF_SUCCESS) return
deallocate(indexList)
locs = ESMF_LocStreamCreate(distgrid=distgrid, &
indexflag=ESMF_INDEX_DELOCAL, coordSys=ESMF_COORDSYS_CART, &
name="Test LocStream", rc=rc)
if (rc .ne. ESMF_SUCCESS) return
field(1) = ESMF_FieldCreate(locs, tmp_data, &
indexflag=ESMF_INDEX_DELOCAL, name="temperature", rc=rc)
if (rc .ne. ESMF_SUCCESS) return
field(2) = ESMF_FieldCreate(locs, hum_data, &
indexflag=ESMF_INDEX_DELOCAL, name="humidity", rc=rc)
if (rc .ne. ESMF_SUCCESS) return
field(3) = ESMF_FieldCreate(locs, prs_data, &
indexflag=ESMF_INDEX_DELOCAL, name="pressure", rc=rc)
if (rc .ne. ESMF_SUCCESS) return
fieldbundle = ESMF_FieldBundleCreate(fieldList=field, &
name="fieldbundle data", rc=rc)
if (rc/=ESMF_SUCCESS) return
call ESMF_StateAdd(importState, fieldbundleList=(/fieldbundle/), rc=rc)
if (rc/=ESMF_SUCCESS) return
print *, de_id, "User Comp2 Init returning"
end subroutine user_init