Changeset 5099 for LMDZ6/branches/Amaury_dev/libf/dynphy_lonlat/phylmd

Timestamp:

Jul 22, 2024, 9:29:09 PM (4 months ago)

Author:

abarral

Message:

Replace most uses of CPP_DUST by the corresponding logical defined in lmdz_cppkeys_wrapper.F90
Convert several files from .F to .f90 to allow Dust to compile w/o rrtm/ecrad
Create lmdz_yoerad.f90
(lint) Remove "!" on otherwise empty line

Location:

LMDZ6/branches/Amaury_dev/libf/dynphy_lonlat/phylmd

Files:

• callphysiq_mod.F90 (modified) (1 diff)
• ce0l.F90 (modified) (6 diffs)
• etat0dyn_netcdf.F90 (modified) (11 diffs)
• etat0phys_netcdf.F90 (modified) (10 diffs)
• iniphysiq_mod.F90 (modified) (1 diff)
• init_ssrf_m.F90 (modified) (3 diffs)
• limit_netcdf.F90 (modified) (22 diffs)

LMDZ6/branches/Amaury_dev/libf/dynphy_lonlat/phylmd/callphysiq_mod.F90

@@ -1 @@
-!
+
 ! $Id: $
-!
+
 MODULE callphysiq_mod
 

LMDZ6/branches/Amaury_dev/libf/dynphy_lonlat/phylmd/ce0l.F90

@@ -1 +1 @@
 PROGRAM ce0l
-!
+
 !-------------------------------------------------------------------------------
 ! Purpose: Initial states and boundary conditions files creation:
@@ -21 +21 @@
   USE etat0phys,      ONLY: etat0phys_netcdf
   USE limit,          ONLY: limit_netcdf
-  USE netcdf,         ONLY: NF90_OPEN, NF90_NOWRITE, NF90_CLOSE, NF90_NOERR,    &
-         NF90_INQUIRE_DIMENSION, nf90_inq_dimid, NF90_INQ_VARID, NF90_GET_VAR
+  USE netcdf,         ONLY: NF90_OPEN, NF90_NOWRITE, NF90_CLOSE, nf90_noerr,    &
+         nf90_inquire_dimension, nf90_inq_dimid, NF90_INQ_VARID, nf90_get_var
   USE infotrac,       ONLY: init_infotrac
   USE dimphy,         ONLY: klon
@@ -154 +154 @@
 ! weights to ensure ocean fractions are the same for atmosphere and ocean.
 !*******************************************************************************
-  IF(NF90_OPEN("o2a.nc", NF90_NOWRITE, nid_o2a)==NF90_NOERR) THEN
+  IF(NF90_OPEN("o2a.nc", NF90_NOWRITE, nid_o2a)==nf90_noerr) THEN
     iret=NF90_CLOSE(nid_o2a)
     WRITE(lunout,*)'BEWARE !! Ocean mask "o2a.nc" file found'
@@ -188 +188 @@
     masque(iip1 ,:)=masque(1,:)
     DEALLOCATE(ocemask)
-  ELSE IF(NF90_OPEN("startphy0.nc", NF90_NOWRITE, nid_sta)==NF90_NOERR) THEN
+  ELSE IF(NF90_OPEN("startphy0.nc", NF90_NOWRITE, nid_sta)==nf90_noerr) THEN
     WRITE(lunout,*)'BEWARE !! File "startphy0.nc" found.'
     WRITE(lunout,*)'Getting the land mask from a previous run.'
     iret=nf90_inq_dimid(nid_sta,'points_physiques',nid_nph)
-    iret=NF90_INQUIRE_DIMENSION(nid_sta,nid_nph,len=nphys)
+    iret=nf90_inquire_dimension(nid_sta,nid_nph,len=nphys)
     IF(nphys/=klon) THEN
       WRITE(lunout,*)'Mismatching dimensions for land mask'
@@ -201 +201 @@
     ALLOCATE(masktmp(klon))
     iret=NF90_INQ_VARID(nid_sta,'masque',nid_msk)
-    iret=NF90_GET_VAR(nid_sta,nid_msk,masktmp)
+    iret=nf90_get_var(nid_sta,nid_msk,masktmp)
     iret=NF90_CLOSE(nid_sta)
     CALL gr_fi_dyn(1,klon,iip1,jjp1,masktmp,masque)
@@ -253 +253 @@
 
 END PROGRAM ce0l
-!
-!-------------------------------------------------------------------------------
+
+!-------------------------------------------------------------------------------

