[1765] | 1 | module find_coord_m |
---|
| 2 | |
---|
| 3 | implicit none |
---|
| 4 | |
---|
| 5 | contains |
---|
| 6 | |
---|
| 7 | subroutine find_coord(ncid, name, dimid, varid, std_name) |
---|
| 8 | |
---|
| 9 | ! This procedure returns the name, dimension id or variable id of |
---|
| 10 | ! the NetCDF coordinate with standard name "std_name", if such a |
---|
| 11 | ! coordinate exists. The standard name is only used to know what |
---|
| 12 | ! to search, it is not used for the search itself. The search |
---|
| 13 | ! itself is done via a string match on the attribute "units". |
---|
| 14 | |
---|
| 15 | use netcdf, only: NF90_MAX_NAME, NF90_NOERR |
---|
| 16 | use nf95_get_att_m, only: nf95_get_att |
---|
| 17 | use nf95_inq_varid_m, only: nf95_inq_varid |
---|
| 18 | use nf95_inquire_dimension_m, only: nf95_inquire_dimension |
---|
| 19 | use nf95_inquire_m, only: nf95_inquire |
---|
| 20 | use nf95_inquire_variable_m, only: nf95_inquire_variable |
---|
| 21 | |
---|
| 22 | integer, intent(in):: ncid |
---|
| 23 | |
---|
| 24 | character(len=*), intent(out), optional:: name ! blanks if not found |
---|
| 25 | ! The actual character argument should normally have the length |
---|
| 26 | ! "NF90_MAX_NAME". |
---|
| 27 | |
---|
| 28 | integer, intent(out), optional:: dimid ! 0 if not found |
---|
| 29 | integer, intent(out), optional:: varid ! 0 if not found |
---|
| 30 | |
---|
| 31 | character(len=*), intent(in):: std_name |
---|
| 32 | ! standard name : "latitude", "longitude" or "time" |
---|
| 33 | |
---|
| 34 | ! Variables local to the procedure: |
---|
| 35 | |
---|
| 36 | character(len=13) units |
---|
| 37 | logical exact ! "units" must be matched exactly |
---|
| 38 | |
---|
| 39 | integer ncerr, nDimensions, dimid_local, varid_local |
---|
| 40 | character(len=NF90_MAX_NAME) name_local |
---|
| 41 | integer, pointer:: dimids(:) |
---|
| 42 | character(len=80) values |
---|
| 43 | logical found |
---|
| 44 | |
---|
| 45 | !---------------------------------------------- |
---|
| 46 | |
---|
| 47 | select case (std_name) |
---|
| 48 | case("longitude") |
---|
| 49 | units="degrees_east" |
---|
| 50 | exact=.true. |
---|
| 51 | case("latitude") |
---|
| 52 | units="degrees_north" |
---|
| 53 | exact=.true. |
---|
| 54 | case("time") |
---|
| 55 | units=" since" |
---|
| 56 | exact=.false. |
---|
| 57 | case default |
---|
| 58 | print *, "find_coord: bad value of std_name" |
---|
| 59 | print *, "std_name = ", std_name |
---|
| 60 | stop 1 |
---|
| 61 | end select |
---|
| 62 | |
---|
| 63 | call nf95_inquire(ncid, nDimensions) |
---|
| 64 | dimid_local = 0 |
---|
| 65 | found = .false. |
---|
| 66 | |
---|
| 67 | ! Loop on dimensions: |
---|
| 68 | do while (.not. found .and. dimid_local < nDimensions) |
---|
| 69 | dimid_local = dimid_local + 1 |
---|
| 70 | call nf95_inquire_dimension(ncid, dimid_local, name_local) |
---|
| 71 | call nf95_inq_varid(ncid, name_local, varid_local, ncerr) |
---|
| 72 | if (ncerr == NF90_NOERR) then |
---|
| 73 | call nf95_inquire_variable(ncid, varid_local, dimids=dimids) |
---|
| 74 | if (size(dimids) == 1) then |
---|
| 75 | if (dimids(1) == dimid_local) then |
---|
| 76 | ! We have found a coordinate |
---|
| 77 | call nf95_get_att(ncid, varid_local, "units", values, ncerr) |
---|
| 78 | if (ncerr == NF90_NOERR)then |
---|
| 79 | if (exact) then |
---|
| 80 | found = values == units |
---|
| 81 | else |
---|
| 82 | found = index(values, trim(units)) /= 0 |
---|
| 83 | end if |
---|
| 84 | end if |
---|
| 85 | end if |
---|
| 86 | end if |
---|
| 87 | deallocate(dimids) ! pointer |
---|
| 88 | end if |
---|
| 89 | end do |
---|
| 90 | |
---|
| 91 | if (found) then |
---|
| 92 | if (present(name)) name = name_local |
---|
| 93 | if (present(dimid)) dimid = dimid_local |
---|
| 94 | if (present(varid)) varid = varid_local |
---|
| 95 | else |
---|
| 96 | if (present(name)) name = "" |
---|
| 97 | if (present(dimid)) dimid = 0 |
---|
| 98 | if (present(varid)) varid = 0 |
---|
| 99 | end if |
---|
| 100 | |
---|
| 101 | end subroutine find_coord |
---|
| 102 | |
---|
| 103 | end module find_coord_m |
---|