Changeset 2567 for trunk/LMDZ.MARS/util/localtime.F90

Timestamp:

Oct 5, 2021, 8:34:09 AM (3 years ago)

Author:

emillour

Message:

Mars GCM utilities:
Minor fixes to run with picky gfortran 10.3.0 which requires one element arrays
(rather than scalars) when calling NetCDF routines, andf that stop statements
should not be followed by strings. While at it replaced tabs with spaces.
EM

File:

• trunk/LMDZ.MARS/util/localtime.F90 (modified) (41 diffs)

trunk/LMDZ.MARS/util/localtime.F90

@@ -12 +12 @@
 include "netcdf.inc" ! NetCDF definitions
 
-character (len=50)  file
+character (len=100)  file
 ! file(): input file(s) names(s)
 character (len=30), dimension(16) :: notprocessed
 ! notprocessed(): names of the (16) variables that won't be processed
-character (len=50), dimension(:), allocatable :: var
+character (len=100), dimension(:), allocatable :: var
 ! var(): name(s) of variable(s) that will be processed
 character (len=100) :: tmpvar,long_name,units
@@ -121 +121 @@
 if (len_trim(file).eq.0) then
    write(*,*) "no file... game over"
-   stop ""
+   stop
 endif
 
@@ -132 +132 @@
    write(*,*) 'ERROR: Pb opening file '//trim(file)
    write(*,*) NF_STRERROR(ierr)
-   stop ""
+   stop
 endif
 
@@ -140 +140 @@
    write(*,*) 'ERROR: Pb with NF_INQ_NVARS'
    write(*,*) NF_STRERROR(ierr)
-   stop ""
+   stop
 endif
 
@@ -199 +199 @@
 write(*,*) "(an empty line , i.e: just <Enter>, implies end of list)"
 nbvar=0
-read(*,'(a50)') tmpvar
+read(*,'(a100)') tmpvar
 do while ((tmpvar/=' ').AND.(trim(tmpvar)/='all'))
    nbvar=nbvar+1
    var(nbvar)=tmpvar
-   read(*,'(a50)') tmpvar
+   read(*,'(a100)') tmpvar
 enddo
 
@@ -219 +219 @@
 else if(nbvar==0) then
    write(*,*) "no variable... game over"
-   stop ""
+   stop
 endif ! of if (tmpvar=="all")
 
@@ -238 +238 @@
       write(*,*) 'ERROR: Pb opening file '//trim(file)
       write(*,*) NF_STRERROR(ierr)
-      stop ""
+      stop
    endif
 
@@ -251 +251 @@
          if (ierr.NE.NF_NOERR) then
             write(*,*) 'ERROR: Dimension <latitude> and <lat> is missing in file '//trim(file)
-            stop ""  
+            stop  
          endif
    endif
@@ -260 +260 @@
       if (ierr.NE.NF_NOERR) then
         write(*,*) 'ERROR: Field <latitude> and <lat> is missing in file '//trim(file)
-        stop ""  
+        stop  
       endif
    endif
@@ -272 +272 @@
          if (ierr.NE.NF_NOERR) then
             write(*,*) 'ERROR: Dimension <longitude> and <lon> is missing in file '//trim(file)
-            stop ""  
+            stop  
          endif 
    endif
@@ -281 +281 @@
       if (ierr.NE.NF_NOERR) then
         write(*,*) 'ERROR: Field <longitude> and <lon> is missing in file '//trim(file)
-        stop ""  
+        stop  
       endif
    endif
@@ -290 +290 @@
    if (ierr.NE.NF_NOERR) then
       write(*,*) 'ERROR: Dimension <altitude> is missing in file '//trim(file)
-      stop ""  
+      stop  
    endif
    ierr=NF_INQ_VARID(nid,"altitude",altvar)
    if (ierr.NE.NF_NOERR) then
       write(*,*) 'ERROR: Field <altitude> is missing in file '//trim(file)
