Changeset 5067 for LMDZ6/trunk/libf

Timestamp:

Jul 18, 2024, 12:11:41 PM (5 months ago)

Author:

abarral

Message:

(lint) replace obsolete operators

File:

• LMDZ6/trunk/libf/dyn3d/guide_mod.F90 (modified) (44 diffs)

LMDZ6/trunk/libf/dyn3d/guide_mod.F90

@@ -123 +123 @@
     CALL getpar('iguide_sav',4,iguide_sav,'freq. sauvegarde guidage')
     ! frequences f>0: fx/jour; f<0: tous les f jours; f=0: 1 seule fois.
-    IF (iguide_sav.GT.0) THEN
+    IF (iguide_sav>0) THEN
        iguide_sav=day_step/iguide_sav
     ELSE if (iguide_sav == 0) then
@@ -143 +143 @@
     CALL getpar('iguide_read',4,iguide_read,'freq. lecture guidage')
     CALL getpar('iguide_int',4,iguide_int,'freq. interpolation vert')
-    IF (iguide_int.EQ.0) THEN
+    IF (iguide_int==0) THEN
         iguide_int=1
-    ELSEIF (iguide_int.GT.0) THEN
+    ELSEIF (iguide_int>0) THEN
         iguide_int=day_step/iguide_int
     ELSE
@@ -171 +171 @@
 ! ---------------------------------------------
     ncidpl=-99
-    if (guide_plevs.EQ.1) then
-       if (ncidpl.eq.-99) then
+    if (guide_plevs==1) then
+       if (ncidpl==-99) then
           rcod=nf90_open('apbp.nc',Nf90_NOWRITe, ncidpl)
-          if (rcod.NE.NF90_NOERR) THEN
+          if (rcod/=NF90_NOERR) THEN
              abort_message=' Nudging error -> no file apbp.nc'
              CALL abort_gcm(modname,abort_message,1)
           endif
        endif
-    elseif (guide_plevs.EQ.2) then
-       if (ncidpl.EQ.-99) then
+    elseif (guide_plevs==2) then
+       if (ncidpl==-99) then
           rcod=nf90_open('P.nc',Nf90_NOWRITe,ncidpl)
-          if (rcod.NE.NF90_NOERR) THEN
+          if (rcod/=NF90_NOERR) THEN
              abort_message=' Nudging error -> no file P.nc'
              CALL abort_gcm(modname,abort_message,1)
@@ -189 +189 @@
 
     elseif (guide_u) then
-           if (ncidpl.eq.-99) then
+           if (ncidpl==-99) then
                rcod=nf90_open('u.nc',Nf90_NOWRITe,ncidpl)
-               if (rcod.NE.NF90_NOERR) THEN
+               if (rcod/=NF90_NOERR) THEN
                   CALL abort_gcm(modname, &
                        ' Nudging error -> no file u.nc',1)
@@ -198 +198 @@
 
     elseif (guide_v) then
-           if (ncidpl.eq.-99) then
+           if (ncidpl==-99) then
                rcod=nf90_open('v.nc',nf90_nowrite,ncidpl)
-               if (rcod.NE.NF90_NOERR) THEN
+               if (rcod/=NF90_NOERR) THEN
                   CALL abort_gcm(modname, &
                        ' Nudging error -> no file v.nc',1)
@@ -206 +206 @@
            endif
     elseif (guide_T) then
-           if (ncidpl.eq.-99) then
+           if (ncidpl==-99) then
                rcod=nf90_open('T.nc',nf90_nowrite,ncidpl)
-               if (rcod.NE.NF90_NOERR) THEN
+               if (rcod/=NF90_NOERR) THEN
                   CALL abort_gcm(modname, &
                        ' Nudging error -> no file T.nc',1)
@@ -214 +214 @@
            endif
     elseif (guide_Q) then
-           if (ncidpl.eq.-99) then
+           if (ncidpl==-99) then
                rcod=nf90_open('hur.nc',nf90_nowrite, ncidpl)
-               if (rcod.NE.NF90_NOERR) THEN
+               if (rcod/=NF90_NOERR) THEN
                   CALL abort_gcm(modname, &
                        ' Nudging error -> no file hur.nc',1)
