atmos_init Subroutine

public subroutine atmos_init(comp, importState, exportState, clock, rc)

Arguments

Type IntentOptional Attributes Name
type(ESMF_GridComp) :: comp
type(ESMF_State) :: importState
type(ESMF_State) :: exportState
type(ESMF_Clock) :: clock
integer, intent(out) :: rc

Calls

proc~~atmos_init~2~~CallsGraph proc~atmos_init~2 atmos_init esmf_fieldcreate esmf_fieldcreate proc~atmos_init~2->esmf_fieldcreate esmf_fieldget esmf_fieldget proc~atmos_init~2->esmf_fieldget esmf_stateadd esmf_stateadd proc~atmos_init~2->esmf_stateadd interface~esmf_vmget ESMF_VMGet proc~atmos_init~2->interface~esmf_vmget proc~esmf_arrayspecset ESMF_ArraySpecSet proc~atmos_init~2->proc~esmf_arrayspecset proc~esmf_gridcompget ESMF_GridCompGet proc~atmos_init~2->proc~esmf_gridcompget proc~make_grid_sph make_grid_sph proc~atmos_init~2->proc~make_grid_sph proc~esmf_vmgetdefault ESMF_VMGetDefault interface~esmf_vmget->proc~esmf_vmgetdefault proc~esmf_vmgetpetspecific ESMF_VMGetPetSpecific interface~esmf_vmget->proc~esmf_vmgetpetspecific proc~esmf_logseterror ESMF_LogSetError proc~esmf_arrayspecset->proc~esmf_logseterror proc~esmf_compget ESMF_CompGet proc~esmf_gridcompget->proc~esmf_compget proc~esmf_compstatusget ESMF_CompStatusGet proc~esmf_gridcompget->proc~esmf_compstatusget proc~esmf_gridcompgetinit ESMF_GridCompGetInit proc~esmf_gridcompget->proc~esmf_gridcompgetinit proc~esmf_imerr ESMF_IMErr proc~esmf_gridcompget->proc~esmf_imerr proc~esmf_logfounderror ESMF_LogFoundError proc~esmf_gridcompget->proc~esmf_logfounderror interface~esmf_gridaddcoord ESMF_GridAddCoord proc~make_grid_sph->interface~esmf_gridaddcoord interface~esmf_gridadditem ESMF_GridAddItem proc~make_grid_sph->interface~esmf_gridadditem interface~esmf_gridcreate1peridim ESMF_GridCreate1PeriDim proc~make_grid_sph->interface~esmf_gridcreate1peridim interface~esmf_gridcreatenoperidim ESMF_GridCreateNoPeriDim proc~make_grid_sph->interface~esmf_gridcreatenoperidim interface~esmf_gridgetcoord ESMF_GridGetCoord proc~make_grid_sph->interface~esmf_gridgetcoord interface~esmf_gridgetitem ESMF_GridGetItem proc~make_grid_sph->interface~esmf_gridgetitem

Source Code

  subroutine atmos_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_ArraySpec)  :: arrayspec
    type(ESMF_DistGrid)   :: distgrid
    type(ESMF_Grid)       :: grid
    type(ESMF_Field)      :: field
    type(ESMF_VM)         :: vm
    integer               :: petCount, localPet
    real(ESMF_KIND_R8), pointer :: farrayPtr(:,:)   ! matching F90 array pointer
    
    ! Initialize return code
    rc = ESMF_SUCCESS

    ! Determine petCount
    call ESMF_GridCompGet(comp, vm=vm, rc=rc)
    if (rc/=ESMF_SUCCESS) return ! bail out
    call ESMF_VMGet(vm, petCount=petCount, localPet=localPet, rc=rc)
    if (rc/=ESMF_SUCCESS) return ! bail out
    print *, "Atmosphere Init starting, localPet =", localPet
    
    ! Create the source Field and add it to the export State
    call ESMF_ArraySpecSet(arrayspec, typekind=ESMF_TYPEKIND_R8, rank=2, rc=rc)
    if (rc/=ESMF_SUCCESS) return ! bail out
    ! 
    grid = make_grid_sph(120,120,3.,1.5,0.,-90.,rc=rc)
    if (rc/=ESMF_SUCCESS) return ! bail out
    !
    field = ESMF_FieldCreate(arrayspec=arrayspec, grid=grid, &
      indexflag=ESMF_INDEX_GLOBAL, name="F_atm", rc=rc)
    if (rc/=ESMF_SUCCESS) return ! bail out
    call ESMF_StateAdd(importState, (/field/), rc=rc)
    if (rc/=ESMF_SUCCESS) return ! bail out
   
    ! Gain access to actual data via F90 array pointer
    call ESMF_FieldGet(field, localDe=0, farrayPtr=farrayPtr, rc=rc)
    if (rc/=ESMF_SUCCESS) return ! bail out

    ! Fill destination Field with invalid data
    farrayPtr = -999.0
    print *, "Atmos Init returning"

  end subroutine atmos_init