subroutine user_run(comp, importState, exportState, clock, rc)
type(ESMF_GridComp) :: comp
type(ESMF_State) :: importState, exportState
type(ESMF_Clock) :: clock
integer, intent(out) :: rc
type(ESMF_Field) :: humidity
integer, save :: onetime=1 ! static variable
type(mydata), pointer :: mydatablock
type(wrapper) :: mywrapper
rc = ESMF_SUCCESS
print *, "User Comp Run starting"
! In a real application, the coupler would move the export from
! one component to the import before this call. For now, copy the
! field from the export state to import state by hand.
if (onetime .gt. 0) then
call ESMF_StateGet(exportState, "humidity", humidity, rc=rc)
if (rc/=ESMF_SUCCESS) return ! bail on error
call ESMF_StateAdd(importState, (/humidity/), rc=rc)
if (rc/=ESMF_SUCCESS) return ! bail on error
onetime = 0
endif
! Get private data block
print *, "in run, ready to get data block"
nullify(mydatablock)
mywrapper%wrap => mydatablock
call ESMF_GridCompGetInternalState(comp, mywrapper, rc)
if (rc/=ESMF_SUCCESS) return ! bail on error
mydatablock => mywrapper%wrap
print *, "run, local data =", mydatablock%index, &
mydatablock%scale_factor, mydatablock%flag
call ESMF_StatePrint(importState, rc=rc)
if (rc/=ESMF_SUCCESS) return ! bail on error
call ESMF_StateGet(importState, "humidity", humidity, rc=rc)
if (rc/=ESMF_SUCCESS) return ! bail on error
call ESMF_FieldPrint(humidity, rc=rc)
if (rc/=ESMF_SUCCESS) return ! bail on error
! This is where the model specific computation goes.
! Here is where the output state is updated.
!call ESMF_StateAdd(exportState, humidity, rc=status)
call ESMF_StatePrint(exportState, rc=rc)
if (rc/=ESMF_SUCCESS) return ! bail on error
print *, "User Comp Run returning"
end subroutine user_run