@@ -225 +225 @@
     endif 
     error=NF_INQ_DIMID(ncidpl,'LEVEL',rid)
-    IF (error.NE.NF90_NOERR) error=NF_INQ_DIMID(ncidpl,'PRESSURE',rid)
-    IF (error.NE.NF90_NOERR) THEN
+    IF (error/=NF90_NOERR) error=NF_INQ_DIMID(ncidpl,'PRESSURE',rid)
+    IF (error/=NF90_NOERR) THEN
         CALL abort_gcm(modname,'Nudging: error reading pressure levels',1)
     ENDIF
@@ -306 +306 @@
     ENDIF
 
-    IF (guide_plevs.EQ.2) THEN
+    IF (guide_plevs==2) THEN
         ALLOCATE(pnat1(iip1,jjp1,nlevnc), stat = error)
         IF (error /= 0) CALL abort_gcm(modname,abort_message,1)
@@ -314 +314 @@
     ENDIF
 
-    IF (guide_P.OR.guide_plevs.EQ.1) THEN
+    IF (guide_P.OR.guide_plevs==1) THEN
         ALLOCATE(psnat1(iip1,jjp1), stat = error)
         IF (error /= 0) CALL abort_gcm(modname,abort_message,1)
@@ -341 +341 @@
     IF (guide_T) tnat1=tnat2
     IF (guide_Q) qnat1=qnat2
-    IF (guide_plevs.EQ.2) pnat1=pnat2
-    IF (guide_P.OR.guide_plevs.EQ.1) psnat1=psnat2
+    IF (guide_plevs==2) pnat1=pnat2
+    IF (guide_P.OR.guide_plevs==1) psnat1=psnat2
 
   END SUBROUTINE guide_init
@@ -440 +440 @@
 ! Lecture des fichiers de guidage ?
 !-----------------------------------------------------------------------
-    IF (iguide_read.NE.0) THEN
+    IF (iguide_read/=0) THEN
       ditau=real(itau)
       dday_step=real(day_step)
-      IF (iguide_read.LT.0) THEN
+      IF (iguide_read<0) THEN
           tau=ditau/dday_step/REAL(iguide_read)
       ELSE
@@ -449 +449 @@
       ENDIF
       reste=tau-AINT(tau)
-      IF (reste.EQ.0.) THEN
-          IF (itau_test.EQ.itau) THEN
+      IF (reste==0.) THEN
+          IF (itau_test==itau) THEN
             write(lunout,*)trim(modname)//' second pass in advreel at itau=',&
             itau
@@ -460 +460 @@
               IF (guide_T) tnat1=tnat2
               IF (guide_Q) qnat1=qnat2
-              IF (guide_plevs.EQ.2) pnat1=pnat2
-              IF (guide_P.OR.guide_plevs.EQ.1) psnat1=psnat2
+              IF (guide_plevs==2) pnat1=pnat2
+              IF (guide_P.OR.guide_plevs==1) psnat1=psnat2
               step_rea=step_rea+1
               itau_test=itau
@@ -482 +482 @@
 ! Interpolation et conversion des champs de guidage
 !-----------------------------------------------------------------------
-    IF (MOD(itau,iguide_int).EQ.0) THEN
+    IF (MOD(itau,iguide_int)==0) THEN
         CALL guide_interp(ps,teta)
     ENDIF
 ! Repartition entre 2 etats de guidage
-    IF (iguide_read.NE.0) THEN
+    IF (iguide_read/=0) THEN
         tau=reste
     ELSE
@@ -496 +496 @@
 !-----------------------------------------------------------------------
 ! Sauvegarde du guidage?
-    f_out=((MOD(itau,iguide_sav).EQ.0).AND.guide_sav)  
+    f_out=((MOD(itau,iguide_sav)==0).AND.guide_sav)
     IF (f_out) THEN
       ! compute pressures at layer interfaces
@@ -633 +633 @@
         IF (guide_reg) THEN
             DO i=1,iim