LMDZ6/branches/Amaury_dev/libf/dynphy_lonlat/phylmd/etat0dyn_netcdf.F90

@@ -1 +1 @@
 MODULE etat0dyn
-!
+
 !*******************************************************************************
 ! Purpose: Create dynamical initial state using atmospheric fields from a
@@ -6 +6 @@
 !-------------------------------------------------------------------------------
 ! Comments:
-!
+
 !    *  This module is designed to work for Earth (and with ioipsl)
-!
+
 !    *  etat0dyn_netcdf routine can access to NetCDF data through the following
 !  routine (to be called after restget):
 !    CALL startget_dyn3d(varname, lon_in,  lat_in, pls, workvar,&
 !                          champ, lon_in2, lat_in2)
-!
+
 !    *  Variables should have the following names in the NetCDF files:
 !            'U'      : East ward wind              (in "ECDYN.nc")
@@ -20 +20 @@
 !            'R'      : Relative humidity           (in "ECDYN.nc")
 !            'RELIEF' : High resolution orography   (in "Relief.nc") 
-!
+
 !    * The land mask and corresponding weights can be:
 !      1) already known (in particular if etat0dyn has been called before) ;
@@ -28 +28 @@
 !         File name: "o2a.nc"  ;  variable name: "OceMask"
 !      3) computed from topography file "Relief.nc" for forced runs.
-!
+
 !   *   There is a big mess with the longitude size. Should it be iml or iml+1 ?
 !  I have chosen to use the iml+1 as an argument to this routine and we declare
@@ -59 +59 @@
 
 !-------------------------------------------------------------------------------
-!
+
 SUBROUTINE etat0dyn_netcdf(masque, phis)
-!
+
 !-------------------------------------------------------------------------------
 ! Purpose: Create dynamical initial states.
@@ -194 +194 @@
 
 END SUBROUTINE etat0dyn_netcdf
-!
-!-------------------------------------------------------------------------------
-
-
-!-------------------------------------------------------------------------------
-!
+
+!-------------------------------------------------------------------------------
+
+
+!-------------------------------------------------------------------------------
+
 SUBROUTINE startget_dyn3d(var, lon_in,  lat_in,  pls,  workvar,&
                         champ, lon_in2, lat_in2)
@@ -279 +279 @@
 
 END SUBROUTINE startget_dyn3d
-!
-!-------------------------------------------------------------------------------
-
-
-!-------------------------------------------------------------------------------
-!
+
+!-------------------------------------------------------------------------------
+
+
+!-------------------------------------------------------------------------------
+
 SUBROUTINE start_init_dyn(lon_in,lat_in,lon_in2,lat_in2,zs,psol)
-!
+
 !-------------------------------------------------------------------------------
   IMPLICIT NONE
@@ -343 +343 @@
 
 !-------------------------------------------------------------------------------
-!
+
 SUBROUTINE get_var_dyn(title,field)
-!
+
 !-------------------------------------------------------------------------------
   USE conf_dat_m, ONLY: conf_dat2d
@@ -375 +375 @@
 
 END SUBROUTINE get_var_dyn
-!
+
 !-------------------------------------------------------------------------------
 
 END SUBROUTINE start_init_dyn
-!
-!-------------------------------------------------------------------------------
-
-
-!-------------------------------------------------------------------------------
-!
+
+!-------------------------------------------------------------------------------
+
+
+!-------------------------------------------------------------------------------
+
 SUBROUTINE start_inter_3d(var,lon_in,lat_in,lon_in2,lat_in2,pls_in,var3d)
-!
+
 !-------------------------------------------------------------------------------
   USE conf_dat_m, ONLY: conf_dat3d
@@ -461 +461 @@
 
 END SUBROUTINE start_inter_3d
-!
-!-------------------------------------------------------------------------------
-
-
-!-------------------------------------------------------------------------------
-!
+
+!-------------------------------------------------------------------------------
+
+
+!-------------------------------------------------------------------------------
+
 SUBROUTINE interp_startvar(nam,ibeg,lon,lat,vari,lon1,lat1,lon2,lat2,varo)
-!
+
 !-------------------------------------------------------------------------------
   USE inter_barxy_m, ONLY: inter_barxy
@@ -502 +502 @@
 
 END SUBROUTINE interp_startvar
-!
+
 !-------------------------------------------------------------------------------
 
 END MODULE etat0dyn
-!
-!*******************************************************************************
+
+!*******************************************************************************