-      stop ""
+      stop
    endif
    ierr=NF_INQ_DIMLEN(nid,altdim,altlen)
@@ -321 +321 @@
          ctllen=0
       endif
-      !stop ""  
+      !stop  
    endif
 !   ierr=NF_INQ_DIMID(nid,"controle_axe",ctldim)
 !   if (ierr.NE.NF_NOERR) then
-!      !stop ""  
+!      !stop  
 !   endif
    ierr=NF_INQ_VARID(nid,"controle",ctlvar)
@@ -331 +331 @@
       write(*,*) 'Field <controle> is missing in file '//trim(file)
       ctllen=0
-      !stop ""
+      !stop
    else
       ierr=NF_INQ_DIMLEN(nid,ctldim,ctllen)
@@ -414 +414 @@
      if (ierr.NE.NF_NOERR) then
         write(*,*) 'ERROR: Dimension <Time> or is missing in file'//trim(file)
-        stop ""
+        stop
      endif
    endif
@@ -423 +423 @@
      if (ierr.NE.NF_NOERR) then
         write(*,*) 'ERROR: Field <Time> is missing in file '//trim(file)
-        stop ""
+        stop
      endif
    endif
@@ -556 +556 @@
 
      do it=1,timelen_lt
-        ierr= NF_PUT_VARA_REAL(nout,timevarout,it,1,lt_hour(it))
+        ierr= NF_PUT_VARA_REAL(nout,timevarout,(/it/),(/1/),lt_hour(it))
      enddo
    else
      do it=1,timelen_lt
-        ierr= NF_PUT_VARA_REAL(nout,timevarout,it,1,lt_out(it))
+        ierr= NF_PUT_VARA_REAL(nout,timevarout,(/it/),(/1/),lt_out(it))
         if (ierr.NE.NF_NOERR) then
           write(*,*) "Error , failed to write Time"
@@ -582 +582 @@
       if (ierr.NE.NF_NOERR) then
          write(*,*) 'ERROR: Field <',var(j),'> not found in file'//file
-         stop ""
+         stop
       endif
       ierr=nf_inq_varndims(nid,varid,ndim)
@@ -718 +718 @@
       if (ierr.ne.NF_NOERR) then
          write(*,*) 'PUT_VAR ERROR: ',NF_STRERROR(ierr)
-         stop ""
+         stop
       endif
 
@@ -820 +820 @@
    WRITE(*,*)'ERROR: Impossible to create the file ',trim(filename)
    write(*,*) NF_STRERROR(ierr)
-   stop ""
+   stop
 endif
 
@@ -831 +831 @@
    WRITE(*,*)'initiate: error failed to define dimension <latitude>'
    write(*,*) NF_STRERROR(ierr)
-   stop ""
+   stop
 endif
 ierr = NF_DEF_DIM(nout, "longitude", lonlen, londimout)
@@ -837 +837 @@
    WRITE(*,*)'initiate: error failed to define dimension <longitude>'
    write(*,*) NF_STRERROR(ierr)
-   stop ""
+   stop
 endif
 ierr = NF_DEF_DIM(nout, "altitude", altlen, altdimout)
@@ -843 +843 @@
    WRITE(*,*)'initiate: error failed to define dimension <altitude>'
    write(*,*) NF_STRERROR(ierr)
-   stop ""
+   stop
 endif
 if (size(ctl).ne.0) then
@@ -850 +850 @@
     WRITE(*,*)'initiate: error failed to define dimension <index>'
     write(*,*) NF_STRERROR(ierr)
-    stop ""
+    stop
   endif
 endif
@@ -857 +857 @@
    WRITE(*,*)'initiate: error failed to define dimension <Time>'
    write(*,*) NF_STRERROR(ierr)
-   stop ""
+   stop
 endif
 
@@ -867 +867 @@
    WRITE(*,*)'initiate: error failed to switch out of define mode'
    write(*,*) NF_STRERROR(ierr)