-                IF (lond(i).LT.lon_min_g) imin(1)=i
-                IF (lond(i).LE.lon_max_g) imax(1)=i
+                IF (lond(i)<lon_min_g) imin(1)=i
+                IF (lond(i)<=lon_max_g) imax(1)=i
             ENDDO
             lond=rlonv*180./pi
             DO i=1,iim
-                IF (lond(i).LT.lon_min_g) imin(2)=i
-                IF (lond(i).LE.lon_max_g) imax(2)=i
+                IF (lond(i)<lon_min_g) imin(2)=i
+                IF (lond(i)<=lon_max_g) imax(2)=i
             ENDDO
         ENDIF
@@ -960 +960 @@
             do j=1,pjm
                 do i=1,pim
-                    if (typ.eq.2) then
+                    if (typ==2) then
                        zlat=rlatu(j)*180./pi
                        zlon=rlonu(i)*180./pi
-                    elseif (typ.eq.1) then
+                    elseif (typ==1) then
                        zlat=rlatu(j)*180./pi
                        zlon=rlonv(i)*180./pi
-                    elseif (typ.eq.3) then
+                    elseif (typ==3) then
                        zlat=rlatv(j)*180./pi
                        zlon=rlonv(i)*180./pi
@@ -1005 +1005 @@
             enddo
         enddo
-        IF (typ.EQ.2) THEN
+        IF (typ==2) THEN
             do j=1,jjp1
                 do i=1,iim
@@ -1013 +1013 @@
             enddo
         ENDIF
-        IF (typ.EQ.3) THEN
+        IF (typ==3) THEN
             do j=1,jjm
                 do i=1,iip1
@@ -1035 +1035 @@
             enddo
             ! Calcul de gamma
-            if (abs(grossismx-1.).lt.0.1.or.abs(grossismy-1.).lt.0.1) then
+            if (abs(grossismx-1.)<0.1.or.abs(grossismy-1.)<0.1) then
               write(*,*)trim(modname)//' ATTENTION modele peu zoome'
               write(*,*)trim(modname)//' ATTENTION on prend une constante de guidage cste'
@@ -1042 +1042 @@
               gamma=(dxdy_max-2.*dxdy_min)/(dxdy_max-dxdy_min)
               write(*,*)trim(modname)//' gamma=',gamma
-              if (gamma.lt.1.e-5) then
+              if (gamma<1.e-5) then
                 write(*,*)trim(modname)//' gamma =',gamma,'<1e-5'
                 abort_message='stopped'
@@ -1057 +1057 @@
         do j=1,pjm
             do i=1,pim
-                if (typ.eq.1) then
+                if (typ==1) then
                    dxdy_=dxdys(i,j)
                    zlat=rlatu(j)*180./pi
-                elseif (typ.eq.2) then
+                elseif (typ==2) then
                    dxdy_=dxdyu(i,j)
                    zlat=rlatu(j)*180./pi
-                elseif (typ.eq.3) then
+                elseif (typ==3) then
                    dxdy_=dxdyv(i,j)
                    zlat=rlatv(j)*180./pi
                 endif
-                if (abs(grossismx-1.).lt.0.1.or.abs(grossismy-1.).lt.0.1) then
+                if (abs(grossismx-1.)<0.1.or.abs(grossismy-1.)<0.1) then
                 ! pour une grille reguliere, xi=xxx**0=1 -> alpha=alphamin
                     alpha(i,j)=alphamin
@@ -1073 +1073 @@
                     xi=((dxdy_max-dxdy_)/(dxdy_max-dxdy_min))**gamma
                     xi=min(xi,1.)
-                    if(lat_min_g.le.zlat .and. zlat.le.lat_max_g) then
+                    if(lat_min_g<=zlat .and. zlat<=lat_max_g) then
                         alpha(i,j)=xi*alphamin+(1.-xi)*alphamax
                     else
@@ -1118 +1118 @@
          write(*,*) trim(modname)//': opening nudging files '
 ! Niveaux de pression si non constants
-         if (guide_plevs.EQ.1) then
+         if (guide_plevs==1) then
              write(*,*) trim(modname)//' Reading nudging on model levels'
              rcode = nf90_open('apbp.nc', nf90_nowrite, ncidpl)