LMDZ6/branches/Amaury_dev/libf/dynphy_lonlat/phylmd/etat0phys_netcdf.F90

@@ -1 @@
-!
+
 ! $Id$
-!
+
 MODULE etat0phys
-!
+
 !*******************************************************************************
 ! Purpose: Create physical initial state using atmospheric fields from a
@@ -9 +9 @@
 !-------------------------------------------------------------------------------
 ! Comments:
-!
+
 !    *  This module is designed to work for Earth (and with ioipsl)
-!
+
 !    *  etat0phys_netcdf routine can access to NetCDF data through subroutines:
 !         "start_init_phys" for variables contained in file "ECPHY.nc":
@@ -18 +18 @@
 !         "start_init_orog" for variables contained in file "Relief.nc":
 !            'RELIEF' : High resolution orography 
-!
+
 !    * The land mask and corresponding weights can be:
 !      1) computed using the ocean mask from the ocean model (to ensure ocean
@@ -24 +24 @@
 !         File name: "o2a.nc"  ;  variable name: "OceMask"
 !      2) computed from topography file "Relief.nc" for forced runs.
-!
+
 !    * Allowed values for read_climoz flag are 0, 1 and 2:
 !      0: do not read an ozone climatology
@@ -76 +76 @@
 
 !-------------------------------------------------------------------------------
-!
+
 SUBROUTINE etat0phys_netcdf(masque, phis)
-!
+
 !-------------------------------------------------------------------------------
 ! Purpose: Creates initial states
@@ -295 +295 @@
 
 END SUBROUTINE etat0phys_netcdf
-!
-!-------------------------------------------------------------------------------
-
-
-!-------------------------------------------------------------------------------
-!
+
+!-------------------------------------------------------------------------------
+
+
+!-------------------------------------------------------------------------------
+
 SUBROUTINE start_init_orog(lon_in,lat_in,phis,masque)
-!
+
 !===============================================================================
 ! Comment:
@@ -391 +391 @@
 
 END SUBROUTINE start_init_orog
-!
-!-------------------------------------------------------------------------------
-
-
-!-------------------------------------------------------------------------------
-!
+
+!-------------------------------------------------------------------------------
+
+
+!-------------------------------------------------------------------------------
+
 SUBROUTINE start_init_phys(lon_in,lat_in,phis)
-!
+
 !===============================================================================
 ! Purpose:   Compute tsol and qsol, knowing phis.
@@ -449 +449 @@
 
 !-------------------------------------------------------------------------------
-!
+
 SUBROUTINE get_var_phys(title,field)
-!
+
 !-------------------------------------------------------------------------------
   IMPLICIT NONE
@@ -474 +474 @@
 
 END SUBROUTINE get_var_phys
-!
-!-------------------------------------------------------------------------------
-!
+
+!-------------------------------------------------------------------------------
+
 END SUBROUTINE start_init_phys
-!
-!-------------------------------------------------------------------------------
-
-
-!-------------------------------------------------------------------------------
-!
+
+!-------------------------------------------------------------------------------
+
+
+!-------------------------------------------------------------------------------
+
 SUBROUTINE interp_startvar(nam,ibeg,lon,lat,vari,lon2,lat2,varo)
-!
+
 !-------------------------------------------------------------------------------
   USE inter_barxy_m, ONLY: inter_barxy
@@ -518 +518 @@
 
 END SUBROUTINE interp_startvar
-!
-!-------------------------------------------------------------------------------
-!
+
+!-------------------------------------------------------------------------------
+
 !*******************************************************************************
 

LMDZ6/branches/Amaury_dev/libf/dynphy_lonlat/phylmd/iniphysiq_mod.F90

@@ -1 @@
-!
+
 ! $Id$
-!
+
 MODULE iniphysiq_mod
 

LMDZ6/branches/Amaury_dev/libf/dynphy_lonlat/phylmd/init_ssrf_m.F90

@@ -1 +1 @@
 MODULE init_ssrf_m
-!
+
 !*******************************************************************************
 
@@ -24 +24 @@
 
 !-------------------------------------------------------------------------------
-!
+
 SUBROUTINE start_init_subsurf(known_mask)
-!
+
 !-------------------------------------------------------------------------------
 ! Purpose: Subsurfaces initialization.
@@ -141 +141 @@
 
 END SUBROUTINE start_init_subsurf
-!
+
 !-------------------------------------------------------------------------------
 
 END MODULE init_ssrf_m
