source: LMDZ6/branches/contrails/libf/phylmd/grid_index.f90 @ 5449

Last change on this file since 5449 was 5390, checked in by yann meurdesoif, 3 weeks ago
  • Remove UTF8 character that inihibit fortran parsing with GPU morphosis
  • Add missing END SUBROUTINE instead of simple END, that inhibit correct parsing with regulat expression parser (quick and dirty parsing)

YM

File size: 737 bytes
Line 
1function grid_index(lon_deg,lat_deg)
2!--------------------------------------------------------------------------------
3! Get local index of grid point of longitude,latitude lon_deg,lat_deg
4! Author : FH 2024/07/18
5! Please do not put this function in a m*odule not to complexify the replay script
6!--------------------------------------------------------------------------------
7USE dimphy, only : klon
8USE geometry_mod, ONLY: latitude_deg, longitude_deg
9implicit none
10real, intent(in) :: lon_deg,lat_deg
11integer :: grid_index
12integer i
13grid_index=0
14do i=1,klon
15   if ( abs(lon_deg-longitude_deg(i)) < 0.02 .and.  abs(lat_deg-latitude_deg(i)) < 0.02 ) then
16        grid_index=i
17        exit
18   endif
19enddo
20return
21end function grid_index
Note: See TracBrowser for help on using the repository browser.