-             IF (rcode.NE.NF90_NOERR) THEN
+             IF (rcode/=NF90_NOERR) THEN
               abort_message='Nudging: error -> no file apbp.nc'
               CALL abort_gcm(modname,abort_message,1)
              ENDIF
              rcode = nf90_inq_varid(ncidpl, 'AP', varidap)
-             IF (rcode.NE.NF90_NOERR) THEN
+             IF (rcode/=NF90_NOERR) THEN
               abort_message='Nudging: error -> no AP variable in file apbp.nc'
               CALL abort_gcm(modname,abort_message,1)
              ENDIF
              rcode = nf90_inq_varid(ncidpl, 'BP', varidbp)
-             IF (rcode.NE.NF90_NOERR) THEN
+             IF (rcode/=NF90_NOERR) THEN
               abort_message='Nudging: error -> no BP variable in file apbp.nc'
               CALL abort_gcm(modname,abort_message,1)
@@ -1139 +1139 @@
 
 ! Pression si guidage sur niveaux P variables
-         if (guide_plevs.EQ.2) then
+         if (guide_plevs==2) then
              rcode = nf90_open('P.nc', nf90_nowrite, ncidp)
-             IF (rcode.NE.NF90_NOERR) THEN
+             IF (rcode/=NF90_NOERR) THEN
               abort_message='Nudging: error -> no file P.nc'
               CALL abort_gcm(modname,abort_message,1)
              ENDIF
              rcode = nf90_inq_varid(ncidp, 'PRES', varidp)
-             IF (rcode.NE.NF90_NOERR) THEN 
+             IF (rcode/=NF90_NOERR) THEN
               abort_message='Nudging: error -> no PRES variable in file P.nc'
               CALL abort_gcm(modname,abort_message,1)
              ENDIF
              write(*,*) trim(modname)//' ncidp,varidp',ncidp,varidp
-             if (ncidpl.eq.-99) ncidpl=ncidp
+             if (ncidpl==-99) ncidpl=ncidp
          endif
 
@@ -1157 +1157 @@
          if (guide_u) then
              rcode = nf90_open('u.nc', nf90_nowrite, ncidu)
-             IF (rcode.NE.NF90_NOERR) THEN
+             IF (rcode/=NF90_NOERR) THEN
               abort_message='Nudging: error -> no file u.nc'
               CALL abort_gcm(modname,abort_message,1)
              ENDIF
              rcode = nf90_inq_varid(ncidu, 'UWND', varidu)
-             IF (rcode.NE.NF90_NOERR) THEN
+             IF (rcode/=NF90_NOERR) THEN
               abort_message='Nudging: error -> no UWND variable in file u.nc'
               CALL abort_gcm(modname,abort_message,1)
              ENDIF
              write(*,*) trim(modname)//' ncidu,varidu',ncidu,varidu
-             if (ncidpl.eq.-99) ncidpl=ncidu
+             if (ncidpl==-99) ncidpl=ncidu
 
              status=NF90_INQ_DIMID(ncidu, "LONU", dimid)
              status=NF90_INQUIRE_DIMENSION(ncidu,dimid,namedim,lendim)
-             IF (lendim .NE. iip1) THEN
+             IF (lendim /= iip1) THEN
                 abort_message='dimension LONU different from iip1 in u.nc'
                 CALL abort_gcm(modname,abort_message,1)
@@ -1178 +1178 @@
              status=NF90_INQ_DIMID(ncidu, "LATU", dimid)
              status=NF90_INQUIRE_DIMENSION(ncidu,dimid,namedim,lendim)
-             IF (lendim .NE. jjp1) THEN
+             IF (lendim /= jjp1) THEN
                 abort_message='dimension LATU different from jjp1 in u.nc'
                 CALL abort_gcm(modname,abort_message,1)
@@ -1188 +1188 @@
          if (guide_v) then
              rcode = nf90_open('v.nc', nf90_nowrite, ncidv)
-             IF (rcode.NE.NF90_NOERR) THEN
+             IF (rcode/=NF90_NOERR) THEN
               abort_message='Nudging: error -> no file v.nc'
               CALL abort_gcm(modname,abort_message,1)
              ENDIF
              rcode = nf90_inq_varid(ncidv, 'VWND', varidv)