-!
+
 !*******************************************************************************

LMDZ6/branches/Amaury_dev/libf/dynphy_lonlat/phylmd/limit_netcdf.F90

@@ -1 +1 @@
 MODULE limit
-!
+
 !*******************************************************************************
 ! Author : L. Fairhead, 27/01/94
@@ -49 +49 @@
 
 !-------------------------------------------------------------------------------
-!
+
 SUBROUTINE limit_netcdf(masque, phis, extrap)
-!
+
 !-------------------------------------------------------------------------------
 ! Author : L. Fairhead, 27/01/94
@@ -73 +73 @@
   USE netcdf,             ONLY: NF90_OPEN,    NF90_CREATE,  NF90_CLOSE,        &
                   NF90_DEF_DIM, NF90_DEF_VAR, NF90_PUT_VAR, NF90_PUT_ATT,      &
-                  NF90_NOERR,   NF90_NOWRITE,  NF90_GLOBAL,       &
+                  nf90_noerr,   NF90_NOWRITE,  NF90_GLOBAL,       &
                   NF90_CLOBBER, NF90_ENDDEF,  NF90_UNLIMITED, NF90_FLOAT,      &
                   NF90_64BIT_OFFSET
@@ -138 +138 @@
 ! Open file only to test if available
   DO ix_sic=1,SIZE(fsic)
-     IF ( NF90_OPEN(TRIM(fsic(ix_sic)),NF90_NOWRITE,nid)==NF90_NOERR ) THEN
+     IF ( NF90_OPEN(TRIM(fsic(ix_sic)),NF90_NOWRITE,nid)==nf90_noerr ) THEN
         icefile=fsic(ix_sic); varname=vsic(ix_sic); EXIT
      END IF
@@ -197 +197 @@
 ! Open file only to test if available
   DO ix_sst=1,SIZE(fsst)
-     IF ( NF90_OPEN(TRIM(fsst(ix_sst)),NF90_NOWRITE,nid)==NF90_NOERR ) THEN
+     IF ( NF90_OPEN(TRIM(fsst(ix_sst)),NF90_NOWRITE,nid)==nf90_noerr ) THEN
        sstfile=fsst(ix_sst); varname=vsst(ix_sst); EXIT
      END IF
@@ -305 +305 @@
 
 !===============================================================================
-!
+
   CONTAINS
-!
+
 !===============================================================================
 
 
 !-------------------------------------------------------------------------------
-!
+
 SUBROUTINE get_2Dfield(fnam, varname, mode, ndays, champo, flag, mask)
-!
+
 !-----------------------------------------------------------------------------
 ! Comments:
@@ -323 +323 @@
 !-----------------------------------------------------------------------------
   USE netcdf, ONLY: NF90_OPEN, NF90_INQ_VARID, NF90_INQUIRE_VARIABLE, &
-       NF90_CLOSE, nf90_inq_dimid, NF90_INQUIRE_DIMENSION, NF90_GET_VAR, &
+       NF90_CLOSE, nf90_inq_dimid, nf90_inquire_dimension, nf90_get_var, &
        NF90_GET_ATT
   USE pchsp_95_m, only: pchsp_95
@@ -399 +399 @@
 
 !--- Longitude
-  CALL ncerr(NF90_INQUIRE_DIMENSION(ncid, dids(1), name=dnam, len=imdep),fnam)
+  CALL ncerr(nf90_inquire_dimension(ncid, dids(1), name=dnam, len=imdep),fnam)
   ALLOCATE(dlon_ini(imdep), dlon(imdep))
   CALL ncerr(NF90_INQ_VARID(ncid, dnam, varid), fnam)