-   stop ""
+   stop
 endif
 
@@ -888 +888 @@
    WRITE(*,*)'initiate: error failed writing variable <latitude>'
    write(*,*) NF_STRERROR(ierr)
-   stop ""
+   stop
 endif
 
@@ -903 +903 @@
    WRITE(*,*)'initiate: error failed writing variable <longitude>'
    write(*,*) NF_STRERROR(ierr)
-   stop ""
+   stop
 endif
 
@@ -926 +926 @@
    WRITE(*,*)'initiate: error failed writing variable <altitude>'
    write(*,*) NF_STRERROR(ierr)
-   stop ""
+   stop
 endif
 
@@ -948 +948 @@
       WRITE(*,*)'initiate: error failed writing variable <controle>'
       write(*,*) NF_STRERROR(ierr)
-      stop ""
+      stop
    endif
 endif
@@ -1014 +1014 @@
   hybrid=.true.
   if (ierr.ne.NF_NOERR) then
-    stop "init2 Error: Failed reading aps"
+    write(*,*) "init2 Error: Failed reading aps"
+    stop
   endif
 
@@ -1031 +1032 @@
     ierr=NF_GET_VAR_REAL(infid,tmpvarid,bps)
     if (ierr.ne.NF_NOERR) then
-      stop "init2 Error: Failed reading bps"
+      write(*,*) "init2 Error: Failed reading bps"
+      stop
     endif
   endif
@@ -1050 +1052 @@
   ierr=NF_GET_VAR_REAL(infid,tmpvarid,phisinit)
   if (ierr.ne.NF_NOERR) then
-    stop "init2 Error: Failed reading phisinit"
+    write(*,*) "init2 Error: Failed reading phisinit"
+    stop
   endif
   phis = .true.
@@ -1069 +1072 @@
                (/layerdimout/),apsid,ierr)
   if (ierr.ne.NF_NOERR) then
-    stop "Error: Failed to def_var aps"
+    write(*,*) "Error: Failed to def_var aps"
+    stop
   endif
 
@@ -1075 +1079 @@
   ierr=NF_PUT_VAR_REAL(outfid,apsid,aps)
   if (ierr.ne.NF_NOERR) then
-    stop "Error: Failed to write aps"
+    write(*,*) "Error: Failed to write aps"
+    stop
   endif
   
@@ -1082 +1087 @@
                (/layerdimout/),bpsid,ierr)
   if (ierr.ne.NF_NOERR) then
-    stop "Error: Failed to def_var bps"
+    write(*,*) "Error: Failed to def_var bps"
+    stop
   endif
 
@@ -1088 +1094 @@
   ierr=NF_PUT_VAR_REAL(outfid,bpsid,bps)
   if (ierr.ne.NF_NOERR) then
-    stop "Error: Failed to write bps"
+    write(*,*) "Error: Failed to write bps"
+    stop
   endif
   
@@ -1105 +1112 @@
               (/londimout,latdimout/),phisinitid,ierr)
   if (ierr.ne.NF_NOERR) then
-     stop "Error: Failed to def_var phisinit"
+     write(*,*) "Error: Failed to def_var phisinit"
+     stop
   endif
 
@@ -1111 +1119 @@
   ierr=NF_PUT_VAR_REAL(outfid,phisinitid,phisinit)
   if (ierr.ne.NF_NOERR) then
-    stop "Error: Failed to write phisinit"
+    write(*,*) "Error: Failed to write phisinit"
+    stop
   endif
 
@@ -1226 +1235 @@
      !write(*,*) 'NF_NOERR', NF_NOERR
      !CALL abort
-     stop ""
+     stop
   endif
 endif
@@ -1238 +1247 @@
   !    WRITE(*,*) 'NF_NOERR', NF_NOERR
   !    CALL abort
-     stop ""
+     stop
   endif
 endif