-             IF (rcode.NE.NF90_NOERR) THEN
+             IF (rcode/=NF90_NOERR) THEN
               abort_message='Nudging: error -> no VWND variable in file v.nc'
               CALL abort_gcm(modname,abort_message,1)
              ENDIF
              write(*,*) trim(modname)//' ncidv,varidv',ncidv,varidv
-             if (ncidpl.eq.-99) ncidpl=ncidv
+             if (ncidpl==-99) ncidpl=ncidv
              
              status=NF90_INQ_DIMID(ncidv, "LONV", dimid)
              status=NF90_INQUIRE_DIMENSION(ncidv,dimid,namedim,lendim)
              
-                IF (lendim .NE. iip1) THEN
+                IF (lendim /= iip1) THEN
                 abort_message='dimension LONV different from iip1 in v.nc'
                 CALL abort_gcm(modname,abort_message,1)
@@ -1211 +1211 @@
              status=NF90_INQ_DIMID(ncidv, "LATV", dimid)
              status=NF90_INQUIRE_DIMENSION(ncidv,dimid,namedim,lendim)
-             IF (lendim .NE. jjm) THEN
+             IF (lendim /= jjm) THEN
                 abort_message='dimension LATV different from jjm in v.nc'
                 CALL abort_gcm(modname,abort_message,1)
@@ -1221 +1221 @@
          if (guide_T) then
              rcode = nf90_open('T.nc', nf90_nowrite, ncidt)
-             IF (rcode.NE.NF90_NOERR) THEN
+             IF (rcode/=NF90_NOERR) THEN
               abort_message='Nudging: error -> no file T.nc'
               CALL abort_gcm(modname,abort_message,1)
              ENDIF
              rcode = nf90_inq_varid(ncidt, 'AIR', varidt)
-             IF (rcode.NE.NF90_NOERR) THEN
+             IF (rcode/=NF90_NOERR) THEN
               abort_message='Nudging: error -> no AIR variable in file T.nc'
               CALL abort_gcm(modname,abort_message,1)
              ENDIF
              write(*,*) trim(modname)//' ncidT,varidT',ncidt,varidt
-             if (ncidpl.eq.-99) ncidpl=ncidt
+             if (ncidpl==-99) ncidpl=ncidt
 
              status=NF90_INQ_DIMID(ncidt, "LONV", dimid)
              status=NF90_INQUIRE_DIMENSION(ncidt,dimid,namedim,lendim)
-             IF (lendim .NE. iip1) THEN
+             IF (lendim /= iip1) THEN
                 abort_message='dimension LONV different from iip1 in T.nc'
                 CALL abort_gcm(modname,abort_message,1)
@@ -1242 +1242 @@
              status=NF90_INQ_DIMID(ncidt, "LATU", dimid)
              status=NF90_INQUIRE_DIMENSION(ncidt,dimid,namedim,lendim)
-             IF (lendim .NE. jjp1) THEN
+             IF (lendim /= jjp1) THEN
                 abort_message='dimension LATU different from jjp1 in T.nc'
                 CALL abort_gcm(modname,abort_message,1)
@@ -1252 +1252 @@
          if (guide_Q) then
              rcode = nf90_open('hur.nc', nf90_nowrite, ncidQ)
-             IF (rcode.NE.NF90_NOERR) THEN
+             IF (rcode/=NF90_NOERR) THEN
               abort_message='Nudging: error -> no file hur.nc'
               CALL abort_gcm(modname,abort_message,1)
              ENDIF
              rcode = nf90_inq_varid(ncidQ, 'RH', varidQ)
-             IF (rcode.NE.NF90_NOERR) THEN
+             IF (rcode/=NF90_NOERR) THEN
               abort_message='Nudging: error -> no RH variable in file hur.nc'
               CALL abort_gcm(modname,abort_message,1)
              ENDIF
              write(*,*) trim(modname)//' ncidQ,varidQ',ncidQ,varidQ