-  CALL ncerr(NF90_GET_VAR(ncid, varid, dlon_ini), fnam)
+  CALL ncerr(nf90_get_var(ncid, varid, dlon_ini), fnam)
   CALL msg(5,'variable '//TRIM(dnam)//' dimension ', imdep)
 
 !--- Latitude
-  CALL ncerr(NF90_INQUIRE_DIMENSION(ncid, dids(2), name=dnam, len=jmdep),fnam)
+  CALL ncerr(nf90_inquire_dimension(ncid, dids(2), name=dnam, len=jmdep),fnam)
   ALLOCATE(dlat_ini(jmdep), dlat(jmdep))
   CALL ncerr(NF90_INQ_VARID(ncid, dnam, varid), fnam)
-  CALL ncerr(NF90_GET_VAR(ncid, varid, dlat_ini), fnam)
+  CALL ncerr(nf90_get_var(ncid, varid, dlat_ini), fnam)
   CALL msg(5,'variable '//TRIM(dnam)//' dimension ', jmdep)
 
 !--- Time (variable is not needed - it is rebuilt - but calendar is)
-  CALL ncerr(NF90_INQUIRE_DIMENSION(ncid, dids(3), name=dnam, len=lmdep), fnam)
+  CALL ncerr(nf90_inquire_dimension(ncid, dids(3), name=dnam, len=lmdep), fnam)
   ALLOCATE(timeyear(lmdep+2))
   CALL ncerr(NF90_INQ_VARID(ncid, dnam, varid), fnam)
   cal_in=' '
-  IF(NF90_GET_ATT(ncid, varid, 'calendar', cal_in)/=NF90_NOERR) THEN
+  IF(NF90_GET_ATT(ncid, varid, 'calendar', cal_in)/=nf90_noerr) THEN
     SELECT CASE(mode)
       CASE('RUG', 'ALB'); cal_in='360_day'
@@ -451 +451 @@
   CALL ncerr(NF90_INQ_VARID(ncid, varname, varid), fnam)
   DO l=1, lmdep
-    CALL ncerr(NF90_GET_VAR(ncid,varid,champ,[1,1,l],[imdep,jmdep,1]),fnam)
+    CALL ncerr(nf90_get_var(ncid,varid,champ,[1,1,l],[imdep,jmdep,1]),fnam)
     CALL conf_dat2d(title, dlon_ini, dlat_ini, dlon, dlat, champ, .TRUE.)
 
@@ -459 +459 @@
       !--- DETERMINE THE UNIT: READ FROM FILE OR ASSUMED USING FIELD VALUES
       ierr=NF90_GET_ATT(ncid, varid, 'units', units)
-      IF(ierr==NF90_NOERR) THEN !--- ATTRIBUTE "units" FOUND IN THE FILE
+      IF(ierr==nf90_noerr) THEN !--- ATTRIBUTE "units" FOUND IN THE FILE
         CALL strclean(units)
         IF(mode=='SIC'.AND.is_in(units,Perc)) units="%"
@@ -472 +472 @@
       END IF
       CALL msg(0,'INPUT FILE '//TRIM(title)//' UNIT IS: "'//TRIM(units)//'".')
-      IF(ierr/=NF90_NOERR) CALL msg(0,'WARNING ! UNIT TO BE CHECKED ! '      &
+      IF(ierr/=nf90_noerr) CALL msg(0,'WARNING ! UNIT TO BE CHECKED ! '      &
         //'No "units" attribute, so only based on the fields values.')
 
@@ -509 +509 @@
 !--- FIRST RECORD: LAST ONE OF PREVIOUS YEAR (CURRENT YEAR IF UNAVAILABLE)
   fnam_m=fnam(1:idx)//'_m.nc'
-  IF(NF90_OPEN(fnam_m,NF90_NOWRITE,ncid)==NF90_NOERR) THEN
+  IF(NF90_OPEN(fnam_m,NF90_NOWRITE,ncid)==nf90_noerr) THEN
     CALL msg(0,'Reading previous year file ("'//TRIM(fnam_m)//'") last record for '//TRIM(title))
     CALL ncerr(NF90_INQ_VARID(ncid, varname, varid),fnam_m)
     CALL ncerr(NF90_INQUIRE_VARIABLE(ncid, varid, dimids=dids),fnam_m)
-    CALL ncerr(NF90_INQUIRE_DIMENSION(ncid, dids(3), len=l), fnam_m)
-    CALL ncerr(NF90_GET_VAR(ncid,varid,champ,[1,1,l],[imdep,jmdep,1]),fnam_m)
+    CALL ncerr(nf90_inquire_dimension(ncid, dids(3), len=l), fnam_m)
+    CALL ncerr(nf90_get_var(ncid,varid,champ,[1,1,l],[imdep,jmdep,1]),fnam_m)
     CALL ncerr(NF90_CLOSE(ncid), fnam_m)
     CALL conf_dat2d(title, dlon_ini, dlat_ini, dlon, dlat, champ, .TRUE.)
@@ -532 +532 @@
 !--- LAST RECORD: FIRST ONE OF NEXT YEAR (CURRENT YEAR IF UNAVAILABLE)
   fnam_p=fnam(1:idx)//'_p.nc'
-  IF(NF90_OPEN(fnam_p,NF90_NOWRITE,ncid)==NF90_NOERR) THEN
+  IF(NF90_OPEN(fnam_p,NF90_NOWRITE,ncid)==nf90_noerr) THEN
     CALL msg(0,'Reading next year file ("'//TRIM(fnam_p)//'") first record for '//TRIM(title))
     CALL ncerr(NF90_INQ_VARID(ncid, varname, varid),fnam_p)
-    CALL ncerr(NF90_GET_VAR(ncid,varid,champ,[1,1,1],[imdep,jmdep,1]),fnam_p)
+    CALL ncerr(nf90_get_var(ncid,varid,champ,[1,1,1],[imdep,jmdep,1]),fnam_p)
     CALL ncerr(NF90_CLOSE(ncid), fnam_p)
     CALL conf_dat2d(title, dlon_ini, dlat_ini, dlon, dlat, champ, .TRUE.)
@@ -641 +641 @@
 
 END SUBROUTINE get_2Dfield
-!
-!-------------------------------------------------------------------------------
-
-
-!-------------------------------------------------------------------------------
-!
+
+!-------------------------------------------------------------------------------
+
+
+!-------------------------------------------------------------------------------
+
 SUBROUTINE start_init_orog0(lon_in,lat_in,phis,masque)
-!
+
 !-------------------------------------------------------------------------------
   USE grid_noro_m, ONLY: grid_noro0
@@ -704 +704 @@
 
 END SUBROUTINE start_init_orog0
-!
-!-------------------------------------------------------------------------------
-
-
-!-------------------------------------------------------------------------------
-!
+
+!-------------------------------------------------------------------------------
+
+
+!-------------------------------------------------------------------------------
+
 SUBROUTINE msg(lev,str1,i,str2)
-!
+
 !-------------------------------------------------------------------------------
 ! Arguments:
@@ -730 +730 @@
 
 END SUBROUTINE msg
-!
-!-------------------------------------------------------------------------------
-
-
-!-------------------------------------------------------------------------------
-!
+
+!-------------------------------------------------------------------------------
+
+
+!-------------------------------------------------------------------------------
+
 SUBROUTINE ncerr(ncres,fnam)
-!
+
 !-------------------------------------------------------------------------------
 ! Purpose: NetCDF errors handling.
 !-------------------------------------------------------------------------------
-  USE netcdf, ONLY : NF90_NOERR, NF90_STRERROR
+  USE netcdf, ONLY : nf90_noerr, NF90_STRERROR
   IMPLICIT NONE
 !-------------------------------------------------------------------------------
@@ -748 +748 @@
   CHARACTER(LEN=*), INTENT(IN) :: fnam
 !-------------------------------------------------------------------------------
-  IF(ncres/=NF90_NOERR) THEN
+  IF(ncres/=nf90_noerr) THEN
     WRITE(lunout,*)'Problem with file '//TRIM(fnam)//' in routine limit_netcdf.'
     CALL abort_physic('limit_netcdf',NF90_STRERROR(ncres),1)
@@ -754 +754 @@
 
 END SUBROUTINE ncerr
-!
-!-------------------------------------------------------------------------------
-
-
-!-------------------------------------------------------------------------------
-!
+
+!-------------------------------------------------------------------------------
+
+
+!-------------------------------------------------------------------------------
+
 SUBROUTINE strclean(s)
-!
+
 !-------------------------------------------------------------------------------
   IMPLICIT NONE
@@ -776 +776 @@
 
 END SUBROUTINE strclean
-!
-!-------------------------------------------------------------------------------
-
-
-!-------------------------------------------------------------------------------
-!
+
+!-------------------------------------------------------------------------------
+
+
+!-------------------------------------------------------------------------------
+
 FUNCTION is_in(s1,s2) RESULT(res)
-!
+
 !-------------------------------------------------------------------------------
   IMPLICIT NONE
@@ -796 +796 @@
 
 END FUNCTION is_in
-!
-!-------------------------------------------------------------------------------
-
-
-!-------------------------------------------------------------------------------
-!
+
+!-------------------------------------------------------------------------------
+
+
+!-------------------------------------------------------------------------------
+
 ELEMENTAL FUNCTION strLow(s) RESULT(res)
-!
+
 !-------------------------------------------------------------------------------
   IMPLICIT NONE
@@ -822 +822 @@
 
 END FUNCTION strLow
-!
+
 !-------------------------------------------------------------------------------
 
@@ -830 +830 @@
 
 END MODULE limit
-!
+
 !*******************************************************************************