IMPORTANT: Global and EAST longitudes are now allowed in regrid.F

[arango]

This one has been in my TODO lists for a long time. Finally, I have a global forcing data with longitudes ranging between 0-360 for an application.

• Updated regrid.F to allow global grids and EAST longitudes. Introduced a new work variable GRID(:)%MyLon for the interpolation of data to ROMS grid, and passed as argument:

       SUBROUTINE regrid (ng, model, ncname, ncid,                       &
      &                   ncvname, ncvarid, gtype, iflag,                &
      &                   Nx, Ny, Finp, Amin, Amax,                      &
      &                   LBi, UBi, LBj, UBj,                            &
      &                   Imin, Imax, Jmin, Jmax,                        &
      &                   MyXout, Xout, Yout, Fout)
  
  

  as MyXout dummy argument. The Following code is added for processing application grid longitude. Notice that input data longitude is not modified:

  !
  !  Copy longitude coordinate Xout to MyXout. If the longitude of the
  !  data is from a global grid [0-360] or in degrees_east, convert Xout
  !  to east longitudes (MyXout) to facilitate regridding. In such case,
  !  positive multiples of 360 map to 360 and negative multiples of 360
  !  map to zero using the MODULO intrinsic Fortran function.
  !
        IF ((Xmin.ge.0.0_r8).and.(Xmax.gt.0.0_r8).and.                    &
       &    ((Xmax-Xmin).gt.315.0_r8)) THEN
          EastLon=.TRUE.
          MyLonMin=MODULO(LonMin(ng), 360.0_r8)
          IF ((MyLonMin.eq.0.0_r8).and.                                   &
       &          (LonMin(ng).gt.0.0_r8)) MyLonMin=360.0_r8
          MyLonMax=MODULO(LonMax(ng), 360.0_r8)
          IF ((MyLonMax.eq.0.0_r8).and.                                   &
       &          (LonMax(ng).gt.0.0_r8)) MyLonMax=360.0_r8
        ELSE
          EastLon=.FALSE.
          MyLonMin=LonMin(ng)
          MyLonMax=LonMax(ng)
        END IF
        IF (EastLon) THEN
          DO j=Jmin,Jmax
            DO i=Imin,Imax
              MyXout(i,j)=MODULO(Xout(i,j), 360.0_r8)   ! range [0 360]
              IF ((MyXout(i,j).eq.0.0_r8).and.                            &
       &          (Xout(i,j).gt.0.0_r8)) MyXout(i,j)=360.0_r8
            END DO
          END DO
        ELSE
          DO j=Jmin,Jmax
            DO i=Imin,Imax
              MyXout(i,j)=Xout(i,j)                     ! range [-180 180]
            END DO
          END DO
        END IF
  

  The routine nf_fread2d.F was modified to add the extra argument to the regrid call. the mod_grid.F now includes the MyLon variable in the structure.

• The routines in check_multifile.F and inquiry.F are more robust by allowing lowercase, uppercase or a combination when searching for the time coordinste associated with a NetCDF variable

        IF ((INDEX(TRIM(lowercase(var_name(i))),'time').ne.0).and.      &
     &            (var_ndim(i).eq.1)) THEN
          TvarName=TRIM(var_name(i))
          foundit=.TRUE.
          EXIT
          ...
        END IF

Notice that we convert to lowercase to simplify the comparison in the conditional.

• Corrected argument to load_l in read_stapar.F:

             CASE ('Sout(idTsur)')
               Npts=load_l(Nval, Cval, MT, Ngrids, Ltracer)
               DO ng=1,Ngrids
                 DO itrc=1,NAT
                   Sout(idTsur(itrc),ng)=Ltracer(itrc,ng)
                 END DO
               END DO
  

  Many thanks to John Warner for reporting this issue.