-             if (ncidpl.eq.-99) ncidpl=ncidQ
+             if (ncidpl==-99) ncidpl=ncidQ
 
              status=NF90_INQ_DIMID(ncidQ, "LONV", dimid)
              status=NF90_INQUIRE_DIMENSION(ncidQ,dimid,namedim,lendim)
-             IF (lendim .NE. iip1) THEN
+             IF (lendim /= iip1) THEN
                 abort_message='dimension LONV different from iip1 in hur.nc'
                 CALL abort_gcm(modname,abort_message,1)
@@ -1273 +1273 @@
              status=NF90_INQ_DIMID(ncidQ, "LATU", dimid)
              status=NF90_INQUIRE_DIMENSION(ncidQ,dimid,namedim,lendim)
-             IF (lendim .NE. jjp1) THEN
+             IF (lendim /= jjp1) THEN
                 abort_message='dimension LATU different from jjp1 in hur.nc'
                 CALL abort_gcm(modname,abort_message,1)
@@ -1283 +1283 @@
          if ((guide_P).OR.(guide_modele)) then
              rcode = nf90_open('ps.nc', nf90_nowrite, ncidps)
-             IF (rcode.NE.NF90_NOERR) THEN
+             IF (rcode/=NF90_NOERR) THEN
               abort_message='Nudging: error -> no file ps.nc'
               CALL abort_gcm(modname,abort_message,1)
              ENDIF
              rcode = nf90_inq_varid(ncidps, 'SP', varidps)
-             IF (rcode.NE.NF90_NOERR) THEN
+             IF (rcode/=NF90_NOERR) THEN
               abort_message='Nudging: error -> no SP variable in file ps.nc'
               CALL abort_gcm(modname,abort_message,1)
@@ -1295 +1295 @@
          endif
 ! Coordonnee verticale
-         if (guide_plevs.EQ.0) then
+         if (guide_plevs==0) then
               rcode = nf90_inq_varid(ncidpl, 'LEVEL', varidpl)
-              IF (rcode.NE.0) rcode = nf90_inq_varid(ncidpl, 'PRESSURE', varidpl)
+              IF (rcode/=0) rcode = nf90_inq_varid(ncidpl, 'PRESSURE', varidpl)
               write(*,*) trim(modname)//' ncidpl,varidpl',ncidpl,varidpl
          endif
 ! Coefs ap, bp pour calcul de la pression aux differents niveaux
-         if (guide_plevs.EQ.1) then
+         if (guide_plevs==1) then
              status=NF90_GET_VAR(ncidpl,varidap,apnc,[1],[nlevnc])
              status=NF90_GET_VAR(ncidpl,varidbp,bpnc,[1],[nlevnc])
-         ELSEIF (guide_plevs.EQ.0) THEN
+         ELSEIF (guide_plevs==0) THEN
              status=NF90_GET_VAR(ncidpl,varidpl,apnc,[1],[nlevnc])
 !FC Pour les corrections la pression est deja en Pascals on commente la ligne ci-dessous
@@ -1329 +1329 @@
 
 ! Pression 
-     if (guide_plevs.EQ.2) then
+     if (guide_plevs==2) then
          status=NF90_GET_VAR(ncidp,varidp,pnat2,start,count)
          IF (invert_y) THEN
@@ -1421 +1421 @@
          write(*,*)trim(modname)//' : opening nudging files '
 ! Ap et Bp si niveaux de pression hybrides
-         if (guide_plevs.EQ.1) then
+         if (guide_plevs==1) then
            write(*,*)trim(modname)//' Reading nudging on model levels'
            rcode = nf90_open('apbp.nc', nf90_nowrite, ncidpl)
-           IF (rcode.NE.NF90_NOERR) THEN
+           IF (rcode/=NF90_NOERR) THEN
              abort_message='Nudging: error -> no file apbp.nc'
            CALL abort_gcm(modname,abort_message,1)
            ENDIF
            rcode = nf90_inq_varid(ncidpl, 'AP', varidap)
-           IF (rcode.NE.NF90_NOERR) THEN
+           IF (rcode/=NF90_NOERR) THEN
              abort_message='Nudging: error -> no AP variable in file apbp.nc'
            CALL abort_gcm(modname,abort_message,1)
            ENDIF
            rcode = nf90_inq_varid(ncidpl, 'BP', varidbp)
-           IF (rcode.NE.NF90_NOERR) THEN
+           IF (rcode/=NF90_NOERR) THEN
              abort_message='Nudging: error -> no BP variable in file apbp.nc'
              CALL abort_gcm(modname,abort_message,1)
@@ -1441 +1441 @@
          endif
 ! Pression
-         if (guide_plevs.EQ.2) then
+         if (guide_plevs==2) then
            rcode = nf90_open('P.nc', nf90_nowrite, ncidp)
-           IF (rcode.NE.NF90_NOERR) THEN
+           IF (rcode/=NF90_NOERR) THEN
              abort_message='Nudging: error -> no file P.nc'
              CALL abort_gcm(modname,abort_message,1)
            ENDIF
            rcode = nf90_inq_varid(ncidp, 'PRES', varidp)
-           IF (rcode.NE.NF90_NOERR) THEN
+           IF (rcode/=NF90_NOERR) THEN
              abort_message='Nudging: error -> no PRES variable in file P.nc'
              CALL abort_gcm(modname,abort_message,1)
            ENDIF
            write(*,*)trim(modname)//' ncidp,varidp',ncidp,varidp
-           if (ncidpl.eq.-99) ncidpl=ncidp
+           if (ncidpl==-99) ncidpl=ncidp
          endif
 ! Vent zonal
          if (guide_u) then
            rcode = nf90_open('u.nc', nf90_nowrite, ncidu)
-           IF (rcode.NE.NF90_NOERR) THEN
+           IF (rcode/=NF90_NOERR) THEN
              abort_message='Nudging: error -> no file u.nc'
              CALL abort_gcm(modname,abort_message,1)
            ENDIF
            rcode = nf90_inq_varid(ncidu, 'UWND', varidu)
-           IF (rcode.NE.NF90_NOERR) THEN
+           IF (rcode/=NF90_NOERR) THEN
              abort_message='Nudging: error -> no UWND variable in file u.nc'
              CALL abort_gcm(modname,abort_message,1)
            ENDIF
            write(*,*)trim(modname)//' ncidu,varidu',ncidu,varidu
-           if (ncidpl.eq.-99) ncidpl=ncidu
+           if (ncidpl==-99) ncidpl=ncidu
          endif
 ! Vent meridien
          if (guide_v) then
            rcode = nf90_open('v.nc', nf90_nowrite, ncidv)
-           IF (rcode.NE.NF90_NOERR) THEN
+           IF (rcode/=NF90_NOERR) THEN
              abort_message='Nudging: error -> no file v.nc'
              CALL abort_gcm(modname,abort_message,1)
            ENDIF
            rcode = nf90_inq_varid(ncidv, 'VWND', varidv)
-           IF (rcode.NE.NF90_NOERR) THEN
+           IF (rcode/=NF90_NOERR) THEN
              abort_message='Nudging: error -> no VWND variable in file v.nc'
              CALL abort_gcm(modname,abort_message,1)
            ENDIF
            write(*,*)trim(modname)//' ncidv,varidv',ncidv,varidv
-           if (ncidpl.eq.-99) ncidpl=ncidv
+           if (ncidpl==-99) ncidpl=ncidv
          endif
 ! Temperature
          if (guide_T) then
            rcode = nf90_open('T.nc', nf90_nowrite, ncidt)
-           IF (rcode.NE.NF90_NOERR) THEN
+           IF (rcode/=NF90_NOERR) THEN
              abort_message='Nudging: error -> no file T.nc'
              CALL abort_gcm(modname,abort_message,1)
            ENDIF
            rcode = nf90_inq_varid(ncidt, 'AIR', varidt)
-           IF (rcode.NE.NF90_NOERR) THEN
+           IF (rcode/=NF90_NOERR) THEN
              abort_message='Nudging: error -> no AIR variable in file T.nc'
              CALL abort_gcm(modname,abort_message,1)
            ENDIF
            write(*,*)trim(modname)//' ncidT,varidT',ncidt,varidt
-           if (ncidpl.eq.-99) ncidpl=ncidt
+           if (ncidpl==-99) ncidpl=ncidt
          endif
 ! Humidite
          if (guide_Q) then
            rcode = nf90_open('hur.nc', nf90_nowrite, ncidQ)
-           IF (rcode.NE.NF90_NOERR) THEN
+           IF (rcode/=NF90_NOERR) THEN
              abort_message='Nudging: error -> no file hur.nc'
              CALL abort_gcm(modname,abort_message,1)
            ENDIF
            rcode = nf90_inq_varid(ncidQ, 'RH', varidQ)
-           IF (rcode.NE.NF90_NOERR) THEN
+           IF (rcode/=NF90_NOERR) THEN
              abort_message='Nudging: error -> no RH,variable in file hur.nc'
              CALL abort_gcm(modname,abort_message,1)
            ENDIF
            write(*,*)trim(modname)//' ncidQ,varidQ',ncidQ,varidQ
-           if (ncidpl.eq.-99) ncidpl=ncidQ
+           if (ncidpl==-99) ncidpl=ncidQ
          endif
 ! Pression de surface
          if ((guide_P).OR.(guide_modele)) then
            rcode = nf90_open('ps.nc', nf90_nowrite, ncidps)
-           IF (rcode.NE.NF90_NOERR) THEN
+           IF (rcode/=NF90_NOERR) THEN
              abort_message='Nudging: error -> no file ps.nc'
              CALL abort_gcm(modname,abort_message,1)
            ENDIF
            rcode = nf90_inq_varid(ncidps, 'SP', varidps)
-           IF (rcode.NE.NF90_NOERR) THEN
+           IF (rcode/=NF90_NOERR) THEN
              abort_message='Nudging: error -> no SP variable in file ps.nc'
              CALL abort_gcm(modname,abort_message,1)
@@ -1530 +1530 @@
          endif
 ! Coordonnee verticale
-         if (guide_plevs.EQ.0) then
+         if (guide_plevs==0) then
            rcode = nf90_inq_varid(ncidpl, 'LEVEL', varidpl)
-           IF (rcode.NE.0) rcode = nf90_inq_varid(ncidpl, 'PRESSURE', varidpl)
+           IF (rcode/=0) rcode = nf90_inq_varid(ncidpl, 'PRESSURE', varidpl)
            write(*,*)trim(modname)//' ncidpl,varidpl',ncidpl,varidpl
          endif
 ! Coefs ap, bp pour calcul de la pression aux differents niveaux
-         if (guide_plevs.EQ.1) then
+         if (guide_plevs==1) then
              status=NF90_GET_VAR(ncidpl,varidap,apnc,[1],[nlevnc])
              status=NF90_GET_VAR(ncidpl,varidbp,bpnc,[1],[nlevnc])
-         elseif (guide_plevs.EQ.0) THEN
+         elseif (guide_plevs==0) THEN
              status=NF90_GET_VAR(ncidpl,varidpl,apnc,[1],[nlevnc])
              apnc=apnc*100.! conversion en Pascals
@@ -1563 +1563 @@
 
 !  Pression
-     if (guide_plevs.EQ.2) then
+     if (guide_plevs==2) then
          status=NF90_GET_VAR(ncidp,varidp,zu,start,count)
          DO i=1,iip1
@@ -1629 +1629 @@
 
 !  Pression de surface
-     if ((guide_P).OR.(guide_plevs.EQ.1))  then
+     if ((guide_P).OR.(guide_plevs==1))  then
          start(3)=timestep
          start(4)=0
@@ -1681 +1681 @@
 
     write(*,*)trim(modname)//': output timestep',timestep,'var ',varname
-    IF (timestep.EQ.0) THEN 
+    IF (timestep==0) THEN
 ! ----------------------------------------------
 ! initialisation fichier de sortie
@@ -1828 +1828 @@
     do l=1,nl
         do i=2,iim-1
-            if(abs(x(i,l)).gt.1.e10) then
+            if(abs(x(i,l))>1.e10) then
                zz=0.5*(x(i-1,l)+x(i+1,l))
               print*,'correction ',i,l,x(i,l),zz